Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,539,906 A→C 5.0% R376R (CGT→CGG)  PP_4872 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,539,9060AC 5.0% 229.5 / 4.4 80R376R (CGT→CGG) PP_4872hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (42/34);  new base C (2/2);  total (44/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.83e-01

TCGTGCGAGCGGTTCTGGATGATGCGGCCTTCGTGCTCGGTGATCGAGCAGAAGGTACAGCCACCGAAGCAGCCACGCATGATGTTCACCGAGAAACGGATCATCTCGTAGGCCGGGATGCGCTCCTTGCCATAGGCCGG‑GTGCGGCACACGGGCGTATGGCATGCCGAACACGTAGTCCATTTCCTCGGTGCTCATGGGGATGGGTGGCGGGTTGAACCACACATCCACTTCACCGTGCTTCTGCACCAGCGCGCGAGCATTGCCCGGGTTGGTCTCCAGGTGCAGCACGCGGTT  >  NC_002947/5539757‑5540052
                                                                                                                                                      |                                                                                                                                                  
tCGTGCGAGCGGTTCTGGATGATGCGGCCTTCGTGCTCGGTGATCGAGCAGAAGGTACAGCCACCGAAGCAGCCACGCATGATGTTCACCGAGAAACGGATCATCTCGTAGGCCGGGATGCGCTCCTTGCCATAGGCCGG‑GTGCGGcaca                                                                                                                                                    >  1:1935585/1‑150 (MQ=255)
        gCGGTTCTGGATGATGCGGCCTTCGTGCTCGGTGATCGAGCAGAAGGTACAGCCACCGAAGCAGCCACGCATGATGTTCACCGAGAAACGGATCATCTCGTAGGCCGGGATGCGCTCCTTGCCATAGGCCGG‑GTGCGGCACACGGGCGta                                                                                                                                            >  2:70694‑M2/1‑148 (MQ=255)
            ttCTGGATGATGCGGCCTTCGTGCTCGGTGATCGAGCAGAAGGTACAGCCACCGAAGCAGCCACGCATGATGTTCACCGAGAAACGGATCATCTCGTAGGCCGGGATGCGATCCTTCACATATGCCGG‑GTGCGGCACACGGGCGtatggc                                                                                                                                        <  2:1610017‑M2/150‑7 (MQ=255)
               tGGATGATGCGGCCTTCGTGCTCGGTGATCGAGCAGAAGGTACAGCCACCGAAGCAGCCACGCATGATGTTCACCGAGAAACGGATCATCTCGTAGGCCGGGATGCGCTCCTTGCCATAGGCCGG‑GTGCGGCACACGGGCGTATGGCATg                                                                                                                                     >  1:1490614/1‑150 (MQ=255)
                    gatgCGGCCTTCGTGCTCGGTGATCGAGCAGAAGGTACAGCCACCGAAGCAGCCACGCATGATGTTCACCGAGAAACGGATCATCTCGTAGGCCGGGATGCGCTCCTTGCCATAGGCCGG‑GTGCGGCACACGGGCGTATGGCATGCCGaa                                                                                                                                >  1:1099626/1‑150 (MQ=255)
                      tgCGGCCTTCGTGCTCGGTGATCGAGCAGAAGGTACAGCCACCGAAGCAGCCACGCATGATGTTCACCGAGAAACGGATCATCTCGTAGGCCGGGATGCGCTCCTTGCCATAGGCCGG‑GTGCGGCACACGGGCGTATGGCATGCCGAaca                                                                                                                              >  1:430693/1‑150 (MQ=255)
                      tgCGGCCTTCGTGCTCGGTGATCGAGCAGAAGGTACAGCCACCGAAGCAGCCACGCATGATGTTCACCGAGAAACGGATCATCTCGTAGGCCGGGATGCGCTCCTTGCCATAGGCCGG‑GTGCGGCACACGGGCGTATGGCATGCCGAaca                                                                                                                              >  1:655192/1‑150 (MQ=255)
                      tgCGGCCTTCGTGCTCGGTGATCGAGCAGAAGGTACAGCCACCGAAGCAGCCACGCATGATGTTCACCGAGAAACGGATCATCTCGTAGGCCGGGATGCGCTCCTTGCCATAGACCGG‑TAGCGGCCCACGGGCGTATGGCATGCCGAgca                                                                                                                              >  2:694644/1‑150 (MQ=255)
                       gCGGCCTTCGTGCTCGGTGATCGAGCAGAAGGTACAGCCACCGAAGCAGCCACGCATGATGTTCACCGAGAAACGGATCATCTCGTAGGCCGGGATGCGCTCCTTGCCATAGGCCGG‑GTGCGGCACACGGGCGTATGGCATGCCGAacac                                                                                                                             <  2:1671009/150‑1 (MQ=255)
                           ccTTCGTGCTCGGTGATCGAGCAGAAGGTACAGCCACCGAAGCAGCCACGCATGATGTTCACCGAGAAACGGATCATCTCGTAGGCCGGGATGCGCTCCTTGCCATAGGCCGG‑GTGCGGCACACGGGCGTATGGCATGCCGAACACGTAg                                                                                                                         >  1:186298/1‑150 (MQ=255)
                            cTTCGTGCTCGGTGATCGAGCAGAAGGTACAGCCACCGAAGCAGCCACGCATGATGTTCACCGAGAAACGGATCATCTCGTAGGCCGGGATGCGCTCCTTGCCATACGCCCGTGTGCCGCACCCGGCCatggcaaggagcgcatcccg                                                                                                                           <  1:1689692‑M2/149‑21 (MQ=255)
                            cTTCGTGCTCGGTGATCGAGCAGAAGGTACAGCCACCGAAGCAGCCACGCATGATGTTCACCGAGAAACGGATCATCTCGTAGGCCGGGATGCGCTCCTTGCCATACGCCCGTGTGCCGCACCCGGCCatggcaaggagcgcatcccg                                                                                                                           >  2:1689692‑M2/1‑129 (MQ=255)
                                gTGCTCGGTGATCGAGCAGAAGGTACAGCCACCGAAGCAGCCACGCATGATGTTCACCGAGAAACGGATCATCTCGTAGGCCGGGATGCGCTCCTTGCCATAGGCCGG‑GTGCGGCACACGGGCGTATGGCATGCCGAACACGTAGTCCAt                                                                                                                    >  1:1898529/1‑150 (MQ=255)
                                  gCTCGGTGATCGAGCAGAAGGTACAGCCACCGAAGCAGCCACGCATGATGTTCACCGAGAAACGGATCATCTCGTAGGCCGGGATGCGCTCCTTGCCATAGGCCGG‑GTGCGGCACACGGGCGTATGGCATGCCGAACACGTAGTCCAttt                                                                                                                  >  1:1091023/1‑150 (MQ=255)
                                    tCGGTGATCGAGCAGAAGGTACAGCCACCAACGCAGCCACGCATGCTGTACAGCGAGGTCCGGATCATCTCGTATGCCGGCAGTGGCACGTTGCCTTATGGCTG‑GGGGGGCTCAAGGGCGTATGGCATGCCGGACACGTAGTCCATTTct                                                                                                                >  2:219809/1‑149 (MQ=255)
                                          aTCGAGCAGAAGGTACAGCCACCGAAGCAGCCACGCATGATGTTCACCGAGAAACGGATCATCTCGTAGGCCGGGATGCGCTCCTTGCCATAGGCCGG‑GTGCGGCACACGGGCGTATGGCATGCCGAACACGTAGTCCATTTCCTCGGTg                                                                                                          <  2:186298/150‑1 (MQ=255)
                                              aGCAGAAGGTACAGCCACCGAAGCAGCCACGCATGATGTTCACCGAGAAACGGATCATCTCGTAGGCCGGGATGCGCTCCTTGCCATAGGCCGG‑GTGCGGCACACGGGCGTATGGCATGCCGAACACGTAGTCCATTTCCTCGGTGCTCa                                                                                                      >  2:980680/1‑150 (MQ=255)
                                               gCAGAAGGTACAGCCACCGAAGCAGCCACGCATGATGTTCACCGAGAAACGGATCATCTCGTAGGCCGGGATGCGCTCCTTGCCATAGGCCGG‑GTGCGGCACACGGGCGTATGGCATGCCGAACACGTAGTCCATTTCCTCGGTGCTCAt                                                                                                     <  2:655192/150‑1 (MQ=255)
                                                 aGAAGGTACAGCCACCGAAGCAGCCACGCATGATGTTCACCGAGAAACGGATCATCTCGTAGGCCGGGATGCGCTCCTTGCCATAGGCCGG‑GTGCGGCACACGGGCGTATGGCATGCCGAACACGTAGTCCATTTCCTCGGTGCTCATgg                                                                                                   <  2:430693/150‑1 (MQ=255)
                                                  gAAGGTACAGCCACCGAAGCAGCCACGCATGATGTTCACCGAGAAACGGATCATCTCGTAGGCCGGGATGCGCTCCTTGCCATAGGCCGG‑GTGCGGCACACGGGCGTATGGCATGCCGAACACGTAGTCCGTTTCCCCGGTGCTCATggg                                                                                                  >  2:1281225/1‑150 (MQ=255)
                                                  gAAGGTACAGCCACCGAAGCAGCCACGCATGATGTTCACCGAGAAACGGATCATCTCGTAGGCCGGGATGCGCTCCTTGCCATAGGCCGG‑GTGCGGCACACGGGCGTATGGCATGCCGAACACGTAGTCCATTTCCt                                                                                                               <  1:1463233/137‑1 (MQ=255)
                                                  gAAGGTACAGCCACCGAAGCAGCCACGCATGATGTTCACCGAGAAACGGATCATCTCGTAGGCCGGGATGCGCTCCTTGCCATAGGCCGG‑GTGCGGCACACGGGCGTATGGCATGCCGAACACGTAGTCCATTTCCt                                                                                                               >  2:1463233/1‑137 (MQ=255)
                                                      gTACAGCCACCGAAGCAGCCACGCATGATGTTCACCGAGAAACGGATCATCTCGTAGGCCGGGATGCGCTCCTTGCCCCGGCCTTT‑GTGCGGCACACGGGCGTATGGCATGCCGAACACGTAGTCCATTTCCTCGGTGCTCATGGGGATg                                                                                              >  1:1575828/1‑150 (MQ=255)
                                                        aCAGCCACCGAAGGAGCCACGCATGATGTTCACCGAGAAACGGATCATCTCGTAGGCCGGGATGCGCTTCTTGCCATAGGCCGG‑GTGCGGCACACGGGCGTATGGCATGCCGCCCTGGTAGTCCATTTCCTCGGTGCTCATGGGGATggg                                                                                            <  1:1610017/150‑1 (MQ=255)
                                                          aGCCACCGAAGCAGCCACGCATGATGTTCACCGAGAAACGGATCATCTCGTAGGCCGGGATGCGCTCCTTGCCATAGGCCGG‑GTGCGGCACACGGGCGTATGGCATGCCGAACACGTAGTCCATTTCCTCGGTGCTCATGGGGATGGGTg                                                                                          >  1:1550864/1‑150 (MQ=255)
                                                           gCCACCGAAGCAGCCACGCATGATGTTCACCGAGAAACGGATCATCTCGTAGGCCGGGATGCGCTCCTTGCCATAGGCCGG‑GTGCGGCACACGGGCGTATGGCATGCCGAACACGTAGTCCATTTCCTCGGTGCTCATGGGGATGGGTgg                                                                                         >  2:1932516/1‑150 (MQ=255)
                                                            ccACCGAAGCAGCCACGCATGATGTTCACCGAGAAACGGATCATCTCGTAGGCCGGGATGCGCTCCTTGCCATAGGCCGG‑GTGCGGCACACGGGCGTATGGCATGCCGAACACGTAGTCCATTTCCTCGGTGCTCATGGGGATGGGTGGc                                                                                        >  2:79720/1‑150 (MQ=255)
                                                              aCCGAAGCAGCCACGCATGATGTTCACCGAGAAACGGATCATCTCGTAGGCCGGGATGCGCTCCTTGCCATAGGCCGG‑GTGCGGCACACGGGCGTATGGCATGCCGAACACGTAGTCCATTTCCTCGGTGCTCATGGGGATGGGTGGCgg                                                                                      <  2:1101546/150‑1 (MQ=255)
                                                              aCCGAAGCAGCCACGCATGATGTTCACCGAGAAACGGATCATCTCGTAGGCCGGGATGCGCTCCTTGCCATACGCCCGTGTGCCGCACCCGGCCatggcaagga                                                                                                                                     >  1:939761‑M2/1‑95 (MQ=255)
                                                              aCCGAAGCAGCCACGCATGATGTTCACCGAGAAACGGATCATCTCGTAGGCCGGGATGCGCTCCTTGCCATACGCCCGTGTGCCGCACCCGGCCatggcaagga                                                                                                                                     <  2:939761‑M2/105‑11 (MQ=255)
                                                                 gAAGCAGCCACGCATGATGTTCACCGAGAAACGGATCATCTCGTAGGCCGGGATGCGCTCCTTGCCATAGGCCGG‑GTGCGGCACACGGGCGTATGGCATGCCGAACACGTAGTCCATTTCCTCGGTGCTCATGGGGATGGGTGGCGGGtt                                                                                   <  2:1550864/150‑1 (MQ=255)
                                                                    gcagcCACGCATGATGTTCACCGAGAAACGGATCATCTCGTAGGCCGGGATGCGCTCCTTGCCATAGGCCGG‑GTGCGGCACACGGGCGTATGGCATGCCGAACACGTAGTCCATTTCCTCGGTGCTCATGGGGATGGGTGGCGGGTTGaa                                                                                <  1:79720/150‑1 (MQ=255)
                                                                     cagcCACGCATGATGTTCACCGAGAAACGGATCATCTCGTAGGCCGGGATGCGCTCCTTGCCATAGGCCGG‑GTGCGGCACACGGGCGTATGGCATGCCGAACACGTAGTCCATTTCCTCGGTGCTCATGGGGATGGGTGGCGGGTTGAAc                                                                               <  2:1898529/150‑1 (MQ=255)
                                                                       gcCACGCATGATGTTCACCGAGAAACGGATCATCTCGTAGGCCGGGATGCGCTCCTTGCCATAGGCCGG‑GTGCGGCACACGGGCGTATGGCATGCCGAACACGTAGTCCATTTCCTCGGTGCTCATGGGGATGGGTGGCGGGTTGAACca                                                                             >  2:1702122/1‑150 (MQ=255)
                                                                         cACGCATGATGTTCACCGAGAAACGGATCATCTCGTAGGCCGGGATGCGCTCCTTGCCATAGGCCGG‑GTGCGGCACACGGGCGTATGGCATg                                                                                                                                     <  2:1364411/92‑1 (MQ=255)
                                                                         cACGCATGATGTTCACCGAGAAACGGATCATCTCGTAGGCCGGGATGCGCTCCTTGCCATAGGCCGG‑GTGCGGCACACGGGCGTATGGCATg                                                                                                                                     >  1:1364411/1‑92 (MQ=255)
                                                                         cACGCATGATGTTCACCGAGAAACGGATCATCTCGTAGGCCGGGATGCGCTCCTTGCCATAGGCCGG‑GTGCGGCACACGGGCGTATGGCATGCCGAACACGTAGTCCATTTCCTCGGTGCTCATGGGGATGGGTGGCg                                                                                       <  2:1604339/138‑1 (MQ=255)
                                                                         cACGCATGATGTTCACCGAGAAACGGATCATCTCGTAGGCCGGGATGCGCTCCTTGCCATAGGCCGG‑GTGCGGCACACGGGCGTATGGCATGCCGAACACGTAGTCCATTTCCTCGGTGCTCATGGGGATGGGTGGCg                                                                                       >  1:1604339/1‑138 (MQ=255)
                                                                             cATGATGTTCACCGAGAAACGGATCATCTCGTAGGCCGGGATGCGCTCCTTGCCATAGGCCGG‑GTGCGGCACACGGGCGTATGGCATGCCGAACACGTAGTCCATTTCCTCGGTGCTCATgg                                                                                                   <  2:718908/122‑1 (MQ=255)
                                                                             cATGATGTTCACCGAGAAACGGATCATCTCGTAGGCCGGGATGCGCTCCTTGCCATAGGCCGG‑GTGCGGCACACGGGCGTATGGCATGCCGAACACGTAGTCCATTTCCTCGGTGCTCATgg                                                                                                   >  1:718908/1‑122 (MQ=255)
                                                                             cATGATGTTCACCGAGAAACGGATCATCTCGTAGGCCGGGATGCGCTCCTTGCCATAGGCCGG‑GTGCGGCACACGGGCGTATGGCATGCCGAACACGTAGTCCATTTCCTCGGTGCTCATGGGGATGGGTGGCGGGTTGAACCACACATc                                                                       <  1:980680/150‑1 (MQ=255)
                                                                             cATGATGTTCACCGAGAAACGGATCATCTCGTAGGCCGGGATGCGCTCCTTGCCATAGGCCGG‑GTGCGGCACACGGGCGTATGGCATGCCGAACACGTAGTCCATTTCCTCGGTGCTCATGGGGATGGGTGGCGGGTTGAACCACACATc                                                                       <  2:1490614/150‑1 (MQ=255)
                                                                                gatgTTCACCGAGAAACGGATCATCTCGTAGGCCGGGATGCGCTCCTTGCCATAGGCCGG‑GTGCGGCACACGGGCGTATGGCATGCCGAACACGTAGTCCATTTCCTCGGTGCTCATGGGGAt                                                                                               >  2:1031367/1‑123 (MQ=255)
                                                                                gatgTTCACCGAGAAACGGATCATCTCGTAGGCCGGGATGCGCTCCTTGCCATAGGCCGG‑GTGCGGCACACGGGCGTATGGCATGCCGAACACGTAGTCCATTTCCTCGGTGCTCATGGGGAt                                                                                               <  1:1031367/123‑1 (MQ=255)
                                                                                gatgCTCACCGAGAAACGGATCATCTCGTAGGCCGGGATGCGCTCCTTGCCATAGGCCGG‑GTGCGGCACACGGGCGTATGGCATGCCGAACACGTAGTCCATTTCCTCGGTGCTCATGGGGATGGGTGGCGGGTTGAACCACACATCCAc                                                                    <  1:1287553/150‑1 (MQ=255)
                                                                                 atgTTCACCGAGAAACGGATCATTTCGGAGGCCAGGATGAGCTCCGTAAAATATACCGT‑AGGTGGCAGATGGACGTATGGCATGCAGAATACGTAGTCCATTGCATAGGTGATCATGGGGATGAGCTGCGTGTTGAACCACACATCCACt                                                                   >  2:400612/1‑150 (MQ=255)
                                                                                   gTTCACCGAGAAACGGATCATCTCGTAGGCCGGGATGCGCTCCTTGCCATAGGCCGG‑GTGCGGCACACGGGCGTATGGCATGCCGAACACGTAGTCCATTTCCTCGGTGCTCATg                                                                                                    <  2:1085768/115‑1 (MQ=255)
                                                                                   gTTCACCGAGAAACGGATCATCTCGTAGGCCGGGATGCGCTCCTTGCCATAGGCCGG‑GTGCGGCACACGGGCGTATGGCATGCCGAACACGTAGTCCATTTCCTCGGTGCTCATg                                                                                                    >  1:1085768/1‑115 (MQ=255)
                                                                                      cACCGAGAAACGGATCATCTCGTAGGCCGGGATGCGCTCCTTGCCATAGGCCGG‑GTGCGGCACACGGGCGTATGGCATGCCGAACACGTAGTCCATTTCCTCGGTGCTCATGGGGATGGGTGGCGGGTCGAACCACACATCCACTTCAcc                                                              <  2:1575828/150‑1 (MQ=255)
                                                                                       aCCGAGAAACGGATCATCTCGTAGGCCGGGATGCGCTCCTTGCCATAGGCCGG‑GTGCGGCACACGGGCGTATGGCATGCCGAACACGTAGTCCATTTCCTCGGTGCTCATGGGGATGGGTGGCGGGTTGAACCACACATCCACTTCACCg                                                             >  1:334585/1‑150 (MQ=255)
                                                                                         cGAGAAACGGATCATCTCGTAGGCCGGGATGCGCTCCTTGCCATAGGCCGG‑GTGCGGCACACGGGCGTATGGCATGCCGAACACGTAGTCCATTTCCTCGGTGCTCATGGGGATGGGTGGCGGGTTGAACCACACATCCACTTCACCGTg                                                           >  1:446983/1‑150 (MQ=255)
                                                                                          gagaAACGGATCATCTCGTAGGCCGGGATGCGCTCCTTGCCATAGGCCGG‑GTGCGGCACACGGGCGTATGGCATGCCGAACACGTAGTCCATTTCCTCGGTGCTCATGGGGATGGGTGGCGGGTTGAACCACACATCCACTTCACCGTGc                                                          >  1:706456/1‑150 (MQ=255)
                                                                                          gagaAACGGATCATCTCGTAGGCCGGGATGCGCTCCTTGCCATAGGCCGG‑GTGCGGCACACGGGCGTATGGCATGCCGAACACGTAGTCCATTTCCTCGGTGCTCATGGGGATGGGTGGCGGGTTGAACCACACATCCACTTCACCGTGc                                                          <  1:1702122/150‑1 (MQ=255)
                                                                                          gagaAACGGATCATCTCGTAGGCCGGGATGCGCTCCTTGCCATAGGCCGG‑GTGCGGCACACGGGCGTATGGCATGCCGAACACGTAGTCCATTTCCTCGGTGCTCATGGGGATGGGTGGCGGGTTGAACCACACACCCACTTCACCGTGc                                                          <  2:1091023/150‑1 (MQ=255)
                                                                                               aCGGATCATCTCGTAGGCCGGGATGCGCTCCTTGCCATAGGCCGG‑GTGCGGCACACGGGCGTATGGCATGCCGAACACGTAGTCCATTTCCTCGGTGCTCATGGGGATGGGTGGCGGGTTGAACCACACATCCACTTCACCGTGCTTCTg                                                     <  1:1932516/150‑1 (MQ=255)
                                                                                                  gATCATCTCGTAGGCCGGGATGCGCTCCTTGCCATAGGCCGG‑GTGCGGCACACGGGCGTATGGCATGCCGAACACGTAGTCCATTTCCTCGGTGCTCATGGGGATGGGTGGCGGGTTGAACCACACATCCACTTCACCGTGCTTCTGTAc                                                  <  2:1935585/150‑1 (MQ=255)
                                                                                                  gATCATCTCGTAGGCCGGGATGCGCTCCTTGCCATAGGCCGG‑GTGCGGCACACGGGCGTATGGCATGCCGAACACGTAGTCCATTTCCTCGGTGCTCATGGGGATGGGTGGCGGGTTGAACCACACATCCACTTCACCGTGCTTCTGCAc                                                  >  2:1693008/1‑150 (MQ=255)
                                                                                                    tcatcTCGTAGGCCGGGATGCGCTCCTTGCCATAGGCCGG‑GTGCGGCACACGGGCGTATGGCATGCCGAACACGTAGTCCATTTCCTCGGTGCTCATGGGGATGGGTGGCGGGTTGAACCACACATCCACTTCACCGTGCTTCTGCACCa                                                <  2:260612/150‑1 (MQ=255)
                                                                                                       tctcGTAGGCCGGGATGCGCTCCTTGCCATAGGCCGG‑GTGCGGCACACGGGCGTATGGCATGCCGAACACGTAGTCCATTTCCTCGGTGCTCATGGGGATGGGTGGCGGGTTGAACCACACATCCACTTCACCGTGCTTCTGCACCAgcg                                             >  1:362546/1‑150 (MQ=255)
                                                                                                              ggCCGGGATGCGCTCCTTGCCATAGGCCGG‑GTGCGGCACACGGGCGTATGGCATGCCGAACACGTAGTCCATTTCCTCGGTGCTCATGGGGATGGGTGGCGGGTTGAACCACACATCCACTTCACCGTGCTTCTGCACCAGCGCGCGAGc                                      >  1:1689221/1‑150 (MQ=255)
                                                                                                               gCCGGGATGCGCTCCTTGCCATAGGCCGG‑GTGCGGCACACGGGCGTATGGCATGCCGAACACGTAGTCCATTTCCTCGGTGCTCATGGGGATGGGTGGCGGGTTGACCCACACATCCACTTCAc                                                               <  2:989476/124‑1 (MQ=255)
                                                                                                               gCCGGGATGCGCTCCTTGCCATAGGCCGG‑GTGCGGCACACGGGCGTATGGCATGCCGAACACGTAGTCCATTTCCTCGGTGCTCATGGGGATGGGTGGCGGGTTGAACCACACATCCACTTCAc                                                               >  1:989476/1‑124 (MQ=255)
                                                                                                               gCCGGGATGCGCTCCTTGCCATAGGCCGG‑GTGCGACACACGGGCGTATGGCATGCCAAACACGTAGTCCATTTCCTCGGTGCTCATGGGGATGGGGGGCGGGTTGAACCACACATCCACTTCACCGTACTTCTGCACCAGCGCGCGAGCa                                     >  2:526811/1‑150 (MQ=255)
                                                                                                                   ggATGCGCTCCTTGCCATAGGCCGG‑GTGCGGCACACGGGCGTATGGCATGCCGAACACGTAGTCCATTTCCTCGGTGCTCATGGGGATGGGTGGCGGGTTGAACCACACATCCACTTCACCGTGCTTCTGCACCAGCGCGCGAGCATTGc                                 >  1:779184/1‑150 (MQ=255)
                                                                                                                       gcgcTCCTTGCCATAGGCCGG‑GTGCGGCACACGGGCGTATGGCATGCCGAACACGTAGTCCATTTCCTCGGTGCTCATGGGGATGGGTGGCGGGTTGAACCACACATCCACTTCACCGTGCTTCTGCAc                                                  >  2:1799649/1‑129 (MQ=255)
                                                                                                                       gcgcTCCTTGCCATAGGCCGG‑GTGCGGCACACGGGCGTATGGCATGCCGAACACGTAGTCCATTTCCTCGGTGCTCATGGGGATGGGTGGCGGGTTGAACCACACATCCACTTCACCGTGCTTCTGCAc                                                  <  1:1799649/129‑1 (MQ=255)
                                                                                                                         gcTCCTTGCCATAGGCCGG‑GTGCGGCACACGGGCGTATGGCATGCCGAACACGTAGTCCATTTCCTCGGTGCTCATGGGGATGGGTGGCGGGTTGAACCACACATCCACTTCACCGTGCTTCTGCACCAGCGCGCGAGCATTGCCCGGGt                           >  2:1694414/1‑150 (MQ=255)
                                                                                                                         gcTCCTTGCCATAGGCCGG‑GTGCGGCACACGGGCGTATGGCATGCCGAACACGTAGTCCATTTCCTCGGTGCTCATGGGGATGGGTGGCGGGTTGAACCACACATCCACTTCACCGTGCTTCTGCACCAGCGCGCGAGCATTGCCCGGGt                           <  1:219809/150‑1 (MQ=255)
                                                                                                                          cTCCTTGCCATAGGCCGG‑GTGCGGCACACGGGCGTATGGCATGCCGAACACGTAGTCCATTTCCTCGGTGCTCATGGGGATGGGTGGCGGGTTGAACCACACATCCACTTCACCGTGCTTCTGCACCAGCGCGCGAGCATTGCCCGGGtt                          <  2:1099626/150‑1 (MQ=255)
                                                                                                                           tCCTTGCCATAGGCCGG‑GTGCGGCACACGGGCGTATGGCATGCCGAACACGTAGTCCATTTCCTCGGTGCTCATGGGGATGGGTGGCGGGTTGAACCACACATCCACTTCACCGTGCTTCTGCACCAGCGCGCGAGCATTGCCCGGGTTg                         >  1:953091/1‑150 (MQ=255)
                                                                                                                             cTTGCCATAGGCCGG‑GTGCGGCACACGGGCGTATGGCATGCCGAACACGTAGTCCATTTCCTCGGTGCTCATGGGGATGGGTGGCGGGTTGAACCACACATCCACTTCACCGTGCTTCTGCACCAGCGCGCGAGCATTGCCCGGGTTGGt                       >  2:1247340/1‑150 (MQ=255)
                                                                                                                             cTTGCCATAGGCCGG‑GTGCGGCACACGGGCGTATGGCATGCCGAACACGTAGTCCATTTCCTCGGTGCTCATGGGGATGGGTGGCGGGTTGAACCACACATCCACTTCACCGTGCTTCTGCACCAGCGCGCGAGCATTGCCCGGGTTGGt                       >  2:329884/1‑150 (MQ=255)
                                                                                                                             cTTGCCATAGGCCGG‑GTGCGGCACACGGGCGTATGGCATGCCGAACACGTAGTCCATTTCCTCGGTGCTCATGGGGATGGGTGGCGGGTTGAACCACACATCCACTTCACCGTGCTTCTGCACCAGCGCGCGAGCATTGCCCGGGTTGGt                       >  2:1324447/1‑150 (MQ=255)
                                                                                                                              ttGCCATAGGCCGG‑GTGCGGCACACGGGCGTATGGCATGCCGAACACGTAGTCCATTTCCTCGGTGCTCATGGGGATGGGTGGCGGGTTGAACCACACATCCACTTCACCGTGCTTCTGCACCAGCGCGCGAGCATTGCCCGGGTTg                         >  2:261859/1‑147 (MQ=255)
                                                                                                                              ttGCCATAGGCCGG‑GTGCGGCACACGGGCGTATGGCATGCCGAACACGTAGTCCATTTCCTCGGTGCTCATGGGGATGGGTGGCGGGTTGAACCACACATCCACTTCACCGTGCTTCTGCACCAGCGCGCGAGCATTGCCCGGGTTg                         <  1:261859/147‑1 (MQ=255)
                                                                                                                              ttGCCATAGGCCGG‑GTGCGGCACACGGGCGTATGGCATGCCGAACACGTAGTCCATTTCCTCGGTGCTCATGGGGATGGGTGGCGGGTTGAACCACACATCCACTTCACCGTGCTTCTGCACCAGCGCGCGAGCATTGCCCGGGTTGGtc                      >  2:466162/1‑150 (MQ=255)
                                                                                                                              ttGCCATAGGCCGG‑GTGCGGCACACGGGCGTATGGCATGCCGAACACGTAGTCCATTTCCTCGGTGCTCATGGGGATGGGTGGCGGGTTGAACCACACATCCACTTCACCGTGCTTCTGCACCAGCGCGCGAGCATTGCCCGGGTTGGtc                      <  2:779184/150‑1 (MQ=255)
                                                                                                                                   aTAGGCCGG‑GTGCGGCACACGGGCGTATGGCATGCCGAACACGTAGTCCATTTCCTCGGTGCTCATGGGGATGGGTGGCGGGTTGAACCACACATCCACTTCACCGTGCTTCTGCACCAGCGCGCGAGCATTGCCCGGGTTGGTCTCCAg                 <  2:1209336/150‑1 (MQ=255)
                                                                                                                                     aGGCCGG‑GTGCGGCACACGGGCGTATGGCATGCCGAACACGTAGTCCATTTCCTCGGTGCTCATGGGGATGGGTGGCGGGTTGAACCACACATCCACTTCACCGTGCTTCTGCACCAGCGCGCGAGCATTGCCCGGGTTGGTCTCCAggc               <  2:953091/150‑2 (MQ=255)
                                                                                                                                         cGG‑GTGCGGCACACGGGCGTATGGCATGCCGAACACGTAGTCCATTTCCTCGGTGCTCATGGGGATGGGTGGCGGGTTGAACCACACATCCACTTCACCGTGCTTCTGCACCAGCGCGCGAGCATTGCCCGGGTTGGTCTCCAGGTgcag           >  2:296640/1‑150 (MQ=255)
                                                                                                                                                cGGCACACGGGCGTATGGCATGCCGAACACGTAGTCCATTTCCTCGGTGCTCATGGGGATGGGTGGCGGGTTGCTCATTACATCCACTTCGTGGTACGCCTGCACCAGCGCGCGAGCATTGCCCGGGTTGGTCTCCAGGTGCAGCAcgcg     <  1:1693008/150‑1 (MQ=255)
                                                                                                                                                   cacaCGGGCGTATGGCATGCCGAACACGTAGTCCATTTCCTCGGTGCTCATGGGGATGGGTGGCGGGTTGAACCACACATCCACTTCACCGTGCTTCTGCACCAGCGCGCGAGCATTGCCCGGGTTGGTCTCCAGGTGCAGCACGCGGtt  >  2:1261704/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                  
TCGTGCGAGCGGTTCTGGATGATGCGGCCTTCGTGCTCGGTGATCGAGCAGAAGGTACAGCCACCGAAGCAGCCACGCATGATGTTCACCGAGAAACGGATCATCTCGTAGGCCGGGATGCGCTCCTTGCCATAGGCCGG‑GTGCGGCACACGGGCGTATGGCATGCCGAACACGTAGTCCATTTCCTCGGTGCTCATGGGGATGGGTGGCGGGTTGAACCACACATCCACTTCACCGTGCTTCTGCACCAGCGCGCGAGCATTGCCCGGGTTGGTCTCCAGGTGCAGCACGCGGTT  >  NC_002947/5539757‑5540052

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: ‑