Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,296,544 T→C 7.8% Q125R (CAG→CGG) ‡ PP_2898 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,296,5440TC 7.8% 262.9 / 9.1 91Q125R (CAG→CGG) ‡PP_2898hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (43/40);  new base C (3/4);  total (46/45)
Fisher's exact test for biased strand distribution p-value = 7.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.38e-01

TCGCGATGATGTGCAGTGGGTACAGCACCCCTCAATCGCACGGGTAGCTGCTCAGCAAACTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCGGCCTGCGCCAGCGTTTGTCGGGCCAGGTCGGGCAAGCGGGCTGTGGGCCCGCTGAAAT  >  NC_002947/3296395‑3296685
                                                                                                                                                     |                                                                                                                                             
tCGCGATGATGTGCAGTGGGTACAGCACCCCTCAATCGCACGGGTAGCTGCTCAGCAAACTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCt                                                                                                                                               >  1:1567891/1‑150 (MQ=255)
   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)
                                                                                                                                                     |                                                                                                                                             
TCGCGATGATGTGCAGTGGGTACAGCACCCCTCAATCGCACGGGTAGCTGCTCAGCAAACTGGCGCTACCCAGATAGGCTTGAGGGGTGGCGCGCATGCTCGGCGGCAACGTGCGCCAATCGACCAGCAGGCACAGCGCGTGCAACGGCTGGCCGTTGCGCTGGCCGATGACCTGGTCGGCAAAGTGGCGGCGGTACAGTACCTGGGTCTGATTCAGTGGCTGTGCCTGGTCGGCCTGCGCCAGCGTTTGTCGGGCCAGGTCGGGCAAGCGGGCTGTGGGCCCGCTGAAAT  >  NC_002947/3296395‑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: ‑