Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,786,024 T→A 9.8% intergenic (+40/‑146) rpsB → / → tsf 30S ribosomal protein S2/elongation factor Ts

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,786,0240TA 9.8% 201.4 / 10.9 71intergenic (+40/‑146)rpsB/tsf30S ribosomal protein S2/elongation factor Ts
Reads supporting (aligned to +/- strand):  ref base T (32/32);  new base A (4/3);  total (36/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

GAGCTGTACATGACTTCGATGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGTTGAAAAACTTTCCAAGAGGATTTTGAAAT  >  NC_002947/1785877‑1786171
                                                                                                                                                   |                                                                                                                                                   
gAGCTGTACATGACTTCGATGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCttt                                                                                                                                                     <  4:166357/148‑1 (MQ=255)
   cTGTACATGACTTCGATGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGccc                                                                                                                                                 >  6:23421/1‑149 (MQ=255)
     gTACATGACTTCGATGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCCCGCAAAAAGGGGGCTTTGccccc                                                                                                                                               >  7:85641/1‑149 (MQ=255)
     gTACATGACTTCGATGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGccccc                                                                                                                                               >  4:370871/1‑149 (MQ=255)
        cATGACTTCGATGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCttt                                                                                                                                            >  1:133450/1‑149 (MQ=255)
        cATGACTTCGATGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCAAAGCCCCCttt                                                                                                                                            <  2:113446/149‑1 (MQ=255)
          tGACTTCGATGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAATCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCAAAGCCCCCttttt                                                                                                                                          <  5:273849/149‑1 (MQ=255)
          tGACTTCGATGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCttttt                                                                                                                                          <  5:155097/149‑1 (MQ=255)
                cGATGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCAAAGCCCCCTTTTTGcgt                                                                                                                                      <  7:100113/147‑3 (MQ=255)
                 gATGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCAcc                                                                                                                                    <  5:340699/148‑1 (MQ=255)
                  aTGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCtt                                                                                                                                  >  5:323917/1‑149 (MQ=255)
                  aTGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCtt                                                                                                                                  <  4:9237/149‑1 (MQ=255)
                   tGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTg                                                                                                                                 <  6:339405/149‑1 (MQ=255)
                      cTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGaaa                                                                                                                              <  8:249434/149‑1 (MQ=255)
                        gACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATc                                                                                                                            >  3:264759/1‑149 (MQ=255)
                             aGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCctgct                                                                                                                       <  4:225562/149‑1 (MQ=255)
                             aGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCctgct                                                                                                                       >  2:192609/1‑149 (MQ=255)
                                 aTCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCa                                                                                                                   <  5:23421/149‑1 (MQ=255)
                                    cgcgGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCa                                                                                                                <  7:154286/149‑1 (MQ=255)
                                      cgGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGc                                                                                                                                                   <  2:60926/112‑1 (MQ=255)
                                      cgGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGc                                                                                                                                                   >  1:60926/1‑112 (MQ=255)
                                          cGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGcc                                                                                                          >  1:330581/1‑149 (MQ=255)
                                           gCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGccc                                                                                                         >  1:373618/1‑149 (MQ=255)
                                            caacaaCGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGcccc                                                                                                        <  7:342348/149‑1 (MQ=255)
                                             aacaacGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCg                                                                                                       >  1:29586/1‑149 (MQ=255)
                                              acaacGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTAGAAATCCTGCTGTCAGCATCGGCCCCGc                                                                                                      <  8:1326/149‑1 (MQ=255)
                                                    ttGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATg                                                                                                >  6:187361/1‑149 (MQ=255)
                                                       ccGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGc                                                                                             >  7:337671/1‑149 (MQ=255)
                                                           cggcACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCctgc                                                                                         <  2:330581/149‑1 (MQ=255)
                                                            ggcACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCctgct                                                                                        <  4:56878/149‑1 (MQ=255)
                                                            ggcACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCctgct                                                                                        <  6:302906/149‑1 (MQ=255)
                                                            ggcACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCctgct                                                                                        >  7:184557/1‑149 (MQ=255)
                                                            ggcACCGAAGTGTACGCTGAAGAAGCGGCTGCAGCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCGCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCctgct                                                                                        <  4:148989/149‑1 (MQ=255)
                                                                 cGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCtgaaat                                                                                   <  1:192609/149‑1 (MQ=255)
                                                                          cGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCtgaaa                                                                                    >  8:138574/1‑139 (MQ=255)
                                                                          cGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCtgaaa                                                                                    <  7:138574/139‑1 (MQ=255)
                                                                          cGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAg                                                                           <  7:131767/148‑1 (MQ=255)
                                                                              gaagaaGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGAt                                                                      >  7:123626/1‑149 (MQ=255)
                                                                                agaagCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGAttt                                                                    <  6:311188/149‑1 (MQ=255)
                                                                                agaagCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGAttt                                                                    <  2:29586/149‑1 (MQ=255)
                                                                                 gaagCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTg                                                                   <  8:184557/149‑1 (MQ=255)
                                                                                    gCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTg                                                                   >  7:315667/1‑146 (MQ=255)
                                                                                    gCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTg                                                                   >  5:67640/1‑146 (MQ=255)
                                                                                       gCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAAt                                                             >  5:51564/1‑149 (MQ=255)
                                                                                        cTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAAtt                                                            <  1:10177/149‑1 (MQ=255)
                                                                                           cACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGa                                                          <  5:187361/148‑1 (MQ=255)
                                                                                                gctgctgAGTAATAGGCGCTAGCGTCTACTTGGCCCGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCGCGGATCATGGGCCTGCTGAAATGAAAACAGCGGCGTTGCGGAATGGCACGcgc                                                    >  4:245713/1‑147 (MQ=255)
                                                                                                    ctgAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGc                                                                                             >  2:211241/1‑104 (MQ=255)
                                                                                                    ctgAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTCCTGTCAGCATCGGCCCCGCATCATGGGc                                                                                             <  1:211241/104‑1 (MQ=255)
                                                                                                             aGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCtgaaa                                                                                    >  2:359689/1‑104 (MQ=255)
                                                                                                             aGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCAGCATGGGCCTGCtgaaa                                                                                    <  1:359689/104‑1 (MQ=255)
                                                                                                               gcgcTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTg                                      >  3:140004/1‑148 (MQ=255)
                                                                                                                 gcTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGaa                                   <  6:323917/149‑1 (MQ=255)
                                                                                                                  cTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAAt                                  >  2:320599/1‑149 (MQ=255)
                                                                                                                   tAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATg                                 >  6:251614/1‑149 (MQ=255)
                                                                                                                     gCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATgg                                <  3:370871/148‑1 (MQ=255)
                                                                                                                     gCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGt                               <  2:133450/149‑1 (MQ=255)
                                                                                                                             ttGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGTTGAAAAAc                       >  1:362529/1‑149 (MQ=255)
                                                                                                                               ggCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGTTGAAAAACtt                     >  1:36865/1‑149 (MQ=255)
                                                                                                                                ggtggCAAAAAGGGGGCAAAGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCCCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGTTGAAAAACtt                     >  8:298021/5‑148 (MQ=255)
                                                                                                                                  aCGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATgg                                <  7:10530/135‑1 (MQ=255)
                                                                                                                                  aCGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATgg                                >  8:10530/1‑135 (MQ=255)
                                                                                                                                  aCGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGTTGAAAAACTTTc                   >  5:44293/1‑148 (MQ=255)
                                                                                                                                   ggCAAAAAGGGGGCAAAGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGTTGAAAAACTTTcc                  >  2:134207/2‑148 (MQ=255)
                                                                                                                                    gCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGTTGAAAAACTTTCCaa                <  2:373618/149‑1 (MQ=255)
                                                                                                                                      aaaaaGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGTTGAAAAACTTTCCaa                >  7:250279/1‑147 (MQ=255)
                                                                                                                                      aaaaaGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGTTGAAAAACTTTCCAaga              <  2:18926/149‑1 (MQ=255)
                                                                                                                                      aaaaaGGGGGCAAAGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGTTGAAAAACTTTCCaa                >  4:227359/1‑147 (MQ=255)
                                                                                                                                       aaaaGGGGGCAAAGCCCCCTTTTTGCCACCTTGAAATCCTGCAGTCAGCATCGGCCCCGCATCATGGGCATGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGACCGTGATTGACGGGTGGAATGGTTGAAAAACTTTCCAagag             >  5:267959/1‑149 (MQ=255)
                                                                                                                                        aaaGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAAGAGCGGATTTGCAGAATTGAACGCGCGTGAGTGACGGGTGGAATGGTTGAATAAGTTTCCATGAgg            >  6:40675/1‑149 (MQ=255)
                                                                                                                                             gggCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGTTGAAAAACTTTCCAAGAGGAtttt       <  8:315667/149‑1 (MQ=255)
                                                                                                                                                cTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGTTGAAAAACTTTCCAAGAGGATTTTGa     >  2:107065/1‑148 (MQ=255)
                                                                                                                                                   tGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGTTGAAAAACTTTCCAAGAGGATTTTGAAAt  <  6:67640/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GAGCTGTACATGACTTCGATGGCTGACGCAGTCATCCGCGGCCGCAACAACGTTGCCGGCGGCACCGAAGTTTACGCTGAAGAAGCGGCTGCACCTGCTGCTGAGTAATAGGCGCTAGCGTCTACTTGGCACGCAAAAAGGGGGCTTTGCCCCCTTTTTGCCACCTTGAAATCCTGCTGTCAGCATCGGCCCCGCATCATGGGCCTGCTGAAATGAAAACAGCGGATTTGCAGAATTGAACGCCCGTGATTGACGGGTGGAATGGTTGAAAAACTTTCCAAGAGGATTTTGAAAT  >  NC_002947/1785877‑1786171

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