Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,134,668 T→A 6.9% intergenic (+111/‑61) flgF → / → flgG flagellar component of cell‑proximal portion of basal‑body rod/flagellar component of cell‑distal portion of basal‑body rod

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,134,6680TA 6.9% 238.8 / 5.8 86intergenic (+111/‑61)flgF/flgGflagellar component of cell‑proximal portion of basal‑body rod/flagellar component of cell‑distal portion of basal‑body rod
Reads supporting (aligned to +/- strand):  ref base T (38/42);  new base A (3/3);  total (41/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.94e-01

CGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAGCAAACCAATATGGACGTCATTGCCAACAACCTGGCAAACGT  >  NC_000913/1134521‑1134814
                                                                                                                                                   |                                                                                                                                                  
cGCAGGCCGTGCCAACCACCTGCTGTCGATGAGTTAATTGAAAGGAGACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTTTCGATTTTGCGAGCAATTGTCGGCCGGATAAGGCGttt                                                                                                                                                   >  2:1253713/1‑149 (MQ=255)
  cAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTAc                                                                                                                                                 >  1:572628/1‑149 (MQ=255)
   cggCCGTGCCAACCAACTGCTGGCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTAcg                                                                                                                                                <  1:537467/148‑1 (MQ=255)
     gCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTATAcgcc                                                                                                                                              >  2:429944/1‑149 (MQ=255)
             aaccaacTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATGCGGCGTAAACGCCTTATCCgg                                                                                                                                      <  1:934959/149‑1 (MQ=255)
             aaccaacTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCgg                                                                                                                                      <  2:780257/149‑1 (MQ=255)
              cccaccTGCTGTCGATGAGTTAATTGCCCGGATACATGACCAGTATACGTTGCCCGATGCGCAAGGTTATCGGGGCTATGGGGGCAATCGCAAGTTATCGATTTTGCGAGCACTTGTAGGCCGGATACGGCGTTTACGCCGCATCCGGc                                                                                                                                     <  1:260483/148‑1 (MQ=38)
              accaacTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTAATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGc                                                                                                                                     <  2:1230000/149‑1 (MQ=255)
               ccaacTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCATTCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCa                                                                                                                                    >  1:936563/1‑149 (MQ=255)
               ccaacTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCa                                                                                                                                    <  2:1425438/149‑1 (MQ=255)
                 aacTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCaag                                                                                                                                  >  1:147304/1‑149 (MQ=255)
                   ctgctgTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCaagaa                                                                                                                                <  2:802478/149‑1 (MQ=255)
                    tgctgTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAAATTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCaagaag                                                                                                                               >  1:335252/1‑149 (MQ=255)
                     gctgTCGATGAGTAAATTGAAAGGATCAATGTCAAGTATAAGTTTCCCGATGCGCAATTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAagaaga                                                                                                                              <  2:155409/149‑1 (MQ=38)
                      ctgTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCAGATGCGCAATTTTATCGGGTCTATGGGGGCAATCGCAATTTAACGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGAc                                                                                                                             >  2:500131/1‑149 (MQ=255)
                          cGATTAGTTAAATGAAACCACACATGACAACTATAATGTGCCCGATGCTCAAGGATATCGGTGCTATGGGGGCAATCGCAATATATCGATTATGCGATCACTTGTAGTCCGGATAAGGAGTTTACGCCGCATCCGGCAAGAAGACatat                                                                                                                         <  2:1013451/149‑1 (MQ=255)
                           gATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATg                                                                                                                        >  1:545331/1‑149 (MQ=255)
                           gATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATg                                                                                                                        <  1:574143/149‑1 (MQ=255)
                                     ttcaaaGGAAACATGACAAGTATAAGTTTCCCTCTGCGCAAGTTTATGGGTGCTATGGGGGCATTCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCTGCAAGAAGACATATGCACTTTGTCa                                                                                                              <  2:1061810/146‑1 (MQ=37)
                                      tGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCAc                                                                                                             >  2:262101/1‑149 (MQ=255)
                                       gAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACt                                                                                                            >  1:1318024/1‑149 (MQ=255)
                                       gAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCAGCATCCGGCAAGAAGACATATGCACTTTGTCACt                                                                                                            >  2:1483279/1‑149 (MQ=255)
                                        aaaGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGAGAGCAGTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTa                                                                                                           >  2:373208/1‑149 (MQ=255)
                                         aaGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTaa                                                                                                          >  1:739596/1‑149 (MQ=255)
                                           ggATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGTGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATc                                                                                                        <  1:303161/149‑1 (MQ=255)
                                           ggATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATc                                                                                                        <  1:247556/149‑1 (MQ=255)
                                              aaCATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCAc                                                                                                     <  1:145351/148‑1 (MQ=255)
                                                  tGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACa                                                                                                 >  1:809566/1‑149 (MQ=255)
                                                    aCAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAgg                                                                                               >  2:376384/1‑149 (MQ=255)
                                                    aCAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAgg                                                                                               <  1:166424/149‑1 (MQ=255)
                                                         tataAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGTTAAGGCGTTTACGCCGCATCCGGCAAGATGACGTATGCACTTAGTACCTAATCCACTACAGGACAtt                                                                                          >  2:1499384/1‑149 (MQ=39)
                                                         tataAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACAtt                                                                                          >  2:1418630/1‑149 (MQ=255)
                                                         tataAGGTGCCCGATGCGCAAGTTTATCGGGTCTATGGGTGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCCTCCGGCAAGAAGACATATGCACTTAGTCACTAAACCACTACAGGACAtt                                                                                          >  2:1412861/1‑149 (MQ=39)
                                                          ataAGTTGCCCGATGCGCCAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACAttt                                                                                         <  2:739441/149‑1 (MQ=255)
                                                          ataAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACAttt                                                                                         >  1:758717/1‑149 (MQ=255)
                                                          ataAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTATACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACAttt                                                                                         <  1:429944/149‑1 (MQ=255)
                                                           taAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACAtttt                                                                                        <  1:573394/149‑1 (MQ=255)
                                                           taAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACAtttt                                                                                        <  1:579822/149‑1 (MQ=255)
                                                              gTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATg                                                                                     <  1:947703/149‑1 (MQ=255)
                                                               ttGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGa                                                                                    >  2:241803/1‑149 (MQ=255)
                                                                 gCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATc                                                                                  <  2:951029/149‑1 (MQ=255)
                                                                 gCCCGAGGCGCAAGTTTATCGGGTCTATGGCGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGagc                                                                                  <  2:1410520/149‑3 (MQ=255)
                                                                    cGATGCTCATTTTTATCGGGTCTTTGTGGGTAATCTCAATTTATCGATTTTGCGAGCACTTGTAGTCCGGATAAGGAGTTTAGGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGt                                                                               <  2:128347/149‑1 (MQ=255)
                                                                     gATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACTTATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGtt                                                                              >  1:1087714/1‑149 (MQ=255)
                                                                     gATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGtt                                                                              >  2:1482417/1‑149 (MQ=255)
                                                                       tGCGCAAGTTTATCGGGTCTAGGGGGGCAATCGCAATTTAGCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCACTTCa                                                                            <  1:861014/149‑1 (MQ=255)
                                                                          gcAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTa                                                                         >  1:1221468/1‑149 (MQ=255)
                                                                               tttATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGAt                                                                    <  1:735261/149‑1 (MQ=255)
                                                                                tgcgcGGCTCCATGGGGTCATTCGCTTTTTTTCTATTTTTCGAGCCCTTGCCCGCCGCCTCTGGCGTATACGCCGCAACCGGCAAGAAGACATCTGCACTTTGTCACTCATCCACTACAGGACATTTTATGATCCGTTCATAATGGATc                                                                   <  2:1009117/145‑1 (MQ=255)
                                                                                   tCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCAGTTGTAGGCCGGATTAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGcc                                                                <  1:1365138/149‑1 (MQ=255)
                                                                                   tCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGcc                                                                <  2:835831/149‑1 (MQ=255)
                                                                                          tATGGGGGCAATCGCAATTTATCGATTTTGCGAGCTCTTGTAGGACGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACgg                                                         <  2:1432672/149‑1 (MQ=255)
                                                                                          tATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACgg                                                         <  2:1544644/149‑1 (MQ=255)
                                                                                           aTGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATTCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACggg                                                        <  1:880981/149‑1 (MQ=255)
                                                                                           aTGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACggg                                                        >  1:57705/1‑149 (MQ=255)
                                                                                           aTGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACggg                                                        >  2:312194/1‑149 (MQ=255)
                                                                                            tGGGGGCAATCGCTATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGc                                                       <  2:194135/149‑1 (MQ=255)
                                                                                               gggCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGTCCGGATAAGGCGTTTACGCCGCATCCTGCAAGAAGACATATGCACTTTGTCAATAATACAATACAGGAAATTTAATGATCAGTTCATTATGGAACGCCAAAACTGTCCtt                                                    >  2:1434758/1‑149 (MQ=38)
                                                                                               gggCAATCGCAATTAATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCtt                                                    <  2:649745/149‑1 (MQ=255)
                                                                                                 gCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGa                                                  >  2:1434517/1‑149 (MQ=255)
                                                                                                 gCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGa                                                  >  2:29398/1‑149 (MQ=255)
                                                                                                   aaTCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACg                                                <  1:376384/149‑1 (MQ=255)
                                                                                                      cGCAATTTATCGATTTTGCGAGCACTTGTCGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTACTCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGccc                                             <  2:604405/149‑1 (MQ=255)
                                                                                                      cGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGccc                                             >  2:282449/1‑149 (MQ=255)
                                                                                                         aaTTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCACAACGCGCATTGACGCCCAGc                                          >  1:739147/1‑149 (MQ=255)
                                                                                                            ttATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAGCaaa                                       <  1:373208/149‑1 (MQ=255)
                                                                                                             tATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAGCAAAc                                      <  2:956630/149‑1 (MQ=255)
                                                                                                                  aTTTTGCGAGCAATTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAGCAAACCAata                                 <  2:1394155/149‑1 (MQ=255)
                                                                                                                      tGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAGCAAACCAATATGGa                             <  2:147304/149‑1 (MQ=255)
                                                                                                                       gCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGAACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAGCAAACCAATATGGAc                            <  2:1186993/149‑1 (MQ=255)
                                                                                                                         gAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAGCAAACCAATATGGACGt                          >  1:696073/1‑149 (MQ=255)
                                                                                                                         gAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAGCAAACCAATATGGACGt                          >  2:694246/1‑149 (MQ=255)
                                                                                                                          aGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAGCAAACCAATATGGACGTc                         >  2:415098/1‑149 (MQ=255)
                                                                                                                            cACTTGTAGGCCGGATAAGGCGGTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAGCAAACCAATATGGACGTCAt                       <  1:157206/149‑1 (MQ=255)
                                                                                                                              cTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAGCAAACCAATATGGACGTCATTg                     >  1:89160/1‑149 (MQ=255)
                                                                                                                               ttGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAGCAAACCAATATGGACGTCATTGc                    >  2:59373/1‑149 (MQ=255)
                                                                                                                               ttGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAGCAAACCAATATGGACGTCATTGc                    <  2:936563/149‑1 (MQ=255)
                                                                                                                                   ttgCCGGATGCGGCGTAAACGCCTTATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAGCAAACCAATATGGACGTCATTGCcaac                >  1:1260233/3‑149 (MQ=255)
                                                                                                                                   ttgCCGGATGCGGCGTAAACGCCTTATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCATTGACGCCCAGCAAACCAATATGGAAGTCATTGAcaac                >  2:982981/3‑149 (MQ=255)
                                                                                                                                     gCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAGCAAACCAATATGGACGTCATTGCcaacac              >  2:153647/1‑148 (MQ=255)
                                                                                                                                     gCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTAATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAGCAAACCAATATGGACGTCATTGCcaacaa              >  2:15433/1‑149 (MQ=255)
                                                                                                                                      ccGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAGCAAACCAATATGGACGTCATTGCCaacaac             >  2:1326529/1‑149 (MQ=255)
                                                                                                                                       cGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAGCAACCCAATATGGACGTCATTGCCACCAccc            >  1:121066/1‑149 (MQ=255)
                                                                                                                                          aTAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAGCAAACCAATATGGACGTCATTGCCAACAACCTgg         <  1:1418630/149‑1 (MQ=255)
                                                                                                                                          aTAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATACACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAGCAAACCAATATGGACGTCATTGCCAACAACCTgg         <  1:1412861/149‑1 (MQ=255)
                                                                                                                                               gCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGATCATTATGGATCGCCAAAACGGGCCTTGACGCCCAGCAAACCAATATGTACGTCATTGCCAACAACCTGGCAAAc    >  2:1264150/1‑149 (MQ=255)
                                                                                                                                                 gTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAGCAAACCAATATGGACGTCATTGCCAACAACCTGGCAAACGt  <  2:739596/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CGCAGGCCGTGCCAACCAACTGCTGTCGATGAGTTAATTGAAAGGATACATGACAAGTATAAGTTGCCCGATGCGCAAGTTTATCGGGTCTATGGGGGCAATCGCAATTTATCGATTTTGCGAGCACTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAAGAAGACATATGCACTTTGTCACTAATCCACTACAGGACATTTTATGATCAGTTCATTATGGATCGCCAAAACGGGCCTTGACGCCCAGCAAACCAATATGGACGTCATTGCCAACAACCTGGCAAACGT  >  NC_000913/1134521‑1134814

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