Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 151,062 T→C 5.7% intergenic (‑113/+19) PP_0142 ← / ← PP_0143 ABC transporter permease/membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947151,0620TC 5.7% 264.6 / 4.4 87intergenic (‑113/+19)PP_0142/PP_0143ABC transporter permease/membrane protein
Reads supporting (aligned to +/- strand):  ref base T (44/38);  new base C (3/2);  total (47/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.39e-01

TGTCCAGGGTGGCGCTGGGGTTGGTCATAGGCTCCATGCAATTGACTCGGCGCTAAGGCTTAGGGGCTGATCATAGCGCAGGGGGCATGGGCTGTTTGGTGATTATGTGCTGTGTGTACGCTTGGCAACAGTGACCTGCAGCCGCCTTGTGCTGCGCAAGGCAGGGGGTCAACTACCCGCGGGGGCATCATTGACCACCTTGAACCGCAGCACGCCGATCACTTGCCCATCCTCGGTCAGCACCCTCACCTGCCACTTGCCTACCGGGTTGGGCGGGAAGTTCTGCTTGTGGGTCC  >  NC_002947/150913‑151208
                                                                                                                                                     |                                                                                                                                                  
tgtCCAGGGTGGCGCTGGGGTTGGTCATAGGCTCCATGCAATTGACTCGGCGCTAAGGCTTAGGGGCTGATCATAGCGCAGGGGGCATGGGCTGTTTGGTGATTATGTGCTGTGTGTACGCTTGGCAACAGTGACCTGCAGCCGCCTtgt                                                                                                                                                    <  2:1292794/150‑1 (MQ=255)
    cAGGGTGGCGCTGGGGTTGGTCATAGGCTCCATGCAATTGACTCGGCGCTAAGGCTTAGGGGCTGATCATAGCGCAGGGGGCATGGGCTGTTTGGTGATTATGTGCTGTGTGTACGCTTGGCAACAGTGACCTGCAGCCGCCTTGtgctg                                                                                                                                                >  1:2254689/1‑150 (MQ=255)
    cAGGGTGGCGCTGGGGTTGGTCATAGGCTCCATGCAATTGACTCGGCGCTAAGGCTTAGGGGCTGATCATAGCGCACGGGGCATGGGCTGTTTGGTGATTATGTGCTGTGTGTACTCTGGGCAACAGTGACCTGCAGCCGCCTTGtgctg                                                                                                                                                >  1:1738150/1‑150 (MQ=255)
        gTGGCGCTGGGGTTGGTCATAGGCTCCATGCAATTGACTCGGCGCTAAGGCTTAGGGGCTGATCATAGCGCAGGGGGCATGGGCTGTTTGGTGATTATGTGCTGTGTGTACGCTTGGCAACAGTGACCTGCAGCCGCCTTGTGCTGCGCa                                                                                                                                            >  2:227519/1‑150 (MQ=255)
            cgcTGGGGTTGGTCATAGGCTCCATGCAATTGACTCGGCGCTAAGGCTTAGGGGCTGATCATAGCGCAGGGGGCATGGGCTGTTTGGTGATTATGTGCTGTGTGTACGCTTGGCAACAGTGACCTGCAGCCGCCTTGTGCTGCGCAAGGc                                                                                                                                        <  1:822330/150‑1 (MQ=255)
             gcTGGGGTTGGTCATAGGCTCCATGCAATTGACTCGGCGCTAAGGCTTAGGGGCTGATCATAGCGCAGGGGGCATGGGCTGTTTGGTGATTATGTGCTGTGTGTACGCTTGGCAACAGTGACCTGCAGCCGCCTTGTGCTGCGCAAGGCa                                                                                                                                       >  2:2047301/1‑150 (MQ=255)
             gcTGGGGTTGGTCATAGGCTCCATGCAATTGACTCGGCGCTAAGGCTTAGGGGCTGATCATAGCGCAGGGGGCATGGGCTGTTTGGTGATTATGTGCTGTGTGTACGCTTGGCAACAGTGACCTGCAGCCGCCTTGTGCTGCGCAAGGCa                                                                                                                                       >  2:2033181/1‑150 (MQ=255)
             gcTGGGGTTGGTCATAGGCTCCATGCAATTGACTCGGCGCTAAGGCTTAGGGGCTGATCATAGCGCAGGGGGCATGGGCTGTTTGGTGATTATGTGCTGTGTGTACGCTTGGCAACAGTGACCTGCAGCCGCCTTGCGCAGCACAAGGcg                                                                                                                                       <  2:370023/150‑2 (MQ=255)
              cTGGGTTTGGTTCTAGGCTCTATGCAATTGTCTCGGCGCTAGGGGTTAGGGGCTGATCATAGCGCAGGGGGCATGGGCTGTTTGTTGATTATGTGCTGTGTGGACGCTTGGCAACAGTGACCTGCAGCTGCCTTGTGCTGCGCAAGGCAg                                                                                                                                      <  2:545625/150‑1 (MQ=255)
               tGGGGTTGGTCATAGGCTCCATGCAATTGACTCGGCGCTAAGGCTTAGGGGCTGATCATAGCGCAGGGGGCATGGGCTGTTTGGTGATTATGTGCTGTGTGTACGCTTGGCAACAGTGACCTGCAGCCGCCTTGTGCTGCGCAAGGCAgg                                                                                                                                     >  2:2375068/1‑150 (MQ=255)
                        tCATAGGCTCCATGCAATTGACTCGGCGCTAAGGCTTAGGGGCTGATCATAGCGCAGGGGGCATGGGCTGTTTGGTGATTATGTGCTGTGTGTACGCTTGGCAACAGTGACCTGCAGCCGCCTTGTGCTGCGCAAGGCAGGGGGTCAACt                                                                                                                            <  1:851425/150‑1 (MQ=255)
                           tAGGCTCCATGCAATTGACTCGGCGCTAAGGCTTAGGGGCTGATCATAGCGCAGGGGGCATGGGCTGTTTGGTGATTATGTGCTGTGTGTACGCTTGGCAACAGTGACCTGCAGCCGCCTTGTGCTGCGCAAGGCAGGGGGTCAACTAcc                                                                                                                         >  1:829195/1‑150 (MQ=255)
                            aGGCTCCATGCAATTGACTCGGCGCTAAGGCTTAGGGGCTGATCATAGCGCAGGGGGCATGGGCTGTTTGGTGATTATGTGCTGTGTGTACGCTTGGCAACAGTGACCTGCAGCCGCCTTGTGCTGCGCAAGGCAGGGGGTCAACTAccc                                                                                                                        <  1:2047301/150‑1 (MQ=255)
                              gCTCCATGCAATTGACTCGGCGCTAAGGCTTAGGGGCTGATCATAGCGCAGGGGGCATGGGCTGTTTGGTGATTATGTGCTGTGTGTACGCTTGGCAACAGTGACCTGCAGCCGCCTTGTGCTGCGCAAGGCAGGGGGTCAACTACCcgc                                                                                                                      >  1:1205510/1‑150 (MQ=255)
                               cTCCATGCAATTGACTCGGCGCTAAGGCTTAGGGGCTGATCATAGCGCAGGGGGCATGGGCTGTTTGGTGATTATGTGCTGTGTGTACGCTTGGCAACAGTGACCTGCAGCCGCCTTGTGCTGCGCAAGGCAGGGGGTCaa                                                                                                                              >  2:610309/1‑141 (MQ=255)
                               cTCCATGCAATTGACTCGGCGCTAAGGCTTAGGGGCTGATCATAGCGCAGGGGGCATGGGCTGTTTGGTGATTATGTGCTGTGTGTACGCTTGGCAACAGTGACCTGCAGCCGCCTTGTGCTGCGCAAGGCAGGGGGTCaa                                                                                                                              <  1:610309/141‑1 (MQ=255)
                               cTCCATGCAATTGACTCGGCGCTAAGGCTTAGGGGCTGATCATAGCGCAGGGGGCATGGGCTGTTTGGTGATTATGTGCTGTGTGTACGCTTGGCAACAGTGACCTGCAGCCGCCTTGTGCTGCGCAAGGCAGGGGGTCAACTACCcgcg                                                                                                                     <  1:2375068/150‑1 (MQ=255)
                                tCCATGCAATTGACTCGGCGCTAAGGCTTAGGGGCTGATCATAGCGCAGGGGGCATGGGCTGTTTGGTGATTATGTGCTGTGTGTACGCTTGGCAACAGTGACCTGCAGCCGCCTTGTGCTGCGCAAGGCAGTGGGTCAACTACCCGCgg                                                                                                                    >  2:1346375/1‑150 (MQ=255)
                                       aaTTGACTCGGCGCTAAGGCTTAGGGGCTGATCATAGCGCAGGTGGCATGGGCTGTTTGGTGATTATGTGCTGTGTGTACGCTTGGCAACAGTGACCTGCAGCCGCCTTGtgctgc                                                                                                                                               <  1:2097369/116‑1 (MQ=255)
                                       aaTTGACTCGGCGCTAAGGCTTAGGGGCTGATCATAGCGCAGGGGGCATGGGCTGTTTGGTGATTATGTGCTGTGTGTACGCTTGGCAACAGTGACCTGCAGCCGCCTTGtgctgc                                                                                                                                               >  2:2097369/1‑116 (MQ=255)
                                       aaTTGACTCGGCGCTAAGGCTTAGGGGCTGATCATAGCGCAGGGGGCATGGGCTGTTTGGTGATTATGTGCTGTGTGTACGCTTGGCAACAGTGACCTGCAGCCGCCTTGTGCTGCGCAAGGCAGGGGGt                                                                                                                                 >  2:608497/1‑130 (MQ=255)
                                       aaTTGACTCGGCGCTAAGGCTTAGGGGCTGATCATAGCGCAGGGGGCATGGGCTGTTTGGTGATTATGTGCTGTGTGTACGCTTGGCAACAGTGACCTGCAGCCGCCTTGTGCTGCGCAAGGCAGGGGGt                                                                                                                                 <  1:608497/130‑1 (MQ=255)
                                          tGACTCGGCGCTAAGGCTTAGGGGCTGATCATAGCGCAGGGGGCATGGGCTGTTTGGTGATTATGTGCTGTGTGTACGCTTGGCAACAGTGACCTGCAGCCGCCTTGTGCTGCGCAAGGCAGGGGGTCAACTACCCGCGGGGGCATCAtt                                                                                                          <  2:1608535/150‑1 (MQ=255)
                                             cTCGGCGCTAAGGCTTAGGGGCTGATCATAGCGCAGGGGGCATGGGCTGTTTGGTGATTATGTGCTGTGTGTACGCTTGGCAACAGTGACCTGCAGCCGCCTTGTGCTGCGCAAGGCAGGGGGTCAACTACCCGCGGGGGCATCATTGac                                                                                                       <  1:226720/150‑1 (MQ=255)
                                              tCGGCGCTAAGGCTTAGGGGCTGATCATAGCGCAGGGGGCATGGGCTGTTTGGTGATTATGTGCTGTGTGTACGTTTGGCAACAGTGACCTGCAGCCGCCTTGTGCTGCGCAAGGCAGGGGGTCAACTACCCGCGGGGGCATCATTGacc                                                                                                      >  2:2308547/1‑150 (MQ=255)
                                              tCGGCGCTAAGGCTTAGGGGCTGATCATAGCGCAGGGGGCATGGGCTGTTCGGTGATTATGTGCTGTGTGTACGCTTGGCAACAGTGACCTGCAGCCGCCTTGTGCTGCGCAAGGCAGGGGGTCAACTACCCGCGGGGGCATCATTGacc                                                                                                      >  2:1542179/1‑150 (MQ=255)
                                                ggCGCTAAGGCTTAGGGGCTGATCATAGCGCAGGGGGCATGGGCTGTTTGGTGATTATGTGCTGTGTGTACGCTTGGCAACAGTGACCTGCAGCCGCCTTGTGCTGCGCAAGGCAGGGGGTCAACTACCCGCGGGGGCATCATTGaccac                                                                                                    >  2:844842/1‑150 (MQ=255)
                                                ggCGCTAAGGCTTAGGGGCTGATCATAGCGCAGGGGGCATGGGCTGTTTGGTGATTATGTGCTGTGTGTACGCTTGGCAACAGTGACCTGCAGCCGCCTTGTGCTGCGCAAGGCAGCGGGTCAACTACCCGCGGCGGCATCATTGaccac                                                                                                    <  2:1971212/150‑1 (MQ=255)
                                                 gcgcTAAGGCTTAGGGGCTGATCATAGCGCAGGGGGCATGGGCTGTTTGGTGATTATGTGCTGTGTGTACGCTTGGCAACAGTGACCTGCAGCCGCCTTGTGCTGCGCAAGGCAgg                                                                                                                                     >  1:1344102/1‑116 (MQ=255)
                                                 gcgcTAAGGCTTAGGGGCTGATCATAGCGCAGGGGGCATGGGCTGTTTGGTGATTATGTGCTGTGTGTACGCTTGGCAACAGTGACCTGCAGCCGCCTTGTGCTGCGCAAGGCAgg                                                                                                                                     <  2:1344102/116‑1 (MQ=255)
                                                   gcTAAGGCTTAGGGGCTGATCATAGCGCAGGGGGCATGGGCTGTTTGGTGATTATGTGCTGTGTGTACGCTTGGCAACAGTGACCTGCAGCCGCCTTGTGCTGCGCaa                                                                                                                                           <  2:1820157/108‑1 (MQ=255)
                                                   gcTAAGGCTTAGGGGCTGATCATAGCGCAGGGGGCATGGGCTGTTTGGTGATTATGTGCTGTGTGTACGCTTGGCAACAGTGACCTGCAGCCGCCTTGTGCTGCGCaa                                                                                                                                           >  1:1820157/1‑108 (MQ=255)
                                                         gCTTAGGGGCTGATCATAGCGCAGGGGGCATGGGCTGTTTGGTGATTATGTGCTGTGTGTACGCTTGGCAACAGTGACCTGCAGCCGCCTTGTGCTGCGCAAGGCAGGGGGTCAACTACCCGCGGGGGCATCATTGACCACCTTGAACCg                                                                                           >  2:2024413/1‑150 (MQ=255)
                                                           ttAGGGGCTGATCATAGCGCAGGGGGCATGGGCTGTTTGGTGATTATGTGCTGTGTGTACGTTTGGCAACAGTGACCTGCAGCCGCCTTGTGCTGCGCAAGGCAGGGGGTCAACTACCCGCGGGGGCATCATTGACCACCTTGAACCgca                                                                                         <  1:2308547/150‑1 (MQ=255)
                                                           ttAGGGGCTGATCATAGCGCAGGGGGCATGGGCTGTTTGGTGATTATGTGCTGTGTGTACGCTTGGCAACAGTGACCTGCAGCCGCCTTGTGCTGCGCAAGGCAGGGGGTCAACTACCCGCGGGGGCATCATTGACCACCTTGAACCgca                                                                                         >  1:2381543/1‑150 (MQ=255)
                                                                 gCTGATCATAGCGCAGGGGGCATGGGCTGTTTGGTGATTATGTGCTGTGTGTACGCTTGGCAACAGTGACCTGCAGCCGCCTtgtg                                                                                                                                                   <  1:2349543/86‑1 (MQ=255)
                                                                 gCTGATCATAGCGCAGGGGGCATGGGCTGTTTGGTGATTATGTGCTGTGTGTACGCTTGGCAACAGTGACCTGCAGCCGCCTtgtg                                                                                                                                                   >  2:2349543/1‑86 (MQ=255)
                                                                  cTGATCATAGCGCAGGGGGCATGGCCTGTTTGGTGATTATGTGCTGTGTGTACGCTTGGCAACAGTGACCTGCAGCCGCCTTGTGCTGCGCAAGGCAGGGGGTCAACTACCCGCGGGGGCATCATTGACCACCTTGAACCGCAGCACGcc                                                                                  >  2:57681/1‑150 (MQ=255)
                                                                   tGATCATAGCGCAGGGGGCATGGGCTGTTTGGTGATTATGTGCTGTGTGTACGCTTGGCAACAGTGACCTGCAGCCGCCTTGTGCTGCGCAAGGCAGGGGGTCAACTACCCGCGGGGGCATCATTGACCACCTTGAACCGCAGCACGCCg                                                                                 <  2:2254689/150‑1 (MQ=255)
                                                                     aTCATAGCGCAGGGGGCATGGGCTGTTTGGTGATTATGTGCTGTGTGTACGCTTGGCAACAGTGACCTGCAGCCGCCTTGTGCTGCGCAAGGCAGGGGGTCAACTACCCGCGGGGGCATCATTGACCACCTTGAACCGCAGCACGCCGAt                                                                               <  1:227519/150‑1 (MQ=255)
                                                                     aTCATAGCGCAGGGGGCATGGGCTGTTTGGTGATTATGTACTGTGTGTACGCTTGGCAACAGTGACCTGCAGCCGCCTTGTGCTGCGCAAGGCAGGGGGTCAACTACCCGCGGGGGCATCATTGACCACCTTGAACCGCAGCACGCCGAt                                                                               <  1:563101/150‑1 (MQ=255)
                                                                      tCATAGCGCAGGGGGCATGGGCTGTTTGGTGATTATGTGCTGTGTGTACGCTTGGCAACAGTGACCTGCAGCCGCCTTGTGCTGCGCAAGGCAGGGGGTCAACTACCCGCgggg                                                                                                                  <  1:1758071/114‑1 (MQ=255)
                                                                      tCATAGCGCAGGGGGCATGGGCTGTTTGGTGATTATGTGCTGTGTGTACGCTTGGCAACAGTGACCTGCAGCCGCCTTGTGCTGCGCAAGGCAGGGGGTCAACTACCCGCgggg                                                                                                                  >  2:1758071/1‑114 (MQ=255)
                                                                           gcgcAGGGGGCATGGGCTGATTGGTGATTATGTGCTGTGTGTACGCTTGGCAACAGTGACCTGCAGCCGCCTTGTGCTGCGCAAGGCAGCGGGGCACACATCCGCGGGGGCATCATTGACCCCCTTGAACCGCAGCACGCCGATCACTTg                                                                         >  2:1271265/1‑150 (MQ=255)
                                                                            cgcAGGGGGCATGGGCTGTTTGGTGATTATGTGCTGTGTGTACGCTTGGCAACAGTGACCTGCAGCCGCCTTGTGCTGCGCAAGGCAGGGGGTCAACTACCCGCGGGGGCATCATTGACCACCTTGAACCGCAGCACGCCGATCACTTGc                                                                        >  2:1546388/1‑150 (MQ=255)
                                                                            cgcAGGGGGCATGGGCTGTTTGGTGATTATGTGCTGTGAGGACGCTTGGCAACAGTGACCTGCAGCCGCCTTGTGCTGCGTAAGGCAGGGGGTCAACTACCCGCGGGGGCATCATTGACCACCTTGAACCGCAGCACGCCGATCACTTGc                                                                        <  1:738606/150‑1 (MQ=255)
                                                                                   ggCATGGGCTGTTTGGTGATTATGTGCTGTGTGTACGCTTGGCAACAGTGACCTGCAGCCGCCTTGTGCTGCGCAAGGCAGGGGGTCAACTACCCGCGGGGGCATCATTGACCACCTTGAACCGCAGCACGCCGATCACTTGCCCATCCt                                                                 >  2:875842/1‑150 (MQ=255)
                                                                                      aTGGGCTGTTTGGTGATTATGTGCTGTGTGTACGCTTGGCAACAGTGACCTGCAGCCGCCTTGTGCTGCGCAAGGCAGGGGGTCAACTACCCGCGGGGGCATCATTGACCACCTTGAACCGCAGCACGCCGATCACTTGCCCATCCTCgg                                                              <  1:2033181/150‑1 (MQ=255)
                                                                                       tGGGCTGTTTGGTGATTATGTGCTGTGTGTACGCTTGGCAACAGTGACCTGCAGCCGCCTTGTGCTGCGCAAGGCAGGGGGTCAACTACCCGCGGGGGCATCATTGACCACCTTGAACCGCAGCACGCCGATCACTTGCCCATCCTCGGt                                                             <  2:1205510/150‑1 (MQ=255)
                                                                                          gCTGTTTGGTGATTATGTGCTGTGTGTACGCTTGGCAACAGTGACCTGCAGCCGCCTTGTGCTGCGCAAGGCAGGGGGTCAACTACCCGCgg                                                                                                                    >  1:476303/1‑92 (MQ=255)
                                                                                          gCTGTTTGGTGATTATGTGCTGTGTGTACGCTTGGCAACAGTGACCTGCAGCCGCCTTGTGCTGCGCAAGGCAGGGGGTCAACTACCCGCgg                                                                                                                    <  2:476303/92‑1 (MQ=255)
                                                                                           cTGTTTGGTGATTATGTGCTGTGTGTACGCTTGGCAACAGTGACCTGCAGCCGCCTTGTGCTGCGCAAGGCAGGGGGTCAACTACCcgcg                                                                                                                     <  2:485212/90‑1 (MQ=255)
                                                                                           cTGTTTGGTGATTATGTGCTGTGTGTACGCTTGGCAACAGTGACCTGCAGCCGCCTTGTGCTGCGCAAGGCAGGGGGTCAACTACCcgcg                                                                                                                     >  1:485212/1‑90 (MQ=255)
                                                                                              tttGGTGATTATGTGCTGTGTGTACGCTTGGCAACAGTGACCTGCAGCCGCCTTGTGCTGCGCAAGGCAGGGGGTCAACTACCCGCGGGGGCATCATTGACCACCTTGAACCGCAGCACGCCGATCACTTGCCCATCCTCGGTCAGCAcc                                                      >  2:2140952/1‑150 (MQ=255)
                                                                                              tttGGTGATTATGTGCTGTGTGTACGCTTGGCAACAGTGACCTGCAGCCGCCTTGTGCTGCGCAAGGCAGGGGGTCAACTACCCGCGGGGGCATCATTGACCACCTTGAACCGCAGCACGCCGATCACTTGCCCATCCTCGGTCAGCAcc                                                      >  1:1412871/1‑150 (MQ=255)
                                                                                               ttGGTGATTATGTGCTGTGTGTACGCTTGGCAACAGTGACCTGCAGCCGCCTTGTGCTGCGCAAGGCAGGGGGTCAACTACCCGCGGGGGCATCATTGACCACCTTGAACCgcag                                                                                        >  2:2181080/1‑115 (MQ=255)
                                                                                               ttGGTGATTATGTGCTGTGTGTACGCTTGGCAACAGTGACCTGCAGCCGCCTTGTGCTGCGCAAGGCAGGGGGTCAACTACCCGCGGGGGCATCATTGACCACCTTGAACCgcag                                                                                        <  1:2181080/115‑1 (MQ=255)
                                                                                                  gTGATTATGTGCTGTGTGTACGCTTGGCAACAGTGACCTGCAGCCGCCTTGTGCTGCGCAAGGCAGGGGGTCAACTACCCGCGGGGGCATCATTGACCACCTTGAACCGCAGCACGCCGATCACTTGCCCATCCTCGGTCAGCACCCTCa                                                  >  1:1304778/1‑150 (MQ=255)
                                                                                                   tGATTATGTGCTGTGTGTACGCTTGGCAACAGTGACCTGCAGCCGCCTTGTGCTGCGCAAGGCAGGGGGTCAACTACCCGCGGGGGCATCATTGACCACCTTGAACCGCAGCACGCCGATCACTTGCCCATCCTCGGTCAGCACCCTCAc                                                 >  2:2026511/1‑150 (MQ=255)
                                                                                                     aTTATGTGCTGTGTGTACGCTTGGCAACAGTGACCTGCAGCCGCCTTGTGCTGCGCAAGGCAGGGGGTCAACTACCCGCGGGGGCATCATTGACCACCTTGAACCGCAGCACGCCGATCACTTGCCCATCCTCGGTCAGCACCCTCACCt                                               >  2:743755/1‑150 (MQ=255)
                                                                                                          gtgCTGTGTGTACGCTTGGCAACAGTGACCTGCAGCCGCCTTGTGCTGCGCAAGGCAGGGGGTCAACTACCCGCGGGGGCATCATTGACCACCTTGAACCGCAGCACGCCGATCACTTGCCCATCCTCGGTCAGCACCCTCACCTGCCAc                                          >  1:1288758/1‑150 (MQ=255)
                                                                                                             cTGTGTGTACGCTTGGCAACAGTGACCTGCAGCCGCCTTGTGCTGCGCAAGGCAGGGGGTCAACTACCCGCGGGGGCATCATTGACCACCTTGAACCGCAGCACGCCGATCACTTGCCCATCCTCGGTCAGCACCCTCACCTGCCACTTg                                       >  1:51003/1‑150 (MQ=255)
                                                                                                               gtgtgtACGCTTGGCAACAGTGACCTGCAGCCGCCTTGTGCTGCGCAGGGCAGGGGGTCAACTACCCGCGGGGGCATCATTGACCACCTTGAACCGCAGCACGCCGATCACTTGCCCATCCTCGGTCAGCACCCTCACCTGCCACTTGcc                                     >  2:1615215/1‑150 (MQ=255)
                                                                                                                tgtgtACGCTTGGCAACAGTGACCTGCAGCCGCCTTGTGCTGCGCAAGGCAGGGGGTCAACTACCCGCGGGGGCATCATTGACCACCTTGAACCGCAGCACGCCGATCACTTGCCCATCCTCGGTCAGCACCCTCACCTGCCACTTGCCt                                    >  1:1691639/1‑150 (MQ=255)
                                                                                                                   gtACGCTTGGCAACAGTGACCTGCAGCCGCCTTGTGCTGCGCAAGGCAGGGGGTCAACTACCCGCGGGGGCATCATTGACCACCTTGAACCGCAGCACGCCGATCACTTGCCCATcc                                                                  <  1:1064362/117‑1 (MQ=255)
                                                                                                                   gtACGCTTGGCAACAGTGACCTGCAGCCGCCTTGTGCTGCGCAAGGCAGGGGGTCAACTACCCGCGGGGGCATCATTGACCACCTTGAACCGCAGCACGCCGATCACTTGCCCATcc                                                                  >  2:1064362/1‑117 (MQ=255)
                                                                                                                    tACGCTTGGCAACAGTGACCTGCAGCCGCCTTGTGCTTCGCAAGGCAGGGGGTCAACTACCCGCGGGGGCATCATTGACCACCTTGAACCGCAGCACGCCGATCACTTGCCCATCCTCGGTCAGCACCCTCACCTGCCACTTGCCTACCg                                >  1:192403/1‑150 (MQ=255)
                                                                                                                    tACGCTTGGCAACAGTGACCTGCAGCCGCCTTGTGCTGCGCAAGGCAGGGGGTCAACTACCCGCGGGGGCATCATTGACCACCTTGAACCGCAGCACGCCGATCACTTGCCCATCCTCGGTCAGCACCCTCACCTGCCACTTGCCTACCg                                <  1:684979/150‑1 (MQ=255)
                                                                                                                      cGCTTGGCAACAGTGACCTGCAGCCGCCTTGTGCTGCGCAAgg                                                                                                                                         <  2:361735/43‑1 (MQ=255)
                                                                                                                      cGCTTGGCAACAGTGACCTGCAGCCGCCTTGTGCTGCGCAAgg                                                                                                                                         >  1:361735/1‑43 (MQ=255)
                                                                                                                        cTTGGCAACAGTGACCTGCAGCCGCCTTGTGCTGCGCAAGGCAGGGGGTCAACTACCCGCGGGGGCATCATTGACCACCTTGAACCGCAGCACGCCGATCACTTGCCCATCCTCGGTCAGCACCCTCACCTGCCACTTGCCTACCGGGtt                            <  1:2024413/150‑1 (MQ=255)
                                                                                                                         ttGGCAACAGTGACCTGCAGCCGCCTTGTGCTGCGCAAGGCAGGGGGTCAACTACCCGCGGGGGCATCATTGACCACCtt                                                                                                 <  1:801971/80‑1 (MQ=255)
                                                                                                                         ttGGCAACAGTGACCTGCAGCCGCCTTGTGCTGCGCAAGGCAGGGGGTCAACTACCCGCGGGGGCATCATTGACCACCtt                                                                                                 >  2:801971/1‑80 (MQ=255)
                                                                                                                           ggCAACAGTGACCTGCAGCCGCCTTGTGCTGCGCAAGGCAGGGGGTCAACTACCCGCGGGGGCATCATTGACCACCTTGACCCGCAGCACGCCGATCACTTGCCCATCCTCGGTCAGCACCCTCACCTGCCACTTGCCTACCGGGTTggg                         >  2:1669606/1‑150 (MQ=255)
                                                                                                                            gCAACAGTGACCTGCAGCCGCCTTGTGCTGCGCAAGGCAGGGGGTCAACTACCCGCGGGGGCATCATTGACCACCTTGAACCGCAGCACGCCGATCACTTGCCCATCCTCGGTCAGCACCCTCACCTGCCACTTGCCTACCGGGTTGGGc                        <  2:2057608/150‑1 (MQ=255)
                                                                                                                              aaCAGTGACCTGCAGCCGCCTTGTGATGCGCAAGGCAGGGGGTCAACTACCCGCGGGGGCATCATTGACCACCTTGAACCGCAGCACGCCGATCACTTGCCCATCCTCGGTCAGCACCCTCACCTGCCACTTGCCTACCGGGTTGGGCgg                      <  2:2111472/150‑1 (MQ=255)
                                                                                                                               aCAGTGACCTGCAGCCGCCTTGTGCTGCGCAAGGCAGGGGGTCAACTACCCGCGGGGGCATCATTGACCACCTTGAACCGCAGCACGCCGATCACTTGCCCATCCTCGGTCAGCACCCTCACCTGCCACTTGCCTACCGGGTTGGGCggg                     >  1:399421/1‑150 (MQ=255)
                                                                                                                                cAGTGACCTGCAGCCGCCTTGTGCTGCGCAAGGCAGGGGGTCAACTACCCGCGGGGGCATCATTGACCACCTTGAACCGCAGCACGCCGATCACTTGCCCATCCTCGGTCAGCACCCTCACCTGCCACTTg                                       <  2:1979776/131‑1 (MQ=255)
                                                                                                                                cAGTGACCTGCAGCCGCCTTGTGCTGCGCAAGGCAGGGGGTCAACTACCCGCGGGGGCATCATTGACCACCTTGAACCGCAGCACGCCGATCACTTGCCCATCCTCGGTCAGCACCCTCACCTGCCACTTg                                       >  1:1979776/1‑131 (MQ=255)
                                                                                                                                     aCCTGCAGCCGCCTTGTGCTGCGCAAGGCAGGGGGTCAACTACCCGCGGGGGCATCATTGACCACCTTGAACCGCAGCACGCCGATCACTTGCCCATCCTCGGTCAGCACCCTCACCTGCCACTTGCCTACCGGGTTGGGCGGGAAGTTc               <  1:743755/150‑1 (MQ=255)
                                                                                                                                       cTGCAGCCGCCTTGTGCTGCGCAAGGCAGGGGGTGAACGACCCGCGGGGGCATCATTGACCACCTTGAACCGCAGAACGCCGAGCACTTGCCCATCCTCGGTCAGCACCCTCACCTGCCACTTGCCTACCGGGTTGGGCGGGAAGTTCTg             >  1:804301/1‑150 (MQ=255)
                                                                                                                                          cccctgccTTGCGCAGCACAAGGCAGGGGGTCAACTACCCGCGGGGGCATCATTGACCACCTTGAACCGCAGCACGCCGATCACTTGCCCATCCTCGGTCAGCACCCTCACCTGCCACTTGCCTACCGGGTTGGGCGGGAAGtt                <  2:2095654/139‑1 (MQ=255)
                                                                                                                                          cccctgccTTGCGCAGCACAAGGCAGGGGGTCAACTACCCGCGGGGGCATCATTGACCACCTTGAACCGCAGCACGCCGATCACTTGCCCATCCTCGGTCAGCACCCTCACCTGCCACTTGCCTACCGGGTTGGGCGGGAAGtt                >  1:2095654/6‑144 (MQ=255)
                                                                                                                                          cccctgccTTGCGCAGCACAAGGCAGGGGGTCAACTACCCGCGGGGGCATCATTGACCACCTTGAACCGCAGCACGCCGATCACTTGCCCATCCTCGGTCAGCACCCTCACCTGCCACTTGCCTACCGGGTTGGGCGGGAAGTTCTGCtt          >  1:1803801/6‑150 (MQ=255)
                                                                                                                                          cccctgccTTGCGCAGCACAAGGCAGGGGGTCAACTACCCGCGGGGGCATCATTGACCACCTTGAACCGCAGCACGCCGATCACTTGCCCATCCTCGGTCAGCACCCTCACCTGCCACTTGCCTACCGGGTTGGGCGGGAAGTTCTGCtt          >  2:865944/6‑150 (MQ=255)
                                                                                                                                            gccgccTTGTGCTGCGCAAGGCAGGGGGTCAACTACCCGCGGGGTCATCATTGACCACCTTGAACCGCAGCACGCCGATCACTTGCCCATCCTCGGTCAGCACCCTCACCTGCCACTTGCCTACCGGGTTGGGCGGGAAGTTCTGCTtgt        <  2:1738150/150‑1 (MQ=255)
                                                                                                                                            gccgccTTGTGCTGCGCAAGGCAGGGGGTCAACTACCCGCGGGGGCATCATTGACCACCTTGAACCGCAGCACGCCGATCACTTGCCCATCCTCGGTCAGCACCCTCACCTGCCACTTGCCTACCGGGTTGGGCGGGAAGTTCTGCTtgt        <  1:875842/150‑1 (MQ=255)
                                                                                                                                            gccgccTTGTGCTGCGCAAGGCAGGGGGTCAACTACCCGCGGGGGCATCATTGACCACCTTGAACCGCAGCACGCCGATCACTTGCCCATCCTCGGTCAGCACCCTCACCTGCCACTTGCCTACCGGGTTGGGCGGGAAGTTCTGCTtgt        <  1:57681/150‑1 (MQ=255)
                                                                                                                                                 cTTGTGCTGCGCAAGGCAGGGGGTCAACTACCCGCGGGGGCATCATTGACCACCTTGAACCGCAGCACGCCGATCACTTGCCCATCCTCGGTCAGCACCCTCACCTGCCACTTGCCTACCGGGTTggg                         <  1:1920542/128‑1 (MQ=255)
                                                                                                                                                 cTTGTGCTGCGCAAGGCAGGGGGTCAACTACCCGCGGGGGCATCATTGACCACCTTGAACCGCAGCACGCCGATCACTTGCCCATCCTCGGTCAGCACCCTCACCTGCCACTTGCCTACCGGGTTggg                         >  2:1920542/1‑128 (MQ=255)
                                                                                                                                                 cTTGTGCTGCGCAAGGCAGGGGGTCAACTACCCGCGGGGGCATCATTGACCACCTTGAACCGCAGCACGCCGATCACTTGCCCATCCTCGGTCAGCACCCTCACCTGCCACTTGCCTACCGGGTTGGGCGGGAAGTTCTGCTTGTGGGTc   >  1:4311/1‑150 (MQ=255)
                                                                                                                                                  ttGTGCTGCGCAAGGCAGGGGGTCAACTACCCGCGGGGGCATCATTGACCACCTTGAACCGCAGCACGCCGATCACTTGCCCATCCTCGGTCAGCACCCTCACCTGCCACTTGCCTACCGGGTTGGGCGGGAAGTTCTGCTTGTGGGTcc  <  1:2140952/150‑1 (MQ=255)
                                                                                                                                                  ttGTGCTGCGCAAGGCAGGGGGTCAACTACCCGCGGGGGCATCATTGACCACCTTGAACCGCAGCACGCCGATCACTTGCCCATCCTCGGTCAGCACCCTCACCTGCCACTTGCCTACCGGGTTGGGCGGGAAGTTCTGCTTGTGGGTcc  >  1:548913/1‑150 (MQ=255)
                                                                                                                                                  ttGTGCTGCGCAAGGCAGGGGGTCAACTACCCGCGGGGGCATCATTGACCACCTTGAACCGCAGCACGCCGATCACTTGCCCATCCTCGGTCAGCACCCTCACCTGCCACTTGCCTACCGGGTTGGGCGGGAAGTTCTGCTTGTGGGTcc  <  1:844842/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
TGTCCAGGGTGGCGCTGGGGTTGGTCATAGGCTCCATGCAATTGACTCGGCGCTAAGGCTTAGGGGCTGATCATAGCGCAGGGGGCATGGGCTGTTTGGTGATTATGTGCTGTGTGTACGCTTGGCAACAGTGACCTGCAGCCGCCTTGTGCTGCGCAAGGCAGGGGGTCAACTACCCGCGGGGGCATCATTGACCACCTTGAACCGCAGCACGCCGATCACTTGCCCATCCTCGGTCAGCACCCTCACCTGCCACTTGCCTACCGGGTTGGGCGGGAAGTTCTGCTTGTGGGTCC  >  NC_002947/150913‑151208

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