Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,076,155 A→C 5.3% S155A (TCC→GCC)  PP_5328 ← phosphate ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,076,1550AC 5.3% 280.3 / 5.8 94S155A (TCC→GCC) PP_5328phosphate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (45/44);  new base C (2/3);  total (47/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.42e-01

GCTGACGTGGTCCAGCGCACGCCCCTGCTCGTCGAGCACGAACGATTGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCGTTACCGGTCTTGGTATCGAAGAACAATGCCTGGCCCTTGTCCGAAACACGC  >  NC_002947/6076009‑6076303
                                                                                                                                                  |                                                                                                                                                    
gCTGACGTGGTGCAGCGCACGCCCCTGCTCGTCGAGCACGAACGATTGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAg                                                                                                                                                   <  2:1618009/150‑1 (MQ=255)
gCTGACGTGGTCCAGCACACGCCCCTGCTCGTCGAGCACGAACGATTGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAg                                                                                                                                                   >  2:2050172/1‑150 (MQ=255)
 cTGACGTGGTCCAGCGCACGCCCCTGCTCGTCGAGCACGAACGATTGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAgg                                                                                                                                                  <  2:2310967/150‑1 (MQ=255)
  tGACGTGGTCCAGCGCACGCCCCTGCTCGTCGAGCACGAACGATTGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGTCCAATGGCCAGGc                                                                                                                                                 <  1:397307/150‑1 (MQ=255)
  tGACGTGGTCCAGCGCACGCCCCTGCTCGTCGAGCACGAACGATTGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGTCCAATGGCCAGGc                                                                                                                                                 <  2:326874/150‑1 (MQ=255)
    aCGTGGTCCAGCGCACGCCCCTGCTCGTCGAGCACGAACGATTGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGccc                                                                                                                                               <  1:1567986/150‑1 (MQ=255)
      gTGGTCCAGCGCACGCCCCTGCTCGTCGAGCACGAACGATTGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGTCCAATGGCCAGGcgctg                                                                                                                                             <  1:1546264/150‑5 (MQ=255)
      gTGGTCCAGCGCACGCCCCTGCTCGTCGAGCACGAACGATTGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAg                                                                                                                                             >  2:988043/1‑150 (MQ=255)
      gTGGTCCAGCGCACGCCCCTGCTCGTCGAGCACGAACGATTGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAg                                                                                                                                             <  1:1788519/150‑1 (MQ=255)
          tCCAGCGCACGCCCCTGCTCGTCGAGCACGAACGATTGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGCAAAACCAACGCCTGCCCATTGGACAGGCCCAGcacc                                                                                                                                         >  2:755359/1‑150 (MQ=255)
           ccAGCGCACGCCCCTGCTCGTCGAGCACGAACGATTGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGcacca                                                                                                                                        <  1:116578/150‑1 (MQ=255)
            cAGCGCACGCCCCTGCTCGTCGAGCACGAACGATTGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGcaccac                                                                                                                                       >  2:1726281/1‑150 (MQ=255)
              tctcACGCCCCTGCTCGTCGAGCACGAACGATTGCTTTCCGGACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCAccacca                                                                                                                                     <  2:1769898/147‑1 (MQ=255)
              gcgcACGCCCCTGCTCGTCGAGCACGACAGATTGCTCGTCATACATCGAGTCGATGGCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCAccacca                                                                                                                                     <  1:540869/150‑1 (MQ=255)
                  aCGCCCCTGCTCGTCGAGCACGAACGATTGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAgaga                                                                                                                                 >  1:1769898/1‑147 (MQ=255)
                      cccTGCTCGTCGAGCACGAACGATTGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTg                                                                                                                             >  2:1399841/1‑150 (MQ=255)
                       ccTGCTCGTCGAGCACGAACGATTGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGc                                                                                                                            >  2:76575/1‑150 (MQ=255)
                        cTGCTCGTCGAGCACGAACGATTGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGcc                                                                                                                           >  2:1493535/1‑150 (MQ=255)
                             cgtcgAGCACGAACGATTGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCt                                                                                                                      <  1:755359/150‑1 (MQ=255)
                                cgAGCACGAACGATTGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGt                                                                                                                   <  1:2050172/150‑1 (MQ=255)
                                     acgaacgaTTGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGcc                                                                                                                                                <  2:762398/116‑1 (MQ=255)
                                     acgaacgaTTGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGcc                                                                                                                                                >  1:762398/1‑116 (MQ=255)
                                           gaTTGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATg                                                                                                        >  1:679505/1‑150 (MQ=255)
                                             ttGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGa                                                                                                      <  1:1493535/150‑1 (MQ=255)
                                             ttGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGa                                                                                                      >  2:2283772/1‑150 (MQ=255)
                                              tGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGAc                                                                                                     <  1:76575/150‑1 (MQ=255)
                                                 tcgacGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGAt                                                                                                         <  2:2305507/139‑1 (MQ=255)
                                                 tCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTCCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTg                                                                                                  >  1:2305604/1‑150 (MQ=255)
                                                 tCGACGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGAt                                                                                                         >  1:2305507/1‑143 (MQ=255)
                                                   gCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCt                                                                                                                      <  2:1410766/128‑1 (MQ=255)
                                                   gCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCt                                                                                                                      >  1:1410766/1‑128 (MQ=255)
                                                      gTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCa                                                                                                                 >  1:678677/1‑130 (MQ=255)
                                                      gTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCa                                                                                                                 <  2:678677/130‑1 (MQ=255)
                                                       tACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTgg                                                                                            <  1:1726281/150‑1 (MQ=255)
                                                         cGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTATGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGgt                                                                                          <  1:521108/150‑1 (MQ=255)
                                                         cGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGgt                                                                                          <  1:1399841/150‑1 (MQ=255)
                                                               gTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGc                                                                                    <  1:701954/150‑1 (MQ=255)
                                                                  gATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGt                                                                                                                   >  1:835991/1‑116 (MQ=255)
                                                                  gATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGt                                                                                                                   <  2:835991/116‑1 (MQ=255)
                                                                  gATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGa                                                                                                 >  2:1913173/1‑134 (MQ=255)
                                                                  gATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGa                                                                                                 <  1:1913173/134‑1 (MQ=255)
                                                                           gCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCACAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAg                                                                        >  2:232508/1‑150 (MQ=255)
                                                                            cTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGt                                                                       >  1:1942129/1‑150 (MQ=255)
                                                                            cTGATGGTTTTCTTGTTGTCCGGGCAGGTGATCTTGTAGGGGTGCTGGAAAACCAGCGCCTGGCTATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATAGACGTGACCTGGGTGTCTGCCGGCACGGGTAGt                                                                       >  2:400820/1‑150 (MQ=255)
                                                                               aTGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTgg                                                                                            <  1:2032129/126‑1 (MQ=255)
                                                                               aTGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTgg                                                                                            >  2:2032129/1‑126 (MQ=255)
                                                                                tGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGa                                                                   >  2:1012227/1‑150 (MQ=255)
                                                                                 ggTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGcc                                                                                                                           <  2:1148169/93‑1 (MQ=255)
                                                                                 ggTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGcc                                                                                                                           >  1:1148169/1‑93 (MQ=255)
                                                                                   ttttCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGcc                                                                                   >  2:2094652/1‑131 (MQ=255)
                                                                                   ttttCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGcc                                                                                   <  1:2094652/131‑1 (MQ=255)
                                                                                    tttCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACg                                                               >  1:196630/1‑150 (MQ=255)
                                                                                     ttCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGcaccac                                                                                                                                       >  2:633906/1‑77 (MQ=255)
                                                                                     ttCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGcaccac                                                                                                                                       <  1:633906/77‑1 (MQ=255)
                                                                                         tgttgtCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGa                                                          <  1:988043/150‑1 (MQ=255)
                                                                                          gttgtCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAg                                                         >  2:1379951/1‑150 (MQ=255)
                                                                                           ttgtCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGt                                                        >  1:59844/1‑150 (MQ=255)
                                                                                             gtCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTc                                                      <  2:1942129/150‑1 (MQ=255)
                                                                                              tCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCg                                                     >  1:2056320/1‑150 (MQ=255)
                                                                                                cGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGACATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCGtt                                                   >  1:1758036/1‑150 (MQ=255)
                                                                                                 gggTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCGTTa                                                  >  1:1556458/1‑150 (MQ=255)
                                                                                                       gTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCGTTAc                                                 >  2:1702914/1‑145 (MQ=255)
                                                                                                       gTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCGTTAc                                                 <  1:1702914/145‑1 (MQ=255)
                                                                                                          aTCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCGTTACCGGTCTTg                                         <  1:232508/150‑1 (MQ=255)
                                                                                                             ttGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCGTTACCGGTCt                                           >  2:509799/1‑145 (MQ=255)
                                                                                                             ttGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCGTTACCGGTCTTGGTa                                      <  2:1811274/150‑1 (MQ=255)
                                                                                                             ctGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCGTTACCGGTCt                                           <  1:509799/144‑1 (MQ=255)
                                                                                                              tGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCGTTACCGGTCTTGGTAt                                     <  2:2056320/150‑1 (MQ=255)
                                                                                                                      tgCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCGTTACCGGTCTTGATATCGAAGAAc                             >  2:364050/1‑150 (MQ=255)
                                                                                                                        cTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCGTTACCGGTCTTGGTATCGAAGAACaa                           >  1:2396579/1‑150 (MQ=255)
                                                                                                                        cTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCGTTACCGGTCTTGGTATCGAAGAACaa                           <  1:2283772/150‑1 (MQ=255)
                                                                                                                         tGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGc                                                                                    >  2:1806892/1‑92 (MQ=255)
                                                                                                                         tGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGc                                                                                    <  1:1806892/92‑1 (MQ=255)
                                                                                                                         tGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCGTTACCGGTCtt                                          <  1:1774124/134‑1 (MQ=255)
                                                                                                                         tGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCGTTACCGGTCtt                                          >  2:1774124/1‑134 (MQ=255)
                                                                                                                         tGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCGTTACCGGTCTTGGTATCGAAGAACAAt                          >  1:511801/1‑150 (MQ=255)
                                                                                                                              aaCCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCGTTACCGGTCTTGGTATCGAAGAACAATGCCTg                     >  2:1963215/1‑150 (MQ=255)
                                                                                                                                 cAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCGTTACCGGTCTTGGTATCGAAGAACAATGCCTGGc                   >  2:2270812/1‑149 (MQ=255)
                                                                                                                                 cAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCGTTACCGGTCTTGGTATCGAAGAACAATGCCTGGc                   <  1:2270812/149‑1 (MQ=255)
                                                                                                                                  aGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACgg                                                                            <  2:170717/91‑1 (MQ=255)
                                                                                                                                  aGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACgg                                                                            >  1:170717/1‑91 (MQ=255)
                                                                                                                                   gggcCTGTCCAATGGCCAGGCGCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCGTTACCGGTCTTGGTATCGAAGAACAATGCCTGGCCCt                >  2:701728/3‑150 (MQ=255)
                                                                                                                                     gcCTGTCCAATGGCCAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCGTTACCGGTCTTGGTATCGAAGAACAATGCCTGGCCCTTg              >  1:466451/1‑150 (MQ=255)
                                                                                                                                      ccTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCGTTACCGGTCTTGGTATCGAAGAACAATGCCTAGCCCTTGt             <  2:511801/150‑1 (MQ=255)
                                                                                                                                       cTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCGTTACCGGTCTTGGTATCGAAGAACAATGCCTGGCCCTTGTc            >  1:834671/1‑150 (MQ=255)
                                                                                                                                       cTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCGTTACCGGTCTTGGTATCGAAGAACAATGCCTGGCCCTTGTc            >  1:800546/1‑150 (MQ=255)
                                                                                                                                          gCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCGTTACCGGTCTTGGTATCGAAGAACAATGCCTGGCCCTTGTCCGa         <  2:1758036/150‑1 (MQ=255)
                                                                                                                                            cATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCGTTACCGGTCTTGGTATCGAAGAACAATGCCTGGCCCTTGTCCGaaa       <  2:1556458/150‑1 (MQ=255)
                                                                                                                                            cATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCGTTACCGGTCTTGGTATCGAAGAACAATGCCTGGCCCTTGTCCGaaa       >  2:1181679/1‑150 (MQ=255)
                                                                                                                                             aTTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGc                                                                                    <  1:2201154/72‑1 (MQ=255)
                                                                                                                                             aTTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGc                                                                                    >  2:2201154/1‑72 (MQ=255)
                                                                                                                                               tGGACAGGTCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCGTTACCGGTCTTGGTATCGAAGAACAATGCCTGGCCCTTGTCCGAAacac    <  1:1012227/150‑1 (MQ=255)
                                                                                                                                               tGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCGTTACCGGTCTTGGTATCGAAGAACAATGCCTGGCCCTTGTCCGAAacac    <  1:400820/150‑1 (MQ=255)
                                                                                                                                                 gACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCGTTACCGGTCTTGGTATCGAAGAACAATGCCTGGCCCTTGTCCGAAACACGc  <  2:800546/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GCTGACGTGGTCCAGCGCACGCCCCTGCTCGTCGAGCACGAACGATTGCTCGCCGTACGGGTAGTCGATGCCAGGGCTGATGGTTTTCTTGTTGTCCGGGTAGGTGATCTTGTAGGTGTGCTGGAAAACCAGCGCCTGGCCATTGGACAGGCCCAGCACCACCAGCGCGCTGCCCGGCTGGTCAGTACTGATGGACGTGACCTGGGTGTCTGCCGGCACGGGTAGTTCGACACGCTTGAGTTCGTTACCGGTCTTGGTATCGAAGAACAATGCCTGGCCCTTGTCCGAAACACGC  >  NC_002947/6076009‑6076303

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: