Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,235,581 T→A 5.1% intergenic (+174/‑325) pgi → / → yjbE glucosephosphate isomerase/extracellular polysaccharide production threonine‑rich protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,235,5810TA 5.1% 478.6 / 12.3 137intergenic (+174/‑325)pgi/yjbEglucosephosphate isomerase/extracellular polysaccharide production threonine‑rich protein
Reads supporting (aligned to +/- strand):  ref base T (66/64);  new base A (3/4);  total (69/68)
Fisher's exact test for biased strand distribution p-value = 7.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.28e-01

Only 100 of 138 total aligned reads displayed.

GCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTAT  >  NC_000913/4235435‑4235728
                                                                                                                                                  |                                                                                                                                                   
gCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTAt                                                                                                                                                   >  4:790214/1‑149 (MQ=255)
gCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTAt                                                                                                                                                   <  2:1394281/149‑1 (MQ=255)
gCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTAt                                                                                                                                                   >  4:1339972/1‑149 (MQ=255)
   ttCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTc                                                                                                                                                >  1:1487962/1‑149 (MQ=255)
    tCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCg                                                                                                                                               <  1:902472/149‑1 (MQ=255)
       cgccgcATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACATAACGTCGGCa                                                                                                                                            <  2:205036/149‑1 (MQ=255)
        gccgcATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACATAACGTCGGCAt                                                                                                                                           <  1:1553744/149‑1 (MQ=255)
         ccgcATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACATAACGTCGGCAtt                                                                                                                                          <  2:1151673/149‑1 (MQ=255)
         ccgcATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACATAACGTCGGCAtt                                                                                                                                          <  1:1575024/149‑1 (MQ=255)
              tCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCAttttttt                                                                                                                                     <  4:333393/149‑1 (MQ=255)
              tCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCAttttttt                                                                                                                                     <  2:1572547/149‑1 (MQ=255)
              tCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCAttttttt                                                                                                                                     <  2:241053/149‑1 (MQ=255)
                cGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTAt                                                                                                                                   <  3:538562/149‑1 (MQ=255)
                 ggCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATc                                                                                                                                  >  1:130053/1‑149 (MQ=255)
                 ggCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGAATAAATAACCTTAATAACAATGCCGACGTTATGTCGGCAtttttt                                                                                                                                      >  3:588986/1‑145 (MQ=255)
                          tGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATc                                                                                                                         <  3:284872/149‑1 (MQ=255)
                           gCCTGATAAGACGCGACCGCGTCGCATCAGGCCTACAGGTCGATGCCGATATGAACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATcc                                                                                                                        >  4:1083425/1‑149 (MQ=255)
                            ccTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATccc                                                                                                                       >  4:871219/1‑149 (MQ=255)
                              tgataagACGCGACCGCGTCGCATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCAtttttt                                                                                                                                      >  3:769631/7‑132 (MQ=255)
                              tGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCt                                                                                                                     >  1:324241/1‑149 (MQ=255)
                              tGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCt                                                                                                                     <  3:956461/149‑1 (MQ=255)
                               gataagACGCGACCGCGTCGCATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATcccc                                                                                                                      >  1:1185153/6‑147 (MQ=255)
                               gATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCTATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCtt                                                                                                                    <  3:1021658/149‑1 (MQ=255)
                                 tGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGt                                                                                                                  <  1:1176483/149‑1 (MQ=255)
                                 tGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGt                                                                                                                  <  3:196619/149‑1 (MQ=255)
                                   cGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCt                                                                                                                >  3:737779/1‑149 (MQ=255)
                                    gACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTg                                                                                                               >  1:270564/1‑149 (MQ=255)
                                       gcgGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATACCAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCGTGTCTGTaa                                                                                                            <  4:476030/149‑1 (MQ=255)
                                       gcgGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTaa                                                                                                            <  1:290397/149‑1 (MQ=255)
                                         ggTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAAtt                                                                                                          <  4:1273862/149‑1 (MQ=255)
                                            cgcgTCGCATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTaa                                                                                                       <  3:1083425/149‑1 (MQ=255)
                                             gcgtcgcaTCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAAc                                                                                                      <  4:1555489/142‑1 (MQ=255)
                                               gTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACgg                                                                                                    >  2:730639/1‑149 (MQ=255)
                                                  ttATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGaaa                                                                                                 <  3:360330/149‑1 (MQ=255)
                                                   tATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAAt                                                                                                <  1:961077/149‑1 (MQ=255)
                                                       aGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATa                                                                                            >  2:983233/1‑149 (MQ=255)
                                                          ccTACAGGTCGATGCCGATATGTACATCGTATTCGGCGATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCg                                                                                         <  2:1567386/149‑1 (MQ=255)
                                                             aCAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGt                                                                                                                                                 >  1:1484839/1‑90 (MQ=255)
                                                             aCAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGt                                                                                                                                                 <  2:1484839/90‑1 (MQ=255)
                                                             aCAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGa                                                                                                                               >  3:1150824/1‑108 (MQ=255)
                                                              cAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAg                                                                                     <  1:755779/149‑1 (MQ=255)
                                                              cAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAg                                                                                     <  4:277098/149‑1 (MQ=255)
                                                                  tCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTa                                                                                 >  4:1463149/1‑149 (MQ=255)
                                                                         cgttatGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTaa                                                                          >  1:618890/4‑149 (MQ=255)
                                                                         cGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTaa                                                                          >  3:602005/1‑149 (MQ=255)
                                                                          gATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTaaa                                                                         >  4:1104764/1‑149 (MQ=255)
                                                                            tatGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAAt                                                                       <  2:1112519/149‑1 (MQ=255)
                                                                                   aTCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGAtttt                                                                <  4:671684/149‑1 (MQ=255)
                                                                                    tCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGAttttt                                                               >  1:478772/1‑149 (MQ=255)
                                                                                      gTATTCGGCAATTAATACATGGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTaa                                                             >  2:203339/1‑149 (MQ=255)
                                                                                      gTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATa                                                                                            >  1:497854/1‑118 (MQ=255)
                                                                                      gTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATa                                                                                            <  2:497854/118‑1 (MQ=255)
                                                                                      gTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTaa                                                             <  2:1523508/149‑1 (MQ=255)
                                                                                       tATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAAt                                                            <  3:551982/149‑1 (MQ=255)
                                                                                        aTTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGt                                                                                                                                                 >  3:948732/1‑63 (MQ=255)
                                                                                         ttCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCAtttt                                                                                                                                        <  4:414934/71‑1 (MQ=255)
                                                                                          tCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCgt                                                         >  3:471676/1‑149 (MQ=255)
                                                                                           cGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCgtt                                                        <  2:298303/149‑1 (MQ=255)
                                                                                               aaTTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGttt                                                    >  4:583058/1‑149 (MQ=255)
                                                                                                      aCATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCgg                                             >  1:1365337/1‑149 (MQ=255)
                                                                                                       cATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGa                                            <  2:768939/149‑1 (MQ=255)
                                                                                                       cATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGa                                            >  1:1553289/1‑149 (MQ=255)
                                                                                                        aTAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAAt                                                                               >  4:1514716/1‑113 (MQ=255)
                                                                                                        aTAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAAt                                                                               <  3:1514716/113‑1 (MQ=255)
                                                                                                        aTAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGaa                                           >  4:1524252/1‑149 (MQ=255)
                                                                                                           gCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGa                                                                                      >  3:747305/1‑103 (MQ=255)
                                                                                                           gCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAAt                                        >  4:650783/1‑149 (MQ=255)
                                                                                                            cACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATa                                       <  4:471676/149‑1 (MQ=255)
                                                                                                               gattgattAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGc                                    >  3:200753/1‑149 (MQ=255)
                                                                                                                attgattAAATAACCTTAATAATAATGCCGACGTTATGTCGGCATTTTTTTATCAGATATATCCCCTTGTCTGAAATTTAACGGAAATCATACCGTGAGGTTAATCATAAAATAGTTTTTTAATTGTTGTTTATTTCGGAAAATagtca                                   >  4:1361367/1‑145 (MQ=255)
                                                                                                                attgattAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCa                                   >  1:623568/1‑149 (MQ=255)
                                                                                                                attgattAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCa                                   >  2:244811/1‑149 (MQ=255)
                                                                                                                attgattAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCa                                   <  2:324241/149‑1 (MQ=255)
                                                                                                                 ttgattAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAg                                  >  1:175905/1‑149 (MQ=255)
                                                                                                                    attAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGa                                                                                      >  2:936620/1‑94 (MQ=255)
                                                                                                                    attAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGa                                                                                      <  1:936620/94‑1 (MQ=255)
                                                                                                                    attAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGAtt                               >  1:749685/1‑149 (MQ=255)
                                                                                                                    attAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGAtt                               >  1:726927/1‑149 (MQ=255)
                                                                                                                    attAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGAtt                               >  2:789442/1‑149 (MQ=255)
                                                                                                                       aaaTAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAAt                            >  3:40224/1‑149 (MQ=255)
                                                                                                                        aaTAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAAtt                           >  4:1400277/1‑149 (MQ=255)
                                                                                                                         aTAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTg                          >  1:361055/1‑149 (MQ=255)
                                                                                                                          tAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGc                         <  3:391913/149‑1 (MQ=255)
                                                                                                                           aaCCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCt                        <  2:362513/149‑1 (MQ=255)
                                                                                                                             ccTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCttt                      >  2:822456/1‑149 (MQ=255)
                                                                                                                             ccTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCttt                      >  3:511547/1‑149 (MQ=255)
                                                                                                                              cTTACTAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCtttt                     <  2:670638/149‑1 (MQ=255)
                                                                                                                              cTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCtttt                     <  3:1104764/149‑1 (MQ=255)
                                                                                                                              cTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCtttt                     >  1:1191989/1‑149 (MQ=255)
                                                                                                                               ttAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTg                    >  4:1431902/1‑149 (MQ=255)
                                                                                                                               ttAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTg                    >  1:563184/1‑149 (MQ=255)
                                                                                                                                  aaaaaaaTGCCGACATAACGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGttt                 >  4:5271/6‑149 (MQ=255)
                                                                                                                                   aaaaaaTGCCGACATAACGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGtttt                >  3:851291/5‑149 (MQ=255)
                                                                                                                                     aCAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGTtttta              <  2:255642/149‑1 (MQ=255)
                                                                                                                                       aaTGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGTttttatt            >  3:937846/1‑149 (MQ=255)
                                                                                                                                        aTGCCGACATAACGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGTttttattt           >  3:1451251/1‑149 (MQ=255)
                                                                                                                                         tGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGTttttacttt          <  4:596503/149‑1 (MQ=255)
                                                                                                                                           ccGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGa                                            <  2:1399637/113‑1 (MQ=255)
                                                                                                                                           ccGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTaa        >  1:1348667/1‑149 (MQ=255)
                                                                                                                                                 ttATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTAt  <  2:468912/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTGATTAAATAACCTTAATAACAATGCCGACGTTATGTCGGCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTAT  >  NC_000913/4235435‑4235728

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