Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,563,870 A→C 5.1% intergenic (+94/+28) PP_1372 → / ← pitB hypothetical protein/phosphate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,563,8700AC 5.1% 230.8 / 4.1 81intergenic (+94/+28)PP_1372/pitBhypothetical protein/phosphate transporter
Reads supporting (aligned to +/- strand):  ref base A (39/37);  new base C (2/2);  total (41/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 = 4.11e-01

GCGGCAGTTGGTGCGCGGGTTGATGGGGTCGTTGTTAGGGGGCAAGAAAAGGTAGTCTTCCTGCGCCGCGCCATACCCTTGGTAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGGACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCATCGAGGCAGGCAGGGTCAGGACCCATGCCAGCAGGATGGTTTTCACAGTGCCGCCTTGCA  >  NC_002947/1563722‑1564007
                                                                                                                                                     |                                                                                                                                          
gcggcAGTTGGTGCGCGGGTTGATGGGGTCGTTGTTAGGGGGCAAGAAAAGGTAGTCTTCCTGCGCCGCGCCATACCCTCCATAAGATCGGAAGAGCACACGTATGAACAAGGTGGCGTAGGCATAACGGCTATAAAAAAGGCGCCCGAg                                                                                                                                           >  1:878463/1‑150 (MQ=255)
 cggcAGTTGGTGCGCGGGTTGATGGGGTCGTTGTTAGGGGGCAAGAAAAGGTAGTCTTCCTGCGCCGCGCCATATACTTGGTAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAgg                                                                                                                                          <  2:298366/150‑1 (MQ=255)
   gcAGTTGGTGCGCGGGTTGATGGGGTCGTTGTTAGGGGGCAAGAAAAGGTAGTCTTCCTGCGCCGCGCCATACCCTTGGTAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGgc                                                                                                                                        >  1:1268721/1‑150 (MQ=255)
      gTTGGTGCGCGGGTTGATGGGGTCGTTGTTAGGGGGCAAGAAAAGGTAGTCTTCCTGTGCCGCGCCATACCCTTGGTAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCTCGGGCGcc                                                                                                                                     <  2:1314957/150‑1 (MQ=255)
      gTTGGTGCGCGGGTTGATGGGGTCGTTGTTAGGGGGCAAGAAAAGGTAGTCTTCCTGCGCCGCGCCATACCCTTGGTCGAACTAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGcc                                                                                                                                     >  1:1157827/1‑150 (MQ=255)
                  gTTGATGGGGTCGTTGTTAGGGGGCAAGAAAAGGTAGTCTTCCTGCGCCGCGCCATACCCTTGGTAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGc                                                                                                                         <  2:1139485/150‑1 (MQ=255)
                  gTTGATGGGGTCGTTGTTAGGGGGCAAGAAAAGGTAGTCTTCCTGCGCCGCGCCATACCCTTGGTAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGc                                                                                                                         <  2:1092320/150‑1 (MQ=255)
                    tGATGGGGTCGTTGTTAGGGGGCAAGAAAAGGTAGTCTTCCTGCGCCGCGCCATACCCTTGGTAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCgg                                                                                                                       <  1:369444/150‑1 (MQ=255)
                           gTCGTTGTTAGGGGGCAAGAAAAGGTAGTCTTCCTGCGCCGCGCCATACCCTTGGTAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCttc                                                                                                                                  <  2:821936/132‑1 (MQ=255)
                           gTCGTTGTTAGGGGGCAAGAAAAGGTAGTCTTCCTGCGCCGCGCCATACCCTTGGTAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCttc                                                                                                                                  >  1:821936/1‑132 (MQ=255)
                           gTCGTTGTTAGGGGGCAAGAAAAGGTAGTCTTCCTGCGCCGCGCCATACCCTTGGTAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACt                                                                                                                >  2:1148022/1‑150 (MQ=255)
                             cgcttgtAGGGGGCGAAAAAAGGTAGTCGTACGGCGCCGCGCCATACCCTTGGTAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGTTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCa                                                                                                             <  2:112680/144‑1 (MQ=255)
                             cGTTGTTAGGGGGCAAGAAAAGGTAGTCTTCCTGCGCCGCGCCATACCCTTGGTAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCa                                                                                                             <  1:1216412/150‑1 (MQ=255)
                             cGTTGTTAGGGGGCAAGAAAAGGTAGTCTTCCTGCGCCGCGCCATACCCTTGGTAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCa                                                                                                             <  2:73201/150‑1 (MQ=255)
                             cGTTGTTAGGGGGCAAGAAAAGGTAGTCTTCCTGCGCCGCGCCATACCCTTGGTAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCa                                                                                                             <  1:48820/150‑1 (MQ=255)
                                  ttAGGGGGCGTGAAAAGGTAGTCTTCCTGCGCCGCGCCATACCCTTGGTAGCACAAGCCCGCTCCGGCAAAGAACACGGTGGCGTAGGCATAACGGCTAGGAAAAAGGCGCCCGAGGGCGCCGTCATTCATTGCGGTGTCACGCAGCCaa                                                                                                        <  2:257827/150‑1 (MQ=255)
                                  ttAGGGGGCAAGAAAAGGTAGTCTTCCTGCGCCGCGCCATACCCTTGGTAGCACAAGCCCGCTCCGGCAAAGAACATGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCaa                                                                                                        >  1:786857/1‑150 (MQ=255)
                                         gCAAGAAAAGGTAGTCTTCCTGCGCCGCGCCATACCCTTGGTAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCttt                                                                                                 <  1:1148022/150‑1 (MQ=255)
                                             gAAAAGGTAGTCTTCCTGCGCCGCGCCATACCCTTGGTAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCttt                                                                                                 <  1:622763/146‑1 (MQ=255)
                                             gAAAAGGTAGTCTTCCTGCGCCGCGCCATACCCTTGGTAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCttt                                                                                                 >  2:622763/1‑146 (MQ=255)
                                             gAAAAGGTAGTCTTCCTGCGCCGCGCCATACCCTTGGTAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGGAc                                                                                             >  1:223606/1‑150 (MQ=255)
                                             gAAAAGGTAGTCTTCCTGCGCCGCGCCATACCCTTGGTAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGGAc                                                                                             >  1:609938/1‑150 (MQ=255)
                                             gAAAAGGTAGTCTTCCTGCGCCGCGCCATACCCTTGGTAGCACAAGCCCGCTCCGGAAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGGAc                                                                                             <  2:297510/150‑1 (MQ=255)
                                               aaaGGTAGTCTTCCTGCGCCGCGCCATACCCTTGGTAGCACAAGCCCGCTCCGGCAAAGAACATGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGGACgc                                                                                           <  2:786857/150‑1 (MQ=255)
                                                    tAGTCTTCCTGCGCCGCGCCATACCCTTGGTAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGGACgccagcc                                                                                      >  2:853838/1‑150 (MQ=255)
                                                         ttCCTGCGCCGCGCCATACCCTTGGTAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGGACgccagcca                                                                                     >  1:140300/1‑146 (MQ=255)
                                                         ttCCTGCGCCGCGCCATACCCTTGGTAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGGACgccagcca                                                                                     <  2:140300/146‑1 (MQ=255)
                                                                gccgcgccATACCCTTGGTAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGgcgc                                                                                                                                      >  2:768768/1‑91 (MQ=255)
                                                                gccgcgccATACCCTTGGTAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGgcgc                                                                                                                                      <  1:768768/91‑1 (MQ=255)
                                                                  cgcgccATACCCTTGGTAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGGACGCCAGCCAGAACAggccggcc                                                                        >  1:763990/1‑150 (MQ=255)
                                                                      ccATACCCTTGGTAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGGACGCCAGCCAGAACAGGCCGGCCGCCa                                                                    >  1:1252533/1‑150 (MQ=255)
                                                                       cATACCCTTGGTAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGGACgccagc                                                                                       >  2:440412/1‑130 (MQ=255)
                                                                       cATACCCTTGGTAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGGACgccagc                                                                                       <  1:440412/130‑1 (MQ=255)
                                                                       cATACCCTTGGTAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGGACGCCAGCCAGAACAGGCCGGCCGCCAg                                                                   >  1:590932/1‑150 (MQ=255)
                                                                         tACCCTTGGTAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCGGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGGACGCCAGCCAGAACAGGCCGGCCGCCAGGc                                                                 >  2:984487/1‑150 (MQ=255)
                                                                            ccTTGGTAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGGACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCa                                                              >  2:414431/1‑150 (MQ=255)
                                                                            ccTTGGTAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGGACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCa                                                              >  2:335115/1‑150 (MQ=255)
                                                                              ttGGTAGCACAAGCCTTCGCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGGACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCATc                                                            >  1:49242/1‑150 (MQ=255)
                                                                              ttGGTAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGGACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCATc                                                            <  2:1268721/150‑1 (MQ=255)
                                                                              ttGGTAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTCGTCGTTAGTCAGAACAGGCCGGCCGCCAGGCCCATc                                                            <  1:984487/150‑1 (MQ=255)
                                                                              ttGGTAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTGCCGGCATTGCGGTGTCACTCAGCCAATGGCTTTGGACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCATc                                                            >  1:738038/1‑150 (MQ=255)
                                                                                  tAGCACAAGCCCGATCGGGCAAAAAATAACTTTGCGGATGCATAACGGCTATGCAAAAGGTGTCCGTAGAATGCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGGACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCATCGagg                                                        <  1:386112/150‑1 (MQ=255)
                                                                                    gCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAAt                                                                                                       >  1:682588/1‑101 (MQ=255)
                                                                                    gCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAAt                                                                                                       <  2:682588/101‑1 (MQ=255)
                                                                                    gCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGGACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCATCGaggca                                                      <  1:853838/150‑1 (MQ=255)
                                                                                    gCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCAGATGCAGAACGGCGAGGAAAAAAGCACACGAGGGCGCCCCAGTCACTTCCAAGGATCTCGTATGAAGGGCTTTGGACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCATCGaggca                                                      >  1:1391797/1‑150 (MQ=255)
                                                                                     cacaAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGGACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCATCGaggcag                                                     <  2:1391797/150‑1 (MQ=255)
                                                                                      acaAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGGACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCATCGaggcagg                                                    >  1:1163055/1‑150 (MQ=255)
                                                                                          gCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTATTTTGCATGTCACTCAGCCAATGGCTTTGGACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCATCGAGGcaggcagg                                                <  2:49242/150‑1 (MQ=255)
                                                                                               gCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCt                                                                                                                                    <  2:753057/63‑1 (MQ=255)
                                                                                               gCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCt                                                                                                                                    >  1:753057/1‑63 (MQ=255)
                                                                                                cTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGGACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCATCGAGGCAGGCAGGGTCAg                                           <  1:414431/150‑1 (MQ=255)
                                                                                                     gCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGGACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCATCGaggcagg                                                    <  2:1399470/136‑1 (MQ=255)
                                                                                                     gCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGGACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCATCGaggcagg                                                    >  1:1399470/1‑136 (MQ=255)
                                                                                                     gCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGGACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCATCGAGGCAGGCAGGGTCAGGa                                         >  1:353593/1‑147 (MQ=255)
                                                                                                     gCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGGACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCATCGAGGCAGGCAGGGTCAGGa                                         <  2:353593/147‑1 (MQ=255)
                                                                                                     gCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGGACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCATCGAGGCAGGCAGGGTCAGGAccc                                      >  2:351359/1‑150 (MQ=255)
                                                                                                      cAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGGACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCATCGAGGCAGGCAGGGTCAGGACCCa                                     >  2:403825/1‑150 (MQ=255)
                                                                                                       aaaGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCaa                                                                                                        <  2:1364718/82‑1 (MQ=255)
                                                                                                       aaaGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCaa                                                                                                        >  1:1364718/1‑82 (MQ=255)
                                                                                                       aaaGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGGACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCATCGAGGCAGGCAGGGTCAGGACCCAt                                    <  2:763990/150‑1 (MQ=255)
                                                                                                             cAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGGACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCATCGaggcaggc                                                   <  1:1310055/129‑1 (MQ=255)
                                                                                                             cAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGGACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCATCGaggcaggc                                                   >  2:1310055/1‑129 (MQ=255)
                                                                                                             cAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGGACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCATCGAGGCAGGCAGGGTCAGGACCCAt                                    >  1:179411/1‑144 (MQ=255)
                                                                                                             cAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGGACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCATCGAGGCAGGCAGGGTCAGGACCCAt                                    <  2:179411/144‑1 (MQ=255)
                                                                                                              aaGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGGACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCATCGAGGCAGGCAGGGTCAGGACCCATGCcagca                             <  2:223606/150‑1 (MQ=255)
                                                                                                              aaGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGGACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCATCGAGGCAGGCAGGGTCAGGACCCATGCcagca                             <  2:1252533/150‑1 (MQ=255)
                                                                                                                    gCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGGACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCATCGAGGCAGGCAGGGTCAGGACCCATGCCAGCAGGATgg                       >  2:806218/1‑150 (MQ=255)
                                                                                                                     cGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCGTTCATTCATTGCGGTGTCACTCATCCAATGGCTTTGGACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCATCGAGGCAAGCAGGGTCAGGACCCATGCCAGCAGGAtaga                      >  2:576264/1‑147 (MQ=255)
                                                                                                                      gTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGGACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCATCGAGGCAGGCAGGGTCAGGACCCATGCCAGCAGGATGGtt                     >  2:1025687/1‑150 (MQ=255)
                                                                                                                          gCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGGACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCATCGAGGCAGGCAGGGTCAGGACCCATGCCAGCAGGATGGTTTTca                 >  2:1092308/1‑150 (MQ=255)
                                                                                                                            aTAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGGACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCATCGAGGCAGGCAGGGTCAGGACCCATGCCAGCAGGATGGTTTTcaca               <  1:403825/150‑1 (MQ=255)
                                                                                                                            aTAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGGACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCATCGAGGCAGGCAGGGTCAGGACCCATGCCAGCAGGATGGTTTTcaca               <  2:590932/150‑1 (MQ=255)
                                                                                                                               aCGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGGACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCATCGAGGCAGGCAGGGTCAGGACCCATGCCAGCAGGATGGTTTTCACAGTg            >  1:860978/1‑150 (MQ=255)
                                                                                                                                cGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGGACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCATCGAGGCAGGCAGGGTCAGGACCCATGCCAGCAGGATGGTTTTCACAGTgc           <  2:860978/150‑1 (MQ=255)
                                                                                                                                 ggCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGGACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCATCGAGGCAGGCAGGGTCAGGACCCATGCCAGCAGGATGGTTTTCACAGTgcc          >  2:723549/1‑150 (MQ=255)
                                                                                                                                      tGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGGACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCATCGAGGCAGGCAGGGTCAGGACCCATGCCAGCAGGATGGTTTTCACAGTGCCGCCtt     <  2:878463/150‑1 (MQ=255)
                                                                                                                                       aatgaaGGCGCCCTCGGGCGCCTTTTTCATTGCGGTGTCACTCAGCCAATGGCTTTGGACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCATCGAGGCAGGCAGGGTCAGGACCCATGCCAGCAGGATGGTTTTCACAGTGCCGCCTTGc   >  2:1330712/5‑150 (MQ=255)
                                                                                                                                         aaaaGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGGACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCATCGAGGCAGGCAGGGTCAGGACCCATGCCAGCAGGATGGTTTTCACAGTGCCGCCTTGCa  >  1:338280/1‑150 (MQ=255)
                                                                                                                                            aGGCGCCCTCGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGGACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCATCGagg                                                        <  2:1211597/93‑1 (MQ=255)
                                                                                                                                            aGGCGCCCTCGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGGACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCATCGagg                                                        >  1:1211597/1‑93 (MQ=255)
                                                                                                                                                     |                                                                                                                                          
GCGGCAGTTGGTGCGCGGGTTGATGGGGTCGTTGTTAGGGGGCAAGAAAAGGTAGTCTTCCTGCGCCGCGCCATACCCTTGGTAGCACAAGCCCGCTCCGGCAAAGAACAAGGTGGCGTAGGCATAACGGCTATGAAAAAGGCGCCCGAGGGCGCCTTCATTCATTGCGGTGTCACTCAGCCAATGGCTTTGGACGCCAGCCAGAACAGGCCGGCCGCCAGGCCCATCGAGGCAGGCAGGGTCAGGACCCATGCCAGCAGGATGGTTTTCACAGTGCCGCCTTGCA  >  NC_002947/1563722‑1564007

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