Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,764,701 G→A 5.7% intergenic (‑315/+12) sufA ← / ← rydB Fe‑S cluster assembly protein/novel sRNA, function unknown

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,764,7010GA 5.7% 219.2 / 3.3 70intergenic (‑315/+12)sufA/rydBFe‑S cluster assembly protein/novel sRNA, function unknown
Reads supporting (aligned to +/- strand):  ref base G (39/27);  new base A (2/2);  total (41/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.70e-01

AAAACATGCTTTTTGCCTGCTTTTAATAAGCATAGACCCTGATGTGTGGGTTAACAGGCACGCTAAGAATTAGGTATCTCATTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAA‑GGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCGGCATTAAGTCAGGATCAATGCACGCCCAGACGCCAGGCAAAGTAG  >  NC_000913/1764479‑1764964
                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                        
aaaaCATGCTTTTTGCCTGCTTTTAATAAGCATAGACCCTGATGTGTGGGTTAACAGGCACGCTAAGAATTAGGTATCTCATTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAA‑GGGTTTCTACCCCATCCGGCGCTTATCTCCg                                                                                                                                                                                                                                     >  4:268524/1‑259 (MQ=255)
    cATGCTTTTTGCCTGCTTTTAATAAGCATAGACCCTGATGTGTGGGTTAACAGGCACGCTAAGAATTAGGTATCTCATTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAA‑GGGTTTCTACCCCATCCGGCGCTTATCTCCGGCAc                                                                                                                                                                                                                                 <  2:268528/259‑1 (MQ=255)
                          taaGCATAGACCCTGATGGGTGGGTTAACAGGCACGCTAAGAATTAGGTAACTCATTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAA‑GGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGaataa                                                                                                                                                                                           <  2:169067/275‑1 (MQ=255)
                                        gATGTGTGGGTTAACAGGGACGCTAAGAATTAGGTATCTCATTGGTAGATGATGGGTATTATCTAAGAAGCATTCAGGCTGTGTTGAAAAAGACAAAGTTTTATTGAAAATGACTATTTAAGACATAGGTAAAAAAGTCCAGCGTTCAGAAATACGAAAACCCTTAAGTCTGTGCGACACAGGCTTAA‑AGGTTTCTACTCAATCCGGCTCTTATCCACGGCACCCTCAGTGGCTTGCATCTTGAAGGGGCGGAACTAATAATCCCAGAAAGCaaa                                                                                                                                                                             >  4:190518/1‑275 (MQ=255)
                                                             nttAAGAATTAGGTATCTTATTGTTAGATAATGGTTATTATCTAAGGAGCATTANGGCTNTGATGACACAGAAAAAGTTGTATNGAAAATGACTATTTAAGAGATAGGTACAAAAAGGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGGGCGACACAGGCTTAA‑GGGTTTCTACCCCCTCCGGCGCTTATCTCCGGCACTCTCa                                                                                                                                                                                                                            <  4:508080/205‑1 (MQ=255)
                                                                             ctcATTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAA‑GGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTAttt                                                                                                                                        >  2:202752/1‑275 (MQ=255)
                                                                               cATTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGCCTGTGTCGCACAGACTTaa                                                                                                                                                                                                                                                                     <  8:330053/149‑1 (MQ=255)
                                                                                  tGTTAGATAATGGTAATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTAATGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTCCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAA‑GGGTTTCTACCCCAACCGGAGCGTATCTCCGGCACTCTCAGTGCCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGAGTAGTATTATTTTTGaa                                                                                                                                   <  2:190520/275‑1 (MQ=255)
                                                                                    ttAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAA‑g                                                                                                                                                                                                                                                                   <  7:357297/145‑1 (MQ=255)
                                                                                    ttAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAA‑g                                                                                                                                                                                                                                                                   >  8:357297/1‑145 (MQ=255)
                                                                                          aaTGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAA‑GGGTTTCTAcc                                                                                                                                                                                                                                                         >  5:277828/1‑149 (MQ=255)
                                                                                           aTGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAA‑GGGTTTCTAccc                                                                                                                                                                                                                                                        >  8:259989/1‑149 (MQ=255)
                                                                                                tattatCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGc                                                                                                                                                                                                                                                                         <  8:65289/128‑1 (MQ=255)
                                                                                                tattatCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGc                                                                                                                                                                                                                                                                         >  7:65289/1‑128 (MQ=255)
                                                                                                  ttatCTAAGGAGCATTAAGGGTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGAGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAA‑GGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGtttt                                                                                                                                                                  <  1:85789/228‑1 (MQ=255)
                                                                                                         aGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAAGGGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGc                                                                                                             >  2:421084/1‑275 (MQ=255)
                                                                                                         aGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAA‑GGGTTTCTACCCCATCCGGCGCTTAt                                                                                                                                                                                                                                          >  8:237691/1‑149 (MQ=255)
                                                                                                               aTTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAA‑GGGTTTCTACCCCATCCGGCGCTTATCTCCgg                                                                                                                                                                                                                                    >  2:533131/1‑149 (MQ=255)
                                                                                                               aTTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAA‑GGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAaca                                                                                                                                                              >  4:263573/1‑219 (MQ=255)
                                                                                                                   aGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAA‑GGGTTTCTACCCCATCCGGCGCTTATCTCCGGCAct                                                                                                                                                                                                                                <  4:533126/149‑1 (MQ=255)
                                                                                                                   aGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAA‑GGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACaaa                                                                                                                                                          <  2:263577/219‑1 (MQ=255)
                                                                                                                        gtgATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAA‑GGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCa                                                                                             >  2:320823/1‑275 (MQ=255)
                                                                                                                         tgatgaAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAA‑GGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTAtttt                                                                                                                                       >  2:150496/1‑232 (MQ=255)
                                                                                                                            tgaAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAA‑GGGtt                                                                                                                                                                                                                                                               <  5:218507/109‑1 (MQ=255)
                                                                                                                            tgaAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAA‑GGGtt                                                                                                                                                                                                                                                               >  6:218507/1‑109 (MQ=255)
                                                                                                                             gaaaaagaaaaaGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAA‑GGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGaa                                                                                                                                   <  4:150495/232‑1 (MQ=255)
                                                                                                                              aaaaagaaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGCCTGTGTCGCACAGACTTaa                                                                                                                                                                                                                                                                     >  6:119842/1‑102 (MQ=38)
                                                                                                                                 aagaaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAA‑GGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGc                                                                                                                                                                                                                  <  6:277828/149‑1 (MQ=255)
                                                                                                                                   gaaaaagTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAA‑GGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGctc                                                                                                                                                                                                                >  6:156941/1‑149 (MQ=255)
                                                                                                                                          ttGTATTGAAAATGACTATTTCAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAA‑GGGTTTCTACCCCATCCGGCGCTTATCTCCGGCAc                                                                                                                                                                                                                                 <  6:350638/125‑1 (MQ=255)
                                                                                                                                          ttGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAA‑GGGTTTCTACCCCATCCGGCGCTTATCTCCGGCAc                                                                                                                                                                                                                                 >  5:350638/1‑125 (MQ=255)
                                                                                                                                                   aaaTGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAA‑GGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAg                                                                                                                                                                                                >  8:241919/1‑149 (MQ=255)
                                                                                                                                                         cTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAA‑GGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGaataat                                                                                                                                                                                          >  7:227148/1‑149 (MQ=255)
                                                                                                                                                               aaGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAA‑GGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATa                                                                                                                                                                                    >  7:34521/1‑149 (MQ=255)
                                                                                                                                                                  agaTAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAA‑GGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTAtt                                                                                                       >  2:327424/1‑223 (MQ=255)
                                                                                                                                                                    aTAGGTAAAAAAGTGCAGCGTTCAGACATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAA‑GGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAg                                                                                                                                                                                                                           <  8:342650/105‑1 (MQ=255)
                                                                                                                                                                    aTAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAA‑GGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAg                                                                                                                                                                                                                           >  7:342650/1‑105 (MQ=255)
                                                                                                                                                                    aTAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAA‑GGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCt                                                                                                                                                                               >  8:359505/1‑149 (MQ=255)
                                                                                                                                                                      aGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAA‑GGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGctct                                                                                                                                                                                                               >  2:364629/1‑115 (MQ=255)
                                                                                                                                                                       ggTAAAAAAGTGCAGCGTTCAGCCATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAA‑GTGTTTCTACCCCATCCGTCGCTTACCTCCGGCACTCTCACTGGCTTAGCTCTTGAAGTGGCGATAAGAATAATCTCATAATGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGTCGTTAGATCAGCAGCATTAAGCAGGTAATCGAACCTGGACGGAGTTGAATACAGGTCGATTTGGCGAACTTGCg                                               >  4:90911/1‑274 (MQ=255)
                                                                                                                                                                         tAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAA‑GGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGaataat                                                                                                                                                                                          >  7:105913/1‑133 (MQ=255)
                                                                                                                                                                         tAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAA‑GGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGaataat                                                                                                                                                                                          <  8:105913/133‑1 (MQ=255)
                                                                                                                                                                         tAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAA‑GGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAAccc                                                                                                                                                                          >  5:449585/1‑149 (MQ=255)
                                                                                                                                                                          aaaaaaGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAA‑GGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGaa                                                                                                                                                                                                           <  4:364624/115‑1 (MQ=255)
                                                                                                                                                                                     gCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAA‑GGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAaca                                                                                                                                                              <  5:156941/149‑1 (MQ=255)
                                                                                                                                                                                     gCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAA‑GGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAaca                                                                                                                                                              <  7:241919/149‑1 (MQ=255)
                                                                                                                                                                                      cGTTCATAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAA‑GGGTTTCTACCCCATCTGGCGCTTATCTCCGGCACTCTCAGTGGCt                                                                                                                                                                                                                      <  2:382731/92‑1 (MQ=255)
                                                                                                                                                                                        ttCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAA‑GGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACaa                                                                                                                                                           >  5:140827/1‑149 (MQ=255)
                                                                                                                                                                                        ttCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAA‑GGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTAtt                                                                                                       >  4:76067/1‑201 (MQ=255)
                                                                                                                                                                                            gAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAA‑GGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTaaga                                                                                                   <  2:76068/201‑1 (MQ=255)
                                                                                                                                                                                            acaaTAAGATAACCCTGAAGTCTGTGCGACACAGGTTTAA‑GTGGTTCTCCCCCCTCCGTCGTTTTGCTCCGGCACTCTCAGTTGTTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCGGCATTAAGTCAGGATCAATGCa                         <  4:202749/273‑1 (MQ=255)
                                                                                                                                                                                             aaaTAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAA‑GGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGc                                                                                                                                                      >  6:277682/1‑149 (MQ=255)
                                                                                                                                                                                             aaaTAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAA‑GGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGc                                                                                                                                                      <  7:237691/149‑1 (MQ=255)
                                                                                                                                                                                                   gAAAACCCTTAAGTCTGTGCGACACAGGCTTAA‑GGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAg                                                                                                                                                >  5:381079/1‑149 (MQ=255)
                                                                                                                                                                                                           ttAAGCCTGTGTCGCACAGACTTAA‑GGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACccgcc                                                                                                                                                                       <  7:192434/118‑1 (MQ=38)
                                                                                                                                                                                                           ttAAGCCTGTGTCGCACAGACTTAA‑GGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACccgcc                                                                                                                                                                       >  8:192434/1‑118 (MQ=39)
                                                                                                                                                                                                            taggTCTGTGCGACACAGGCTTAA‑GGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTAtttt                                                                                                                                       <  8:34521/146‑1 (MQ=255)
                                                                                                                                                                                                             aaGTCTGTGCGACACAGGCTTAA‑GGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTGGCAGGGGCGATAAGATTAATCTCATAAAGCTAACACGCCGTTTTAACACAAACGGCGATTAGTATTATTTTTGAACAATAACAGGCGGTAGATAAGCAGTATTAAGAAGGTAATTGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCCAACTTGCCGCATTAATTCAGGAGCAATGCACGCCCAAACGACaaggc        >  4:428628/1‑271 (MQ=255)
                                                                                                                                                                                                                tCTGTGCGACACAGGCTTAA‑GGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGaa                                                                                                                                   <  6:381079/149‑1 (MQ=255)
                                                                                                                                                                                                                tCTGTGCGACACAGGCTTAA‑GGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGaa                                                                                                                                   >  5:207384/1‑149 (MQ=255)
                                                                                                                                                                                                                aggcacagACCCAGGTGTAA‑GGGTTCNNNCCCCGTCCGGGCTTTTTCTCCGGCACTCTCTTTGTCTTGACTCTTGCAGGGGCTATCAGAATAATCTCATAAATCTCACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAAAATCCGGCGGTAGATAAGCAGTATTAAGAAAGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCGGCATTAAGTCAGGATCAATGCACGCCCAGACGCCAGGCa        <  3:154004/265‑1 (MQ=255)
                                                                                                                                                                                                                   gtgCGACACAGGCTTAA‑GGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAAc                                                                                                                                                                            >  1:33930/1‑105 (MQ=255)
                                                                                                                                                                                                                   gtgCGACACAGGCTTAA‑GGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTaa                                                                                                                                                                >  2:308621/1‑117 (MQ=255)
                                                                                                                                                                                                                   gtgCGACACAGGCTTAA‑GGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACACTATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCGGCATTAATTCAGGATCAATGCCCGCCAAGACTCCAGCCAAAGTAg  >  2:4249/1‑275 (MQ=255)
                                                                                                                                                                                                                   gtgCGACACAGGCTTAA‑GGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCt                                                                                     >  1:38582/1‑192 (MQ=255)
                                                                                                                                                                                                                   gtgCGACACAGGCTTAA‑GGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCgg                                              >  4:329982/1‑231 (MQ=255)
                                                                                                                                                                                                                   gtgCGACACAGGCTTAA‑GGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTANTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTANATTTGGCGAATTTGCGGCATTAAGTTAGGANNAATGCACGCACAGCCCCCAGGCAAAATAg  >  3:57771/1‑275 (MQ=255)
                                                                                                                                                                                                                       gACACAGGCTTAA‑GGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACccgc                                                                                                                                                                        <  3:33930/105‑1 (MQ=255)
                                                                                                                                                                                                                       gACACAGGCTTAA‑GGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCGGCAtt                                          <  2:329986/231‑1 (MQ=255)
                                                                                                                                                                                                                           caGGCTTAA‑GGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGc                                                                                                                        >  8:101641/1‑149 (MQ=255)
                                                                                                                                                                                                                              gCTTAA‑GGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGt                                                                                                                     <  5:277682/149‑1 (MQ=255)
                                                                                                                                                                                                                              gCTTAA‑GGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCCCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGCACACNATCAGGCGGTAGa                                                                                                                  >  3:67127/1‑152 (MQ=255)
                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                        
AAAACATGCTTTTTGCCTGCTTTTAATAAGCATAGACCCTGATGTGTGGGTTAACAGGCACGCTAAGAATTAGGTATCTCATTGTTAGATAATGGTTATTATCTAAGGAGCATTAAGGCTGTGATGAAAAAGAAAAAGTTGTATTGAAAATGACTATTTAAGAGATAGGTAAAAAAGTGCAGCGTTCAGAAATAAGAAAACCCTTAAGTCTGTGCGACACAGGCTTAA‑GGGTTTCTACCCCATCCGGCGCTTATCTCCGGCACTCTCAGTGGCTTAGCTCTTGAAGGGGCGATAAGAATAATCTCATAAAGCTAACCCGCCGTTTTAACACAAACTGCGATTAGTATTATTTTTGAACAATATCAGGCGGTAGATAAGCAGTATTAAGAAGGTCATCGAACCTGGACGGAGGTTAATCCAGGTCGATTTGGCGAACTTGCGGCATTAAGTCAGGATCAATGCACGCCCAGACGCCAGGCAAAGTAG  >  NC_000913/1764479‑1764964

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑