Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,805,475 T→C N159S (AAT→AGT)  waaP ← kinase that phosphorylates core heptose of lipopolysaccharide

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,805,4750TC99.2% 473.0 / NA 134N159S (AAT→AGT) waaPkinase that phosphorylates core heptose of lipopolysaccharide
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (58/75);  minor base A (1/0);  total (59/75)
Fisher's exact test for biased strand distribution p-value = 4.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 134 total aligned reads displayed.

CCCGATGCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGATTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGTATATCAGGCGGGTTGACTGCCCA  >  NC_000913/3805380‑3805566
                                                                                               |                                                                                           
cccGATGCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTAt                                                                                         >  2:2894899/1‑100 (MQ=255)
  cGATGCAGATCGATCACTGAAATTTTTAATTCATCTTCCCGGCCAGTACATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATc                                                                                        >  1:188786/1‑99 (MQ=255)
  cGATGCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATcc                                                                                       >  2:1776578/1‑100 (MQ=255)
  cGATGCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATcc                                                                                       <  1:93496/100‑1 (MQ=255)
   gATGCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATccc                                                                                      <  1:427309/100‑1 (MQ=255)
   gATGCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATccc                                                                                      <  1:1734480/100‑1 (MQ=255)
     tGCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTg                                                                                    <  2:960770/100‑1 (MQ=255)
      gCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGc                                                                                   <  1:3109287/100‑1 (MQ=255)
      gCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGc                                                                                   <  2:2766800/100‑1 (MQ=255)
       cAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCa                                                                                  <  1:1909577/100‑1 (MQ=255)
        aGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAg                                                                                 <  2:293695/100‑1 (MQ=255)
         gATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGt                                                                                >  2:193818/1‑100 (MQ=255)
            cGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATg                                                                             >  2:1533765/1‑100 (MQ=255)
            cGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATg                                                                             <  2:2287178/100‑1 (MQ=255)
            cGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATATAACAGTCGCGGTGACTTATCCCTGCAGTATg                                                                             <  2:3291820/100‑1 (MQ=255)
              aTAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCa                                                                           <  2:2825126/100‑1 (MQ=255)
                aaCTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCAtt                                                                         >  2:2488486/1‑100 (MQ=255)
                aaCTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCAtt                                                                         <  2:2316710/100‑1 (MQ=255)
                aaCTGAAATTTTTAATTCATCTTCCCGACCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCAtt                                                                         <  2:1146272/100‑1 (MQ=255)
                 aCTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGAATTAAACCTGCAGTATGCCttt                                                                        >  1:1505972/1‑100 (MQ=255)
                  cTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACATTCGCGGTGACTTATCCCTGCAGTATGCAtttt                                                                       >  2:510537/1‑100 (MQ=255)
                  cTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCAtttt                                                                       >  2:2273905/1‑100 (MQ=255)
                  cTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCAtttt                                                                       >  1:1267471/1‑100 (MQ=255)
                   ttaaaTTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTa                                                                      <  1:1150274/98‑1 (MQ=255)
                    gAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTAc                                                                     <  2:2799806/100‑1 (MQ=255)
                    gAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTAc                                                                     >  2:3078471/1‑100 (MQ=255)
                      aaTTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGc                                                                   <  1:600020/100‑1 (MQ=255)
                      aaTTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGc                                                                   <  2:2526478/100‑1 (MQ=255)
                      aaTTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGc                                                                   <  1:2202071/100‑1 (MQ=255)
                      aaTTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGc                                                                   >  2:2059692/1‑100 (MQ=255)
                        tttttAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCAc                                                                 >  2:1174151/1‑100 (MQ=255)
                         ttttAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCAcc                                                                >  2:2029755/1‑100 (MQ=255)
                          tttAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCa                                                               <  1:235488/100‑1 (MQ=255)
                             aaTTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAg                                                            >  1:2348316/1‑100 (MQ=255)
                             aaTTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAg                                                            >  1:2786888/1‑100 (MQ=255)
                              aTTCAGCTGCCCGGACAGTAAATGGCAAATGAAGCAACAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCa                                                                  <  1:2879811/93‑1 (MQ=255)
                                tCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTg                                                         >  2:1637737/1‑100 (MQ=255)
                                tCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTg                                                         <  1:1422734/100‑1 (MQ=255)
                                 cATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGc                                                        >  2:3179485/1‑100 (MQ=255)
                                 cATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGc                                                        >  1:590665/1‑100 (MQ=255)
                                 cATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGc                                                        <  2:782694/100‑1 (MQ=255)
                                 cATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGc                                                        <  2:3456640/100‑1 (MQ=255)
                                  aTCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCt                                                       >  1:3356853/1‑100 (MQ=255)
                                  aTCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCt                                                       <  1:1388388/100‑1 (MQ=255)
                                  aTCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCt                                                       <  2:629603/100‑1 (MQ=255)
                                  aTCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCt                                                       >  2:944577/1‑100 (MQ=255)
                                    cTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTac                                                     <  2:791967/100‑1 (MQ=255)
                                    cTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTac                                                     >  1:1771385/1‑100 (MQ=255)
                                      tCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTacac                                                   >  1:1516255/1‑100 (MQ=255)
                                       cccGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACg                                                  <  2:1613338/100‑1 (MQ=255)
                                        ccGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACga                                                 >  2:593330/1‑99 (MQ=255)
                                           gccatttagtGGCAAATCAAGCAAAAATTGGCAAATTTAACAGTCGCTGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTg                                                <  1:3137667/89‑1 (MQ=255)
                                            ccAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCg                                              <  1:22032/99‑1 (MQ=255)
                                             cAGTAAATGGCAAATGAAGCAAAAAGTGGCAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATc                                           >  1:1149192/1‑100 (MQ=255)
                                               gTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCACTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCa                                          >  1:724105/1‑100 (MQ=255)
                                                tAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAg                                         >  1:1520391/1‑100 (MQ=255)
                                                  aaTGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCa                                       <  1:1237033/100‑1 (MQ=255)
                                                  aaTGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCa                                       <  2:570416/100‑1 (MQ=255)
                                                     ggCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATAc                                    >  2:133192/1‑100 (MQ=255)
                                                      gCAAATTAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATAc                                    <  1:2277686/99‑1 (MQ=255)
                                                      gCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACg                                   >  1:1949276/1‑100 (MQ=255)
                                                       cAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGc                                  <  2:802577/100‑1 (MQ=255)
                                                        aaaTGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCt                                 >  1:1755506/1‑100 (MQ=255)
                                                         aaTGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCtt                                >  1:2482784/1‑100 (MQ=255)
                                                          aTGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTa                               <  2:1068048/100‑1 (MQ=255)
                                                          aTGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTa                               <  2:1047851/100‑1 (MQ=255)
                                                            gAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAac                             <  2:1656447/100‑1 (MQ=255)
                                                            gAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAac                             >  1:1109535/1‑100 (MQ=255)
                                                            gAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAac                             <  2:2861390/100‑1 (MQ=255)
                                                             aaGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAaca                            <  1:972059/100‑1 (MQ=255)
                                                              aGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAacac                           >  1:2825315/1‑100 (MQ=255)
                                                              aGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAacac                           >  1:2041435/1‑100 (MQ=255)
                                                                cAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGt                         >  1:3256720/1‑100 (MQ=255)
                                                                cAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGt                         <  1:2252078/100‑1 (MQ=255)
                                                                 aaaaaGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGta                        <  1:133190/100‑1 (MQ=255)
                                                                      gTGGCAAATGTAACAGTGGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGTATATc                    <  2:1507133/99‑1 (MQ=255)
                                                                        ggCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGTATATCAgg                 >  1:1397565/1‑100 (MQ=255)
                                                                        ggCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGTATATCAgg                 >  2:360103/1‑100 (MQ=255)
                                                                        ggCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGTAGATCAgg                 >  1:2374799/1‑100 (MQ=255)
                                                                         gCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGTATATCAGGc                >  1:2886351/1‑100 (MQ=255)
                                                                            aaTGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGTATATCAGGCggg             <  2:1481228/100‑1 (MQ=255)
                                                                            aaTGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGTATATCAGGCggg             <  2:1771410/100‑1 (MQ=255)
                                                                            aaTGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGTATATCAGGCggg             <  2:3468318/100‑1 (MQ=255)
                                                                             aTGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGTATATCAGGCGGGt            <  2:2389594/100‑1 (MQ=255)
                                                                              tGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGTATATCAGGCGGGtt           <  1:570173/100‑1 (MQ=255)
                                                                               gTAACAGTCGGGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGTATATCAGGCGGGTTg          <  1:2275499/100‑1 (MQ=255)
                                                                               gTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGTATATCAGGCGGGTTg          <  2:675077/100‑1 (MQ=255)
                                                                                tAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGTATATCAGGCGGGTTGa         >  2:1070843/1‑100 (MQ=255)
                                                                                 aaCAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGTATATCAGGCGGGTTGAc        <  1:97662/100‑1 (MQ=255)
                                                                                 aaCAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGTATATCAGGCGGGTTGAc        <  1:263951/100‑1 (MQ=255)
                                                                                  aCAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGTATATCAGGCGGGTTGACt       <  2:2382714/100‑1 (MQ=255)
                                                                                   cAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGTATATCAGGCGGGTTGACTg      <  1:3465461/100‑1 (MQ=255)
                                                                                   cAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGTATATCAGGCGGGTTGACTg      <  2:2055156/100‑1 (MQ=255)
                                                                                   cAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGTATATCAGGCGGGTTGACTg      >  1:2032892/1‑100 (MQ=255)
                                                                                   cAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGTATATCAGGCGGGTTGACTg      >  1:382211/1‑100 (MQ=255)
                                                                                   cAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGTATATCAGGCGGGTTGACTg      <  2:3176561/100‑1 (MQ=255)
                                                                                     gTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGTATATCAGGCGGGTTGACTGcc    >  2:678076/1‑100 (MQ=255)
                                                                                     gTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGTATATCAGGCGGGTTGACTGcc    <  2:2960402/100‑1 (MQ=255)
                                                                                      tCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGTATATCAGGCGGGTTGACTGccc   <  1:1637715/100‑1 (MQ=255)
                                                                                       cgcgGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGTATATCAGGCGGGTTGACTGCCCa  >  1:2219945/1‑100 (MQ=255)
                                                                                               |                                                                                           
CCCGATGCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGATTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGTATATCAGGCGGGTTGACTGCCCA  >  NC_000913/3805380‑3805566

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