Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,296,545 G→A 9.9% Q125* (CAG→TAG) ‡ PP_2898 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,296,5450GA 9.9% 238.0 / 15.0 91Q125* (CAG→TAG) ‡PP_2898hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (42/40);  new base A (4/5);  total (46/45)
Fisher's exact test for biased strand distribution p-value = 7.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.39e-01

CGATGATGTGCAGTGGGTACAGCACCCCTCAATCGCACGGGTAGCTGCTCAGCAAACTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCGGCCTGCGCCAGCGTTTGTCGGGCCAGGTCGGGCAAGCGGGCTGTGGGCCCGCTGAAAT  >  NC_002947/3296398‑3296685
                                                                                                                                                   |                                                                                                                                            
cgATGATGTGCAGTGGGTACAGCACCCCTCAATCGCACGGGTAGCTGCTCAGCAAACTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGc                                                                                                                                            <  1:953247/150‑1 (MQ=255)
   tgatgTGCAGTGGGTACAGCACCCCTCAGTCGCACGGGTAGCTGATCAGGAATCTGGCGCTACCCAAATAGGCTTGAGTGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCTGt                                                                                                                                         <  2:662344/150‑1 (MQ=255)
   tgatgTGCAGTGGGTACAGCACCCCTCAATCGCACGGGTAGCTGCTCAGCAAACTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGt                                                                                                                                         >  1:1088715/1‑150 (MQ=255)
    gatgTGCAGTGGGTACAGCACCCCTCAATCGCACGGGTAGCTGCTCAGCAAACTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGAAAAGTCAGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGtt                                                                                                                                        <  2:2201080/150‑1 (MQ=255)
     atgTGCAGTGGGTACAGCACCCCTCAATCGCACGGGTAGCTGCTCAGCAAACTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTg                                                                                                                                       >  1:1395561/1‑150 (MQ=255)
      tgtgCAGTGGGTACAGCACCCCTCAATCGCACGGGTAGCTGCTCAGCAAACTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTgg                                                                                                                                      >  1:228529/1‑149 (MQ=255)
      tgtgCAGTGGGTACAGCACCCCTCAATCGCACGGGTAGCTGCTCAGCAAACTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTgc                                                                                                                                      <  2:103618/150‑1 (MQ=255)
       gtgCAGTGGGTACAGCACCCCTCAATCGCACGGGTAGCTGCTCAGCAAACTGGCGCTACCCAGATAGGCTTGAGGGGTGGGGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTgcg                                                                                                                                     >  1:1312386/1‑150 (MQ=255)
       gtgCAGTGGGTACAGCACCCCTCAATCGCACGGGTAGCTGCTCAGCAAACTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCCAGCCGTTgca                                                                                                                                     <  1:703069/150‑2 (MQ=255)
         gCAGTGGGTACAGCACCCCTCAATCGCACGGGTAGCTGCTCAGCAAACTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCt                                                                                                                                   <  1:2304017/150‑1 (MQ=255)
           aGTGGGTACAGCACCCCTCAATCGCACGGGTAGCTGCTCAGCAAACTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTgg                                                                                                                                 <  2:1852874/150‑1 (MQ=255)
            gTGGGTACAGCACCCCTCAATCGCACGGGTAGCTGCTCAGCAAACTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGc                                                                                                                                >  2:104954/1‑150 (MQ=255)
                gTACAGCACCCCTCAATCGCACGGGTAGCTGCTCAGCAAACTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGAt                                                                                                                            <  1:322572/150‑1 (MQ=255)
                   cAGCACCCCTCAATCGCACGGGTAGCTGCTCAGCAAACTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGAc                                                                                                                         <  2:302334/150‑1 (MQ=255)
                     gCACCCCTCAATCGCACGGGTAGCTGCTCAGCAAACTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCTGAGCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCt                                                                                                                       >  1:23903/1‑150 (MQ=255)
                      cACCCCTCAATCGCACGGGTAGCTGCTCAGCAAACTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTg                                                                                                                      <  1:1621841/150‑1 (MQ=255)
                      cACCCCTCAATCGCACGGGTAGCTGCTCAGCAAACTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCCAGCCGTTgcacgcgc                                                                                                                                <  2:223716/140‑7 (MQ=255)
                      cACCCCTCAATCGCACGGGTAGCTGCTCAGCAAACTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCCAGCCGTTgcacgcgc                                                                                                                                >  1:223716/1‑134 (MQ=255)
                       aCCCCTCAATCGCACGGGTAGCTGCTCAGCAAACTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGTCGTTGCGCTGGCCGATGACCTgg                                                                                                                     >  2:1087759/1‑150 (MQ=255)
                         cccTCAATCGCACGGGTAGCTGCTCAGCAAACTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGAGCTGGCCGATGACCTGGTc                                                                                                                   >  1:1251263/1‑150 (MQ=255)
                          ccTCAATCGTACGGGTAGCTGCTCAGCAAACTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCg                                                                                                                  <  2:120104/150‑1 (MQ=255)
                           ctgaaTCGCAAGGGTAGCTGCTCAGCAAACTGGCGCTACCCAGATAGGCTTGAGGGGGGGCGCGCATGCTCGGCGGCAACGTGCGCCCATTGCCCAGCAGGCACAGCGCGTGCAACGGCTGGCCCTTGCGCTGGCCGATGACCTGGTCgg                                                                                                                 <  2:1496884/147‑1 (MQ=255)
                           cTCAATCGCACGGGTAGCTGCTCAGCAAACTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCgg                                                                                                                 >  2:2114151/1‑150 (MQ=255)
                           cTCAATCGCACGGGTAGCTGCTCAGCAAACTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCgg                                                                                                                 >  2:2167980/1‑150 (MQ=255)
                           cTCAATCGCACGGGTAGCTGCTCAGCAAACTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCgg                                                                                                                 <  2:1395561/150‑1 (MQ=255)
                            tCAATCGCACGGGTAGCTGCTCAGCAAACTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGc                                                                                                                >  1:2009270/1‑150 (MQ=255)
                              aaTCGCACGGGTAGCTGCTCAGCAAACTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCaa                                                                                                              >  1:203683/1‑150 (MQ=255)
                                     cGGGTAGCTGCTCAGCAAACTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTggcg                                                                                                       <  1:2167980/150‑1 (MQ=255)
                                         tAGCTGCTCAGCAAACTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATTACCTGGTCGGCAAAGTGGcggcgg                                                                                                   <  2:1312386/150‑1 (MQ=255)
                                         tAGCTGCTCAGCAAACTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGcggcgg                                                                                                   >  1:2085228/1‑150 (MQ=255)
                                          aGCTGCTCAGCAAACTGGCGCTAGCCAGATAGGCTTGAGGGGTGGCCCGCATGCCCGGCGGCAACAGGCGCCAATCGACCAGCCGGCACAGAGCGTGCAACGGCTGCGTAGGGCGCTTTCCGCTTACGCGGTCGGCAAAGTGgcggcgtg                                                                                                  >  2:123738/1‑148 (MQ=255)
                                           gctgctCAGCAAACTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTa                                                                                                 <  2:2119114/150‑1 (MQ=255)
                                                  aGCAAACTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTg                                                                                                          >  1:774154/1‑134 (MQ=255)
                                                  aGCAAACTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTg                                                                                                          <  2:774154/134‑1 (MQ=255)
                                                  aGCAAACTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTAcc                                                                                          >  2:1575100/1‑150 (MQ=255)
                                                   gCAAACTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGcc                                                                                                                               >  1:943352/1‑112 (MQ=255)
                                                   gCAAACTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGcc                                                                                                                               <  2:943352/112‑1 (MQ=255)
                                                   gCAAACTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCt                                                                                         >  1:1818845/1‑150 (MQ=255)
                                                      aaCTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTggg                                                                                      >  2:974324/1‑150 (MQ=255)
                                                           gcgcTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGa                                                                                 >  2:1606728/1‑150 (MQ=255)
                                                           gcgcTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGa                                                                                 >  1:1274968/1‑150 (MQ=255)
                                                           gcgcTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGa                                                                                 >  1:297022/1‑150 (MQ=255)
                                                            cgcTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGAt                                                                                <  1:220634/150‑1 (MQ=255)
                                                             gcTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGAtt                                                                               <  2:2009270/150‑1 (MQ=255)
                                                                   cAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTgg                                                                         >  2:2269400/1‑150 (MQ=255)
                                                                   cAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTgg                                                                         >  1:280474/1‑150 (MQ=255)
                                                                     gATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCt                                                                       >  2:1490715/1‑150 (MQ=255)
                                                                           cTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCtgt                                                                     <  2:1746670/146‑1 (MQ=255)
                                                                           cTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCtgt                                                                     >  1:1746670/1‑146 (MQ=255)
                                                                           cTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCt                                                                 <  2:1274968/150‑1 (MQ=255)
                                                                            ttGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGTCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTg                                                                >  2:1717488/1‑150 (MQ=255)
                                                                            ttGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTAGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTg                                                                >  2:703105/1‑150 (MQ=255)
                                                                                ggggTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCg                                                            <  2:2160686/150‑1 (MQ=255)
                                                                                   gTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCGGcc                                                         >  1:558012/1‑150 (MQ=255)
                                                                                   gTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCGGcc                                                         >  1:805201/1‑150 (MQ=255)
                                                                                      gcgcgcATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGTCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCGGCCTgc                                                      <  1:1087759/150‑1 (MQ=255)
                                                                                      gcgcgcATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCGGCCTgc                                                      <  2:1088715/150‑1 (MQ=255)
                                                                                      gcgcgcATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCGGCCTgc                                                      >  2:2187023/1‑150 (MQ=255)
                                                                                        gcgcATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTAGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCGGCCTgcgc                                                    <  1:703105/150‑1 (MQ=255)
                                                                                           cATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCCAGCCGTTg                                                                                                                                       >  2:579420/1‑64 (MQ=255)
                                                                                           cATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCCAGCCGTTg                                                                                                                                       <  1:579420/64‑1 (MQ=255)
                                                                                             tGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCGGCCTGCGCCAGCg                                               <  2:280474/150‑1 (MQ=255)
                                                                                                   gcggcAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCGGCCTGCGCCAGCGTTTGTc                                         <  2:297022/150‑1 (MQ=255)
                                                                                                        aaCGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTATCTGGGTCTGATTCAGTGGCTGTGCCTGGTCGGCCTGCGCCAGCGTTTGTCGGGcc                                    >  1:2190544/1‑150 (MQ=255)
                                                                                                            tGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCGGCCTGCGCTAGCGTTTGTCGGGCCAGGt                                >  2:1539314/1‑150 (MQ=255)
                                                                                                               gcCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCGGCCTGCTCCAGCGTTTGTCGGGCCAGGTCgg                             <  2:203683/150‑1 (MQ=255)
                                                                                                               gcCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCGGCCTGCGCCAGCGTTTGTCGGGCCAGGTCgg                             <  2:170736/150‑1 (MQ=255)
                                                                                                               gcCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCGGCCTGCGCCAGCGTTTGTCGGGCCAGGTCgg                             <  1:104954/150‑1 (MQ=255)
                                                                                                               gcCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCGGCCTGCGCCAGCGTTTGTCGGGCCAGGTCgg                             <  2:2085228/150‑1 (MQ=255)
                                                                                                                   aTCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCGGCCTgc                                                      >  1:1947621/1‑121 (MQ=255)
                                                                                                                   aTCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCGGCCTgc                                                      <  2:1947621/121‑1 (MQ=255)
                                                                                                                   aTCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCGGCCTGCGCCAGCGTTTGTCGGGCCAGGTCGGGCaa                         >  1:581958/1‑150 (MQ=255)
                                                                                                                      gACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCGGCCTGCGCCAGCGTTTGTCGGGCCAGGTCGGGCAAGCg                      <  1:974324/150‑1 (MQ=255)
                                                                                                                       aCCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCGGCCTGCGCCAGCGTTTGTCGGGCCAGGTCGGGCAAGCgg                     <  1:1539314/150‑1 (MQ=255)
                                                                                                                         cagcagGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCGGCCTGCGCCAg                                                 >  1:1440107/1‑120 (MQ=255)
                                                                                                                         cagcagGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCGGCCTGCGCCAGCGTTTGTCGGGCCAGGTCGGGCAAGCgaga                   >  1:2309345/1‑147 (MQ=255)
                                                                                                                         cagcagGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCGGCCTGCGCCAGCGTTTGTCGGGCCAGGTCGGGCAAGCGGGc                   >  1:1822213/1‑150 (MQ=255)
                                                                                                                         cagcagGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCGGCCTGCGCCAGCGTTTGTCGGGCCAGGTCGGGCAAGCGGGc                   <  1:1717488/150‑1 (MQ=255)
                                                                                                                         cagcagGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGACTCAGTGGCTGTGCCTGGTCGGCCTGCGCCAg                                                 <  2:1440107/120‑1 (MQ=255)
                                                                                                                              ggCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCGGCCTGCGCCAGCGTTTGTCGGGCCAGGTCGGGCAAGCGGGCTGTgg              <  2:1822213/150‑1 (MQ=255)
                                                                                                                               gCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGAACTGGTTCAGTGGCTGTGCCTGGTAGGCCTGCGCCAGCGGGTGTCGGGCCAGGTCGGGCAAGCGGGCTGTggg             >  1:521932/1‑150 (MQ=255)
                                                                                                                                 acaGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCGGCCTGCGCCAGc                                                <  2:259831/113‑1 (MQ=255)
                                                                                                                                 acaGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCGGCCTGCGCCAGc                                                >  1:259831/1‑113 (MQ=255)
                                                                                                                                 acaGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCGGCCTGCGCCAGCGTTTGTCGGGCCAGGTCGGGCAAGCGGGCTGTGGGcc           >  2:1347557/1‑150 (MQ=255)
                                                                                                                                    tggcgTGACACCGCAGGCCGTTGCGCTGACCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCGGTGCCTGGTGGGCCTGCGCCAGCGTTTGTCGGGCCAGGTCGGGCAAGCGGGCTGTGGGccctc        <  2:600197/148‑3 (MQ=255)
                                                                                                                                      gcgTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAATAGCTGTGCCTGGTCGGCCTGCGCCAGCGGTTGACGACACAGGGCGGGCCAGCGGGCTGTGGGCCCGCtc      >  1:600197/1‑149 (MQ=255)
                                                                                                                                      cagcgCAACGGCCAGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCGGCCTGCGCCAGCGTTTGTCGGGCCAGGTCGGGCAAGCGGGCTGTGGGCCCGCTg      >  2:338922/5‑150 (MQ=255)
                                                                                                                                         cgCAACGGCCAGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCGGCCTGCGCCAGCGTTTGTCGGGCCAGGTCGGGCAAGCGGGCTGTGGGCCCGCTGaaa   <  1:338922/149‑1 (MQ=255)
                                                                                                                                          gCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCGGCCTGCGCCAGCGTTTGTCGGGCCAGGTCGGGCAAGCGGGCTGTGGGCCCGCTGAAAt  <  2:803885/150‑1 (MQ=255)
                                                                                                                                            aaCGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCGGCCTGCGCCAGCGTTTGTCgg                                       >  2:1103071/1‑111 (MQ=255)
                                                                                                                                            aaCGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCGGCCTGCGCCAGCGTTTGTCgg                                       <  1:1103071/111‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                            
CGATGATGTGCAGTGGGTACAGCACCCCTCAATCGCACGGGTAGCTGCTCAGCAAACTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCGGCCTGCGCCAGCGTTTGTCGGGCCAGGTCGGGCAAGCGGGCTGTGGGCCCGCTGAAAT  >  NC_002947/3296398‑3296685

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