Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,798,659 (T)8→7 coding (654/1074 nt) waaU ← lipopolysaccharide core biosynthesis

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,798,6520T.93.3% 854.5 / 59.0 212coding (661/1074 nt)waaUlipopolysaccharide core biosynthesis
Reads supporting (aligned to +/- strand):  ref base T (12/2);  new base . (74/122);  total (86/126)
Fisher's exact test for biased strand distribution p-value = 9.90e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.35e-01

Only 100 of 242 total aligned reads displayed.

AATCTTGCGGTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACATCTACTGGTAAATGCAAGTCGT  >  NC_000913/3798553‑3798747
                                                                                                   |                                                                                               
 aTCTTGCGGTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATAtt                                                                                                >  1:942389/1‑100 (MQ=255)
 aTCTTGCGGTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATAtt                                                                                                >  2:2696170/1‑100 (MQ=255)
   cTTGCGGTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATAtttt                                                                                              >  1:4391401/1‑100 (MQ=255)
    ttGCGGTAACGCGGTAAATATGATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGGTATATAACTTTTATTTGCTCAAACGTAAGACGGCATAttttt                                                                                             <  2:2947134/100‑1 (MQ=255)
     tGCGGTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATAtttttt                                                                                            <  1:3762791/100‑1 (MQ=255)
      gCGGTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATAttttttt                                                                                           >  2:298543/1‑100 (MQ=255)
      gCGGTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATAttttttt                                                                                           <  1:1757926/100‑1 (MQ=255)
        ggTAACCCGGTAAATATAATCCGATAATTTTCTAAGTGTGTTTTCACTTCTTTATATATAACTTTTATTTGCTCAAACGTTAGACGGCATAtttttttgc                                                                                         <  2:666086/100‑3 (MQ=255)
         gTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATAtttttttgca                                                                                        >  2:3404324/1‑97 (MQ=255)
         gTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATAtttttttgca                                                                                        <  2:3876475/100‑4 (MQ=255)
           aaCCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCAcc                                                                                     <  2:3148812/100‑1 (MQ=255)
            aCCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCt                                                                                    >  2:2620/1‑100 (MQ=255)
             cccGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTa                                                                                   <  2:2792969/100‑1 (MQ=255)
             cccGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTa                                                                                   >  1:2886391/1‑100 (MQ=255)
              ccGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTaa                                                                                  >  1:3850736/1‑100 (MQ=255)
               cGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAAt                                                                                 <  1:1832714/100‑1 (MQ=255)
                  tAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGa                                                                              <  2:2479548/100‑1 (MQ=255)
                    aaTATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGAtt                                                                            >  1:4470609/1‑100 (MQ=255)
                      tataATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATtaa                                                                          <  1:2371417/100‑1 (MQ=255)
                       ataATCCGATGATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCACCAACGTAAGACGGCATATTTTTTTGCACCTAATGGATtaat                                                                         <  2:2236496/100‑1 (MQ=255)
                        taATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATtaata                                                                        <  2:142129/100‑1 (MQ=255)
                         aaTCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATtaataa                                                                       >  1:3491028/1‑100 (MQ=255)
                           tCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATg                                                                     <  1:1479389/100‑1 (MQ=255)
                             cGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTTCTCCTAATTGATTAATAATGAc                                                                   >  2:3533939/1‑100 (MQ=255)
                             cGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGAc                                                                   >  1:4103313/1‑100 (MQ=255)
                             cGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGAc                                                                   <  1:1471468/100‑1 (MQ=255)
                             cGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGAc                                                                   <  1:1537231/100‑1 (MQ=255)
                              gATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACt                                                                  <  1:2542197/100‑1 (MQ=255)
                                tAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCGAAAGTAAGACGGCTTATTTTTTTTCACCTAATGGATTAATAATGACTAt                                                                >  2:4420218/1‑100 (MQ=255)
                                  aTTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTc                                                              <  1:2704660/100‑1 (MQ=255)
                                  aTTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTc                                                              <  1:1390762/100‑1 (MQ=255)
                                   ttttCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCt                                                             <  1:46237/100‑1 (MQ=255)
                                   ttttCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCt                                                             <  2:2210006/100‑1 (MQ=255)
                                    tttCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCtt                                                            >  2:1770020/1‑100 (MQ=255)
                                     ttCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCATCTAATGGATTAATAATGACTATTCTTg                                                           <  2:2381616/100‑1 (MQ=255)
                                     ttCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTg                                                           <  1:3877763/100‑1 (MQ=255)
                                      tCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGt                                                          <  1:242890/100‑1 (MQ=255)
                                       cAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTTTTTTTTTTTTATTAATAATGACTATTCTTGTa                                                         >  2:3743540/1‑100 (MQ=255)
                                       cAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTa                                                         <  1:1945026/100‑1 (MQ=255)
                                       cAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTa                                                         <  1:581431/100‑1 (MQ=255)
                                       cAAAGTGTGTTTTCACTTCTTGATATATAACTCTTATTTGCTCATACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTa                                                         <  1:4190954/100‑1 (MQ=255)
                                        aaaGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTAt                                                        >  1:1913774/1‑100 (MQ=255)
                                        aaaGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTAt                                                        <  2:481715/100‑1 (MQ=255)
                                          aGTGTGTTTTCATTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGAAGGCATATATTTTTGTACCTTATGGGTTAATACTGACCATTCTTGTAtca                                                      >  2:1729407/1‑99 (MQ=255)
                                            tgtgtTTTCACTTCTTGATATATAACTTTTATTTGCTAAAACGTAAGACGGCATATTTTTTTGCCCCTAATGGATTAATAATGACTATTCTTGTATCTcc                                                    <  2:4465696/100‑1 (MQ=255)
                                               gtTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAAt                                                 >  2:3452050/1‑100 (MQ=255)
                                                ttttCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATa                                                <  2:3332592/100‑1 (MQ=255)
                                                 tttCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATaa                                               >  1:3236334/1‑100 (MQ=255)
                                                  ttCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATaaa                                              <  1:2417339/100‑1 (MQ=255)
                                                  ttCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATaaa                                              <  2:1927700/100‑1 (MQ=255)
                                                    cACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAAtt                                            >  1:3342731/1‑100 (MQ=255)
                                                        tcttGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCttt                                        >  1:4548631/1‑100 (MQ=255)
                                                        tcttGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCttt                                        <  1:573917/100‑1 (MQ=255)
                                                        tcttGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCttt                                        >  1:673625/1‑100 (MQ=255)
                                                        tcttGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCttt                                        >  1:2068724/1‑100 (MQ=255)
                                                        tcttGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTCGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAAATCttt                                        >  1:4330722/1‑100 (MQ=255)
                                                        tcttGAGATATAAATTTTATTTGCGCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCttt                                        <  1:3543603/100‑1 (MQ=255)
                                                         cttGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCtttn                                       <  1:3634702/100‑2 (MQ=255)
                                                          ttGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCtttta                                      <  1:776857/100‑1 (MQ=255)
                                                           tGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCttttat                                     >  1:3913182/1‑100 (MQ=255)
                                                            gATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCttttatt                                    >  2:3820290/1‑100 (MQ=255)
                                                             atatatAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCttttattt                                   <  1:1311311/100‑1 (MQ=255)
                                                             atatatAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCttttattt                                   <  2:3636619/100‑1 (MQ=255)
                                                             atatatAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCttttattt                                   <  2:2229287/100‑1 (MQ=255)
                                                               atataACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTTTTTTTTTTTTATTAATAATGACTATTCTTGTATCTCCAATAAATTCttttatttta                                 >  2:1968963/1‑100 (MQ=255)
                                                                tataACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTtttattttat                                >  2:943524/1‑100 (MQ=255)
                                                                tataACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTtttattttat                                >  1:983778/1‑100 (MQ=255)
                                                                  taACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCt                              <  1:1245363/100‑1 (MQ=255)
                                                                  taACTTTTAGTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCAAATAAATTCTTTTATTTTATCt                              <  2:1743176/100‑1 (MQ=255)
                                                                  taAATTTTATTTGCTTAAACGTAAGACGGCATATTTTTTTGCTCCTAATGGTTTAATAATGACTATTCTTGTATAGCCAATAAATTCTTTTATTTTATCt                              >  2:4375645/1‑100 (MQ=255)
                                                                    aCTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTc                            <  1:1620669/100‑1 (MQ=255)
                                                                     cTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCa                           <  2:426340/100‑1 (MQ=255)
                                                                     cTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCa                           <  2:402250/100‑1 (MQ=255)
                                                                     cTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTGTATCTTCa                           <  2:4298005/100‑1 (MQ=255)
                                                                     cTTTTATTTGCTCAAAAGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATTTCCAATAAATTCTTTTATTTTATTTTCa                           >  2:2401210/1‑100 (MQ=255)
                                                                       tttATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAAc                         <  2:536150/100‑1 (MQ=255)
                                                                        ttATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACa                        <  1:1775457/100‑1 (MQ=255)
                                                                         tATTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACAt                       >  2:3832807/1‑100 (MQ=255)
                                                                          aTTTGCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACATc                      >  1:649100/1‑100 (MQ=255)
                                                                          aTTTGCTCAAACGTAAGACGGCATA‑‑TTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACATCt                     >  1:2257925/1‑100 (MQ=255)
                                                                              gCTCAAACGTAAGACGGCATATTTTTTTTTTTTTTTTTTAATAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCATCATCTACt                  >  2:767165/1‑100 (MQ=255)
                                                                              gCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACATCTACt                  <  1:3062676/100‑1 (MQ=255)
                                                                              gCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACATCTACt                  <  1:2944057/100‑1 (MQ=255)
                                                                              gCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACATCTACt                  <  1:3276186/100‑1 (MQ=255)
                                                                              gCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACATCTACt                  <  2:1594762/100‑1 (MQ=255)
                                                                              gCTCAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAAGTATGACTATTCTTGTATCACCAATAAATTCTTTTATTTTCTCTTCAACATCTACt                  >  2:455773/1‑100 (MQ=255)
                                                                                 cAAACGTAAGACGGCATATTTTTTTTTTTTTTTTTAATTAATAATGACTATTCTAGTATCTCCAATAAATTGTTTTATTTTATCTTCAACTTCTatttgt               >  2:926435/1‑95 (MQ=255)
                                                                                 cAAACGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACATCTACTGGt               >  1:1213126/1‑100 (MQ=255)
                                                                                   aaCGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACATCTACTGGTaa             <  2:1713369/100‑1 (MQ=255)
                                                                                    aCGTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACATCTACTGGTaaa            <  1:265814/100‑1 (MQ=255)
                                                                                      gTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACATCTACTGGTAAATg          <  2:1803083/100‑1 (MQ=255)
                                                                                      gTAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACATCTACTGGTAAATg          <  2:3189865/100‑1 (MQ=255)
                                                                                      atAAGACGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACATCTGCTGGTAGATg          <  2:4115490/99‑1 (MQ=255)
                                                                                          gACGGCATATTTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACATCTACTGGTAAATGCaa       <  1:2461667/100‑1 (MQ=255)
                                                                                          gACGGCATATTTTTGTGCACCTAATGCATTAATAATGACTATTCTTGAATCTCCCATAAATTCTTTTATTTTATCTTCAACATCTACTGGTAAATGCAAg      >  2:4193075/1‑100 (MQ=255)
                                                                                           aCGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACATCTACTGGTAAATGCAAGt     <  1:1022516/100‑1 (MQ=255)
                                                                                            cGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACATCTACTGGTAAATGCAAGTc    <  1:261091/100‑1 (MQ=255)
                                                                                            cGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACATCTACTGGTAAATGCAAGTc    >  1:1087732/1‑100 (MQ=255)
                                                                                            cGGCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACATCTACTGGTAAATGCAAGTc    <  2:3011878/100‑1 (MQ=255)
                                                                                             ggCATATTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATGAATTCTTTTATTTTATCTTCAACATCTACTGGTACATGCAAGTCg   >  1:310566/1‑100 (MQ=255)
                                                                                                   |                                                                                               
AATCTTGCGGTAACCCGGTAAATATAATCCGATAATTTTCAAAGTGTGTTTTCACTTCTTGATATATAACTTTTATTTGCTCAAACGTAAGACGGCATATTTTTTTTGCACCTAATGGATTAATAATGACTATTCTTGTATCTCCAATAAATTCTTTTATTTTATCTTCAACATCTACTGGTAAATGCAAGTCGT  >  NC_000913/3798553‑3798747

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