Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 809,217 2 bp→AG coding (40‑41/1290 nt) bioA ← adenosylmethionine‑8‑amino‑7‑oxononanoate aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913809,2170CA88.7% 194.3 / 7.9 72W14M (TGG→CTG) bioAadenosylmethionine‑8‑amino‑7‑oxononanoate aminotransferase
Reads supporting (aligned to +/- strand):  ref base C (5/3);  new base A (22/41);  total (27/45)
Fisher's exact test for biased strand distribution p-value = 2.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_000913809,2180AG90.1% 204.4 / 7.7 72W14M (TGG→CTG) bioAadenosylmethionine‑8‑amino‑7‑oxononanoate aminotransferase
Reads supporting (aligned to +/- strand):  ref base A (3/4);  new base G (24/40);  total (27/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.69e-01

GTTTGCGGCACAGCTCAATGGCTGGCGCATGGGTGATACCGCCAAACATCACATGCGACATGGCATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGA‑TTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCCAG‑ATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCACGCTGGACATTGTCGCAAGTCACAGAATTATTTGAAAAACCGTTGCTGGATCTGCTGTTTGAAGCGCAGCAGGTGCATCGCCAGCATTTCGATCCTCGT  >  NC_000913/808975‑809457
                                                                                                                                                                                                                                                   ||                                                                                                                                                                                                                                                
gTTTGCGGCACAGCTCAATGGCTGGCGCATGGGTGATACCGCCAAACATCACATGCGACATGGCATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGA‑TTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGG‑ATATgg                                                                                                                                                                                                                                          <  1:236518/251‑1 (MQ=255)
   tGCGGCACAGCTCAATGGCTGTCGCATGGGTGATACCGCCAAACATTACATGCGACATGGCATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGA‑TTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCTGCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGG‑ATATGGCg                                                                                                                                                                                                                                        <  2:617872/250‑1 (MQ=255)
      ggCACAGCTCTAGGGCTGGCGCATGGGCGATACCGCCAAACATCACATGCGACATGGCATCAATGTGCGACTTCATCGCAGCATGAAGCTGCGGGTGA‑TTGTAGCCGTGGATCGCCGCCCACCAGGACGACAGACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGG‑ATATGGCGTTgg                                                                                                                                                                                                                                    <  2:124205/251‑1 (MQ=255)
            gCTCAATGGCTGGCGCATGGGTGATACCGCCAAACATCACATGCGACATGGCATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGA‑TTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGG‑ATATGGCGTTGGTCAAAg                                                                                                                                                                                                                              <  1:519309/251‑1 (MQ=255)
             cTCAATGGCTGGCGCATGGGTGATACCGCCAAACATCACATGCGACATGGCATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGA‑TTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGG‑ATATGGCGTTGGTCAAAgg                                                                                                                                                                                                                             >  1:333425/1‑251 (MQ=255)
               cAATGGCTGGCGCATGGGTGATACCGCCAAACATCACATGCGACATGGCATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGA‑TTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGG‑ATATGGCGTTGGTCAAAgg                                                                                                                                                                                                                             <  1:265268/249‑1 (MQ=255)
               cAATGGCTGGCGCATGGGTGATACCGCCAAACATCACATGCGACATGGCATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGA‑TTGTAGCCGTGGAACGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGCGGTCATGGATGTGTATGGGTGCAGG‑ATATGGCGTTGGGCAAAGGc                                                                                                                                                                                                                            >  2:556030/1‑250 (MQ=255)
                      tggcGCATGGGTGATACCGCCAAACATCACATGCGACATGGCATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGA‑TTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGG‑ATCTGGCGTTGGTCAAAGGCAAGAtagt                                                                                                                                                                                                                    >  2:290827/1‑248 (MQ=255)
                       ggcGCATGGGTGATACCGCCAAACATCACATGCGACATGGCATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGA‑TTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGG‑ATATGGCGTTGGTCAAAGGCAAGATCg                                                                                                                                                                                                                     <  1:186243/249‑1 (MQ=255)
                       ggcGCATGGGTGATACCGCCAAACATCACATGCGACATGGCATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGA‑TTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGG‑ATATGGCGTTGGTCAAAGGCAAGATCg                                                                                                                                                                                                                     <  2:237124/249‑1 (MQ=255)
                          gcATGGGTGATACCGCCAAACATCACATGCGACATGGCATCAATTTGCGACTTCATCGCCGCATTAAGCTGTGGGTGA‑TTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCTCCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGG‑ATATGGCGTTGGTCAAAGGCAAGATCGTCCGt                                                                                                                                                                                                                <  1:344119/251‑1 (MQ=255)
                                  gATACCGCCAAACATCACATGCGACATGGCAGCAATTTGCGTCTTCATCGCCGCATTAAGCTGCGGGTGA‑TTGTAGCCGTGGATCGCCGCCCCCCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCATCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGG‑ATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATa                                                                                                                                                                                                         <  1:674155/250‑1 (MQ=255)
                                          caaacatgtaagtcgtattggccgcAATTTGCAACTTCATCGCCGAATCAAGATGCGGTCGA‑TTGTTGCCGTGGATCGCCGCCCACCCGGACGACATACTGTCAATCAGGCGTCTGGCGTCAGACAAAATCAGCTAGCAACTTTCGGCGCTCACCACCGGATAAACCGGCAGAAGGGAGATTAGGGATGGGTATGGGTGCAGG‑ATATGTCGTTTGTCAAAGGCTAGATCGTCCGTTGTCATAATCGACTTg                                                                                                                                                                                                <  2:649944/227‑1 (MQ=255)
                                          cAAACATCACATGCGATATGGCATCACTTTGCGACTTCCTCGCCGCATTAAGCTGCGGGTGA‑TTGTAGACGTGGATCGCCGCCCACCATGACGACATACCGTCAACCAGGCGTATGCCGTCAGACAAAATCCGCTCGCAACATTAGGCGCTCACCACCGCATAACCCCGCAGAGGGGACGTCATGGATGTGTATGGGTGCCAG‑ATATGGCGTTGGGCAAAGGCAAGATCATCAGGTGTCATAATCGACttt                                                                                                                                                                                                >  2:63668/1‑250 (MQ=255)
                                                   caTGCGACATGGCAGCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGA‑TTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACTAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGG‑ATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCaa                                                                                                                                                                                        <  1:686617/250‑1 (MQ=255)
                                                      gCGACATGGCATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGA‑TTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGG‑ATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTg                                                                                                                                                                                    <  1:482530/251‑1 (MQ=255)
                                                      gCGACATGGCATCAATTTGCGACTTCATCGCCGCATGAAGCTGCGGGTGA‑TTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTGCGGCGCTCTCCACCGGATAAACCGGCAGAGGGGAGGGCATGGATGTGTATGGGTGCAGG‑ATAAGGCGTTGGTCAAAGGCAAGATCGTCCGGTGTCATAAGCGACTTGTTAACCAAAtt                                                                                                                                                                                     <  2:105879/250‑1 (MQ=255)
                                                       cGACATGGCATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGA‑TTGTAGCCGTGGATCACCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGGGCCGG‑ATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGAATTGTAAACAAAATTg                                                                                                                                                                                    >  2:132801/1‑250 (MQ=255)
                                                         aggtggactcattgtgcttcttaaTCGCCGGCTTAAGCTGCGGGTGAGTTG‑AGGCGTGGATCGCCGCCCACCATGACGACATAACGGCAACCAGGCGTCTGCCGTCAGACATAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGTTCATGGATGTGTATGGGTGCAGG‑ATATGGCGTTGATCAAAGGCAAGATCGGCCGTTGGCATAATCGACTTGTAAACCAAATTGaaa                                                                                                                                                                                 <  2:248158/228‑1 (MQ=255)
                                                             ggCATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGA‑TTGTAGCTGTGGATCGCCGCCCACCAGGACGACATACCGGCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCCCCACCGGATAAACCGGCAGAGGGGAGGTCATGGCTGGGAATGAGTGCAGG‑ATATAGCGTTGGTCAAAGTCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCacataggaa                                                                                                                                                                                 >  2:178808/1‑239 (MQ=255)
                                                             ggCATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGA‑TTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGG‑ATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGAt                                                                                                                                                                             <  1:600402/251‑1 (MQ=255)
                                                             ggCATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGA‑TTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGG‑ATATGGCGTTGGTCAAAGGAAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAGATTGAAAAgtt                                                                                                                                                                             >  2:544331/1‑249 (MQ=255)
                                                              gCATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGA‑TTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGG‑ATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGAtt                                                                                                                                                                            <  1:178528/251‑1 (MQ=255)
                                                              gCATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGA‑TTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGG‑ATATGGCGTTGGGCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGAtt                                                                                                                                                                            >  1:134304/1‑251 (MQ=255)
                                                              gCATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGA‑GTGAAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGCCAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCCGAGGGGGGGTCATGGATGTGACTGGGTGCAGG‑ATATGGCGTTGGGCACAGGCAAGATCGCCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGAtt                                                                                                                                                                            >  2:98426/1‑251 (MQ=255)
                                                               caacAATTTGCGCCTTCATCGCCGCATTAAGCTGCGGGTGA‑TTGTAGCCGTGGATCGCCGCCCCCCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGG‑ATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGAttt                                                                                                                                                                           <  2:672085/248‑1 (MQ=255)
                                                                aTCAATTTGCGACTTCAGCTCCGCATTAAGCTGCGGGTGA‑TTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGG‑ATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTa                                                                                                                                                                          <  1:637534/251‑1 (MQ=255)
                                                                       tGCGACTTCATCGCCGCATTAAGCTGCGGGTGA‑TTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGG‑ATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTAc                                                                                                                                                                   >  1:581349/1‑251 (MQ=255)
                                                                         cGACTTCATCGCCGCATTAAGCTGCGGGTGA‑TTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGG‑ATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACaa                                                                                                                                                                 <  1:544331/251‑1 (MQ=255)
                                                                           aCTTCATCGCCGCATTAAGCTGCGGGTGA‑TTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGCGTGCAGG‑ATATGGCGTTGGTCAAAGGCACGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGGTTACAAgg                                                                                                                                                               >  2:668764/1‑250 (MQ=255)
                                                                             ttCATCGCCGCATTAAGCTGCGGGTGA‑TTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGG‑ATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGtat                                                                                                                                                             >  2:426190/1‑249 (MQ=255)
                                                                               cATCGCCTCATTAAGCTGCGGGTGA‑TTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCTCCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGG‑ATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGCCTTGTAAACCAAATTGAAAATATTTAGGTTTACAAGTCTa                                                                                                                                                            <  1:266455/250‑1 (MQ=255)
                                                                                 ccgccgcATTAAGCTGCGGGTGA‑TTGGAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGG‑ATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTacac                                                                                                                                                         <  2:196524/250‑1 (MQ=255)
                                                                                  cgccgcATTAAGCTGCGGGTGA‑TTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGGCATGGATGTGTATGGGTGCAGG‑CTATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGGCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTAAAcc                                                                                                                                                        >  2:390457/1‑251 (MQ=255)
                                                                                   gccgcATTAAGCTGCGGGTGA‑TTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGG‑ATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCg                                                                                                                                                       <  1:63668/251‑1 (MQ=255)
                                                                                     cgcATTAAGCTGCGGGTGA‑TTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGG‑ATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGaa                                                                                                                                                     <  1:290827/251‑1 (MQ=255)
                                                                                           aaGCTGCGGGTTA‑TTTTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGGTAACCAGGCGGCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGTGAGGTCATGGATGTGTATGGGTGCAGG‑ATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAAATGAAAAAATTTAGGTTTACAAGTCTACACCGAATTaaca                                                                                                                                               <  2:212272/251‑1 (MQ=255)
                                                                                                 cGGGTGA‑TTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGG‑ATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACaaaa                                                                                                                                         <  1:394631/251‑1 (MQ=255)
                                                                                                         ttgcaGCCGTGGAGCGCCTCCCACCAGGACGACATCCCGTACACCAGGCGTCTGACGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGG‑ATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATCGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACg                                                                                                                                   <  2:124030/246‑1 (MQ=255)
                                                                                                             agtcttggatctccgcctgccatcacgacATAACGTCAACCAGGCGTCCGCCGTCAGACAAAATCAGCTCGCACCCTTCGGCGCTCCTCACCGGATAAACCGGCAGAGTGGAGGTCAGGGATGTGTATGGGTGCATG‑ATATGGCGTTGGTCAAAGGCAAGTTCGTTCGTTTTCATAATCGACTTGTAAACCAAATTGAAAAGATTGAGTTTTACAAGTCTACACCGAATTAACAACCAAAAACACGTATTg                                                                                                                              <  2:507292/227‑1 (MQ=255)
                                                                                                             aGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGG‑ATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTg                                                                                                                              >  1:345558/1‑251 (MQ=255)
                                                                                                                       cgccgcCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGG‑ATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGcccc                                                                                                                    >  1:690214/1‑251 (MQ=255)
                                                                                                                             ccaccaGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGG‑ATATGGCGTTGGTCAAAGGCAAGATCGTCCGCTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGACCCATGGCt                                                                                                              >  2:114680/1‑251 (MQ=255)
                                                                                                                                     acgacATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGG‑ATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGccc                                                                                                      >  1:466375/1‑251 (MQ=255)
                                                                                                                                      cgacATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGG‑ATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGcccc                                                                                                     >  1:198513/1‑250 (MQ=255)
                                                                                                                                          aTACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGG‑ATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCACGCt                                                                                                 >  1:360177/1‑251 (MQ=255)
                                                                                                                                                 cAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGG‑ATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCACGCTGGACAt                                                                                           <  1:426190/250‑1 (MQ=255)
                                                                                                                                                  aaCCAGGCGGCTGACGTCAGACAAAATCAGCGCGCAACCGTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGG‑ATATGGCGTTGGTCAAAGGCAAGATAGTCCGTTGTCATAATCGACGTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACGCGTTTTGGAGAAGCCCCATGGCTCACCGCCCACGCTGGACATTg                                                                                         <  1:47025/251‑1 (MQ=255)
                                                                                                                                                       ggCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGTCAG‑ATATGGCGTTGGTCAAAGGCAAGAGCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCACGCTGGACATTGTCGCa                                                                                    <  1:36920/251‑1 (MQ=255)
                                                                                                                                                       ggCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCCAG‑ATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCACGCTGGACa                                                                                            <  1:282865/243‑1 (MQ=255)
                                                                                                                                                       ggCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCCAG‑ATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCACGCTGGACa                                                                                            >  2:282865/1‑243 (MQ=255)
                                                                                                                                                        gCGTCGGCCGGCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGG‑ATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCACGCTGGACATTGTCGCaa                                                                                   <  2:1421/251‑1 (MQ=255)
                                                                                                                                                         cggcTGCCGTCAGACAAAAGCAGCTCGCAACCTGCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGG‑ATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCACGCTGGACATTGTCGCAAg                                                                                  <  2:466375/248‑1 (MQ=255)
                                                                                                                                                              gccctcggACAAAACCAGGGCGCAACCTTCGGCGATCACCCCCGGAGAAACCGGCAGAGGGGAGGTCATGGATGTGTAGGTGTGCAGG‑ATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGGCTTGTAAACCAAATTGAAAAGATTTATTTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTGTTTGAGAAGCCCCATGGCTCACCGCCCACGCTGGACATTGTCGCAAGTca                                                                               <  1:683320/242‑1 (MQ=255)
                                                                                                                                                                      aCAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGG‑ATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCACGCTGGACATTGTCGCAAGTCACAGAattatt                                                                     >  1:262147/1‑251 (MQ=255)
                                                                                                                                                                                cTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGG‑ATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCACGCTGGACATTTACGCAAGTCACAGAATTATTTGAAAAACCg                                                           >  2:526792/1‑251 (MQ=255)
                                                                                                                                                                                 ccGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGG‑ATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCACGCTGGACATTGTCGCAAGTCACAGAATTATTTGAAAAACCGt                                                          <  1:184656/250‑1 (MQ=255)
                                                                                                                                                                                     aaCCTTCGGCGCTAACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGG‑ATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCTCCGCCCACGCTGGACATTGTCGCAAGTCACAGAAGTATTTGAAAAACCGTTGCt                                                      <  1:629987/251‑1 (MQ=255)
                                                                                                                                                                                       ccTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCCAG‑ATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCACGCTGGACATTGTCGCAAGTCACAGAATTATTTGAAAAACCGTTGCTgg                                                    >  1:588999/1‑251 (MQ=255)
                                                                                                                                                                                             gcgcTCACCACCGTATAAACCGGCAGAGGTGATGTCATGGATGTGTATGGGTGCAGG‑ATATGGCGTTGGTCAAAGTCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGAGTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCCCGCTGGACATTGTCGCAAGTCACAGAATTATTTGAAAAACCGTTGCTGGATctgc                                              <  1:438477/251‑1 (MQ=255)
                                                                                                                                                                                              cgcTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGG‑ATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCACGCTGGACATTGTCGCAAGTCACAGAATTATTTGAAAAACCGTTGCTGGATctgct                                             <  1:122457/251‑1 (MQ=255)
                                                                                                                                                                                               gcTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGG‑ATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCACGCTGGACATTGTCGCAAGTCACAGAATTATTTGAAAAACCGTTGCTGGATctgctg                                            <  1:439134/251‑1 (MQ=255)
                                                                                                                                                                                                  caccacCGGATAAACCGGCAGAGGGGAGGACATGGATGTGTATGGGTGCAGG‑ATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCACGCTGGACATTGTCGTAAGTCACAGAATTATTTGAAAAACCGTTGCTGGATCTGCTGtt                                          <  2:125387/250‑1 (MQ=255)
                                                                                                                                                                                                   accaccGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCCGG‑ATATGGCGTTAGACAGAGGCAAGATGGTCAGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTAAAAGTATACACCGCAGTACCAACCAAAAACACGTTTAGGAGAAGCACCATGGCTCACCGCCAACCCTGGACATTGGCGCACGTCACAGAACTATTTGAAAAACCGCTGCTGGATCTGCTGTTTg                                        >  2:602507/1‑251 (MQ=255)
                                                                                                                                                                                                           aTAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCAGG‑ATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCACGCTGGACATTGTCGCAAGTCACAGAATTATTTGAAAAACCGTTGCTGGATCTGCTGTTTGAAGCgcag                                <  1:278695/251‑1 (MQ=255)
                                                                                                                                                                                                              aaCCGGCAGAGGGGAGGTCATGGATGTGTATGGGAGCAGG‑ATAAGGCGTTGGTCAACGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGATTTGGAGAAGCCCCATGGATCACCGCCCAAGCTGGACATTGTCGCAAGTCACAGACTTATTTGAAAAACCGTTGCTGGATCTGCTGTTTGAAGCgcagccg                             >  2:590587/1‑249 (MQ=255)
                                                                                                                                                                                                                                cATGGATGTGTATGGGTGCAGG‑ATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCACGCTGGACATTGTCGCAAGTCACAGAATTATTTGAAAAACCGTTGCTGGATCTGCTGTTTGAAGCGCAGCAGGTGCATCGCCAGCATTTc           >  1:313977/1‑251 (MQ=255)
                                                                                                                                                                                                                                   ggATGTGTATTGGTGCAGG‑ATATGGCGTTGGGCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAAGTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCACGCTGGACATTGTCGCAAGTCACAGAATTATTTGAAAAACCGTTGCTGGATCTGCTGTTTGAAGCGCAGCAGGTGCATCGCCAGCATTTCGAt        <  2:227750/251‑1 (MQ=255)
                                                                                                                                                                                                                                     aTGTGTATGGGTGCAGG‑ATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGttt                                                                                                                                <  2:62436/129‑1 (MQ=255)
                                                                                                                                                                                                                                     aTGTGTATGGGTGCAGG‑ATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGttt                                                                                                                                >  1:62436/1‑129 (MQ=255)
                                                                                                                                                                                                                                         gtATGGGTGCCAG‑ATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCACGCTGGACATTGTCGCAAGTCACAGAATTATTTGAAAAACCGTTGCTGGATCTGCTGTTTGAAGCGCAGCAGGTGCATCGCCAGCATTTCTATCCTCg   <  2:588999/250‑1 (MQ=255)
                                                                                                                                                                                                                                          tATGGGTGC‑AGCATATGTCTTTGGTCAAAGGCAAGATCGTCCGTTTTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAATTCTACACCGAATTAACAACAAAAAACACGTTTTGGATAAGCCCCATGGCTCACCGCCCACGCTGGACATTGTCGCAAGTCACAGAATTATTTTAAAAACCTTTCCTGTATCTTCTGTTTGAAGCGCAGCAGTTGCATCGCCAGCATTTCAATCCTCGt  <  2:313977/250‑1 (MQ=255)
                                                                                                                                                                                                                                                   ||                                                                                                                                                                                                                                                
GTTTGCGGCACAGCTCAATGGCTGGCGCATGGGTGATACCGCCAAACATCACATGCGACATGGCATCAATTTGCGACTTCATCGCCGCATTAAGCTGCGGGTGA‑TTGTAGCCGTGGATCGCCGCCCACCAGGACGACATACCGTCAACCAGGCGTCTGCCGTCAGACAAAATCAGCTCGCAACCTTCGGCGCTCACCACCGGATAAACCGGCAGAGGGGAGGTCATGGATGTGTATGGGTGCCAG‑ATATGGCGTTGGTCAAAGGCAAGATCGTCCGTTGTCATAATCGACTTGTAAACCAAATTGAAAAGATTTAGGTTTACAAGTCTACACCGAATTAACAACAAAAAACACGTTTTGGAGAAGCCCCATGGCTCACCGCCCACGCTGGACATTGTCGCAAGTCACAGAATTATTTGAAAAACCGTTGCTGGATCTGCTGTTTGAAGCGCAGCAGGTGCATCGCCAGCATTTCGATCCTCGT  >  NC_000913/808975‑809457

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑