Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,968,140 A→G 5.9% intergenic (‑36/+39) PP_4377 ← / ← fliC flagellin FlaG/flagellin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,968,1400AG 5.9% 249.9 / 4.5 85intergenic (‑36/+39)PP_4377/fliCflagellin FlaG/flagellin
Reads supporting (aligned to +/- strand):  ref base A (44/36);  new base G (2/3);  total (46/39)
Fisher's exact test for biased strand distribution p-value = 6.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.80e-01

TTGGTCGGCACGTCTGGTAGCTCAGCATTCGATCTGCTCGATGGCTGTGCACTCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCAGCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTTGCTGCAGGGTTTGCTGCTTGCTCAGCTCCGCAGTTTCGG  >  NC_002947/4967992‑4968285
                                                                                                                                                    |                                                                                                                                                 
ttGGTCGGCACGTCTGGTAGCTCAGCATTCGATCTGCTCGATGGCTGTGCACTCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCAGCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTa                                                                                                                                                   >  1:388717/1‑149 (MQ=255)
      ggCACGTCTGGTAGCTCAGCATTCGATCTGCTCGATGGCTGTGCACTCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCAGCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGAACGTTGc                                                                                                                                             >  6:219350/1‑149 (MQ=255)
      ggCAAGTCTGGTAGCTCAGCATTCGATCTTCTCGATGTCTTTGCACTCCTGTCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCAGCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGc                                                                                                                                             <  5:311059/149‑1 (MQ=255)
        cACGTCTGGTAGCTCAGCATTCGATCTGCTCGATGGCTGTGCACTCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCAGCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCg                                                                                                                                           <  5:50804/149‑1 (MQ=255)
          cGTCTGGTAGCTCAGCATTCGATCTGCTCGATGGCTGTGCACTCCAGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCAGCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCa                                                                                                                                         <  1:194988/149‑1 (MQ=255)
                gTAGCTCAGCATTCGATCTGCTCGATGGCTGTGCACTCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCAGCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCttt                                                                                                                                   <  7:117085/149‑1 (MQ=255)
                gTAGCTCAGCATTCGATCTGCTCGATGGCTGTGCACTCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCAGCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCtt                                                                                                                                    <  2:388717/148‑1 (MQ=255)
                 tAGCTCAGCATTCGATCTGCTCGATGGCTGTGCACTCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCAGCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTc                                                                                                                                  >  3:78977/1‑149 (MQ=255)
                   gCTCAGCATTCGATCTGCTCGATGGCTGTGCACTCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCAGCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGGAAGGGTACGTTGCCGCACCCtt                                                                                                                                    >  1:372760/1‑145 (MQ=255)
                   gCTCAGCATTCGATCTGCTCGATGGCTGTGCACTCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCAGCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTGCGGCAACGTACCCTTTcg                                                                                                                                 <  3:367207/148‑1 (MQ=255)
                   gCTCAGCATTCGATCTGCTCGATGGCTGTGCACTCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCAGCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTcgc                                                                                                                                >  5:346414/1‑149 (MQ=255)
                     tCAGCATTCGATCTGCTCGATGGCTGTGCACTCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCAGCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGa                                                                                                                              >  7:362222/1‑149 (MQ=255)
                       cgCATTCGATCTGCGCGTTGGCTGTGCACTCCTGGCTTCGTGAGTGGCGGGCGCAGAGAGAACAGAGTTCAGCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAACACGAAACGCTACGTTGCCGCACCCTTTCGCGa                                                                                                                              <  6:130606/146‑1 (MQ=255)
                       aGCATTCGATCTGCTCGATGGCTGTGCACTCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCAGCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGAt                                                                                                                             >  8:11401/1‑148 (MQ=255)
                        gCATTCGATCTGCTCGATGGCTGTGCACTCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCAGCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATtc                                                                                                                           <  7:11401/149‑1 (MQ=255)
                          aTTCTATATGCTCGATGGCTGTGCACTCCTGGATTCTTGTGTTGCGTGCGCAGAGAGAACAGAGTTCAGCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCATAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATtcat                                                                                                                         <  2:14593/149‑1 (MQ=255)
                                    cTCGATGGCTGTGCACTCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCAGCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTGCGGCAACGTACCCTTTcgt                                                                                                                                <  7:69727/132‑2 (MQ=255)
                                    cTCGATGGCTGTGCACTCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCAGCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTGCGGCAACGTACCCTTTcgt                                                                                                                                >  8:69727/1‑131 (MQ=255)
                                    cTCGATGGCTGTGCACTCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCAGCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATcaca                                                                                                                     <  7:30914/143‑1 (MQ=255)
                                    cTCGATGGCTGTGCACTCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCAGCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATc                                                                                                                >  4:271522/1‑148 (MQ=255)
                                      cGATGGCTGTGCACTCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCAGCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGaa                                                                                                             >  4:285971/1‑149 (MQ=255)
                                         tGGCTGTGCACTCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCAGCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTa                                                                                                          <  4:78977/149‑1 (MQ=255)
                                         tGGCTGTGCACTCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCAGCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTa                                                                                                          <  3:96444/149‑1 (MQ=255)
                                                gCACTCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCAGCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGag                                                                                                    >  6:247543/1‑148 (MQ=255)
                                                gCACTCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCAGCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGAGATTCATCACAGAATCGAATTAGCCGag                                                                                                    >  6:241526/1‑148 (MQ=255)
                                                  aCTCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCAGCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGagcag                                                                                                 >  3:261531/1‑149 (MQ=255)
                                                  aCTCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCAGCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGagcag                                                                                                 >  6:229101/1‑149 (MQ=255)
                                                       tGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCAGCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCa                                                                                            >  2:44526/1‑149 (MQ=255)
                                                              ttGAGTGGCGGGCGCAGAGAGAACAGAGTTCAGCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCg                                                                                      >  5:10645/1‑148 (MQ=255)
                                                                gAGTGGCTGGCGCAGAGAGAACAGAGTTCAGCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGTGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGaa                                                                                    <  2:311048/148‑1 (MQ=255)
                                                                gAGTGGCGGGCGCAGAGAGAACAGAGTTCAGCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGaa                                                                                    <  3:285971/148‑1 (MQ=255)
                                                                gAGTGGCGGGCGCAGAGAGAACAGAGTTCAGCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGaa                                                                                    <  5:229101/148‑1 (MQ=255)
                                                                 aGTGGCGGGCGCAGAGAGAACAGAGTTCAGCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGa                                                                                  <  4:261531/149‑1 (MQ=255)
                                                                  gTGGCGGGCGCAGAGAGAACAGAGTTCAGCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTTCCGAAGAt                                                                                 <  8:362222/149‑1 (MQ=255)
                                                                  gTGGCGGGCGCAGAGAGAACAGAGTTCAGCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGAt                                                                                 <  6:346414/149‑1 (MQ=255)
                                                                     gcgggcgCAGAGAGAACAGAGTTCAGCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGCAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGc                                                                              <  2:372760/149‑1 (MQ=255)
                                                                            cAGAGAGAACAGAGTTCAGCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGt                                                                       >  7:317698/1‑149 (MQ=255)
                                                                               agagaACAGAGTTCAGCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATTGCATCTGGTTCg                                                                    >  3:51198/1‑149 (MQ=255)
                                                                                 agaACAGAGTTCAGCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGtt                                                                      >  5:345717/1‑145 (MQ=255)
                                                                                 agaACAGAGTTCAGCTTGACACTCATATCCATGATGCTTGCATCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGtt                                                                      <  6:345717/145‑1 (MQ=255)
                                                                                   aaCAGAGTTCAGCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTg                                                                >  7:187797/1‑149 (MQ=255)
                                                                                        agTTCAGCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGa                                                             >  8:81618/1‑147 (MQ=255)
                                                                                         gTTCAGCTTGACACGCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAg                                                          >  7:367731/1‑149 (MQ=255)
                                                                                          ttCAGCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCCAATTAGCAGAGCAGCTTCAGAACCGAAGATCGCAGATTGTTCGCATGCGACac                                                          >  5:292078/1‑147 (MQ=255)
                                                                                               cTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGc                                                    <  2:219395/149‑1 (MQ=255)
                                                                                                 tGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCgg                                                  >  4:135898/1‑149 (MQ=255)
                                                                                                  gACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGGAGCTGGTTCGCCTGCGACAGGATAGCGGt                                                 >  5:385224/1‑149 (MQ=255)
                                                                                                  gACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCgg                                                  <  3:135898/148‑1 (MQ=255)
                                                                                                  gACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCgg                                                  <  3:271522/148‑1 (MQ=255)
                                                                                                   acacTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGt                                                 <  8:317698/148‑1 (MQ=255)
                                                                                                     acTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTgg                                               >  2:402389/1‑148 (MQ=255)
                                                                                                       tCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAg                                            >  1:167676/1‑149 (MQ=255)
                                                                                                       tCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAg                                            >  5:306936/1‑149 (MQ=255)
                                                                                                         atatCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCg                                                   >  8:15601/1‑140 (MQ=255)
                                                                                                         atatCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCg                                                   <  7:15601/140‑1 (MQ=255)
                                                                                                         atatCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCt                                          >  5:156096/1‑149 (MQ=255)
                                                                                                             ccATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTtgct                                      >  5:408638/1‑149 (MQ=255)
                                                                                                               atgatgCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTtgctg                                     >  2:27170/1‑148 (MQ=255)
                                                                                                                   tgCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTTGCTGCAgg                                 <  6:10645/148‑1 (MQ=255)
                                                                                                                   tgCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCAGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTTGCTGCAgg                                 <  6:274841/148‑1 (MQ=255)
                                                                                                                    gCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTTGCTGCAGGGt                               >  7:351793/1‑149 (MQ=255)
                                                                                                                      ttGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTTGCTGCAGGGttt                             >  6:17650/1‑149 (MQ=255)
                                                                                                                      ttGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTATCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTTGCTGCAGGGttt                             >  1:95901/1‑149 (MQ=255)
                                                                                                                         ccTCGTAGAATCAGAAAACGAAAGGGGACGTTGCCGCACCCTTTCGCGCTTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGGCTGCGTCAGGAGAGCGGTGGACGCTTGCTGCCGGGTTtgct                          >  2:266786/1‑149 (MQ=255)
                                                                                                                           tCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTTGCTGCAGGGTTtgctgc                        >  6:10375/1‑149 (MQ=255)
                                                                                                                             gTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTTGCTGCAGGGTTtgctgc                        <  1:402389/147‑1 (MQ=255)
                                                                                                                             gTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTTGCTGCAGGGTTtgctgc                        <  1:44526/147‑1 (MQ=255)
                                                                                                                             gTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTTGCTGCAGGGTTtgctgc                        <  5:247543/147‑1 (MQ=255)
                                                                                                                             gTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTTGCTGCAGGGTTtgctgc                        <  5:241526/147‑1 (MQ=255)
                                                                                                                             gTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTTGCTGCAGGGTTtgctgc                        <  6:385224/147‑1 (MQ=255)
                                                                                                                             gTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTTGCTGCAGGGTTTGCtgctt                      <  2:357533/149‑1 (MQ=255)
                                                                                                                              tAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGaa                                             <  4:306040/125‑1 (MQ=255)
                                                                                                                              tAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGaa                                             >  3:306040/1‑125 (MQ=255)
                                                                                                                              tAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTTGCTGCAGGGTTTGCtgcttg                     >  5:256617/1‑149 (MQ=255)
                                                                                                                                taaTCATAAAACGAAATGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACATGATAGCGGTGGAAGCTTGCTGCAGGGTTTGCtgcttgct                   <  8:367731/148‑1 (MQ=255)
                                                                                                                                 aaTCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTTGCTGCAGGGt                               <  8:117246/136‑1 (MQ=255)
                                                                                                                                 aaTCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTTGCTGCAGGGt                               >  7:117246/1‑136 (MQ=255)
                                                                                                                                   tCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTTGCTGCAGGGTTTGCTGCTTGCTCAg                >  2:10688/1‑149 (MQ=255)
                                                                                                                                      gAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTTGCTGCAGGGTTTGCTGCTTGCTCAGCTc             >  3:125263/1‑149 (MQ=255)
                                                                                                                                      aatcgcGAAAGGGTGCGGCAACGTACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTTGCTGCAGGGTTtgctg                         <  5:219973/132‑1 (MQ=38)
                                                                                                                                      aatcgcGAAAGGGTGCGGCAACGTACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTTGCTGCAGGGTTtgctg                         >  6:219973/6‑137 (MQ=38)
                                                                                                                                         aaCGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTTGCTGCAGTGTTTGCTGCTTGCTCAGCTcag           >  5:116373/1‑146 (MQ=255)
                                                                                                                                           cGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTTGCTGCAGGGTTTGCTGCTTGCTCAGCTCCGCAg        >  4:297872/1‑149 (MQ=255)
                                                                                                                                           cGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTTGCTGCAGGGTTTGCTGCTTGCTCAGCTACGCAg        >  4:292261/1‑149 (MQ=255)
                                                                                                                                               aGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTTGCTGCAGGGTTTGCTGCTTGCTCAGCTCCGCAGTTTc    <  6:256617/149‑1 (MQ=255)
                                                                                                                                                 ggTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTTGCTGCAGGGTTTGCTGCTTGCTCAGCTCCGCAGTTTCgg  >  1:50378/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
TTGGTCGGCACGTCTGGTAGCTCAGCATTCGATCTGCTCGATGGCTGTGCACTCCTGGCTTCTTGAGTGGCGGGCGCAGAGAGAACAGAGTTCAGCTTGACACTCATATCCATGATGCTTGCCTCGTAGAATCAGAAAACGAAAGGGTACGTTGCCGCACCCTTTCGCGATTCATCACAGAATCGAATTAGCCGAGCAGCTTCAGTACCGAAGATGGCAGCTGGTTCGCCTGCGACAGGATAGCGGTGGAAGCTTGCTGCAGGGTTTGCTGCTTGCTCAGCTCCGCAGTTTCGG  >  NC_002947/4967992‑4968285

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