Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,135,548 T→C 5.3% intergenic (+37/‑16) flgG → / → flgH flagellar component of cell‑distal portion of basal‑body rod/flagellar protein of basal‑body outer‑membrane L ring

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,135,5480TC 5.3% 236.9 / 3.2 75intergenic (+37/‑16)flgG/flgHflagellar component of cell‑distal portion of basal‑body rod/flagellar protein of basal‑body outer‑membrane L ring
Reads supporting (aligned to +/- strand):  ref base T (28/43);  new base C (2/2);  total (30/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.61e-01

ACGTGGCGGAAGAACTGGTCAATATGATTCAGGTGCAACGCGCTTACGAAATCAACAGTAAAGCGGTGTCCACCACCGATCAGATGCTGCAAAAACTGACGCAACTCTAAGGCTTAACCGGTGGCAGGTTCACCGGTTTACTGATTTTTGAAGATGATAGCCATGCAAAAAAACGCTGCGCATACTTATGCCATTTCCAGCTTGTTGGTGCTTTCACTAACCGGCTGCGCCTGGATACCCTCCACGCCGCTGGTGCAGGGGGCGACCAGTGCACAACCGGTTCCCGGTCCG  >  NC_000913/1135402‑1135692
                                                                                                                                                  |                                                                                                                                                
aCGTGGCGGAAGAACTGGTCAATATGATTCAGGTGCAACGCGCTTACGAAATCAACAGTAAAGCGGTGTCCACCACCGATCAGATGCTGCAAAAACTGACGCGACTCTAAGGCTTAACCGGTGGCAGGTTCACCGGTTTACTGAtttt                                                                                                                                                 <  1:261637/148‑1 (MQ=255)
aCGTGGCGGAAGAACTGGTCAATATGATTCAGGTGCAACGCGCTTACGAAATCAACAGTAAAGCGGTGTCCACCACCGATCAGATGCTGCAAAAACTGACGCAACTCTAAGGCTTAACCGGTGGCAGGTTCACCGGTTTACTGAtttt                                                                                                                                                 >  2:133790/1‑148 (MQ=255)
 cGTGGCGGAAGAACTGGTCAATATGATTCAGGTGCAACGCGCTTACGAAATCAACAGTAAAGCGGTGTCCACCACCGATCAGATGCTGCAAAAACTGACGCAACTCTAAGGCTTAACCGGTGGCAGGTTCACCGGTTTACTGAttttt                                                                                                                                                >  2:12474/1‑148 (MQ=255)
  gTGGCTGAAGAACTGTTCAATATGATTCAGGTGCAACGCGCTTACGAAATCAACAGTAAAGCGGTGTCCACCACCGATCAGATGCTGCAAAAACTGACGCAACTCTAAGGCTTAACCGGTGGCAGGTTCACCGGTTTACTGATTTTTg                                                                                                                                               <  2:168822/148‑1 (MQ=255)
     gtggAAGAACTGGTCAATATGATTCAGGTGCAACGCGCTTACGAAATCAACAGTAAAGCGGTGTCCACCACCGATCAGATGATGCAAAAACTGACGCAACTCTAAGGCTTAACCGGTGGCAGGTTCACCGGTTTACTGATTTTTGAAg                                                                                                                                            <  3:261613/146‑1 (MQ=255)
     gCGGAAGAACTGGTCAATATGATTCAGGTGCAACGCGCTTACGAAATCAACAGTAAAGCGGTGTCCACCACCGATCAGATGCTGCAAAAACTGACGCAACTCTAAGGCTTAACCGGTGGCAGGTTCACCGGTTTACTGATTTTTGaa                                                                                                                                             <  2:235665/147‑1 (MQ=255)
     gCGGAAGAACTGGTCAATATGATTCAGGTGCAACGCGCTTACGAAATCAACAGTAAAGCGGTGTCCACCACCGATCAGATGCTGCAAAAACTGACGCAACTCTAAGGCTTAACCGGTGGCAGGTTCACCGGTTTACTGATTTTTGAAg                                                                                                                                            <  5:274109/148‑1 (MQ=255)
     gCGGAAGAACTGGTCAATATGATTCAGGTGCAACGCGCTTACGAAATCAACAGTAAAGCGGTGTCCACCACCGATCAGATGCTGCAAAAACTGACGCAACTCTAAGGCTTAACCGGTGGCAGGTTCACCGGTTTACTGATTTTTGAAg                                                                                                                                            <  5:224744/148‑1 (MQ=255)
     gCGGAAGAACTGGTCAATATGATTCAGGTGCAACGCGCTTACGAAATCAACAGTAAAGCGGTGTCCACCACCGATCAGATGCTGCAAAAACTGACGCAACTCTAAGGCTTAACCGGTGGCAGGTTCACCGGTTTACTGATTTTTGAAg                                                                                                                                            <  1:59666/148‑1 (MQ=255)
        gaagaaCTGGTCAATATGATTCAGGTGCAACGCGCTTACGAAATCAACAGTAAAGCGGTGTCCACCACCGATCAGATGCTGCAAAAACTGACGCAACTCTAAGGCTTAACCGGTGGCAGGTTCACCGGTTTACTGATTTTTGAAgat                                                                                                                                          >  8:344727/1‑147 (MQ=255)
        gaagaaCTGGTCAATATGATTCAGGTGCAACGCGCTTACGAAATCAACAGTAAAGCGGTGTCCACCACCGATCAGATGCTGCAAAAACTGACGCAACTCTAAGGCTTAACCGGTGGCAGGTTCACCGGTTACTGATTTTTGAAgatga                                                                                                                                        <  2:117421/148‑1 (MQ=255)
          agaaCTGGTCAATATGATTCAGGTGCAACGCGCTTACGAAATCAACAGTAAAGCGGTGTCCACCACCGATCAGATGCTGCAAAAACTGACGCAACTCTAAGGCTTAACCGGTGGCAGGTTCACCGGTTTACTGATTTTTGAAgatgat                                                                                                                                       <  2:40725/148‑1 (MQ=255)
             aCTGGTCAATATGATTCAGGTGCAACGCGCTTACGAAATCAACAGTAAAGCGGTGTCCACCACCGATCAGATGCTGCAAAAACTGACTCAACTCTAAGTCTTAACCGGTGGCAGGTTCACCGGTTTACTGATTTTTGAAGATGATAGc                                                                                                                                    <  5:327713/148‑1 (MQ=255)
                      tatGATTCAGGTGCAACGCGCTTACGAAATCAACAGTAAAGCGGTGTCCACCACCGATCAGATGCTGCAAAAACTGACGCAACTCTAAGGCTTAACCGGTGGCAGGTTCACCGGTTTACTGATTTTTGAAGATGATAGCCATGCaaa                                                                                                                            >  8:289809/1‑147 (MQ=255)
                        tGATTCAGGTGCAACGCGCTTACGAAATCAACAGTAAAGCGGTGTCCACCACCGATCAGATGCTGCAAAAACTGACGCAACTCTAAGGCTTAACCGGTGGCAGGTTCACCGGTTTACTGATTTTTGa                                                                                                                                              >  4:150300/1‑127 (MQ=255)
                        tGATTCAGGTGCAACGCGCTTACGAAATCAACAGTAAAGCGGTGTCCACCACCGATCAGATGCTGCAAAAACTGACGCAACTCTAAGGCTTAACCGGTGGCAGGTTCACCGGTTTACTGATTTTTGa                                                                                                                                              <  3:150300/127‑1 (MQ=255)
                        tGATTCAGGTGCAACGCGCTTACGAAATCAACAGTAAAGCGGTGTCCACCACCGATCAGATGCTGCAAAAACTGACGCAACTCTAAGGCTTAACCGGTGGCAGGTTCACCGGTTTACTGATTTTTGAAGATGATAGCCATGCaaaaaa                                                                                                                         <  8:226504/148‑1 (MQ=255)
                             cAGGTGCAACGCGCTTACGAAATCAACAGTAAAGCGGTGTCCACCACCGATCAGATGCTGCAAAAACTGACGCAACTCTAAGGCTTAACCGGTGGCAGGTTCACCGGTTTACTGATTTTTGAAGATGATAGCCATGCAAAAAAACGc                                                                                                                     >  6:204068/1‑147 (MQ=255)
                             cAGGTGCAACGCGCTTACGAAATCAACAGTAAAGCGGTGTCCACCACCGATCAGATGCTGCAAAAACTGACGCAACTCTAAGGCTTAACCGGTGGCAGGTTCACCGGTTTACTGATTTTTGAAGATGATAGCCATGCAAAAAAACGc                                                                                                                     <  5:204068/147‑1 (MQ=255)
                                  gCAACGCGCTTACGAAATCAACAGTAAAGCGGTGTCCACCACCGATCAGATGCTGCAAAAACTGACGCAACTCTAAGGCTTAACCGGTGGCAGGTTCACCGGTTTACTGATTTTTGAAGATGATAGCCATGCAAAAAAACGCTGCGCa                                                                                                               >  6:347898/1‑148 (MQ=255)
                                    aaCGCGCTTACGAAATCAACAGTAAAGCGGTGTCCACCACCGATCAGATGCTGCAAAAACTGACGCAACTCTAAGGCTTAACCGGTGGCAGGTTCACCGGTTTACTGATTTTTGAAGATGATAGCCATGCAAAAAAACGCTGCGCATa                                                                                                             <  6:257304/148‑1 (MQ=255)
                                     aCGCGCTTACGAAATCAACAGTAAAGCGGTGTCCACCACCGATCAGATGCTGCAAAAACTGACGCAACTCTAAGGCTTAACCGGTGGCAGGTTCACCGGTTTACTGATTTTTGAAGATGATAGCCATGCAAAAAAACGCTGCGCATAc                                                                                                            >  1:241237/1‑148 (MQ=255)
                                     aCGCGCTTACGAAATCAACAGTAAAGCGGTGTCCACCACCGATCAGATGCTGCAAAAACTGACGCAACTCTAAGGCTTAACCGGTGGCAGGTTCACCGGTTTACTGATTTTTGAAGATGATAGCCATGCAAAAAAACGCTGCGCATAc                                                                                                            >  6:96740/1‑148 (MQ=255)
                                       gcgcTTACGAAATCAACAGTAAAGCGGTGTCCACCACCGATCAGATGCTGCAAAAACTGACGCAACTCTAAGGCTTAACCGGTGGCAGGTTCACCGGTTTACTGATTTTTGAAGATGATAGCCATGCAAAAAAACGCTGCGCATACtt                                                                                                          <  5:96740/148‑1 (MQ=255)
                                          cTTACGAAATCAACAGTAAAGCGGTGTCCACCACCGATCAGATGCTGCAAAAACTGACGCAACTCTAAGGCTTAACCGGTGGCAGGTTCACCGGTTTACTGATTTTTGAAGATGATAGCCATGCAAAAAAACGCTGCGCATACTTATg                                                                                                       >  3:190872/1‑148 (MQ=255)
                                             aCGCAAGCAACAGTAAAGCGGTGTCCACAACCGATCAGCTGCTGCAAAAACTGACGCAACTCTAAGGCTTAACCGGGGGCAGGTTCACCGGTTTACTGATTTTTGAAGATGATAGCGATGCAAAAAAACGCTGCGCATACTTATGcc                                                                                                     <  5:225331/147‑1 (MQ=255)
                                               gAAATCAACAGTAAAGCGGTGTCCACCACCGATCAGATGCTGCAAAAACTGACGCAACTCTAAGGCTTAACCGGTGGCAGGTTCACCGGTTTACTGATTTTTGAAGATGATAGCCATGCAAAAAAACGCTGCGCATACTTATGCCAtt                                                                                                  <  8:318124/148‑1 (MQ=255)
                                                   tCAACAGTAAAGCGGTGTCCACCACCGATCAGATGCTGCAAAAACTGACGCAACTCTAAGGCTTAACCGGTGGCAGGTTCACCGGTTTACTGATTTTTGAAGATGATAGCCATGCAAAAAAACGCTGCGCATACTTATGCCATTTc                                                                                                >  1:154359/1‑146 (MQ=255)
                                                   tCAACAGTAAAGCGGTGTCCACCACCGATCAGATGCTGCAAAAACTGACGCAACTCTAAGGCTTAACCGGTGGCAGGTTCACCGGTTTACTGATTCTTGa                                                                                                                                              >  8:313192/1‑100 (MQ=255)
                                                   tCAACAGTAAAGCGGTGTCCACCACCGATCAGATGCTGCAAAAACTGACGCAACTCTAAGGCTTAACCGGTGGCAGGTTCACCGGTTTACTGATTCTTGa                                                                                                                                              <  7:307842/100‑1 (MQ=255)
                                                   tCAACAGTAAAGCGGTGTCCACCACCGATCAGATGCTGCAAAAACTGACGCAACTCTAAGGCTTAACCGGTGGCAGGTTCACCGGTTTACTGATTCTTGa                                                                                                                                              <  7:313192/100‑1 (MQ=255)
                                                   tCAACAGTAAAGCGGTGTCCACCACCGATCAGATGCTGCAAAAACTGACGCAACTCTAAGGCTTAACCGGTGGCAGGTTCACCGGTTTACTGATTCTTGa                                                                                                                                              >  8:307842/1‑100 (MQ=255)
                                                      aCAGTAAAGCGGTGTCCACCACCGATCAGATGCTGCAAAAACTGACGCAACTCTAAGGCTTAACCGGTGGCAGGTTCACCGGTTTACTGATTTTTGAAGATGATAGCCATGCAAAAAAACGCTGCGCATACTTATGCCATTTCCAGCt                                                                                           >  1:69580/1‑148 (MQ=255)
                                                       cAGTAAAGCGGTGTCCACCACCGATCAGATGCTGCAAAAACTGACGCAACTCTAAGGCTTAACCGGTGGCAGGTTCACCGGTTTACTGATTTTTGAAGATGATAGCCATGCAAAAAAACGCTGCGCATACTTATGCCATTTCCAg                                                                                             >  2:322101/1‑145 (MQ=255)
                                                       cAGTAAAGCGGTGTCCACCACCGATCAGATGCTGCAAAAACTGACGCAACTCTAAGGCTTAACCGGTGGCAGGTTCACCGGTTTACTGATTTTTGAAGATGATAGCCATGCAAAAAAACGCTGCGCATACTTATGCCATTTCCAg                                                                                             <  1:322101/145‑1 (MQ=255)
                                                       cAGTAAAGCGGTGTCCACCACCGATCAGATGCTGCAAAAACTGACGCAACTCTAAGGCTTAACCGGTGGCAGGTTCACCGGTTTACTGATTTTTGAAGATGATAGCCATGCAAAAAAACGCTGCGCATACTTATGCCATTTCCAGCtt                                                                                          >  1:191860/1‑148 (MQ=255)
                                                        cgTAAAGCGGTGTCCACCACCGATCAGATGCTGCAAAAACTGACGCAACTCTAAGGCTTAACCGGTGGCAGGTTCACCGGTTTACTGATTTTTGAAGATGATAGCCATGCAAAAAAACGCTGCGCATACTTATGCCATTTCCAGCttg                                                                                         <  6:230548/147‑1 (MQ=255)
                                                         gTAAAGCGGTGTCCACCACCGATCAGATGCTGCAAAAACTGACGCAACTCTAAGGCTTAACCGGTGGCAGGTTCACCGGTTTACTGATTTTTGAAGATGATAGCCATGCAAAAAAACGCTGCGCATACTTATGCCATTTCCAGCttgt                                                                                        >  7:49845/1‑148 (MQ=255)
                                                          tAAAGCGGTGTCCACCACCGATCAGATGCTGCAAAAACTGACGCAACTCTAAGGCTTAACCGGTGGCAGGTTCACCGGTTTACTGATTTTTGAAGATGATAGCCATGCAAAAAAACGCTGCGCATACTTATGCCATTTCCAGCttgtt                                                                                       <  7:289809/148‑1 (MQ=255)
                                                           aaaGCGGTGTCCACCACCGATCAGATGCTGCAAAAACTGACGCAACTCTAAGGCTTAACCGGTGGCAGGTTCACCGGTTTACTGATTTTTGAAGATGATAGCCATGCAAAAAAACGCTGCGCATACTTATGCCATTTCCAGCttgttg                                                                                      >  8:43288/1‑148 (MQ=255)
                                                            aaGCGGTGTCCACCACCGATCAGATGCTGCAAAAACTGACGCAACTCTAAGGCTTAACCGGTGGCAGGTTCACCGGTTTACTGATTTTTGAAGATGATAGCCATGCAAAAAAACGCTGCGCATACTTATGCCATTTCCAGCttgttg                                                                                      <  1:12474/147‑1 (MQ=255)
                                                                ggTGTCCACCACCGATCAGATGCTGCAAAAACTGACGCAACTCTAAGGCTTAACCGGTGGCAGGTTCACCGGTTTACTGATTTTTGAAGATGATAGCCATGCAAAAAAACGCTGCGCATACTTATGCCATTTCCAGCTTGTTGGTGCt                                                                                 <  1:133790/148‑1 (MQ=255)
                                                                ggTGTCCACCACCGATCAGATGCTGCAAAAACTGACGCAACTCTAAGGCTTAACCGGTGGCAGGTTCACCGGTTTACTGATTTTTGAAGATGATAGCCATGCAAAAAAACGCTGCGCATACTTATGCCATTTCCAGCTTGTTGGTGCt                                                                                 <  2:154359/148‑1 (MQ=255)
                                                                    tgcaccacCGATCAGAGGCAGCAAAAACTGAAGCAACTCTAAGTCGTAACCGGTGGCATGTTCTCCGGTTTACTGATTTTTGAAGATGATAGCCATGCTAATATAAGCTGCGCATACTTATGCCATTTCCAGCTTGTTGGTGCTTTc                                                                              <  4:139057/145‑1 (MQ=255)
                                                                    tCCACCACCGATCAGATGCTGCAAAAACTGACGCAACTCTAAGGCTTAACCGGTGGCAGGTTCACCGGTTTACTGATTTTTGAAGATGATAGCCATGCAAAAAAACGCTGCGCATACTTATGCCATTTCCAGCTTGTTGGTGCTTTc                                                                              <  5:347898/147‑1 (MQ=255)
                                                                    tCCACCACCGATCAGATGCTGCAAAAACTGACGCAACTCTAAGGCTTAACCGGTGGCAGGTTCACCGGTTTACTGATTTTTGAAGATGATAGCCATGCAAAAAAACGCTGCGCATACTTATGCCATTTCCAGCTTGTTGGTGCTTTc                                                                              <  8:215859/147‑1 (MQ=255)
                                                                    tCCACCACCGATCAGATGCTGCAAAAACTGACGCAACTCTAAGGCTTAACCGGTGGCAGGTTCACCGGTTTACTGATTTTTGAAGATGATAGCCATGCAAAAAAACGCTGCGCATACTTATGCCATTTCCAGCTTGTTGGTGCTTTc                                                                              <  8:127121/147‑1 (MQ=255)
                                                                     ccaccaCCGATCAGATGCTGCAAAAACTGACGCAACTCTAAGGCTTAACCGGTGGCAGGTTCACCGGTTTACTGATTTTTGAAGATGATAGCCATGCAAAAAAACGCTGCGCATACTTATGCCATTTCCAGCTTGTTGGTGCTTTCAc                                                                            >  1:46012/1‑148 (MQ=255)
                                                                     ccaccaCCGATCAGATGCTGCAAAAACTGACGCAACTCTAAGGCTTAACCGGTGGCAGGTTCACCGGTTTACTGATTTTTGAAGATGATAGCCATGCAAAAAAACGCTGCGCATACTTATGCCATTTCCAGCTTGTTGGTGCTTTCAc                                                                            >  2:138078/1‑148 (MQ=255)
                                                                        ccaccGATCAGATGCTGCAAAAACTGACGCAACTCTAAGGCTTAACCGGTGGCAGGTTCACCGGTTTACTGATTTTTGAAGATGATAGCCATGCAAAAAAACGCTGCGCATACTTATGCCATTTCCAGCTTGTTGGTGCTTTCACTaa                                                                         >  3:130972/1‑148 (MQ=255)
                                                                         caccGATCAGATGCTGCAAAAACTGACGCAACTCTAAGGCTTAACCGGTGGCAGGTTCACCGGTTTACTGATTTTTGAAGATGATAGCCATGCAAAAAAACGCTGCGCATACTTATGCCATTTCCAGCTTGTTGGTGCTTTCACTAAc                                                                        <  7:344727/148‑1 (MQ=255)
                                                                          accGATCAGATGCTGCAAAAACTGACGCAACTCTAAGGCTTAACCGGTGGCAGGTTCACCGGTTTACTGATTTTTGAAGATGATAGCCATGCAAAAAAACGCTGCGCATACTTATGCCATTTCCAGCTTGTTGGTGCTTTCACTAAcc                                                                       >  6:248725/1‑148 (MQ=255)
                                                                            cGATCAGATGCTGCAAAAACTGACGCAACTCTAAGGCTTAACCGGTGGCAGGTTCACCGGTTTACTGATTTTTGAAGATGATAGCCATGCAAAAAAACGCTgcgc                                                                                                                <  5:108858/105‑1 (MQ=255)
                                                                            cGATCAGATGCTGCAAAAACTGACGCAACTCTAAGGCTTAACCGGTGGCAGGTTCACCGGTTTACTGATTTTTGAAGATGATAGCCATGCAAAAAAACGCTgcgc                                                                                                                >  6:108858/1‑105 (MQ=255)
                                                                            cGATCAGATGCTGCAAAAACTGACGCAACTCTAAGGCTTAACCGGTGGCAGGTTCACCGGTTTACTGATTTTTGAAGATGATAGCCATGCAAAAAAACGCTGCGCATACTTATGCCATTTCCAGCTTGTTGGTGCTTTCACTAACCgg                                                                     >  8:57290/1‑148 (MQ=255)
                                                                              aTCAGATGCTGCAAAAACTGACGCAACTCTAAGGCTTAACCGGTGGCAGGTTCACCGGTTTACTGATTTTTGAAGATGATAGCCATGCAAAAAAACGCTGCGCATACTTAGGCCATTTCCAGCTTGTTGGTGCTTTCACTACCCGGCt                                                                   >  6:187739/1‑148 (MQ=255)
                                                                                  gATGCTGCAAAAACTGAAGCAACTCTAAGGCTTAACCGGTTGAAGGTTCACCGGTTTACTGATTTTTGAAGCTGATAGCCATGCAAAAAAACGCTGCGCATACTTATGCCATTTCCAGCTTGTTGGTGCTTTCACTAACCGGCTgcgc                                                               <  3:186290/148‑1 (MQ=255)
                                                                                   aTGCTGCAAAAACTGACGCAACTCTAAGGCTTAACCGGTGGCAGGTTCACCGGTTTACTGATTTTTGAAGATGATAGCCATGCAAAAAAACGCTGCGCATACTTATGCCATTTCCAGCTTGTTGGTGCTTTCACTAACCGGCTGCGcc                                                              <  5:266788/148‑1 (MQ=255)
                                                                                    tgctccAAAAACTGACGCAACTCTAAGGCTTAACCGGTGGCAGGTTCACCGGTTTACTGATTTTTGAAGATGCTAGCCATGCAAAAAAACGCTGCGCATACTTATGCCATTTCCAGCTTGTTGGTGCTTTCACTAACCGGCTGCGCCt                                                             <  6:308562/148‑1 (MQ=255)
                                                                                      ctgcAAAAACTGACGCAACTCTAAGGCTTAACCGGTGGCAGGTTCACCGGTTTACTGATTTTTGAAGATGATAGCCATGCAAAAAAACGCTGCGCATACTTATGCCATTTCCAGCTTGTTGGTGCTTTCACTAACCGGCTGCGCCTgg                                                           <  8:129916/148‑1 (MQ=255)
                                                                                        gcAAAAACTGACGCAACTCTAAGGCTTAACCGGTGGCAGGTTCACCGGTTTACTGATTTTTGAAGATGATAGCCATGCAAAAAAACGCt                                                                                                                    >  2:278845/1‑89 (MQ=255)
                                                                                        gcAAAAACTGACGCAACTCTAAGGCTTAACCGGTGGCAGGTTCACCGGTTTACTGATTTTTGAAGATGATAGCCATGCAAAAAAACGCt                                                                                                                    <  1:278845/89‑1 (MQ=255)
                                                                                         cAAAAACTGAAGCAACTCTAAGGCTTAACCGGTGTCAGGTTCACCGGTTTACTGATTTTTGAAGATGATAGCCATGCAAAAAAACGCTGCGCATACTTATGCCATTTCCAGCTTGTTGGTGCTTTCACTAACCGGCTGCGCCTGGAt                                                         <  8:153603/147‑1 (MQ=255)
                                                                                          aaaaaCTGACGCAACTCTAAGGCTTAACCGGTGGCAGGTTCACCGGTTTACTGATTTTTGAAGATGATAGCCATGCAAAAAAACGCTGCGCATACTTATGCCATTTCCAGCTTGTTGGTGCTTTCACTAACCGGCTGCGCCTGGATAc                                                       <  4:190872/148‑1 (MQ=255)
                                                                                                    gCAACTCTAAGGCTTAACCGGTGGCAGGTTCACCGGTTTACTGATTTTTGAAGATGATAGCCATGCAAAAAAACGCTGCGCATACTTATGCCATTTCCAGCTTGTTGGTGCTTTCACTAACCGGCTGCGCCTGGATACCCTCCAcgcc                                             <  8:257641/148‑1 (MQ=255)
                                                                                                           tAAGGCTTAACCGGTGGCAGGTTCACCGGTTTACTGATTTTTGAAGATGATAGCCATGCAAAAAAACGCTGCGCATACTTATGCCATTTCCAGCTTGTTGGTGCTTTCACTAACCGGCTGCGCCTGGATACCCTCCACGCCGCTgg                                        <  2:46012/146‑1 (MQ=255)
                                                                                                            aaGGCTTAACCGGTGGCAGGTTCACCGGTTTACTGATTTTTGAAGATGATAGCCATGCAAAAAAACGCTGCGCATACTTCTGCCATTTCCCGCTTGTTGGTGCTTTCACTAACCGGCTGCGCCTGGATACCCTCCACGCCGCTGGTGc                                     <  3:208925/148‑1 (MQ=255)
                                                                                                             aGGCTTAACCGGTGGCAGGTTCACCGGTTTACTGATTTTTGAAGATGATAGCCATGCAAAAAAACGCTGCGCATACTTATGCCATTTCCAGCTTGTTGGTGCTTTCACTAACCGGCTGCGCCTGGATACCCTCCACGCCGCTGGTGCa                                    <  2:241237/148‑1 (MQ=255)
                                                                                                             aGGCTTAACCGGTGGCAGGTTCACCGGTTTACTGATTTTTGAAGATGATAGCCATGCAAAAAAACGCTGCGCATACTTATGCCATTTCCAGCTTGTTGGTGCTTTCACTAACCGGCTGCGCCTGGATACCCTCCACGCCGCTGGTGCa                                    <  5:318095/148‑1 (MQ=255)
                                                                                                             aGGCTTAACCGGTGGCAGGTTCACCGGTTTACTGATTTTTGAAGAAGATAGCCATGCAAAAAAACGCTGCGCATACTTATGCCATTTCCAGCTTGTTGGTGCTTTCACTAACCGGCTGCGCCTGGATACCCTCCACGCCGCTGGTg                                      >  2:133310/1‑146 (MQ=255)
                                                                                                              ggCTTAACCGGTGGCAGGTTCACCGGTTTACTGATTTTTGAAGATGATAGCCATGCAAAAAAACGCTGCGCATACTTATGCCATTTCCAGCTTGTTGGTGCTTTCACTAACCGGCTGCGCCTGGATACCCTCCACGCCGCTGGTGCAg                                   <  8:49845/148‑1 (MQ=255)
                                                                                                              ggCTTAACCGGTGGCAGGTTCACCGGTTTACTGATTTTTGAAGATGATAGCCATGCAAAAAAACGCTGCGCATACTTATGCCATTTCCAGCTTGTTGGTGCTTTCACTAACCGGCTGCGCCTGGATACCCTCCACGCCGCTGGTGCAg                                   <  2:69580/148‑1 (MQ=255)
                                                                                                                          ggCAGGTTCACCGGTTTACTGATTTTTGAAGATGATAGCCATGCAAAAAAACGCTGCGCATACTTATGCCATTTCCAGCTTGTTGGTGCTTTCACTAACCGGCTGCGCCTGGATACCCTCCACGCCGCTGGTGCAGGGGGCGACCAGt                       >  7:64660/1‑148 (MQ=255)
                                                                                                                               gTTCACCGGTTTACTGATTTTTGAAGATGATAGCCATGCAAAAAAACGCTGCGCATACTTATGCCATTTCCAGCTTGTTGGTGCTTTCACTAACCGGCTGCGCCTGGATACCCTCCACGCCGCTGGTGCAGGGGGCGACCAGTGcac                   >  7:148491/1‑147 (MQ=255)
                                                                                                                                ttCACCGGTTTACTGATTTTTGAAGATGATAGCCATGCAAAAAAACGCTGCGCATACTTATGCCATTTCCAGCTTGTTGGTGCTTTCACTAACCGGCTGCGCCTGGATACCCTCCACGCCGCTGGTGCAGGGGGCGACCAGTGCACaa                 >  8:374796/1‑148 (MQ=255)
                                                                                                                                   aCCGGTTTACTGATTTTTGAAGATGATAGCCATGCAAAAAAACGCTGCGCATACTTATGCCATTTCCAGCTTGTTGGTGCTTTCACTAACCGCCTGCGCCTGGATACCCTCCACGCCGCTGGTGCAGGGGGCGACCAGTGCACAACCg              >  3:141995/1‑148 (MQ=255)
                                                                                                                                     cGGTTTACTGATTTTTGAAGATGATAGCCATGCAAAAAAACGCTGCGCATACTTATGCCATTTCCAGCTTGTTGGTGCTTTCACTAACCGGCTGCGCCTGGATACCCTCCACGCCGCTGGTGCAGGGGGCGACCAGTGCACAACCGGt            <  4:130972/148‑1 (MQ=255)
                                                                                                                                               aTTTTTGAAGATGATAGCCATGCAAAAAAACGCTGCGCATACTTATGCCATTTCCAGCTTGTTGGTGCTTTCACTAACCGGCTGCGCCTGGATACCCTCCACGCCGCTGGTGCAGGGGGCGACCAGTGCACAACCGGTTCCCGGTCCg  >  1:96975/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
ACGTGGCGGAAGAACTGGTCAATATGATTCAGGTGCAACGCGCTTACGAAATCAACAGTAAAGCGGTGTCCACCACCGATCAGATGCTGCAAAAACTGACGCAACTCTAAGGCTTAACCGGTGGCAGGTTCACCGGTTTACTGATTTTTGAAGATGATAGCCATGCAAAAAAACGCTGCGCATACTTATGCCATTTCCAGCTTGTTGGTGCTTTCACTAACCGGCTGCGCCTGGATACCCTCCACGCCGCTGGTGCAGGGGGCGACCAGTGCACAACCGGTTCCCGGTCCG  >  NC_000913/1135402‑1135692

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