Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,805,475 T→C 100% 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.1% 427.8 / ‑4.7 127N159S (AAT→AGT) waaPkinase that phosphorylates core heptose of lipopolysaccharide
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (56/69);  total (57/70)
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.89e-01

Only 100 of 128 total aligned reads displayed.

TGTGCCCGATGCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGATTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGTATATCAGGCGGGTTGACTGCCCA  >  NC_000913/3805376‑3805566
                                                                                                   |                                                                                           
tgtgCCCGATGCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGAt                                                                                             <  2:2676748/100‑1 (MQ=255)
   gCCCGATGCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTa                                                                                          >  1:383652/1‑100 (MQ=255)
     ccGATGCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATc                                                                                        >  1:1474782/1‑100 (MQ=255)
       gATGCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTAt                                                                                         >  2:2033409/1‑97 (MQ=255)
       gATGCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATccc                                                                                      <  1:2278483/100‑1 (MQ=255)
         tGCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTg                                                                                    >  1:244545/1‑100 (MQ=255)
         tGCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTg                                                                                    <  1:1260677/100‑1 (MQ=255)
          gCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGc                                                                                   <  1:2708371/100‑1 (MQ=255)
           cAGATCGATAACTGAAATTTTTAATTAATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAACGTGGCAAATGTCACCGTCGCGGTGACTTAtacgt                                                                                     >  2:2876767/1‑93 (MQ=255)
            aGATCGATAACTGAAATTTTTAATTCATCTTCCCTGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAg                                                                                 <  2:399775/100‑1 (MQ=255)
            aGACCGATAACTGAAATTTTGAATTAATCTGCCCGGTCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAg                                                                                 <  2:2887423/100‑1 (MQ=255)
             gATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAg                                                                                 <  2:2674487/99‑1 (MQ=255)
             gATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGt                                                                                >  1:1692493/1‑100 (MQ=255)
             gATCGATAACTGAAAGTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGt                                                                                <  2:2655871/100‑1 (MQ=255)
              aTCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTa                                                                               <  2:199836/100‑1 (MQ=255)
               tCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTAt                                                                              <  1:1483954/100‑1 (MQ=255)
                cGATAACTGGAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATg                                                                             >  1:2020411/1‑100 (MQ=255)
                cGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATg                                                                             >  2:658708/1‑100 (MQ=255)
                  aTAACTGAAATTTTTAATTCATCTTCCCGTCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCa                                                                           <  1:162280/100‑1 (MQ=255)
                   tAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCa                                                                           >  1:1011656/1‑99 (MQ=255)
                   tAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCAt                                                                          <  1:1963034/100‑1 (MQ=255)
                    aaCTTAAATTTTTAATTCATCTTCCCTGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCAtt                                                                         <  1:1475806/100‑1 (MQ=255)
                    aaCTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCAtt                                                                         <  2:1267046/100‑1 (MQ=255)
                      cTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATAAAACAAAAAATAACAAATGTAACAATCGCAGTAACTTATCCCTGCAGTATGCAtttt                                                                       >  1:4262/1‑100 (MQ=255)
                        gAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTAc                                                                     >  1:5176/1‑100 (MQ=255)
                        gAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACg                                                                    <  1:1516010/100‑1 (MQ=255)
                         aaaTTTTTAATTCATCTTCCCGGCTAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACCGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACg                                                                    <  1:2145135/100‑1 (MQ=255)
                         aaaTTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACg                                                                    <  1:2747184/100‑1 (MQ=255)
                           aTTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGc                                                                   <  2:2683808/99‑1 (MQ=255)
                           aTTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCa                                                                  >  2:1249645/1‑100 (MQ=255)
                           aTTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCa                                                                  <  1:2092955/100‑1 (MQ=255)
                            tttttAATTCATCTTCCCGGCCAGTAAATGTCAAATGAAGTAAATAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCAc                                                                 <  1:18682/100‑1 (MQ=255)
                                   ttCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGtt                                                          >  2:248319/1‑100 (MQ=255)
                                    tCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTg                                                         >  2:2306075/1‑100 (MQ=255)
                                      aTCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCt                                                       >  2:589144/1‑100 (MQ=255)
                                       ncTTCCCGGCCAGGAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTa                                                      >  2:1877781/2‑100 (MQ=255)
                                        cTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATATTACGCACCATAGTTGCTa                                                      >  1:597405/1‑99 (MQ=255)
                                         ttCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTaca                                                    >  1:2156266/1‑100 (MQ=255)
                                           cccGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACg                                                  >  1:1939361/1‑100 (MQ=255)
                                           cccGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACg                                                  <  2:1692473/100‑1 (MQ=255)
                                           cccGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACg                                                  >  2:2086862/1‑100 (MQ=255)
                                             cgtccAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTg                                                <  2:1800570/97‑1 (MQ=255)
                                             cGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGt                                                 <  2:1729272/99‑1 (MQ=255)
                                              ggCCAGTTAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTTGTGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGc                                               <  1:134870/100‑1 (MQ=255)
                                              ggCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGc                                               >  1:1240396/1‑100 (MQ=255)
                                              ggCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGc                                               <  2:1040226/100‑1 (MQ=255)
                                              ggCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGc                                               <  2:1287485/100‑1 (MQ=255)
                                               gCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCg                                              >  1:1608863/1‑100 (MQ=255)
                                               gCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCg                                              <  2:333392/100‑1 (MQ=255)
                                               gCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCg                                              >  1:1685938/1‑100 (MQ=255)
                                                 taGTAAATGGCAAATGAAGTAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGa                                             <  1:1255624/98‑1 (MQ=255)
                                                  aGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATc                                           <  1:541115/100‑1 (MQ=255)
                                                  aGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATc                                           >  2:33671/1‑100 (MQ=255)
                                                  aGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATc                                           <  1:236703/100‑1 (MQ=255)
                                                   gTAAATGCCAAATGAAGCAAAAAGTGGCAAATGTAACAGTTGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACAAGTGCg                                              <  2:1032945/96‑1 (MQ=255)
                                                    tAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTTGTGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAg                                         <  1:1024502/100‑1 (MQ=255)
                                                     aaaTGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCTTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGc                                        <  1:857982/100‑1 (MQ=255)
                                                     aaaTGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGc                                        >  2:1601167/1‑100 (MQ=255)
                                                       aTGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCAt                                      <  2:437651/100‑1 (MQ=255)
                                                        tGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATa                                     >  2:623606/1‑100 (MQ=255)
                                                         ggCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATAc                                    >  2:3019277/1‑100 (MQ=255)
                                                         ggCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATAc                                    <  1:2900946/100‑1 (MQ=255)
                                                         ggCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCACTACAGTATGCATAATACGCACCATAATTGCTACACGTGCGATCAGCATAc                                    >  1:303030/1‑100 (MQ=255)
                                                          gCAAATGAAGCTGAAAGTGGCAAATGTAACCGTCGCTGTGACTTTTCTCTACAGTATGCATTTTACGCTCCATTTTTGCTTCATGGGCGATCAGc                                        <  2:2432777/95‑1 (MQ=255)
                                                           cAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGc                                  >  1:1045727/1‑100 (MQ=255)
                                                            aaaTGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCt                                 <  1:2086885/100‑1 (MQ=255)
                                                            aaaTGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCt                                 >  2:2191553/1‑100 (MQ=255)
                                                             aaTGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCtt                                <  2:1073375/100‑1 (MQ=255)
                                                             aaTGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCtt                                <  2:1490071/100‑1 (MQ=255)
                                                                gAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAac                             <  2:2206341/100‑1 (MQ=255)
                                                                  aGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAacac                           <  1:89943/100‑1 (MQ=255)
                                                                   gCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACg                          >  1:132153/1‑100 (MQ=255)
                                                                    catagatgggCGAATGTCACAGTCGCGGTGGCGATTGCCTGCCGTATGCAGTTTAACAACCAGAGTTGCTTAACGTGAGATCAGCATACGCTTa                               <  2:2914382/86‑1 (MQ=255)
                                                                     aaaaaGTGGCAAATGTAACAGTCGCGGTGATTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACAAGTGCGATCAGCATACGCTTAACACGta                        >  1:1250779/1‑100 (MQ=255)
                                                                     aaaaaGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGta                        <  2:329996/100‑1 (MQ=255)
                                                                     aaaaaGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGta                        <  1:68685/100‑1 (MQ=255)
                                                                     aaaaaGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTAAACGTGCGATCAGCATACGCTTAACACGta                        >  1:1341544/1‑100 (MQ=255)
                                                                      aaaaGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGtat                       <  2:1891373/100‑1 (MQ=255)
                                                                         aGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGTATATc                    <  2:1425490/100‑1 (MQ=255)
                                                                           tGGCAAATGTAACAGTCGCGGTGACTTATTCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGTTCAGCATACGCTTAACACGTATATCan                  <  2:1658271/100‑2 (MQ=255)
                                                                           tGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGTATATCAg                  <  2:2114888/100‑1 (MQ=255)
                                                                            ggCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCGACACGTGCGATCAGCATACGCTTAACACGTATATCAgg                 >  2:1647939/1‑100 (MQ=255)
                                                                               aaaTGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGTATATCAGGCgg              >  1:1050512/1‑100 (MQ=255)
                                                                                aaTGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGTATATCAGGCggg             <  2:244542/100‑1 (MQ=255)
                                                                                aaTGTAACAGTCGCGGTGACTTAACAATGCGGTGTGGATTTTACGCAACATATTTGCTACACGTGCGATGAGCATACGCGTAACACGTATATccg                  >  2:3044801/1‑93 (MQ=255)
                                                                                 aTGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGTATATCAGGTGGGt            >  1:1961919/1‑100 (MQ=255)
                                                                                 aTGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGTATATCAGGCggg             >  1:1417714/1‑99 (MQ=255)
                                                                                 aTGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGTATATCAGGCGGGt            >  1:794690/1‑100 (MQ=255)
                                                                                  tGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGTATATCAGGCGGGtt           >  1:274372/1‑100 (MQ=255)
                                                                                   gTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGTATATCAGGCGGGTTg          <  1:2382628/100‑1 (MQ=255)
                                                                                    tAACAGTCGCGGTGACTTATTTTTTTATTATTCTTTTTTTGTACCATAGTTGCTACACGTTTGATAAGCATACGCTTAACACGTATATCAGGCGGGTTg          >  2:459013/1‑99 (MQ=255)
                                                                                    tAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGTATATCAGGCGGGTTGa         >  1:1754991/1‑100 (MQ=255)
                                                                                     aaCAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGTATATCAGGCGGGTTGAc        <  1:2025567/100‑1 (MQ=255)
                                                                                       cAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGTATATCAGGCGGGTTGACTg      <  2:1133492/100‑1 (MQ=255)
                                                                                        aGTCGCAGTAATTTACCCATGAAGTATACTTTTTCCACAACAAAGTTCCTACACCCGCGGTCAGAATACGCTTAACACGTATATCAGGcaagttcacc       <  2:556129/98‑10 (MQ=255)
                                                                                         ntCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGTATATCAGGCGGGTTGACTGcc    >  2:2906765/2‑100 (MQ=255)
                                                                                          tCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCCTACGCTTAACACGTATATCAGGCGGGTTGACTGccc   <  1:2147287/100‑1 (MQ=255)
                                                                                          tCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGTATATCAGGCGGGTTGACTGccc   <  2:2786018/100‑1 (MQ=255)
                                                                                          tCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGTATATCAGGCGGGTTGACTGccc   <  2:2512228/100‑1 (MQ=255)
                                                                                           cgcgGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGTATATCAGGCGGGTTGACTGccc   <  2:1527528/99‑1 (MQ=255)
                                                                                                   |                                                                                           
TGTGCCCGATGCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGATTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGTATATCAGGCGGGTTGACTGCCCA  >  NC_000913/3805376‑3805566

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