Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,056,697 C→G 8.0% intergenic (‑207/+101) PP_2666 ← / ← PP_2667 hypothetical protein/ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,056,6970CG 8.0% 445.1 / 22.7 138intergenic (‑207/+101)PP_2666/PP_2667hypothetical protein/ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base C (20/107);  new base G (9/2);  total (29/109)
Fisher's exact test for biased strand distribution p-value = 1.06e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

Only 100 of 146 total aligned reads displayed.

GACGATAGACAGTAAATGCAGTACTAAAGTAGTAGATGGCCGGGCGCCTGAAGGCATGGAGACGCCCTAGGGAGCAACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTCAGCGCGGCTTGCGCAGCGCTGCATGCTGGGGGTTGAAGGTGAGGATC  >  NC_002947/3056550‑3056846
                                                                                                                                                   |                                                                                                                                                     
gacgaTAGACAGTAAATGCAGTACTAAAGTAGTAGATGGCCGGGCGCCTGAAGGCATGGAGACGCCCTAGGGAGCAACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGGCGGGCGGGGCGGCAGGGgcgcgccc                                                                                                                                                     >  2:71029/1‑148 (MQ=37)
gacgaTAGACAGTAAATGCAGTACTAAAGTAGTAGATGGCCGGGCGCCTGAAGGCATGGAGACGCCCTAGGGAGCAACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAg                                                                                                                                                     <  2:158890/150‑1 (MQ=255)
gacgaTAGAAAGTAAATGCAGTACTAAAGTAGTAGATGGCCGGGCGCCTGAAGGCATGGAGACGCCCTAGGGAGCAACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCat                                                                                                                                                     <  1:4952324/150‑2 (MQ=255)
      aGACAGTAAATGCAGTACTAAAGTAGTAGATGGCCGGGCGCCTGAAGGCATGGAGACGCCCTAGGGAGCAACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGccc                                                                                                                                               <  1:3412630/150‑1 (MQ=255)
      aGACAGTAAAGGCAGTACTAACGTAGTAGAGGGCCGGGCGCCTGAAGGCATGGAGACGCCCTAGGGAGCAACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGccc                                                                                                                                               <  1:238400/150‑1 (MQ=255)
       gACAGTAAATGCAGTACTAAAGTAGTAGATGGCCGGGCGCCTGAAGGCATGGAGACGCCCTAGGGAGCAACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCTGCGCGGccct                                                                                                                                              <  1:7060559/150‑2 (MQ=255)
        aCAGTAAATGCAGTACTAAAGTAGTAGATGGCCGGGCGCCTGAAGGCATGGAGACGCCCTAGGGAGCAACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGGCTTGCGGCGCGGAATGGCCGCGCCGCGGCCCCa                                                                                                                                             >  2:8189364/1‑150 (MQ=38)
        aCAGTAAATGCAGTACTAAAGTAGTAGATGGCCGGGCGCCTGAAGGCATGGAGACGCCCTAGGGAGCAACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCa                                                                                                                                             <  2:2664842/150‑1 (MQ=255)
        aCAGTAAATGCAGTACTAAAGTAGTAGATGGCCGGGCGCCTGAAGGCATGGAGACGCCCTAGGGAGCAACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCa                                                                                                                                             <  2:7321948/150‑1 (MQ=255)
        aCAGTAAATGCAGTACTAAAGTAGTAGATGGCCGGGCGCCTGAAGGCATGGAGACGCCCTAGGGAGCAACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCa                                                                                                                                             <  2:7322236/150‑1 (MQ=255)
        aCAGTAAATGCAGTACTAAAGTAGTAGATGGCCGGGCGCCTGAAGGCATGGAGACGCCCTAGGGAGCAACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCa                                                                                                                                             <  1:2943825/150‑1 (MQ=255)
           gTAAATGCAGTACTAAAGTAGTAGATGGCCGGGCGCCTGAAGGCATGGAGACGCCCTAGGGAGCAACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCTGCGCGGccctttcg                                                                                                                                          <  1:7260093/150‑6 (MQ=255)
                tGCAGTACTAAAGTAGTAGATGGCCGGGCGCCTGAAGGCGTGGAGACGCCCTAGGGAGCAACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCg                                                                                                                                     <  2:4777267/150‑1 (MQ=255)
                tGCAGTACTAAAGTAGTAGATGGCCGGGCGCCTGAAGGCGTGGAGACGCCCTAGGGAGCAACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCg                                                                                                                                     <  2:8264827/150‑1 (MQ=255)
                tGCAGTACTAAAGTAGTAGATGGCCGGGCGCCTGAAGGCGTGGAGACGCCCTAGGGAGCAACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCg                                                                                                                                     <  2:8264699/150‑1 (MQ=255)
                tGCAGTACTAAAGTAGTAGATGGCCGGGCGCCTGAAGGCATGGAGACGCCCTAGGGAGCAACTGTCTTGCTCATTTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCg                                                                                                                                     <  2:6181050/150‑1 (MQ=255)
                   aGTACTAAAGTAGTAGATGGCCGGGCGCCTGAAGGCATGGAGACGCCCTAGGGAGCAACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGca                                                                                                                                  <  2:816223/150‑1 (MQ=255)
                    gTACTAAAGTAGTAGATGGCCGGGCGCCTGAAGGCATGGAGACGCCCTAGGGAGCAACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGcat                                                                                                                                 <  2:5665315/150‑1 (MQ=255)
                     tACTAAAGTAGTAGATGGCCGGGCGCCTGAAGGCATGGAGACGCCCTAGGGAGCAACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGcatc                                                                                                                                >  2:3423867/1‑150 (MQ=255)
                       cTAAAGTAGTAGATGGCCGGGCGCCTGAAGGCATGGAGACGCCCTAGGGAGCAACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGcatcat                                                                                                                              <  2:7904610/150‑1 (MQ=255)
                         aaaGTAGTAGATGGCCGGGCGCCTGAAGGCATGGAGACGCCCTAGGGAGCAACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGc                                                                                                                            <  1:167858/150‑1 (MQ=255)
                           agtagtAGATGGCCGGGCGCCTGAAGGCATGGAGACGCCCTAGGGAGCAACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCAt                                                                                                                          <  1:4682250/150‑1 (MQ=255)
                            gtagtaGATGGCCGGGCGCCTGAAGGCATGGAGACGCCCTAGGGAGCAACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATa                                                                                                                         <  2:4799038/150‑1 (MQ=255)
                            gtagtaGATGGCCGGGCGCCTGAAGGCATGGAGACGCCCTAGGGAGCAACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATa                                                                                                                         <  1:7445164/150‑1 (MQ=255)
                            gtagtaGATGGCCGGGCGCCTGAAGGCATGGAGACGCCCTAGGGAGCAACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATa                                                                                                                         <  1:7445111/150‑1 (MQ=255)
                             tagtagATGGCCGGGCGCCTGAAGGCATGGAGACGCCCTAGGGAGCAACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAg                                                                                                                        <  2:5713484/150‑1 (MQ=255)
                              agtagATGGCCGGGCGCCTGAAGGCATGGAGACGCCCTAGGGAGCAACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGa                                                                                                                       <  2:8555381/150‑1 (MQ=255)
                               gtagATGGCCGGGCGCCTGAAGGCATGGAGACGCCCTAGGGAGCAACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGAt                                                                                                                      <  2:4890278/150‑1 (MQ=255)
                               gtagATGGCCGGGCGCCTGAAGGCATGGAGACGCCCTAGGGAGCAACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGAt                                                                                                                      <  2:850296/150‑1 (MQ=255)
                                   aTGGCCGGGCGCCTGAAGGCAAGGAGACGCCCTAGGGAGCAACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCt                                                                                                                  <  2:5390960/150‑1 (MQ=255)
                                      gCCGTGCGCCTGAAGGCATGGAGACGCCCTAGGGAGCATCTGTCTTGCTCAATTTCTGAAATCTGGCGCGATCCCTGTGGGAGGGGCCTTGAGTCGCGATAGGGCCGCGCTGCGGCCCCAGGATCTCGGCATCATGCTTAGATTGCTggg                                                                                                               <  2:3958174/150‑1 (MQ=38)
                                      gCCGGGCGCCTGAAGGCATGGAGACGCCCTAGGGAGCAACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTggg                                                                                                               <  2:7159885/150‑1 (MQ=255)
                                      gCCGGGCGCCTGAAGGCATGGAGACGCCCTAGGGAGCAACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTggg                                                                                                               <  1:4758481/150‑1 (MQ=255)
                                       ccGGGCGCCTGAAGGCATGGAGACGCCCTAGGGAGCAACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTgggg                                                                                                              <  2:1367019/150‑1 (MQ=255)
                                        cGGGCGCCTGAAGGCATGGAGACGCCCTAGGGAGCAACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGc                                                                                                             <  1:5758174/150‑1 (MQ=255)
                                                tGAAGGCATGGAGACGCCCTAGGGAGCAACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGc                                                                                                     <  1:7170230/150‑1 (MQ=255)
                                                   aGGCATGGAGACGCCCTAGGGAGCAACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGc                                                                                                  <  2:5155203/150‑1 (MQ=255)
                                                         ggAGACGCCCTAGGGAGCAACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATc                                                                                            <  1:3423867/150‑1 (MQ=255)
                                                           agaCGCCCTAGGGAGCAACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATcgc                                                                                          <  1:8725451/150‑1 (MQ=255)
                                                             aCGCCCTAGGGAGCAACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGa                                                                                        <  1:1621947/150‑1 (MQ=255)
                                                             aCGCCCTAGGGAGCAACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGa                                                                                        <  2:8043020/150‑1 (MQ=255)
                                                                cccTAGGGAGCAACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGacac                                                                                     <  2:4550745/150‑1 (MQ=255)
                                                                    aGGGAGCAACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAgg                                                                                 <  1:2698636/150‑1 (MQ=255)
                                                                    aGGGAGCAACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCACAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAgg                                                                                 <  2:1503041/150‑1 (MQ=255)
                                                                          cAACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGACGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTc                                                                           <  2:4178444/150‑1 (MQ=255)
                                                                               gTCGTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCcaca                                                                      <  1:5192224/150‑1 (MQ=255)
                                                                                tCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGcatcat                                                                                                                              >  2:2707848/1‑93 (MQ=255)
                                                                                tCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGcatcat                                                                                                                              <  1:2707848/93‑1 (MQ=255)
                                                                                tCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCacac                                                                     <  2:1818410/150‑1 (MQ=255)
                                                                                  ttGCTCAATTTCTGAACGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACgg                                                                   <  2:8004108/150‑1 (MQ=255)
                                                                                  ttGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACgg                                                                   <  1:8283651/150‑1 (MQ=255)
                                                                                             cTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCa                                                        <  2:3897782/150‑1 (MQ=255)
                                                                                                 aaGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGtt                                                    <  1:7374006/150‑1 (MQ=255)
                                                                                                     tGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTCa                                                <  2:7814883/150‑1 (MQ=255)
                                                                                                     tGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCATTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTCa                                                <  2:6303676/150‑1 (MQ=255)
                                                                                                                  gtgGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGc                                                                                                          <  1:4002415/79‑1 (MQ=255)
                                                                                                                           gCCTTGCGTCGCTAAAGGGCCCCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTCAGCGCGGCTTGCGCAGCGCTGCa                          <  1:2146212/150‑1 (MQ=255)
                                                                                                                            ccTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTCAGCGCGGCTTGTGCAGCGCTGCAt                         >  2:5401322/1‑150 (MQ=255)
                                                                                                                            ccTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTCAGCGCGGCTTGCGCAGCGCTGCAt                         >  1:4735858/1‑150 (MQ=255)
                                                                                                                             cTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTCAGCGCGGCTTGCGCAGCGCTGCATg                        >  1:1211502/1‑150 (MQ=255)
                                                                                                                              ttGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTCAGCGCGGCTTGCGCAGCGCTGCATGc                       >  1:1851050/1‑150 (MQ=255)
                                                                                                                                 cGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTCAGCGCGGCTTGCGCAGCGCTGCATGCTgg                    <  2:1578463/150‑1 (MQ=255)
                                                                                                                                 cGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTCAGCGCGGCTTGCGCAGCGCTGCATGCTgg                    <  1:5916664/150‑1 (MQ=255)
                                                                                                                                 cGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATGGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTCAGCGCGGCTTGCGCAGCGCTGCATGCTgg                    <  1:6034677/150‑1 (MQ=255)
                                                                                                                                  gTCGCGAACGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATTTATTGCGGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTCAGCGCGGCTTGCGCAGCGCTGCATGCTggg                   <  1:235858/150‑1 (MQ=255)
                                                                                                                                  gTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACACGGCCGCTCCCACACGGGTTACACGCCAGGTTTTCAGCGCGGCTTGCGCAGCGCTGCATGCTggg                   <  2:2410398/150‑1 (MQ=255)
                                                                                                                                   tCGCGAATGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGAATAGTTTGCTGGGGCCGCTTTGCGGCCAACTCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTCAGCGCGGCTTGCGCAGCGCTGCATGCTgggg                  <  1:2281433/150‑1 (MQ=39)
                                                                                                                                   tCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTCAGCGCGGCTTGCGCAGCGCTGCATGCTgggg                  <  1:1256479/150‑1 (MQ=255)
                                                                                                                                   tCGCGAAAGGGCCGCGCAGCGGCCCCAGGAACTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTCAGCGCGGCTTGCGCAGCGCTGCATGCTgggg                  <  1:6900655/150‑1 (MQ=255)
                                                                                                                                     gcgAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTc                                                 >  2:7663657/1‑117 (MQ=255)
                                                                                                                                     gcgAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTc                                                 <  1:7663657/117‑1 (MQ=255)
                                                                                                                                       gAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCACTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTCAGCGCGGCTTGCGCAGCGCTGCATGCTGGGGGTTg              >  1:4291042/1‑150 (MQ=255)
                                                                                                                                       cctggggCCGCTGCGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACAAGCCAGGTTTTAAGCGCGGCTTGCGCAGCGCTGCATGCTGGGGGTTg              >  1:2674914/5‑150 (MQ=38)
                                                                                                                                        aaaGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTCAGCGCGGCTTGCGCAGCGCTGCATGCTGGGGGTTGa             <  2:4248591/150‑1 (MQ=255)
                                                                                                                                         tggggCCGCTGCGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTCAGCGCGGCTTGCGCAGCGCTGCATGCTGGGGGTTGac            >  2:4877242/3‑149 (MQ=255)
                                                                                                                                         tggggCCGCTGCGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTCAGCGCGGCTTGCGCAGCGCTGCATGCTGGGGGTTGaa            >  2:4672477/3‑150 (MQ=255)
                                                                                                                                         tggggCCGCTGCGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTCAGCGCGGCTTGCGCAGCGCTGCATGCTGGGGGTTGaa            >  1:611842/3‑150 (MQ=255)
                                                                                                                                         aaGGGCCGCGCAGCGGCCCCTGGATCTCGGCATCATGCATCGATTGCTGGGGCCGCTTTGCGGCCCAAGCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTCAGCGCGGCTTGCGCAGCGCTGCATGCTGGGGGTTGaa            <  1:809417/150‑1 (MQ=255)
                                                                                                                                         aaGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCTGGTTTTCAGCGCGGCTTGCGCAGCGCTGCATGCTGGGGGTTGaa            <  2:8077949/150‑1 (MQ=255)
                                                                                                                                         aaGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTCAGCGCGGCTTGCGCAGCGCTGCATGCTGGGGGTTGaa            <  1:5851756/150‑1 (MQ=255)
                                                                                                                                         aaGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTCAGCGCGGCTTGCGCAGCGCTGCATGCTGGGGGTTGaa            <  2:2648290/150‑1 (MQ=255)
                                                                                                                                          ggggCCGCTGCGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTCAGCGCGGCTTGCGCAGCGCTGCATGCTGGGGGTTGAAg           >  2:1556657/2‑150 (MQ=255)
                                                                                                                                          ggggCCGCTGCGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTCAGCGCGGCTTGCGCAGCGCTGCATGCTGGGGGTTGAAg           >  2:1557214/2‑150 (MQ=255)
                                                                                                                                             gCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATGGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTCAGCGCGGCTTGCGCAGCGCTGCATGCTGGGGGTTGAAGGTg        <  2:2427099/150‑1 (MQ=255)
                                                                                                                                              ccGCGCAGCGGCCCCAGGATCTCGGCATCTTGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTCAGCGCGGCTTGCGCAGCGCTGCATGCTGGGGGTTTAAGGTGa       <  2:1892830/150‑1 (MQ=255)
                                                                                                                                              ccGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTCAGCGCGGCTTGCGCAGCGCTGCATGCTGGGGGTTGAAGGTGa       <  1:5212144/150‑1 (MQ=255)
                                                                                                                                              ccGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTCAGCGCGGCTTGCGCAGCGCTGCATGCTGGGGGTTGAAGGTGa       <  1:5400336/150‑1 (MQ=255)
                                                                                                                                              ccGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTCAGCGCGGCTTGCGCAGCGCTGCATGCTGGGGGTTGAAGGTGa       <  1:6474378/150‑1 (MQ=255)
                                                                                                                                              ccGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTCAGCGCGGCTTGCGCAGCGCTGCATGCTGGGGGTTGAAGGTGa       <  1:7098197/150‑1 (MQ=255)
                                                                                                                                              ccGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTCAGCGCGGCTTGCGCAGCGCTGCATGCTGGGGGTTGAAGGTGa       <  2:2759088/150‑1 (MQ=255)
                                                                                                                                              ccGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTCAGCGCGGCTTGCGCAGCGCTGCATGCTGGGGGTTGAAGGTGa       <  1:110933/150‑1 (MQ=255)
                                                                                                                                              ccGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTCAGCGCGGCTTGCGCAGCGCTGCATGCTGGAGGTTGAAGGTGa       <  2:7683305/150‑1 (MQ=255)
                                                                                                                                              ccGCGCAGCGCCCCCCGGATCTCGGCATCATGCATAGATTGCGGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTCAGCGCGGCTTGCGCAGCGCTGCATGCTGGGGGTTGAAGGTGa       <  1:71029/150‑1 (MQ=255)
                                                                                                                                               cgcgCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTCAGCGCGGCTTGCGCAGCGCTGCATGCTGGGGGTTGAAGGTGAg      <  1:6321087/150‑1 (MQ=255)
                                                                                                                                               cgcgCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTCAGCGCGGCTTGCGCAGCGCTGCATGCTGGGGGTTGAAGGTGAg      >  1:1758423/1‑150 (MQ=255)
                                                                                                                                               cgcgCAGCGCCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTCAGCGCGGCTTGCGCAGCGCTGCATGCTGGGGGTTGAAGGTGAg      <  2:1137044/150‑1 (MQ=255)
                                                                                                                                                gcgcAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTCAGCGCGGCTTGCGCAGCGCTGCATGCTGGGGGTTGAAGGTGAgg     <  1:3615540/150‑1 (MQ=255)
                                                                                                                                                gcgcAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTCAGCGCGGCTTGCGCAGCGCTGCATGCTGGGGGTTGAAGGTGAgg     >  1:3067795/1‑150 (MQ=255)
                                                                                                                                                gcgcAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTCAGCGCGGCTTGCGCAGCGCTGCATGCTGGGGGTTGAAGGTGAgg     <  2:8811829/150‑1 (MQ=255)
                                                                                                                                                  gcAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTCAGCGCGGCTTGCGCAGCGCTGCATGCTGGGGGTTGAAGGTGAGGAt   >  2:8696433/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GACGATAGACAGTAAATGCAGTACTAAAGTAGTAGATGGCCGGGCGCCTGAAGGCATGGAGACGCCCTAGGGAGCAACTGTCTTGCTCAATTTCTGAAAGCTGGCGCGATCCCTGTGGGAGCGGCCTTGCGTCGCGAAAGGGCCGCGCAGCGGCCCCAGGATCTCGGCATCATGCATAGATTGCTGGGGCCGCTTTGCGGCCCAATCGCGACACAAGGCCGCTCCCACACGGGTTACACGCCAGGTTTTCAGCGCGGCTTGCGCAGCGCTGCATGCTGGGGGTTGAAGGTGAGGATC  >  NC_002947/3056550‑3056846

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