Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,121,583 Δ1 bp 6.4% intergenic (‑29/+26) wcaJ ← / ← cpsG colanic biosynthesis UDP‑glucose lipid carrier transferase/phosphomannomutase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,121,5830G. 6.4% 327.5 / 13.5 93intergenic (‑29/+26)wcaJ/cpsGcolanic biosynthesis UDP‑glucose lipid carrier transferase/phosphomannomutase
Reads supporting (aligned to +/- strand):  ref base G (45/42);  new base . (3/3);  total (48/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 = 1.00e+00

CTGACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTAAAGCGCCA  >  NC_000913/2121436‑2121731
                                                                                                                                                   |                                                                                                                                                     
cTGACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCg                                                                                                                                                       >  2:122756/1‑148 (MQ=255)
cTGACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCg                                                                                                                                                       >  2:860676/1‑148 (MQ=255)
 tGACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCGTTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCa                                                                                                                                                     >  3:548586/1‑149 (MQ=255)
 tGACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCa                                                                                                                                                     >  4:441141/1‑149 (MQ=255)
  gACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAg                                                                                                                                                    <  3:422844/149‑1 (MQ=255)
     ttCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTgg                                                                                                                                                >  4:313426/1‑149 (MQ=255)
      tCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTgg                                                                                                                                                <  2:529655/148‑1 (MQ=255)
      tCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTgg                                                                                                                                                >  2:114446/1‑148 (MQ=255)
           aaCCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAg                                                                                                                                          >  3:475027/1‑149 (MQ=255)
            aCCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAgg                                                                                                                                         <  3:701513/149‑1 (MQ=255)
            aCCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAgg                                                                                                                                         >  3:310739/1‑149 (MQ=255)
             ccAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCCACTGCGGGTAAggg                                                                                                                                        <  3:759390/149‑1 (MQ=255)
              cAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCCACTGCGGGTAAggg                                                                                                                                        <  1:252239/148‑1 (MQ=255)
                 ccATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTgg                                                                                                                                                <  4:135441/137‑1 (MQ=255)
                 ccATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTgg                                                                                                                                                >  3:135441/1‑137 (MQ=255)
                 ccATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAg                                                                                                                                     <  1:778649/148‑1 (MQ=255)
                 ccATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGa                                                                                                                                    >  4:652347/1‑149 (MQ=255)
                  cATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGAt                                                                                                                                   >  3:566029/1‑149 (MQ=255)
                   aTAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATc                                                                                                                                  <  3:326597/149‑1 (MQ=255)
                      gTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCTCTAGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCgt                                                                                                                               <  4:784704/149‑2 (MQ=255)
                       tCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGAc                                                                                                                              >  4:852811/1‑149 (MQ=255)
                          ggCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACAtt                                                                                                                           <  4:132191/149‑1 (MQ=255)
                            cAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACGTTa                                                                                                                          <  2:222970/148‑1 (MQ=255)
                              aaCATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACt                                                                                                                        <  1:679623/148‑1 (MQ=255)
                              aaCATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACt                                                                                                                        <  2:859818/148‑1 (MQ=255)
                                 atgatgGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGtt                                                                                                                    <  3:740308/149‑1 (MQ=255)
                                 atgatgGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGt                                                                                                                     <  2:250143/148‑1 (MQ=255)
                                   gatgGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTc                                                                                                                   >  2:446175/1‑148 (MQ=255)
                                       gTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCaa                                                                                                              <  4:144013/149‑1 (MQ=255)
                                             tCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTTTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGGAGATCCGACATTACTCGTTCATCAACGTcag                                                                                                        >  4:535206/1‑149 (MQ=255)
                                             tCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTcag                                                                                                        >  2:110015/1‑148 (MQ=255)
                                                gAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAg                                                                                                                                          >  1:263211/1‑112 (MQ=255)
                                                gAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAg                                                                                                                                          <  2:263211/112‑1 (MQ=255)
                                                gAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGaa                                                                                                                                      <  3:468194/116‑1 (MQ=255)
                                                gAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGaa                                                                                                                                      >  4:468194/1‑116 (MQ=255)
                                                gAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATAACTCGTTCAGCAACGTcagcag                                                                                                     >  3:778284/1‑149 (MQ=255)
                                                 aaaaGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCaga                                                                                                    <  4:564930/149‑1 (MQ=255)
                                                     gCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAgttc                                                                                                <  3:471496/149‑1 (MQ=255)
                                                                  gagaTTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTc                                                                                   <  3:9669/149‑1 (MQ=255)
                                                                      ttAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATc                                                                               <  3:62153/149‑1 (MQ=255)
                                                                       tAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCa                                                                              <  4:566029/149‑1 (MQ=255)
                                                                             tGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGc                                                                         <  2:457945/148‑1 (MQ=255)
                                                                                aTTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACAt                                                                     >  3:734103/1‑149 (MQ=255)
                                                                                aTTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACAt                                                                     >  3:386352/1‑149 (MQ=255)
                                                                                 ttGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAgttcgcgttcgc                                                                                        <  3:888030/129‑1 (MQ=255)
                                                                                 ttGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAgttcgcgttcgc                                                                                        >  4:888030/1‑129 (MQ=255)
                                                                                 ttGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACAt                                                                     >  2:259132/1‑148 (MQ=255)
                                                                                 ttGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATc                                                                    >  3:653989/1‑149 (MQ=255)
                                                                                     ttttCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCTATCAGCGGCACATCACcg                                                                >  3:341429/1‑149 (MQ=255)
                                                                                          gctcgctcgcGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACcgcgcg                                                            <  1:446175/148‑1 (MQ=255)
                                                                                          gctcgctcgcGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGa                                                           >  4:606709/1‑149 (MQ=255)
                                                                                           ctcgctcgcGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGAt                                                          >  3:568346/1‑149 (MQ=255)
                                                                                              gctcgcGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTc                                                        <  2:270262/148‑1 (MQ=255)
                                                                                                 cgcgCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCcaca                                                    <  3:852811/149‑1 (MQ=255)
                                                                                                    gcTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTCCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTc                                                 <  3:487070/149‑1 (MQ=255)
                                                                                                     cTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCa                                                <  1:110015/148‑1 (MQ=255)
                                                                                                      ttttttAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAg                                               <  3:313426/149‑1 (MQ=255)
                                                                                                       tttttAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAg                                               >  1:790059/1‑148 (MQ=255)
                                                                                                          ttAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGgcgc                                           >  3:758271/1‑149 (MQ=255)
                                                                                                             gATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACAt                                                                     >  1:550715/1‑120 (MQ=255)
                                                                                                             gATTTGTCATCGTTGTTCCTGTGATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACAt                                                                     <  2:550715/120‑1 (MQ=255)
                                                                                                                ttGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGcaccac                                      >  2:428737/1‑148 (MQ=255)
                                                                                                                     aTCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGtt                                 <  1:860676/148‑1 (MQ=255)
                                                                                                                     aTCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTc                                >  3:500671/1‑149 (MQ=255)
                                                                                                                     aTCGTTGTTCCTGTTATTAGCCCCTTACCCCACTGCGGGTAAgg                                                                                                                                         >  3:334976/1‑44 (MQ=17)
                                                                                                                     aTCGTTGTTCCTGTTATTAGCCCCTTACCCCACTGCGGGTAAgg                                                                                                                                         <  4:334976/44‑1 (MQ=17)
                                                                                                                      tCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTc                                >  2:404295/1‑148 (MQ=255)
                                                                                                                      tCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCg                               >  3:95973/1‑149 (MQ=255)
                                                                                                                       cGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGGACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCgg                              <  4:707987/149‑1 (MQ=255)
                                                                                                                       cGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCgg                              <  4:386352/149‑1 (MQ=255)
                                                                                                                       cGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCg                               >  2:575324/1‑148 (MQ=255)
                                                                                                                           gttCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTAt                           <  2:790059/148‑1 (MQ=255)
                                                                                                                            ttCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTg                         >  3:5321/1‑149 (MQ=255)
                                                                                                                             tCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTg                         >  1:435047/1‑148 (MQ=255)
                                                                                                                                tGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAgg                     <  3:606709/149‑1 (MQ=255)
                                                                                                                                 gTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCa                                                >  1:745639/1‑121 (MQ=255)
                                                                                                                                 gTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCa                                                <  2:745639/121‑1 (MQ=255)
                                                                                                                                  ttattaGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTgg                        <  1:675067/144‑1 (MQ=255)
                                                                                                                                  ttattaGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTgg                        >  2:675067/1‑144 (MQ=255)
                                                                                                                                   tattaGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTg                   >  2:761184/1‑148 (MQ=255)
                                                                                                                                         tcccTTACCCCACTGCGGGTAAGGGGAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAgg             >  2:144718/2‑148 (MQ=255)
                                                                                                                                         tcccTTACCCCACTGCGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGt            >  3:603080/2‑149 (MQ=255)
                                                                                                                                         ccccTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAgg             >  2:819731/1‑148 (MQ=255)
                                                                                                                                         ccccTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAg              <  3:864866/147‑1 (MQ=255)
                                                                                                                                         ccccTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAg              >  4:864866/1‑147 (MQ=255)
                                                                                                                                         ccccTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGt            >  3:549206/1‑149 (MQ=255)
                                                                                                                                         ccccTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGt            >  4:722715/1‑149 (MQ=255)
                                                                                                                                           ccTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCGGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTcc                                                       <  4:126478/104‑1 (MQ=255)
                                                                                                                                           ccTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCGGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTcc                                                       >  3:126478/1‑104 (MQ=255)
                                                                                                                                           ccTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTa          >  3:379995/1‑149 (MQ=255)
                                                                                                                                            cTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTa          >  1:1918/1‑148 (MQ=255)
                                                                                                                                             ttACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTaaa        <  4:475027/149‑1 (MQ=255)
                                                                                                                                               aCCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTAAAgc      <  4:653989/149‑1 (MQ=255)
                                                                                                                                               aCCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGAGTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTAAAgc      <  4:548586/149‑1 (MQ=255)
                                                                                                                                                 ccGCAGTGGGGAAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACGTTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTAAAgcgc    >  4:81992/1‑149 (MQ=255)
                                                                                                                                                   gCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTAAAGCGCCa  <  4:310739/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CTGACTTCGCAAACCAGCCATAGTCCGGCAAACATGATGGTGATATCTGAAAAGCGTTGCACCATAGAGATTAACGATGCATTGGTTTTCGCTCGCTCGCGCTTTTTTAGATTTGTCATCGTTGTTCCTGTTATTAGCCCCTTACCCGCAGTGGGGTAAGGGAAGATCCGACATTACTCGTTCAGCAACGTCAGCAGAGTTCGCGTTCGCGCTTCCATCAGCGGCACATCACCGCGCGATTCCACATTCAGGCGCACCACCGGTTCGGTATTGGAGGTGCGCAGGTTAAAGCGCCA  >  NC_000913/2121436‑2121731

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: