Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,805,056 (T)8→7 coding (61/936 nt) waaS ← lipopolysaccharide rhamnose:KdoIII transferase; lipopolysaccharide core biosynthesis protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,805,0490T.98.4% 291.2 / ‑3.4 62coding (68/936 nt)waaSlipopolysaccharide rhamnose:KdoIII transferase; lipopolysaccharide core biosynthesis protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base . (31/30);  minor base A (0/1);  total (31/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

AATGATAATTATGTCATCTTTTTTAATATCATTGAGGTAATGCTTGATGAATAATCGCCCCAGCCATGCATTATCCATCCCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCC  >  NC_000913/3804901‑3805195
                                                                                                                                                    |                                                                                                                                                  
aaTGATAATTATGTCATCTTTTTTAATATCATTGAGGTAATGCTTGATGAATAATCGCCCCAGCCATGCATTATCCATCCCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACttt                                                                                                                                                  <  2:168975/151‑1 (MQ=255)
 aTGATAATTATGTCATCTTTTTTAATATCATTGAGGTAATGCTTGATGAATAATCGCCCCAGCCATGCATTATCCATCCCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACtttt                                                                                                                                                 >  2:898662/1‑151 (MQ=255)
  tGATAATTATGTCATCTTTTTTAATATCATTGAGGTAATGCTTGATGAATAATCGCCCCAGCCATGCATTATCCATCCCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACttttt                                                                                                                                                >  1:600942/1‑151 (MQ=255)
       aTTATGTCATCTTTTTTAATATCATTGAGGTAATGCTTGATGAATAATCGCCCCAGCCATGCATTATCCATCCCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACtttttttcag                                                                                                                                           <  1:630054/151‑4 (MQ=255)
             tCATCTTTTTTAATATCATTGAGGTAATGCTTGATGAATAATCGCCCCAGCCATGCATTATCCATCCCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTAC‑TTTTTTTCAGGTACTg                                                                                                                                    >  2:592728/1‑151 (MQ=255)
              cATCTTTTTTAATATCATTGAGGTAATGCTTGATGAATAATCGCCCCAGCCATGCATTATCCATCCCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTAC‑TTTTTTTCAGGTACTg                                                                                                                                    <  1:917676/150‑1 (MQ=255)
               aTCTTTTTTAATATCATTGAGGTAATGCTTGATGAATAATCGCCCCAGCCATGCATTATCCATCCCTATTTTTGATATTTTTTTATTAATCCAGTTATA‑TTTTTTGGTGTTGTTATATTTTTTATCTTGTAC‑TTTTTTTCAGGTACTGAAt                                                                                                                                 >  1:79772/1‑151 (MQ=255)
                   tttttAATATCATTGAGGTAATGCTTGATGAATAATCGCCCCAGCCATGCATTATCCATCCCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTAC‑TTTTTTTCAGGTACTGAATCAt                                                                                                                              >  2:334275/1‑151 (MQ=255)
                        tatatCATTGAGGTAATGCTTGATGAATAATCGTCCCAGCCATTCCTTATCCATCCCTATTTTTGATATTATTTTATTAATCCAGTTATATTTTTTTGATGTTGTTATATTTTTTATCTTGTAC‑TTTTTTTCAGGTACTGAATCAttttt                                                                                                                          >  1:184300/2‑150 (MQ=255)
                            tCATTGAGGTAATGCTTGATGAATAATCGCCCCAGCCATGCATTATCCATCCCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTAC‑TTTTTTTCAGGTACTGAATCATTTTTAATTc                                                                                                                     >  2:911150/1‑151 (MQ=255)
                                  aGGTAATGCTTGATGAATAATCGCCCCAGCCATGCATTATCCATCCCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTAC‑TTTTTTTCAGGTACTGAATCATTTTTAATTCATAATc                                                                                                               >  2:289395/1‑151 (MQ=255)
                                    gTAATGCTTGATGAATAATCGCCCCAGCCATGCATTATCCATCCCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTAC‑TTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTg                                                                                                             >  2:247734/1‑151 (MQ=255)
                                      aaTGCTTGATGAATAATCGCCCCAGCCATGCATTATCCATCCCTATTTTTTATATTTTTTTATTAATCCAGTTATA‑TTTTTTGGTGTTGTTATATTTTTTATCTTGTA‑ATTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCaa                                                                                                          <  1:56242/151‑1 (MQ=255)
                                      aaTGCTTGATGAATAATCGCCCCAGCCATGCATTATCCATCCCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTAC‑TTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGc                                                                                                            <  2:892976/150‑1 (MQ=255)
                                      aaTGCTTGATGAATAATCGCCCCAGCCATGCATTATCCATCCCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTAC‑TTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGc                                                                                                            <  1:861792/150‑1 (MQ=255)
                                           ttGATGAATAATCGCCCCAGCCATGCATTATCCATCCCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTAC‑TTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAAccc                                                                                                       >  2:902967/1‑150 (MQ=255)
                                                gaATAGTCGCCCCAGCCATGCATTATCCATCCCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTAC‑TTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTta                                                                                                 <  2:379332/151‑1 (MQ=255)
                                                 aataatCGCCCCAGCCATGCATTATCCATCCCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTAC‑TTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTtat                                                                                                >  2:466350/1‑151 (MQ=255)
                                                    aatCGCCCCAGCCATGCATTATCCATCCCTATTTTTGATATTTTTTTATTAATCCAGTTATA‑TTTTTTGGTGTTGTTATATTTTTTATCTTGTAC‑TTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTataaa                                                                                             <  2:592189/150‑1 (MQ=255)
                                                       caccccAGCCATGCATTATCCATCCCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTAC‑TTTTTTTCAAGGACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTattaaa                                                                                            >  2:854877/3‑149 (MQ=255)
                                                          cccAGCCATGCATTATCCATCCCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTAC‑TTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAaaataaat                                                                                       >  2:104163/1‑151 (MQ=255)
                                                          cccAGCCATGCATTATCCATCCCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTAC‑TTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAaaataaat                                                                                       >  2:63195/1‑151 (MQ=255)
                                                                 aTGCATTATCCATCCCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTAC‑TTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCata                                                                                <  2:745680/151‑1 (MQ=255)
                                                                 aTGCATTATCCATCCCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTAC‑TTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCat                                                                                 >  2:547180/1‑150 (MQ=255)
                                                                     aTTATCCATCCCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTAC‑TTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCAtaataa                                                                            <  2:343171/151‑1 (MQ=255)
                                                                        atccatccCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTAC‑TTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAagt                                                                         <  1:13182/151‑1 (MQ=255)
                                                                        atccatccCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTAC‑TTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAag                                                                          <  1:96445/150‑1 (MQ=255)
                                                                          ccatccCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTAC‑TTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAagtt                                                                        >  2:227683/1‑150 (MQ=255)
                                                                           catccCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTAC‑TTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAagtt                                                                        >  2:74254/1‑149 (MQ=255)
                                                                               ccgaTTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGGGTTGTTATATTTTTTATCTTGTAC‑TTTTTTTCAGGTAGTGTATCGTTTTTAATTCATAATCTGCACGCCATTTTATAAAATAAATAGTCATAATAAAGTTAGTTcc                                                                  >  1:783/4‑151 (MQ=255)
                                                                               ccTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTAC‑TTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTcc                                                                  >  2:326074/1‑151 (MQ=255)
                                                                                cTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATTTTGTAC‑TTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTcc                                                                  >  1:184892/1‑150 (MQ=255)
                                                                                cTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTAC‑TTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTcc                                                                  >  2:596037/1‑150 (MQ=255)
                                                                                  aTTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTAC‑TTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAg                                                                <  2:198858/150‑1 (MQ=255)
                                                                                   tttttGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTAC‑TTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAg                                                                >  1:576023/1‑149 (MQ=255)
                                                                                    ttttGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTGC‑TTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTa                                                              <  1:158148/150‑1 (MQ=255)
                                                                                     tttGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTAC‑TTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTAc                                                             <  1:226925/150‑1 (MQ=255)
                                                                                       tGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTAC‑TTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACAt                                                           <  1:720495/150‑1 (MQ=255)
                                                                                        gATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTAC‑TTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACATa                                                          <  2:345570/150‑1 (MQ=255)
                                                                                        gATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTAC‑TTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACATAc                                                         >  2:391836/1‑151 (MQ=255)
                                                                                        gATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTAC‑TTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACATAc                                                         >  1:117089/1‑151 (MQ=255)
                                                                                           atTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTAC‑TTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCGAGTACATACta                                                       <  2:350350/150‑1 (MQ=255)
                                                                                           atTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTAC‑TTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACATACta                                                       >  1:236426/1‑150 (MQ=255)
                                                                                             ttttttATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTAC‑TTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACATACtaat                                                     <  1:114132/150‑1 (MQ=255)
                                                                                             ttttttATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTAC‑TTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACATACtaat                                                     >  2:337178/1‑150 (MQ=255)
                                                                                               ttttATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTGTTATCTTGTAC‑TTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACATACtaataa                                                   <  1:891340/150‑1 (MQ=255)
                                                                                                   attaATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTAC‑TTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACATACTAATAAATAtt                                              >  2:470216/1‑151 (MQ=255)
                                                                                                      aaTCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTAC‑TTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACATACTAATAAATAttttt                                           <  2:859935/151‑1 (MQ=255)
                                                                                                      aaTCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTAC‑TTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACATACTAATAAATAttttt                                           >  1:664656/1‑151 (MQ=255)
                                                                                                      aaTCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTAC‑TTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACATACTAATAAATAtttt                                            <  1:547545/150‑1 (MQ=255)
                                                                                                       aTCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTAC‑TTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACATACTAATAAATAttttt                                           >  1:904405/1‑150 (MQ=255)
                                                                                                         ccAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTAC‑TTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTTCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACATACTAATAAATAtttat                                           <  1:194229/148‑3 (MQ=255)
                                                                                                         ccAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTAC‑TTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACATACTAATAAATATTTTtat                                         <  2:690954/150‑1 (MQ=255)
                                                                                                              tataTTTTTTTGGTGTTGTTATATTTTTTATCTTGTAC‑TTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACATACTAATAAATATTTTTATAATCCt                                   >  2:187721/1‑151 (MQ=255)
                                                                                                                 atTTTTTTGGTGTTGTTATATTTTTTATCTTGTAC‑TTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACATACTAATAAATATTTTTATAATCCTTTg                                <  2:318249/151‑1 (MQ=255)
                                                                                                                 atTTTTTTGGTGTTGTTATATTTTTTATCTTGTAC‑TTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACATACTAATAAATATTTTTATAATCCTTTg                                <  2:600976/151‑1 (MQ=255)
                                                                                                                       ttGGTGTTGTTATATTTTTTATCTTGTAC‑TTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTtg                           <  1:203879/150‑1 (MQ=255)
                                                                                                                         ggTGTTGTTATATTTTTTATCTTGTAC‑TTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTgtgt                        <  2:434760/151‑1 (MQ=255)
                                                                                                                          gtgtTGTTATATTTTTTATCTTGTAC‑TTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACATCCTAATACATATTTTTATAATCCTTTGCGTtgtg                         <  1:51360/149‑1 (MQ=255)
                                                                                                                              tgttATATTTTTTATCTTGTAC‑TTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGc                    >  1:89995/1‑150 (MQ=255)
                                                                                                                              tgttATATTTTTTATCTTGTAC‑TTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGc                    <  1:361936/150‑1 (MQ=255)
                                                                                                                                ttATATTTTTTATCTTGTAC‑TTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGctct                 >  1:413823/1‑151 (MQ=255)
                                                                                                                                  atatTTTTTATCTTGTAC‑TTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCttt               <  1:106224/151‑1 (MQ=255)
                                                                                                                                               tGTAC‑TTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGcc  >  2:258156/1‑151 (MQ=255)
                                                                                                                                               tGTAC‑TTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGc   <  1:289375/150‑1 (MQ=255)
                                                                                                                                               tGTAC‑TTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGc   <  1:1711/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
AATGATAATTATGTCATCTTTTTTAATATCATTGAGGTAATGCTTGATGAATAATCGCCCCAGCCATGCATTATCCATCCCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGCC  >  NC_000913/3804901‑3805195

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑