Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,805,475 T→C 35.0% 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,4750TC35.0% 145.5 / 159.6 163N159S (AAT→AGT) waaPkinase that phosphorylates core heptose of lipopolysaccharide
Reads supporting (aligned to +/- strand):  ref base T (45/61);  new base C (31/26);  total (76/87)
Fisher's exact test for biased strand distribution p-value = 1.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.66e-01

Only 100 of 166 total aligned reads displayed.

TGTGCCCGATGCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGATTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGTATATCAGGCGGGTTGACTGCCCAATCGGCA  >  NC_000913/3805376‑3805573
                                                                                                   |                                                                                                  
 gtgCCCGATGCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGAtt                                                                                                   >  1:471759/1‑100 (MQ=255)
 gtgCCCGATGCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGAtt                                                                                                   <  1:3140325/100‑1 (MQ=255)
   gCCCGATGCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGATTTa                                                                                                 >  1:510853/1‑100 (MQ=255)
   gCCCGATGCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTa                                                                                                 >  2:149768/1‑100 (MQ=255)
     ccGATGCAGATCGATAACTAAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGATTTATc                                                                                               <  1:961256/100‑1 (MQ=255)
      cGATGCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGATTTATcc                                                                                              <  1:2521635/100‑1 (MQ=255)
       gATGCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCGAGTAAAGGGCAAATGAAGCAAAAAGTCGCAAATGTAACAGTCGCGGTGATTTATccc                                                                                             >  1:778567/1‑100 (MQ=255)
       gATGCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGATTTATccc                                                                                             <  2:2591719/100‑1 (MQ=255)
        aTGCAGATCGATAACTGAAATTTTTAATTCGTCTTCCCGGCCAGTAAATGGCAAATGAAACAAAATGTGGCAAATGTAACAGTCGAGGTGATTTATCGCt                                                                                            >  1:2884280/1‑100 (MQ=255)
        aTGCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCt                                                                                            <  2:2713511/100‑1 (MQ=255)
          gCAGATCGATAACTGAAATTTTTAATTCATTTTCCCGGCCAGTAAATGGAAAAAAAAACAAAAAGTGGCAAATGTAACAGTCGCGGTGATTTATCCCTGc                                                                                          <  1:1951918/100‑1 (MQ=255)
           cAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGATTTATCCCTGCa                                                                                         <  1:1427181/100‑1 (MQ=255)
           cAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGATTTATCCCTGCa                                                                                         <  2:759036/100‑1 (MQ=255)
           cAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGATTTATCCCTGCa                                                                                         >  1:2668978/1‑100 (MQ=255)
              aTCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGATTTATCCCTGCAGTa                                                                                      <  2:2667195/100‑1 (MQ=255)
                cGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGATTTATCCCTGCAGTATg                                                                                    <  1:723200/100‑1 (MQ=255)
                cGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATg                                                                                    >  2:1103077/1‑100 (MQ=255)
                cGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATg                                                                                    >  2:3284337/1‑100 (MQ=255)
                 gATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGc                                                                                   >  1:1259558/1‑100 (MQ=255)
                   tAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGGGGCAAATGTAAAAGTCGCGGTGATTTATCCCTGCAGTATGCa                                                                                  >  2:1266504/1‑99 (MQ=255)
                     aCTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATAGCAAATGAAACAAAAAGTGGCACGTGTAACAGGCGCGGGGATTTATacc                                                                                             >  1:2328133/1‑86 (MQ=255)
                      cTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCAtttt                                                                              >  2:452836/1‑100 (MQ=255)
                      cTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCAtttt                                                                              >  1:713125/1‑100 (MQ=255)
                       tGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTa                                                                             <  1:452834/100‑1 (MQ=255)
                        gAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTAc                                                                            <  1:1201598/100‑1 (MQ=255)
                         aaaTTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGATTTATCCCTGCAGTATGCATTTTACg                                                                           <  2:2342543/100‑1 (MQ=255)
                          aaTTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATTAAGCAAAAATTGGCAAATGTAACAGTCGAGGTGATGTATCCCTGCAGTATGTATTTTAc                                                                            >  1:2496940/1‑98 (MQ=255)
                          aaTTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTa                                                                             <  2:1300663/97‑1 (MQ=255)
                            tttttAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGATTTATCCCTGCAGTATGCATTTTACGCAc                                                                        <  1:2071842/100‑1 (MQ=255)
                            tttttAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGATTTATCCCTGCAGTATGCATTTTACGCAc                                                                        <  2:2402905/100‑1 (MQ=255)
                            tttttAATTCATCTTCCCAGCCAGTAAATGGTAAATGAATTTAAAAGTGGCAAACGTAGCAGTCGCGGTGATTTATCCCTGCAGTATGCGTTTTACGc                                                                          <  1:1837987/98‑1 (MQ=255)
                             ttttAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCAcc                                                                       >  1:1405713/1‑100 (MQ=255)
                             ttttAATTCATCTGCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCAcc                                                                       >  2:2301974/1‑100 (MQ=255)
                                  aTTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGATTTATCCCTGCAGTATGCATTTTACGCACCATAGt                                                                  >  1:1597571/1‑100 (MQ=255)
                                  aTTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGATTTATCCCTGCAGTATGCATTTTACGCACCATAGt                                                                  >  2:333172/1‑100 (MQ=255)
                                   ttCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGATTTATCCCTGCAGTATGCATTTTACGCACCATAGtt                                                                 >  2:103551/1‑100 (MQ=255)
                                    tCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGATTTATCCCTGCAGTATGCATTTTACGCACCATAGTTg                                                                <  2:576709/100‑1 (MQ=255)
                                     cATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGc                                                               <  2:3371505/100‑1 (MQ=255)
                                     cATCTTCCCGGACAGTAAATGACAAATGAAGCAAAAAGTGGCAAATGTAAAAGTCGCGATGACTTATCACTGAAGTATGAATTTTACGCAAGATGGTTGc                                                               >  1:3075552/1‑100 (MQ=255)
                                      aTCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGATTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCt                                                              <  2:212343/100‑1 (MQ=255)
                                      aTCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGATTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCt                                                              <  1:1558962/100‑1 (MQ=255)
                                      aTCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGATTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCt                                                              <  1:2426689/100‑1 (MQ=255)
                                       tCTTCCCGGACAGTACATGGGAAATTAATCTAATAGTGGCAAATGTAACAGTCGCGGTGATTTATCCCTGCAGTTTACATTTAACGCACCATAGtt                                                                 <  1:748967/96‑1 (MQ=255)
                                        cTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTac                                                            <  2:678525/100‑1 (MQ=255)
                                         ttCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGATTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTaca                                                           <  2:858942/100‑1 (MQ=255)
                                         ttCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTaca                                                           <  2:539642/100‑1 (MQ=255)
                                           cccGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGATTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACg                                                         >  2:2157645/1‑100 (MQ=255)
                                            ccGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGt                                                        >  2:1157873/1‑100 (MQ=255)
                                            ccGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGt                                                        >  2:615034/1‑100 (MQ=255)
                                              ggCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGATTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGc                                                      >  1:237368/1‑100 (MQ=255)
                                               gCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGATTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCg                                                     <  2:2464844/100‑1 (MQ=255)
                                               gCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGATTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCg                                                     <  2:1498470/100‑1 (MQ=255)
                                                 cAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGATTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGAt                                                   >  2:1968428/1‑100 (MQ=255)
                                                  aGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGATTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATc                                                  >  1:520781/1‑100 (MQ=255)
                                                   gTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGATTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCa                                                 <  2:3358890/100‑1 (MQ=255)
                                                   gTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCa                                                 <  2:480171/100‑1 (MQ=255)
                                                    tAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGATTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCGACACGTGCGATCAg                                                <  1:550083/100‑1 (MQ=255)
                                                     aagtaggaaaGGAAGCAGAGGGTGGCAAATGTAACAGTCGCGGTGGTTGATCCGTGCAATCTGGAGTTGACGGATCATGGTTGCTACACGTGCGATCAgg                                               <  1:133780/93‑2 (MQ=255)
                                                      aaTGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGATTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGcn                                              <  2:2386215/100‑2 (MQ=255)
                                                      aaTGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGATTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCa                                              <  2:14598/100‑1 (MQ=255)
                                                       aTGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCAGTGATTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTAAACGTGCGATCAGCAt                                             <  1:2077102/100‑1 (MQ=255)
                                                        tGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGATTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATa                                            >  2:2866914/1‑100 (MQ=255)
                                                         ggCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGATTTATCCCTGCAGAATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATAc                                           >  2:2211071/1‑100 (MQ=255)
                                                          gCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGATTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACg                                          <  2:850448/100‑1 (MQ=255)
                                                           cAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGATTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGc                                         <  2:3113386/100‑1 (MQ=255)
                                                            aaaTGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGATTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCt                                        <  2:1873141/100‑1 (MQ=255)
                                                              aTGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGATTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTa                                      <  1:677004/100‑1 (MQ=255)
                                                              aTGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGATTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTa                                      <  1:1063740/100‑1 (MQ=255)
                                                               tGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGATTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTaa                                     >  1:1599483/1‑100 (MQ=255)
                                                                   gCAAAAAGTGGCAAATGTAACAGTCGCGGTGATTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACg                                 >  1:3320844/1‑100 (MQ=255)
                                                                   gCAAAAAGTGGCAAATGTAACAGTCGCGGTGATTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACg                                 >  1:1902881/1‑100 (MQ=255)
                                                                     aaaaaGTGGCGAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGta                               <  1:2954132/100‑1 (MQ=255)
                                                                     aaaaaGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGta                               >  1:1105128/1‑100 (MQ=255)
                                                                     aaaaaGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGta                               <  2:2182437/100‑1 (MQ=255)
                                                                     aaaaaGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGta                               >  2:1578454/1‑100 (MQ=255)
                                                                       aaaGTGGCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGtata                             <  2:3166191/100‑1 (MQ=255)
                                                                          gTGGCAAATGTAACAGTCGCGGTGATTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGTATATCa                          <  2:2657514/100‑1 (MQ=255)
                                                                            ggCAAATGTAACAGTCGCGGTGATTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGTATATCAgg                        >  2:2832988/1‑100 (MQ=255)
                                                                            ggCAAATGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGTATATCAgg                        <  2:903889/100‑1 (MQ=255)
                                                                            gacAAATGTAACAATCGCGGTGACTTATCCCAGCAGTATGCAGTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGTATATCAgg                        >  2:2285365/3‑100 (MQ=255)
                                                                             gCAAATGTAACAGTCGCGGTGATTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGTATATCAGGc                       <  1:2994935/100‑1 (MQ=255)
                                                                              cAAATGTAACAGTCGCGGTGATTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGTATATCAGGCg                      >  1:2625788/1‑100 (MQ=255)
                                                                              cAAATGTAACAGTCGCGGTGATTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGTATATCAGGCg                      <  2:1599490/100‑1 (MQ=255)
                                                                              cAAATGTAACAGTCGCGGTGATTTATCCCTGCAGTATGCATTTTACGCACCAGAGTTGCTACACGTGCGATCATCATGCGCTTAACACGTATATCAgg                        >  1:325751/1‑98 (MQ=255)
                                                                               aaaTGTAACAGTCGCGGTGATTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGTATATCAGGCgg                     <  1:2140851/100‑1 (MQ=255)
                                                                               aaaTGTAACAGTCGCGGTGATTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGTATATCAGGCgg                     >  2:2140866/1‑100 (MQ=255)
                                                                                aaTGTAACAGTCGCGGTGATTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGTATATCAGGCggg                    <  2:2625808/100‑1 (MQ=255)
                                                                                aaTGTAACAGTCGCGGTGATTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGTATATCAGGCggg                    <  2:2908582/100‑1 (MQ=255)
                                                                                aaTGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGTATATCAGGCggg                    <  1:2699141/100‑1 (MQ=255)
                                                                                aaTGTAACAGTCGCAGTGATTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGTATATCAGGCggg                    <  1:2910697/100‑1 (MQ=255)
                                                                                  tGTAACAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGTATATCAGGCGGGtt                  <  2:1550269/100‑1 (MQ=255)
                                                                                      aCAGTCGCGGTGATTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGTATATCAGGCGGGTTGACt              <  1:817341/100‑1 (MQ=255)
                                                                                       cAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGTATATCAGGCGGGTTGACTg             >  2:1213986/1‑100 (MQ=255)
                                                                                       cAGTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGTATATCAGGCGGGTTGACTg             >  1:3203549/1‑100 (MQ=255)
                                                                                         gTCGCGGTGATTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGTATATCAGGCGGGTTGACTGcc           >  1:2780702/1‑100 (MQ=255)
                                                                                         gTCGCGGTGACTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGTATATCAGGCGGGTTGACTGcc           >  2:1338025/1‑100 (MQ=255)
                                                                                            gcgGTGATTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGTATATCAGGCGGGTTGACTGCCCaa        >  2:3388452/1‑100 (MQ=255)
                                                                                            gcgGTGATTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGTATATCAGGCGGGTTGACTGCCCaa        >  2:3237985/1‑100 (MQ=255)
                                                                                               gTGATTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGTATATCAGGCGGGTTGACTGCCCAATCg     >  2:2537850/1‑100 (MQ=255)
                                                                                                  aTTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGTATATCAGGCGGGTTGACTGCCCAATCGGca  <  1:1968415/100‑1 (MQ=255)
                                                                                                   |                                                                                                  
TGTGCCCGATGCAGATCGATAACTGAAATTTTTAATTCATCTTCCCGGCCAGTAAATGGCAAATGAAGCAAAAAGTGGCAAATGTAACAGTCGCGGTGATTTATCCCTGCAGTATGCATTTTACGCACCATAGTTGCTACACGTGCGATCAGCATACGCTTAACACGTATATCAGGCGGGTTGACTGCCCAATCGGCA  >  NC_000913/3805376‑3805573

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