Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,250,992 G→C 8.8% intergenic (+38/‑32) yqjF → / → yqjG putative quinol oxidase subunit/putative S‑transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,250,9920GC 8.8% 296.1 / 15.8 91intergenic (+38/‑32)yqjF/yqjGputative quinol oxidase subunit/putative S‑transferase
Reads supporting (aligned to +/- strand):  ref base G (50/33);  new base C (3/5);  total (53/38)
Fisher's exact test for biased strand distribution p-value = 2.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

AACTCGCTGATGTTCATGAAAAACCTGACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTCGTAACTGGCTCACTGCCGATGGCGC  >  NC_000913/3250844‑3251139
                                                                                                                                                    |                                                                                                                                                   
aaCTCGCTGATGTTCATGAAACACCTGACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAg                                                                                                                                                     <  2:365022/149‑1 (MQ=255)
   tCGCTGATGTTCATGAAAAACCTGACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTAAAGCATCGACCGCCGGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGcta                                                                                                                                                  >  2:1340737/1‑148 (MQ=255)
       tGATGTTCATGAAAAACCTGACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCg                                                                                                                                              >  2:606892/1‑149 (MQ=255)
        gATGTTCATGAAAAACCTGACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGCTTCCTCCTCGt                                                                                                                                             <  2:1261381/149‑1 (MQ=255)
        gATGTTCATGAAAAACCTGACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACAGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGCTTCCTCCTCGt                                                                                                                                             <  1:29091/149‑1 (MQ=255)
        ccttgtCATGAAAAACCAGACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGt                                                                                                                                             <  1:109498/144‑1 (MQ=255)
         aTGTTCATGAAAAACCTGACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGCTTCCTCCTCGtt                                                                                                                                            <  1:1009214/149‑1 (MQ=255)
         aTGTTCATGAAAAACCTGACAATTTATGGCTGATTCCTGTTGCTGGCAATTACCGGTCCGGGCGCGTATTGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCTCGCTCTATACTTAACTTTAAATAACGAGGAGCTTCCTCCTCGtt                                                                                                                                            <  2:1073660/149‑1 (MQ=255)
           tttCATTAAAAACCTTACAATTTCTGGCGGATTCATGCTTGTGGCAATTAGCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGCTTCCTCCTCGtttt                                                                                                                                          <  2:668308/148‑1 (MQ=255)
            ttCATGAAAAACCTGACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGttttt                                                                                                                                         >  2:1348847/1‑149 (MQ=255)
            ttCATGAAAAACCTGACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGttttt                                                                                                                                         >  1:29693/1‑149 (MQ=255)
              cATGAAAAACCTGACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTACGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGATCCTCGTTTTTGc                                                                                                                                       >  1:28612/1‑149 (MQ=255)
              aaTGAAAAACCTGACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGc                                                                                                                                       <  1:160835/148‑1 (MQ=255)
                 gAAAAACCTGACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCTGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGAACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTGTGCTAt                                                                                                                                    >  1:121082/1‑149 (MQ=255)
                 gAAAAACCTGACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTAt                                                                                                                                    <  2:354454/149‑1 (MQ=255)
                 gAAAAACCTGACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTAt                                                                                                                                    <  2:953107/149‑1 (MQ=255)
                 gAAAAACCTGACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTAt                                                                                                                                    <  2:348688/149‑1 (MQ=255)
                  caaaaCCTGACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGCCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGTTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTAtt                                                                                                                                   <  1:1252494/148‑1 (MQ=255)
                  aaaaaCCTGACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTAtt                                                                                                                                   <  1:952913/149‑1 (MQ=255)
                   aaaaCCTGACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTg                                                                                                                                  <  1:1340986/149‑1 (MQ=255)
                         tGACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGgaggag                                                                                                                            >  1:1523397/1‑149 (MQ=255)
                          gACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCCTCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGgaga                                                                                                                           >  1:95822/1‑149 (MQ=255)
                          gACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGgaga                                                                                                                           >  2:1383546/1‑149 (MQ=255)
                                  cTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATg                                                                                                                   <  2:143829/149‑1 (MQ=255)
                                    ggCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCATGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATggg                                                                                                                 <  1:1389204/149‑1 (MQ=255)
                                     gCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGt                                                                                                                <  2:708263/149‑1 (MQ=255)
                                         aTTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAAc                                                                                                            <  2:1536708/149‑1 (MQ=255)
                                         aTTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAAc                                                                                                            <  1:658769/149‑1 (MQ=255)
                                         aTTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAAc                                                                                                            <  1:618343/149‑1 (MQ=255)
                                         aTTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAAc                                                                                                            <  2:665190/149‑1 (MQ=255)
                                               gctgctgGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCGTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTg                                                                                                      >  1:1432503/1‑149 (MQ=255)
                                               gctgctgGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAAAGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTg                                                                                                      >  2:1125136/1‑149 (MQ=255)
                                                 tgctgGCAATTACAGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACCCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGAc                                                                                                    >  2:1101245/1‑149 (MQ=255)
                                                    tgGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAATCGAGGAGGAAGATCCGCGGTTTTGCTATGGGTGGAGAGTAAAATGGGTCAACTGATTGACGGc                                                                                                 >  1:879266/1‑149 (MQ=255)
                                                         aTTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTg                                                                                            >  1:381957/1‑149 (MQ=255)
                                                         aTTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGGGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTg                                                                                            <  1:723073/149‑1 (MQ=255)
                                                          ttACCGGTCCGGGCGCATATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTAGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAAATGATTGACGGCGTCTgg                                                                                           >  2:579193/1‑149 (MQ=255)
                                                           tACCGGTCCGGGCGCTTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGc                                                                                          <  2:1430673/149‑1 (MQ=255)
                                                           tACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGc                                                                                          <  1:606892/149‑1 (MQ=255)
                                                              cGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATg                                                                                       >  2:1230766/1‑149 (MQ=255)
                                                              cGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATg                                                                                       >  1:174023/1‑149 (MQ=255)
                                                                      gcgcGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTgg                                                                               <  1:342586/149‑1 (MQ=255)
                                                                           tataGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGa                                                                          >  2:1102730/1‑149 (MQ=255)
                                                                              aGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATAc                                                                       <  2:1171282/149‑1 (MQ=255)
                                                                                     aCCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCt                                                                >  2:1167814/1‑149 (MQ=255)
                                                                                     aCCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCt                                                                >  1:599537/1‑149 (MQ=255)
                                                                                       cGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTAc                                                              <  2:1338816/149‑1 (MQ=255)
                                                                                       cGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTAc                                                              >  2:1403478/1‑149 (MQ=255)
                                                                                         ccTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACcg                                                            >  1:829923/1‑149 (MQ=255)
                                                                                         ccTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACcg                                                            >  1:691341/1‑149 (MQ=255)
                                                                                           tgctgAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAATCGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGAt                                                                         <  2:559603/134‑1 (MQ=255)
                                                                                           tgctgAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAATCGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGAt                                                                         >  1:559603/1‑134 (MQ=255)
                                                                                             ctgAATATAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTATTGCTATTGGAGGAGAGAAAAATGGGTCAAATGATTGACGGCGTCTGGCATGACACCTGGTACCATACAAAATCTACcggcgg                                                        >  2:803591/1‑149 (MQ=255)
                                                                                             ctgAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTTTTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACaggcgt                                                        >  2:306941/1‑148 (MQ=255)
                                                                                             ctgAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGGCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACcggcgg                                                        <  1:1348847/149‑1 (MQ=255)
                                                                                              tgAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGt                                                       >  1:1401630/1‑149 (MQ=255)
                                                                                              tgAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGt                                                       >  1:203687/1‑149 (MQ=255)
                                                                                               gAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAAt                                                                  <  2:1534693/137‑1 (MQ=255)
                                                                                               gAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAAt                                                                  >  1:1534693/1‑137 (MQ=255)
                                                                                               gAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTa                                                      >  2:818173/1‑149 (MQ=255)
                                                                                                aaTAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTaa                                                     >  2:1413120/1‑149 (MQ=255)
                                                                                                  tAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAAt                                                   >  2:345224/1‑149 (MQ=255)
                                                                                                     aaaGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTc                                                >  2:777207/1‑149 (MQ=255)
                                                                                                         tGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACg                                            >  2:810942/1‑149 (MQ=255)
                                                                                                           gTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGtt                                          >  2:1427471/1‑149 (MQ=255)
                                                                                                            tAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTc                                         <  2:197857/149‑1 (MQ=255)
                                                                                                             aaGCACGCGCTATACTTAACGATAAAAAACGAGGATGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGAAACCAAATCTACAGACGGTAAATTTCACCGTTCa                                        >  2:1166418/1‑149 (MQ=255)
                                                                                                               gCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGa                                                                                      <  2:132608/101‑1 (MQ=255)
                                                                                                               gCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGa                                                                                      >  1:132608/1‑101 (MQ=255)
                                                                                                                cACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGCCACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCt                                     >  1:313118/1‑149 (MQ=255)
                                                                                                                cACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCt                                     >  2:694943/1‑149 (MQ=255)
                                                                                                                  cgcgCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTc                                   >  2:234342/1‑149 (MQ=255)
                                                                                                                   gcgcTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTcc                                  <  2:95822/149‑1 (MQ=255)
                                                                                                                     gcTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGc                                >  2:762480/1‑149 (MQ=255)
                                                                                                                      cTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCa                               >  1:130523/1‑149 (MQ=255)
                                                                                                                         taCTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCAttt                            <  1:1230766/149‑1 (MQ=255)
                                                                                                                               aCGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGTTACGATACCAAATCTACCTGCGGTAAATTTCAACGTTCAGCTTCCGCATTTCGTAAc                      >  2:551147/1‑149 (MQ=255)
                                                                                                                               aCGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTCGTAAc                      <  2:29693/149‑1 (MQ=255)
                                                                                                                                atagcaaaaaCGAGGAGCTTCCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTCGTAACt                     >  1:612360/6‑149 (MQ=255)
                                                                                                                                 gATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTCGTAACTg                    >  2:1349068/1‑149 (MQ=255)
                                                                                                                                 gATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTAAACGTTCAGCTTCCGCATTTCGTTACTg                    >  1:1591368/1‑149 (MQ=255)
                                                                                                                                  agcaaaaaCGAGGAGCTTCCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTCGTAACTgg                   >  1:1161579/4‑149 (MQ=255)
                                                                                                                                  aTAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCACCGTTCAGCTTCCGCATTTCGTAACtcg                   >  2:624356/1‑147 (MQ=255)
                                                                                                                                  aTAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTCGTAACTgg                   <  1:1101245/149‑1 (MQ=255)
                                                                                                                                       aaaCGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTCGTAACTGGCTCAc              >  1:1209650/1‑149 (MQ=255)
                                                                                                                                        aaCGAGGAGCTTCCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTCGTAACTGGCTCACt             >  2:650652/1‑149 (MQ=255)
                                                                                                                                         aCGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTCGTAACTGGCTCACTg            >  2:1224783/1‑149 (MQ=255)
                                                                                                                                          cGAGGAGGAAGCTCGTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATTCCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTCGTAACTGGCTCACTGc           <  2:121082/149‑1 (MQ=255)
                                                                                                                                           gaggagGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATTGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTCGTAACTGGCTCACTGcc          <  2:28612/149‑1 (MQ=255)
                                                                                                                                              gaggaAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTTAATTTCAACGTTCAGCATCCGCATTTCGTAACTGGCTCACTGCCGAt       >  1:1533368/1‑149 (MQ=255)
                                                                                                                                                 gaAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTCGTAACTGGCTCACTGCCGATGgc    >  2:1469593/1‑149 (MQ=255)
                                                                                                                                                   aGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGATCAGCTTCCGCATTTCGTAACTGGCTCACTGCCGATGgcgc  <  1:818173/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
AACTCGCTGATGTTCATGAAAAACCTGACAATTTCTGGCGGATTCCTGCTGCTGGCAATTACCGGTCCGGGCGCGTATAGCATCGACCGCCTGCTGAATAAAAAGTGGTAAGCACGCGCTATACTTAACGATAAAAAACGAGGAGGAAGCTCCTCGTTTTTGCTATTGGAGGAGAGAAAAATGGGTCAACTGATTGACGGCGTCTGGCATGACACCTGGTACGATACCAAATCTACCGGCGGTAAATTTCAACGTTCAGCTTCCGCATTTCGTAACTGGCTCACTGCCGATGGCGC  >  NC_000913/3250844‑3251139

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