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)
             tCATCTTTTTTAATATCATTGAGGTAATGCTTGATGAATAATCGCCCCAGCCATGCATTATCCATCCCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTCAGGTACTg                                                                                                                                    >  2:592728/1‑151 (MQ=255)
              cATCTTTTTTAATATCATTGAGGTAATGCTTGATGAATAATCGCCCCAGCCATGCATTATCCATCCCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTCAGGTACTg                                                                                                                                    <  1:917676/150‑1 (MQ=255)
               aTCTTTTTTAATATCATTGAGGTAATGCTTGATGAATAATCGCCCCAGCCATGCATTATCCATCCCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTCAGGTACTGAAt                                                                                                                                 >  1:79772/1‑151 (MQ=255)
                   tttttAATATCATTGAGGTAATGCTTGATGAATAATCGCCCCAGCCATGCATTATCCATCCCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTCAGGTACTGAATCAt                                                                                                                              >  2:334275/1‑151 (MQ=255)
                        tatatCATTGAGGTAATGCTTGATGAATAATCGTCCCAGCCATTCCTTATCCATCCCTATTTTTGATATTATTTTATTAATCCAGTTATATTTTTTTGATGTTGTTATATTTTTTATCTTGTACTTTTTTTCAGGTACTGAATCAttttt                                                                                                                          >  1:184300/2‑150 (MQ=255)
                            tCATTGAGGTAATGCTTGATGAATAATCGCCCCAGCCATGCATTATCCATCCCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTCAGGTACTGAATCATTTTTAATTc                                                                                                                     >  2:911150/1‑151 (MQ=255)
                                  aGGTAATGCTTGATGAATAATCGCCCCAGCCATGCATTATCCATCCCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATc                                                                                                               >  2:289395/1‑151 (MQ=255)
                                    gTAATGCTTGATGAATAATCGCCCCAGCCATGCATTATCCATCCCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTg                                                                                                             >  2:247734/1‑151 (MQ=255)
                                      aaTGCTTGATGAATAATCGCCCCAGCCATGCATTATCCATCCCTATTTTTTATATTTTTTTATTAATCCAGTTATATTTTTTGGTGTTGTTATATTTTTTATCTTGTAATTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCaa                                                                                                          <  1:56242/151‑1 (MQ=255)
                                      aaTGCTTGATGAATAATCGCCCCAGCCATGCATTATCCATCCCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGc                                                                                                            <  2:892976/150‑1 (MQ=255)
                                      aaTGCTTGATGAATAATCGCCCCAGCCATGCATTATCCATCCCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGc                                                                                                            <  1:861792/150‑1 (MQ=255)
                                           ttGATGAATAATCGCCCCAGCCATGCATTATCCATCCCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAAccc                                                                                                       >  2:902967/1‑150 (MQ=255)
                                                gaATAGTCGCCCCAGCCATGCATTATCCATCCCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTta                                                                                                 <  2:379332/151‑1 (MQ=255)
                                                 aataatCGCCCCAGCCATGCATTATCCATCCCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTtat                                                                                                >  2:466350/1‑151 (MQ=255)
                                                    aatCGCCCCAGCCATGCATTATCCATCCCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTataaa                                                                                             <  2:592189/150‑1 (MQ=255)
                                                       caccccAGCCATGCATTATCCATCCCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTCAAGGACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTattaaa                                                                                            >  2:854877/3‑149 (MQ=255)
                                                          cccAGCCATGCATTATCCATCCCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAaaataaat                                                                                       >  2:104163/1‑151 (MQ=255)
                                                          cccAGCCATGCATTATCCATCCCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAaaataaat                                                                                       >  2:63195/1‑151 (MQ=255)
                                                                 aTGCATTATCCATCCCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCata                                                                                <  2:745680/151‑1 (MQ=255)
                                                                 aTGCATTATCCATCCCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCat                                                                                 >  2:547180/1‑150 (MQ=255)
                                                                     aTTATCCATCCCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCAtaataa                                                                            <  2:343171/151‑1 (MQ=255)
                                                                        atccatccCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAagt                                                                         <  1:13182/151‑1 (MQ=255)
                                                                        atccatccCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAag                                                                          <  1:96445/150‑1 (MQ=255)
                                                                          ccatccCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAagtt                                                                        >  2:227683/1‑150 (MQ=255)
                                                                           catccCTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAagtt                                                                        >  2:74254/1‑149 (MQ=255)
                                                                               ccgaTTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGGGTTGTTATATTTTTTATCTTGTACTTTTTTTCAGGTAGTGTATCGTTTTTAATTCATAATCTGCACGCCATTTTATAAAATAAATAGTCATAATAAAGTTAGTTcc                                                                  >  1:783/4‑151 (MQ=255)
                                                                               ccTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTcc                                                                  >  2:326074/1‑151 (MQ=255)
                                                                                cTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATTTTGTACTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTcc                                                                  >  1:184892/1‑150 (MQ=255)
                                                                                cTATTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTcc                                                                  >  2:596037/1‑150 (MQ=255)
                                                                                  aTTTTTGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAg                                                                <  2:198858/150‑1 (MQ=255)
                                                                                   tttttGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAg                                                                >  1:576023/1‑149 (MQ=255)
                                                                                    ttttGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTGCTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTa                                                              <  1:158148/150‑1 (MQ=255)
                                                                                     tttGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTAc                                                             <  1:226925/150‑1 (MQ=255)
                                                                                       tGATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACAt                                                           <  1:720495/150‑1 (MQ=255)
                                                                                        gATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACATa                                                          <  2:345570/150‑1 (MQ=255)
                                                                                        gATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACATAc                                                         >  2:391836/1‑151 (MQ=255)
                                                                                        gATATTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACATAc                                                         >  1:117089/1‑151 (MQ=255)
                                                                                           atTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCGAGTACATACta                                                       <  2:350350/150‑1 (MQ=255)
                                                                                           atTTTTTTATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACATACta                                                       >  1:236426/1‑150 (MQ=255)
                                                                                             ttttttATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACATACtaat                                                     <  1:114132/150‑1 (MQ=255)
                                                                                             ttttttATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACATACtaat                                                     >  2:337178/1‑150 (MQ=255)
                                                                                               ttttATTAATCCAGTTATATTTTTTTGGTGTTGTTATATTTGTTATCTTGTACTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACATACtaataa                                                   <  1:891340/150‑1 (MQ=255)
                                                                                                   attaATCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACATACTAATAAATAtt                                              >  2:470216/1‑151 (MQ=255)
                                                                                                      aaTCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACATACTAATAAATAttttt                                           <  2:859935/151‑1 (MQ=255)
                                                                                                      aaTCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACATACTAATAAATAttttt                                           >  1:664656/1‑151 (MQ=255)
                                                                                                      aaTCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACATACTAATAAATAtttt                                            <  1:547545/150‑1 (MQ=255)
                                                                                                       aTCCAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACATACTAATAAATAttttt                                           >  1:904405/1‑150 (MQ=255)
                                                                                                         ccAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTTCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACATACTAATAAATAtttat                                           <  1:194229/148‑3 (MQ=255)
                                                                                                         ccAGTTATATTTTTTTGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACATACTAATAAATATTTTtat                                         <  2:690954/150‑1 (MQ=255)
                                                                                                              tataTTTTTTTGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACATACTAATAAATATTTTTATAATCCt                                   >  2:187721/1‑151 (MQ=255)
                                                                                                                 atTTTTTTGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACATACTAATAAATATTTTTATAATCCTTTg                                <  2:318249/151‑1 (MQ=255)
                                                                                                                 atTTTTTTGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACATACTAATAAATATTTTTATAATCCTTTg                                <  2:600976/151‑1 (MQ=255)
                                                                                                                       ttGGTGTTGTTATATTTTTTATCTTGTACTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTtg                           <  1:203879/150‑1 (MQ=255)
                                                                                                                         ggTGTTGTTATATTTTTTATCTTGTACTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTgtgt                        <  2:434760/151‑1 (MQ=255)
                                                                                                                          gtgtTGTTATATTTTTTATCTTGTACTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACATCCTAATACATATTTTTATAATCCTTTGCGTtgtg                         <  1:51360/149‑1 (MQ=255)
                                                                                                                              tgttATATTTTTTATCTTGTACTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGc                    >  1:89995/1‑150 (MQ=255)
                                                                                                                              tgttATATTTTTTATCTTGTACTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGc                    <  1:361936/150‑1 (MQ=255)
                                                                                                                                ttATATTTTTTATCTTGTACTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGctct                 >  1:413823/1‑151 (MQ=255)
                                                                                                                                  atatTTTTTATCTTGTACTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCttt               <  1:106224/151‑1 (MQ=255)
                                                                                                                                               tGTACTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGcc  >  2:258156/1‑151 (MQ=255)
                                                                                                                                               tGTACTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGc   <  1:289375/150‑1 (MQ=255)
                                                                                                                                               tGTACTTTTTTTCAGGTACTGAATCATTTTTAATTCATAATCTGCAACCCAATTTATAAAATAAATAGTCATAATAAAGTTAGTTCCAGTACATACTAATAAATATTTTTATAATCCTTTGCGTTGTGTTCGCTCTTTAATACGCTCGGc   <  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: