Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,201,754 A→T 5.9% V2022V (GTT→GTA PP_3691 ← DNA helicase‑related protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,201,7540AT 5.9% 363.4 / 10.5 118V2022V (GTT→GTAPP_3691DNA helicase‑related protein
Reads supporting (aligned to +/- strand):  ref base A (52/59);  new base T (5/2);  total (57/61)
Fisher's exact test for biased strand distribution p-value = 2.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.86e-01

Only 100 of 120 total aligned reads displayed.

TCATCTGCTCCAAAGTGGTTGTATAGCTCGCCTCAAAAAATCGGTCTGGATCAACTGCTTCGACAGCCTCGGTAGGGTCTGCTTCTTTGTAAACCAAGTGCTCCACCGGAACATAAGGAATGACCACTTCCTCAATGGCTTGAGCCGTAACGGCTCTGGCGTACTTCAACTCAGGCAATGCACCAGACGGATCAT‑‑TTGCTACGGCAATTGTGGGAGAAGCTACAACAGCAAGCTGATCCACTACCTCAATGACATCAGCCATCGCTGCATCGACCTTCAGGCGCTCCGCTTCATGCGC  >  NC_002947/4201606‑4201903
                                                                                                                                                      |                                                                                                                                                       
tcatcTGCTCCAAAGTGGTTGTATAGCTCGCCTCAAAAAATCGGTCTGGATCAACTGCTTCGACAGCCTCGGTAGGGTCTGCTTCTTTGTAAACCAAGTGCTCCACCGGAACATAAGGAATGACCACTTCCTCAATGGCTTGAGCCGTaa                                                                                                                                                        <  2:250644/150‑1 (MQ=255)
       ccccACAGTGGCTCTATAGCTCGCCTCAAAAAATCGGTCTGGATCAACGACCTCGACCATCCTCCGGT‑‑CATCTGGGGGGTTCTAAACCGAGTGCTCCACCGGAACATAATGAATGACCACTTCCTCAATGGCTTGAGCCGTAACGGctct                                                                                                                                                 <  1:543554/148‑1 (MQ=255)
         ccAAAGTGGTTGTATAGCTCGCCTCAAAAAATCGGTCTGGATCAACTGCTTCGACAGCCTCGGTAGGGTCTGCTTCTTTGTAAACCAAGTGCTCCACCGGAACATAAGGAATGACCACTTCCTCAATGGCTTGAGCCGTAACGGCTCTgg                                                                                                                                               >  1:1958379/1‑150 (MQ=255)
         ccAAAGTGGTTGTATAGCTCGCCTCAAAAAAGCGGTCTGGATCAACTGCTTCGACAGCCTCGGTAGGGTCTGCTTCTTTGTAAACCAAGTGCTCCACCGGAACATAAGGAATGACCACTTCCTCAATGGCTTGAGCCGTTACGGctcaag                                                                                                                                               <  2:138138/150‑4 (MQ=255)
          cAAAGTGGTTGTATAGCTCGCCTCAAAAAATCGGTCTGGATCAACTGCTTCGACAGCCTCGGTAGGGTCTGCTTCTTTGTAAACCAAGTGCTCCACCGGAACATAAGGAATGACCACTTCCTCAATGGCTTGAGCCGTAACGGCTCTGGc                                                                                                                                              >  1:1830377/1‑150 (MQ=255)
          cAAAGTGGTTGTATAGCTCGCCTCAAAAAATCGGTCTGGATCAACTGCTTCGACAGCCTCGGTAGGGTCTGCTTCTTTGTAAACCAAGTGCTCCACCGGAACATAAGGAATGACCACTTCCTCAATGGCTTGAGCCGTAACGGCTCTGGc                                                                                                                                              >  1:613389/1‑150 (MQ=255)
             aGTGGTTGTATAGCTCGCCTCAAAAAATCGGTCTGGATCAACTGCTTCGACAGCCTCGGTAGGGTCTGCTTCTTTGTAAACCAAGTGCTCCACCGGAACATAAGGAATGACCACTTCCTCAATGGCTTGAGCCGTAACGGCTCTGGCGTa                                                                                                                                           <  1:615806/150‑1 (MQ=255)
                 gTTGTATAGCTCGCCTCAAAAAATCGGTCTGGATCAACTGCTTCGACAGCCTCGGTAGGGTCTGCTTCTTTGTAAACCAAGTGCTCCACCGGAACATAAGGAATGACCACTTCCTCAATGGCTTGAGCCGTAACGGCTCTGGCGTACTTc                                                                                                                                       >  2:579607/1‑150 (MQ=255)
                   tGTATAGCTCGCCTAAAAAAATCGGTCTGGATCAACCGATTCGACAGCCTCGGTAGGGTCTGCTTCTTTGTAAACCAAGTGCTCCACCGGAACATAAGGAATGACCACTTCCTCAATGGCTTGAGCCGTAACGGCTCTGGCGTACTTCaa                                                                                                                                     <  1:726067/150‑1 (MQ=255)
                     tataGCTCGCCTCAAAAAATCGGTCTGGATCAACTGCTTCGACAGCCTCGGTAGGGTCTGCTTCTTTGTAAACCAAGTGCTCCACCGGAACATAAGGAATGACCACTTCCTCAATGGCTTGAGCCGTAACGGCTCTGGCGTACTTCAACt                                                                                                                                   >  1:2245238/1‑150 (MQ=255)
                       taGCTCGCCTCAAAAAATCGGTCTGGATCAACTGCTTCGACAGCCTCGGTAGGGTCTGCTTCTTTGTAAACCAAGTGCTCCACCGGAACATAAGGAATGACCACTTCCTCAATGGCTTGAGCCGTAACGGCTCTGGCGTACTTCAACTCa                                                                                                                                 <  2:2203030/150‑1 (MQ=255)
                        aGCTCGCCTCAAAAAATCGGTCTGGATCAACTGCTTCGACAGCCTCGGTAGGGTCTGCTTCTTTGTAAACCAAGTGCTCCACCGGAACATAAGGAATGACCACTTCCTCTATGGCTTGAGCCGTAACGGCTCTGGCGTACTTCAACTCAg                                                                                                                                <  1:1136710/150‑1 (MQ=255)
                             gCCTCAAAAAATCGGTCTGGATCAACTGCTTCGACAGCCTCGGTAGGGTCTGCTTCTTTGTAAACCAAGTGCTCCACCGGAACATAAGGAATGACCACTTCCTCAATGGCTTGAGCCGTAACGGCTCTGGCGTACTTCAACTCAGTCAAt                                                                                                                           >  2:1750357/1‑150 (MQ=255)
                               cTCAAAAAATCGGTCTGGATCAACTGCTTCGACAGCCTCGGTAGGGTCTGCTTCTTTGTAAACCAAGTGCTCCACCGGAACATAAGGAATGACCACTTCCTCAATGGCTTGAGCCGTAACGGCTCTGGCGTACTTCAACTCAGGCAATGc                                                                                                                         >  2:1576510/1‑150 (MQ=255)
                                tCAAAAAATCGGTCTGGATCAACTGCTTCGACAGCCTCGGTAGGGTCTGCTTCTTTGTAAACCAAGTGCTCCACCGGAACATAAGGAATGACCACTTCCTCAATGGCTTGAGCCGTAACGGCTCTg                                                                                                                                                >  2:754142/1‑126 (MQ=255)
                                tCAAAAAATCGGTCTGGATCAACTGCTTCGACAGCCTCGGTAGGGTCTGCTTCTTTGTAAACCAAGTGCTCCACCGGAACATAAGGAATGACCACTTCCTCAATGGCTTGAGCCGTAACGGCTCTGGCGTACTTCAACTCAGGCAATGCa                                                                                                                        >  1:852089/1‑150 (MQ=255)
                                tCAAAAAATCGGTCTGGATCAACTGCTTCGACAGCCTCGGTAGGGTCTGCTTCTTTGTAAACCAAGTGCTCCACCGGAACATAAGGAATGACCACTTCCTCAATGGCTTGAGCCGTAACGGCTCTGGCGTACTTCAACTCAGGCAATGCa                                                                                                                        >  1:2262375/1‑150 (MQ=255)
                                tCAAAAAATCGGTCTGGATCAACTGCTTCGACAGCCTCGGTAGGGTCTGCTTCTTTGTAAACCAAGTGCTCCACCGGAACATAAGGAATGACCACTTCCTCAATGGCTTGAGCCGTAACGGCTCTGGCGTACTTCAACTCAGGCAATGCa                                                                                                                        >  2:1754448/1‑150 (MQ=255)
                                 cAAAAACTCGGTCTGTAGCAACTGCTTCGACAGCCTCGGTAGGGTCTGCTTCTTTGTCAGACGGTTAATCCACCTGAACCTCCGGAATCACCACCTCCTCAATGGCTTGGGAGGTATCGTCTCTGGCGTACTTCAACTCAGTCAATGCAc                                                                                                                       <  1:1750357/150‑1 (MQ=255)
                                    aaaaTCGGTCTGGATCAACTGCTTCGACAGCCTCGGTAGGGTCTGCTTCTTTGTAAACCAAGTGCTCCACCGGAACATAAGGAATGACCACTTCCTCAATGGCTTGAGCCGTAACGGctc                                                                                                                                                  <  2:1560483/120‑1 (MQ=255)
                                    aaaaTCGGTCTGGATCAACTGCTTCGACAGCCTCGGTAGGGTCTGCTTCTTTGTAAACCAAGTGCTCCACCGGAACATAAGGAATGACCACTTCCTCAATGGCTTGAGCCGTAACGGctc                                                                                                                                                  >  1:1560483/1‑120 (MQ=255)
                                          ggTCTGGATCAACTGCTTCGACAGCCTCGGTAGGGTCTGCTTCTTTGTAAAACAAGTGCTCGACCGGAACATAAGGAATGACCACTTCCTCCATGGCTTGTGCCGTAACGGCTCTGGCGTACTTCAACTCAGGCAATGCACCATACgaga                                                                                                              >  2:459412/1‑147 (MQ=255)
                                          ccgatctaTCAACTGCTTCGACAGCCTCGGTAGGGTCTGCTTCTTTGTAAACCAAGTGCTCCACCGGAACATAAGGAATGACCACTTCCTCAATGGCTTGAGCCGTAACGGCTCTGGCGTACTTCAACTCAGGCAATGCACCAGACGGAt                                                                                                              <  2:1993797/143‑1 (MQ=255)
                                           gTCTGGATCAACTGCTTCGACAGCCTCGGTAGGGTCTGCTTCTTTGTAAACCAAGTGCTCCACCGGAACATAAGGAATGACCACTTCCTCAATGGCTTGAGCCGTAACGGCTCTGGCGTACTTCAACTCAGGCAATGCACCAGACGGATc                                                                                                             >  1:2136617/1‑150 (MQ=255)
                                            tCTGGATCAACTGCTTCGACAGCCTCGGTAGGGTCTGCTTCTTTGTAAACCAAGTGCTCCACCGGAACATAAGGAATGACCACTTCCTCAATGGCTTGAGCCGTAACGGCTCTGGCGTACTTCAACTCAGGCAATGCACCAGACGGATCa                                                                                                            >  2:1973556/1‑150 (MQ=255)
                                              tGGATCAACTGCTTCGACAGCCTCGGTAGGGTCTGCTTCTTTGTAAACCAAGTGCTCCACCGGAACATAAGGAATGACCACTTCCTCAATGGCTTGAGCCGTAACGGCTCTGGCGTACTTCAACTCAg                                                                                                                                <  1:951384/128‑1 (MQ=255)
                                               ggATCAACTGCTTCGACAGCCTCGGTAGGGTCTGCTTCTTTGTAAACCAAGTGCTCCACCGGAACATAAGGAATGACCACTTCCTCAATGGCTTGAGCCGTAACGGCTCTGGCGTACTTCAACTCAGGCAATGCACCAGACGGATCAt‑‑tt                                                                                                       <  1:294373/150‑1 (MQ=255)
                                               ggATCAACTGCTTCGACAGCCTCGGTAGGGTCTGCTTCTTTGTAAACCAAGTGCTCCACCGGAACATAAGGAATGACCACTTCCTCAATGGCTTGAGCCGTAACGGCTCTGGCGTACTTCAACTCAGGCAATGCACCAGACGGATCAt‑‑tt                                                                                                       <  2:1161200/150‑1 (MQ=255)
                                                gATCAACTGCTTCGACAGCCTCGGTAGGGTCTGCTTCTTTGTAAACCAAGTGCTCCACCGGAACATAAGGAATGACCACTTCCTCAATGGCTTGAGCCGTAACGGCTCTg                                                                                                                                                <  2:937105/110‑1 (MQ=255)
                                                gATCAACTGCTTCGACAGCCTCGGTAGGGTCTGCTTCTTTGTAAACCAAGTGCTCCACCGGAACATAAGGAATGACCACTTCCTCAATGGCTTGAGCCGTAACGGCTCTg                                                                                                                                                >  1:937105/1‑110 (MQ=255)
                                                   cAACTGCTTCGACAGCCTCGGTAGGGTCTGCTTCTTTGTAAACCAAGTGCTCCACCGGAACATAAGGAATGACCACTTCCTCAATGGCTTGAGCCGTAACGGCTCTGGCGTACTTCAACTCAGGCAATGCACCAGACGGATCAT‑‑TTGCTa                                                                                                   <  2:2245238/150‑1 (MQ=255)
                                                   aaaCTGCTTCGACATCCTCGGTAGGGTCTGCTTCTTTGTAAACCAAGTGCTCCACCGGAACATAAGGAATGACCACTTCCTCAATGGCTTGAGCCGTAACNGCCCNGGCGTACTTCAACTNAGGCAATGCACCAGACGGATCAT‑‑TTGCTa                                                                                                   <  1:459412/149‑1 (MQ=255)
                                                     aCTGCTTCGACAGCCTCGGTAGGGTCTGCTTCTTTGTAAACCAAGTGCTCCACCGGAACATAAGGAATGACCACTTCCTCAATGGCTTGAGCCGTAACGGCTCTGGCGTACTTCAAc                                                                                                                                    >  1:794194/1‑117 (MQ=255)
                                                     aCTGCTTCGACAGCCTCGGTAGGGTCTGCTTCTTTGTAAACCAAGTGCTCCACCGGAACATAAGGAATGACCACTTCCTCAATGGCTTGAGCCGTAACGGCTCTGGCGTACTTCAAc                                                                                                                                    <  2:794194/117‑1 (MQ=255)
                                                      cTGCTTCGACAGCCTCGGTAGGGTCTGCTTCTTTGTAAACCAAGTGCTCCACCGGAACATAAGGAATGACCACTTCCTCAATGGCTTGAGCCGTAACGGCTCTGGCGTACTTCAACTCAGGCAATGCACCAGACGGATCAT‑‑TTGCTACgg                                                                                                <  1:264911/150‑1 (MQ=255)
                                                       tGCTTCGACAGCCTCGGTAGGGTCTGCTTCTTTGTAAACCAAGTGCTCCACCGGAACATAAGGAATGACCACTTCCTCAATGGCTTGAGCCGTAACGGCTCTGGCGTACTTCAACTCAGGCAATGCACCAGACGGATCAT‑‑TTGCTACGGc                                                                                               <  2:1830377/150‑1 (MQ=255)
                                                       tGCTTCGACAGCCTCGGTAGGGTCTGCTTCTTTGTAAACCAAGTGCTCCACCGGAACATAAGGAATGACCACTTCCTCAATGGCTTGAGCCGTAACGGCTCTGGCGTACTTCAACTCAGGCAATGCACCAGACGGATCAT‑‑TTGCTACGGc                                                                                               <  2:613389/150‑1 (MQ=255)
                                                       tGCTTCGACAGCCTCGGTAGGGTCTGATTCTTTGTAAACCAAGTGCTCCACCGGAACATAAGGAATGACCACTTCCTCAATGGCTTGAGCCGTAACGGCTCTGGCGTACTTCAACTCAGGCAATGCACCAGACGGATCAT‑‑TTGCTACGGc                                                                                               >  1:790737/1‑150 (MQ=255)
                                                         cTTCGACAGCCTCGGTAGGGTCTGCTTCTTTGTAAACCAAGTGCTCCACCGGAACATAAGGAATGACCACTTCCTCAATGGCTTGAGCCGTAACGGCTCTGGCGTACTTCAACTCAGGCAATGCACCAGACGGATCAT‑‑TTGCTACGGCaa                                                                                             >  1:899951/1‑150 (MQ=255)
                                                         cTTCGACAGCCTCGGTAGGGTCTGCTTCTTCGTAAACCAAGTGCTCCACCGGAACATAAGGAATGACCACTTCCTCAATGGCTTGAGCCGTAACGGCTCTGGCGTACTTCAACTCAGGCAATGCACCAGACGGATCAT‑‑TTGCTACGGCaa                                                                                             >  2:293905/1‑150 (MQ=255)
                                                          ttCGACAGCCTCGGTAGGGTCTGCTTCTTTGTAAACCAAGTGCTCCACCGGAACATAAGGAATGACCACTTCCTCAATGGCTTGAGCCGTAACGGCTCTGGCGTACTTCAACTCAGGCAATGCACCAGACGGATCAT‑‑TTGCTACGGCAAt                                                                                            <  1:1973556/150‑1 (MQ=255)
                                                                cAGCCTCGGTAGGGTCTGCTTCTTTGTAAACCAAGTGCTCCACCGGAACATAAGGAATGACCACTTCCTCAATGGCTTGAGCCGTAACGGCTCTGGCGTACTTCAACTCAGGCAATGCACCAGACGGATCAT‑‑TTGCTACGGCAATTGTgg                                                                                       >  2:1157235/1‑150 (MQ=255)
                                                                cAGCCTCGGTAGGGTCTGCTTCTTTGTAAACCAAGTGCTCCACCGGAACATAAGGAATGACCACTTCCTCAATGGCTTGAGCCGTAACGGCTCTGGCGTACTTCAACTCAGGCAATGCACCAGACGGATCAT‑‑TTGCTACGGCAATTGTgg                                                                                       <  2:107054/150‑1 (MQ=255)
                                                                  gCCTCGGTAGGGTCTGCTTCTTTGTAAACCAAGTGCTCCACCGGAACATAAGGAATGACCACTTCCTCAATGGCTTGAGCCGTAACGGCTCTGGCGTACTTCAACTCAGGCAATGCACCAGACGGATCAT‑‑TTGCTACGGCAATTGTGGga                                                                                     <  2:899951/150‑1 (MQ=255)
                                                                   ccTCGGTAGGGTCTGCTTCTTTGTAAACCAAGTGCTCCACCGGAACATAAGGAATGACCACTTCCTCAATGGCTTGAGCCGTAACGGCTCTGGCGTACTTCAACTCAGGCAATGCACCAGACGGATCAT‑‑TTGCTACGGCAATTGTGGgag                                                                                    >  1:577883/1‑150 (MQ=255)
                                                                          tAGGGTCTGCTTCTTTGTAAACCAAGTGCTACACCGGAAAATAACGAATGACGACTTCATGAATGACTTATGCGCGAACTGCTGGGTGTTTCTTGAACTCAGGCAATGCACCAGACGGATCAT‑‑TTGCTACGGCAATTGTGGGAGAAGCTa                                                                              >  2:1465327/1‑150 (MQ=255)
                                                                             ggTCTGCTTCTTTGTAAACCAAGTGCTCCACCGGAACATAAGGAATGACCACTTCCTCAATGGCTTGAGCCGTAACGGCTCTGGCGTACTTCAACTCAGGCAATGCACCAGACGGATCAT‑‑TTGCTACGGCAATTGTGGGAGAAGCTacaa                                                                           <  2:1439469/150‑1 (MQ=255)
                                                                               tCTGCTTCTTTGTAACCCAAGTGCTCCACCGGAACATAAGGAATGACCACTTCCTCAATGGCTTGAGCCGTAACGGCTCTGGCGTACTTCAACTCAGGCAATGCACCAGACGGATCAT‑‑TTGCTACGGCAATTGTgg                                                                                       <  1:486519/136‑1 (MQ=255)
                                                                               tCTGCTTCTTTGTAAACCAAGTGCTCCACCGGAACATAAGGAATGACCACTTCCTCAATGGCTTGAGCCGTAACGGCTCTGGCGTACTTCAACTCAGGCAATGCACCAGACGGATCAT‑‑TTGCTACGGCAATTGTGGGAGAAGCTacaa                                                                           <  2:2283975/148‑1 (MQ=255)
                                                                               tCTGCTTCTTTGTAAACCAAGTGCTCCACCGGAACATAAGGAATGACCACTTCCTCAATGGCTTGAGCCGTAACGGCTCTGGCGTACTTCAACTCAGGCAATGCACCAGACGGATCAT‑‑TTGCTACGGCAATTGTGGGAGAAGCTacaa                                                                           >  1:2283975/1‑148 (MQ=255)
                                                                                cTGCTTGTTTGGAAACCAAGTGCTCCACCGGAACATAAGGAATGACCACTTCCTCAATGGCTTGAGCCGTAACGGCTCTGGCGTACTTCAACTCAGGCAATGCACCAGACGGATCAT‑‑TTGCTACGGCAATTGTGGGAGAAGCTACAACAg                                                                        >  2:1674788/1‑150 (MQ=255)
                                                                                    ttcttGGTATCCCTAGTGCTCTGTCGGAACATAAGGAATGACCACTCCCTCAATGGCTTGAGCAGCAACGGCTCTGGCGTACCTCAGCTCAGTCAATGCACCAGCCGGATCAT‑‑TTGCTACGGCAATTGTGGGAGAAGCTACAACAGCAAg                                                                    <  2:734979/150‑1 (MQ=255)
                                                                                     tcttTGTAAACCAAGTGCTCCACCGGAACATAAGGAATGACCACTTCCTCAATGGCTTGAGCCGTAACGGCTCTGGCGTACTTCAACTCAGGCAATGCACCAGACGGATCAT‑‑TTGCTACGGCAATTGTgg                                                                                       >  1:1633880/1‑130 (MQ=255)
                                                                                      cttTGTAAACCAAGTGCTCCACCGGAACATAAGGAATGACCACTTCCTCAATGGCTTGAGCCGTAACGGCTCTGGCGTACTTCAACTCAGGCAATGCACCAGACGGATCAT‑‑TTGCTACGGCAATTGTGGGAGAAGCTACAACAGCAAGCt                                                                  <  1:1157235/150‑1 (MQ=255)
                                                                                              aCCAAGTGCTCCACCGGAACATAAGGAATGACCACTTCCTCAATGGCTTGAGCCGTAACGGCTCTGGCGTACTTCAACTCAGGCAATGCACCAGACGGATCAT‑‑TTGCTACGGCAATTGTGGGAGAAGCTACAACAGCAAGCTGATCCACt                                                          >  2:888475/1‑150 (MQ=255)
                                                                                                cAAGTGCTCCACCGGAACATAAGGAATGACCACTTCCTCAATGGCTTGAGCCGTAACGGCTCTGGCGTACTTCAACTCAGGCAATGCACCAGACGGATCAT‑‑TTGCTACGGCAATTGTGGGAGAAGCTACAACAGCAAGCTGATCCACTAc                                                        <  2:2136617/150‑1 (MQ=255)
                                                                                                cAAGTGCTCCACCGGAACATAAGGAATGACCACTTCCTCAATGGCTTGAGCCGTAACGGCTCTGGCGTACTTCAACTCAGCCAATGCACCAGACGCCTGCTCATTGCT‑‑GACAATTTTGTGGGAAGCTACAACAGCAAGCTGATCCACTAc                                                        <  1:2254292/150‑1 (MQ=255)
                                                                                                 taGTGCTCCACCGGAACATAAGGAATGACCACTTCCTCAATGGCTTGAGCCGTAACGGCTCTGGCGTACTTCAACTCAGGCAATGCACCAGACGGATCAT‑‑TTGCTACGGCAATTGTGGGAGAAGCTACAACAGCAAGCTGATCCACTAcc                                                       <  2:425583/149‑1 (MQ=255)
                                                                                                 aaGTGCTCCACCGGAACATAAGGAATGACCACTTCCTCAATGGCTTGAGCCGTAACGGCTCTGGCGTACTTCAACTCAGGCAATGCACCAGAc                                                                                                                  >  2:611699/1‑93 (MQ=255)
                                                                                                 aaGTGCTCCACCGGAACATAAGGAATGACCACTTCCTCAATGGCTTGAGCCGTAACGGCTCTGGCGTACTTCAACTCAGGCAATGCACCAGAc                                                                                                                  <  1:611699/93‑1 (MQ=255)
                                                                                                           ccGGAACATAAGGAATGACCACTTCCTCAATGGCTTGAGCCGTAACGGCTCTGGCGTACTTCAACTCAGGCAATGCACCAGACGGATCAT‑‑TTGCTACGGCAATTGTGGGAGAAGCTACAACAGCAAGCTGATCCACTAcc                                                       >  2:344766/1‑140 (MQ=255)
                                                                                                           ccGGAACATAAGGAATGACCACTTCCTCAATGGCTTGAGCCGTAACGGCTCTGGCGTACTTCAACTCAGGCAATGCACCAGACGGATCAT‑‑TTGCTACGGCAATTGTGGGAGAAGCTACAACAGCAAGCTGATCCACTAcc                                                       <  1:344766/140‑1 (MQ=255)
                                                                                                            cGGAACATAAGGAATGACCACTTCCTCAATGGCTTGAGCCGTAACGGCTCTGGCGTACTTCAACTCAGGCAATGc                                                                                                                         >  2:1905882/1‑75 (MQ=255)
                                                                                                            cGGAACATAAGGAATGACCACTTCCTCAATGGCTTGAGCCGTAACGGCTCTGGCGTACTTCAACTCAGGCAATGc                                                                                                                         <  1:1905882/75‑1 (MQ=255)
                                                                                                             ggAACATAAGGAATGACCACTTCCTCAATGGCTTGAGCCGTAACGGCTCTGGCGTACTTCAACTCAGGCAATGCACCAGACGGATCAt‑‑tt                                                                                                       >  2:2017271/1‑90 (MQ=255)
                                                                                                             ggAACATAAGGAATGACCACTTCCTCAATGGCTTGAGCCGTAACGGCTCTGGCGTACTTCAACTCAGGCAATGCACCAGACGGATCAt‑‑tt                                                                                                       <  1:2017271/90‑1 (MQ=255)
                                                                                                                   tAAGGAATGACCACTTCCTCAATGGCTTGAGCCGTAACGGCTCTGGCGTACTTCAACTCAGGCAATGCACCAGAc                                                                                                                  <  1:344908/75‑1 (MQ=255)
                                                                                                                   tAAGGAATGACCACTTCCTCAATGGCTTGAGCCGTAACGGCTCTGGCGTACTTCAACTCAGGCAATGCACCAGAc                                                                                                                  >  2:344908/1‑75 (MQ=255)
                                                                                                                   tAAGGAATGACCACTTCCTCAATGGCTTGAGCCGTAACGGCTCTGGCGTACTTCAACTCAGGCAATGCACCAGACGGATCAT‑‑TTGCTACGGCAATTGTGGGAGAAGCTACAACAGCAAGCTGATCCACTACCTCAATGACATCAGCCATc                                     >  1:900301/1‑150 (MQ=255)
                                                                                                                    aaGGAATGACCACTTCCTCAATGGCTTGAGCCGTAACGGCTCTGGCGTACTTCAACTCAgg                                                                                                                               <  1:1194442/61‑1 (MQ=255)
                                                                                                                    aaGGAATGACCACTTCCTCAATGCCTTGAGCCGTAACGGCTCTGGCGTACTTCAACTCAgg                                                                                                                               >  2:1194442/1‑61 (MQ=255)
                                                                                                                      ggAATGACCACTTCCTCAATGGCTTGAGCCGTAACGGCTCTGGCGTACTTCAACTCAGGCAATGCACCAGACGGATCAT‑‑TTGCTACGGCAATTGTGGGAGAAGCTACAACAGCAAg                                                                    >  1:1280563/1‑116 (MQ=255)
                                                                                                                      ggAATGACCACTTCCTCAATGGCTTGAGCCGTAACGGCTCTGGCGTACTTCAACTCAGGCAATGCACCAGACGGATCAT‑‑TTGCTACGGCAATTGTGGGAGAAGCTACAACAGCAAg                                                                    <  2:1280563/116‑1 (MQ=255)
                                                                                                                       gAATGACCACTTCCTCAATGGCTTGAGCCGTAACGGCTCTGGCGTACTTCAACTCAGGCAATGCACCAGACGGATCAT‑‑TTGCTACGGCAATTGTGGGAGAAg                                                                                 <  1:132477/102‑1 (MQ=255)
                                                                                                                         aTGACCACTTCCTCAATGGCTTGAGCCGTAACGGCTCTGGCGTACTTCAACTCAGGCAATGCACCAGACGGAAAAT‑‑TTGCTACGGCAATTGTGGGAGAAGCTACAACAGCAAGCTGATCCACTACCTCAATGACATCAGCCATCGCTGCa                               >  1:90144/1‑150 (MQ=255)
                                                                                                                            aCCACTTCCTCAATGGCTTGAGCCGTAACGGCTCTGGCGTACTTCAACTCAGGCAATGCACCAGACGGATCAT‑‑TTGCTACGGCAATTGTGGGAGAAGCTACAACAGCAAGCTGATCAACTACCTCAATGACATCAGCCATCGCTGCATCg                            >  1:842776/1‑150 (MQ=255)
                                                                                                                             ccACTTCCTCAATGGCTTGAGCCGTAACGGCTCTGGCGTACTTCAACTCAGGCAATGCACCAGACGGATCAT‑‑TTGCTACGGCAATTGTGGGAGAAGCTACAACAGCAAGCTGATCCACTACCTCAATGACATCa                                           >  1:806041/1‑134 (MQ=255)
                                                                                                                             ccACTTCCTCAATGGCTTGAGCCGTAACGGCTCTGGCGTACTTCAACTCAGGCAATGCACCAGACGGATCAT‑‑TTGCTACGGCAATTGTGGGAGAAGCTACAACAGCAAGCTGATCCACTACCTCAATGACATCa                                           <  2:806041/134‑1 (MQ=255)
                                                                                                                              cACTTCCTCAATGGCTTGAGCCGTAACGGCTCTGGCGTACTTCAACTCAGGCAATGCACCAGACGGATCAG‑‑GTGCTACGGCCAGTGCGCGAAAAACTACAACAGCAAGCTGATCCACTACCTCAATGACATCAGCCAGCGCTGCATCGAc                          >  1:168858/1‑150 (MQ=255)
                                                                                                                                   ccTCAATGGCTTGAGCCGTAACGGCTCTGGCGTACTTCAACTCAGGCAATGCACCAGACGGATCAT‑‑TTGCTACGGCAATTGTGGGAGAAGCTACAACAGCAAGCTGATCCACTACCTCAATGACATCAGCCATCGCTGCATCGACCTTCa                     <  2:577883/150‑1 (MQ=255)
                                                                                                                                    cacAATGGCTTGAGCCGTAACGGCTCCGGCGTACGTCAACTCAGGCAATGCACCAGACGGATCAT‑‑TTGCTACGGCCATTGTGGGAGAAGCTACAACAGCAAGCTGATCCACTACCTCAATGACATCAGCCATCGCTGCATCGACCTGCAg                    >  2:91404/3‑150 (MQ=255)
                                                                                                                                    cTCAATGGCTTGAGCCGTAACGGCTCTGGCGTACTTCAACTCAGGCAATGCACCAGACGGATCAT‑‑TTGCTACGGCAATTGTGGGAGAAGCTACAACAGCAAGCTGATCCACTACCTCAATGACATCAGCCATCGCTGCATCGACCTTCAg                    <  1:582203/150‑1 (MQ=255)
                                                                                                                                    cTCAATGGCTTGAGCCGTAACGGCTCTGGCGTACTTCAACTCAGGCAATGCACCAGACGGATCAT‑‑TTGCTACGGCAATTGTGGGAGAAGCTACAACAGCAAGCTGATCCACTACCTCAATGACATCAGCCATCGCTGCATCGACCTTCAg                    >  2:1900173/1‑150 (MQ=255)
                                                                                                                                     tCAATGGCTTGAGCCGTAACGGCTCTGGCGTACTTCAACTCAGGCAATGCACCAGACGGATCAT‑‑TTGCTACg                                                                                                 <  2:2326652/72‑1 (MQ=255)
                                                                                                                                     tCAATGGCTTGAGCCGTAACGGCTCTGGCGTACTTCAACTCAGGCAATGCACCAGACGGATCAT‑‑TTGCTACg                                                                                                 >  1:2326652/1‑72 (MQ=255)
                                                                                                                                     tCAATGGCTTGAGCCGTAACGGCTCTGGCGTACTTCAACTCAGGCAATGCACCAGACGGATCAT‑‑TTGCTACGGCAATTGTGGGAGAAGCTACAACAGc                                                                       <  1:204979/98‑1 (MQ=255)
                                                                                                                                     tCAATGGCTTGAGCCGTAACGGCTCTGGCGTACTTCAACTCAGGCAATGCACCAGACGGATCAT‑‑TTGCTACGGCAATTGTGGGAGAAGCTACAACAGc                                                                       >  2:204979/1‑98 (MQ=255)
                                                                                                                                     tCAATGGCTTGAGCCGTAACGGCTCTGGCGTACTTCAACTCAGGCAATGCACCAGACGGATCAT‑‑TTGCTACGGCAATTGTGGGAGAAGCTACAACAGc                                                                       >  2:394726/1‑98 (MQ=255)
                                                                                                                                          ggCTTGAGCCGTAACGGCTCTGGCGTACTTCAACTCAGGCAATGCACCAGACGGATCAT‑‑TTGCTACGGCAATTGTGGGAGAAGCTACAACAGCAAGCTGATCCACTACCTCAATGACATCAGCCATCGCTGCATCGACCTTCAGGCGCTc              <  2:276367/150‑1 (MQ=255)
                                                                                                                                            ccagAGCCGTTACGGCTCTGGCGTACTTCAACTCAGGCAATGCACCAGACGGATCAT‑‑TTGCTACGGCAATTGTGGGAGAAGCTACAACAGCAAGCTGATCCACTACCTCAATGACATCAGCCATCGCTGCATCGACCTTCAGGCGCTCCg            >  2:1734003/4‑150 (MQ=255)
                                                                                                                                            ccagAGCCGTTACGGCTCTGGCGTACTTCAACTCAGGCAATGCACCAGACGGATCAT‑‑TTGCTACGGCAATTGTGGGAGAAGCTACAACAGCAAGCTGATCCACTACCTCAATGACATCAGCCATCGCTGCATCGACCTTCAGGCGCTCCg            >  1:182140/4‑150 (MQ=255)
                                                                                                                                            cTTGAGCCGTAACGGCTCTGGCGTACTTCAACTCAGGCAATGCACCAGACGGATCAT‑‑TTGCTACGGCAATTGTTGTAGAAGCTACAACAGCAAGCTGATCCACTACCTCAATGACATCAGCCATCGCTGCATCGACCTTCAGGCGCTCCg            <  2:790737/150‑1 (MQ=255)
                                                                                                                                             ttGAGCCGTAACGGCTCTGGCGTACTTCAACTCAGGCAATGCACCAGACGGATCAT‑‑TTGCTACGGCAATTGTGGGAGAAGCTACAACAGCAAGCTGATCCACTACCTCAATGACATCAGCCATCGCTGCATCGa                           <  2:2046115/134‑1 (MQ=255)
                                                                                                                                             ttGAGCCGTAACGGCTCTGGCGTACTTCAACTCAGGCAATGCACCAGACGGATCAT‑‑TTGCTACGGCAATTGTGGGAGAAGCTACAACAGCAAGCTGATCCACTACCTCAATGACATCAGCCATCGCTGCATCGa                           >  1:2046115/1‑134 (MQ=255)
                                                                                                                                             cagAGCCGTTACGGCTCTGGCGTACTTCAACTCAGGCAATGCACCAGACGGATCAT‑‑TTGCTACGGCAATTGTGGGAGAAGCTACAACAGCAAGCTGATCCACTACCTCAATGACATCAGCCATCGCTGCATCGACCTTCAGGCGCTCCGc           >  2:505186/3‑150 (MQ=255)
                                                                                                                                             cagAGCCGTTACGGCTCTGGCGTACTTCAACTCAGGCAATGCACCAGACGGATCAT‑‑TTGCTACGGCAATTGTGGGAGAAGCTACAACAGCAAGCTGATCCACTACCTCAATGACATCAGCCATCGCTGCATCGACCTTCAGGCGCTCCGc           >  2:903899/3‑150 (MQ=255)
                                                                                                                                              agAGCCGTTACGGCTCTGGCGTACTTCAACTCAGGCAATGCACCAGACGGATCAT‑‑TTGCTACGGCAATTGTgg                                                                                       >  1:319872/2‑73 (MQ=255)
                                                                                                                                               gAGCCGTAACGGCTCTGGCGTACTTCAACTCAGGCAATGCACCAGACGGATCAT‑‑TTGCTACGGCAATTGTGGGAGAAGCTACAACAGCAAGCTGATCCACTACCTCAATGACATCAGCCATCGCTGCATCGACCTTCAGGCGCTCCGCtt         <  1:1754448/150‑1 (MQ=255)
                                                                                                                                                 gCCGTAACGGCTCTGGCGTACTTCAACTCAGGCAATGCACCAGACGGATCAT‑‑TTGCTACGGCAATTGTGGGAGAAGCTACAACAGCAAGCTGATCCACTACCTCAATGACATCAGCCATCGCTGCATCGACCTTCAGGCGCTCCGCTTCa       >  1:2182550/1‑150 (MQ=255)
                                                                                                                                                  ccGTAACGGCTCTGGCGTACTTCAACTCAGGCAATGCACCAGACGGATCAT‑‑TTGCTACGGCAATTGTGGGAGAAGCTACAACAGCAAGCTGATCCACTACCTCAATGACATCAGCCATCGCTGCATCGACCTTCAGGCGCTCCGCTTCAt      <  1:395469/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                       
TCATCTGCTCCAAAGTGGTTGTATAGCTCGCCTCAAAAAATCGGTCTGGATCAACTGCTTCGACAGCCTCGGTAGGGTCTGCTTCTTTGTAAACCAAGTGCTCCACCGGAACATAAGGAATGACCACTTCCTCAATGGCTTGAGCCGTAACGGCTCTGGCGTACTTCAACTCAGGCAATGCACCAGACGGATCAT‑‑TTGCTACGGCAATTGTGGGAGAAGCTACAACAGCAAGCTGATCCACTACCTCAATGACATCAGCCATCGCTGCATCGACCTTCAGGCGCTCCGCTTCATGCGC  >  NC_002947/4201606‑4201903

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: