Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,806,158 T→C 35.6% intergenic (‑75/‑212) ydiY ← / → pfkB acid‑inducible putative outer membrane protein/6‑phosphofructokinase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,806,1580TC35.6% 135.9 / 151.0 155intergenic (‑75/‑212)ydiY/pfkBacid‑inducible putative outer membrane protein/6‑phosphofructokinase II
Reads supporting (aligned to +/- strand):  ref base T (57/42);  new base C (24/31);  total (81/74)
Fisher's exact test for biased strand distribution p-value = 1.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.01e-01

Only 100 of 156 total aligned reads displayed.

CAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGATTAGCGTCCCTGGAAAGGTAACGAATTATAAAAAGGCGCGAATAACTTAGCAATGTATTCTTATTTCATTTTTTGAATAAGCATGTGGCGAAAACAGATTTTTAT  >  NC_000913/1806059‑1806256
                                                                                                   |                                                                                                  
cAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGATTAGCGt                                                                                                    >  1:3259740/1‑100 (MQ=255)
   gAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGATTAGCGTccc                                                                                                 >  2:714378/1‑100 (MQ=255)
       tGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGATTAGCGTCCCTGGa                                                                                             >  2:2924696/1‑100 (MQ=255)
        gTCTCCAAACGCTTCATTTATCAAGAGTCCGTACAACAACAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGATTAGCGCCCCTgg                                                                                              <  1:4196429/98‑1 (MQ=255)
          cTTCAAAAGCTTCATTTATCAAGAGTCCGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGATTAGCGTCCCTGGAAAg                                                                                          >  2:4180593/1‑100 (MQ=255)
          cTTCAAAAGCTTCATTTATCAAGAGTCCGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGATTAGCGCCCCTGGAAAg                                                                                          <  2:3167873/100‑1 (MQ=255)
           ttCAAAAGCTTCATTTATCAAGAGTCCGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGATTAGCGTCCCTGGAAAgg                                                                                         >  1:180665/1‑100 (MQ=255)
           ttCAAAAGCTTCATTTATCAAGAGTCCGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGATTAGCGCCCCTGGAAAgg                                                                                         >  2:2777192/1‑100 (MQ=255)
           ttCAAAAGCTTCATTTATCAAGAGTCCGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGATTAGCGCCCCTGGAAAgg                                                                                         <  2:1107211/100‑1 (MQ=255)
              aaaaGCTTCATTTATCAAGAGTCCGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGATTAGCGCCCCTGGAAAGGTaa                                                                                      <  2:94796/100‑1 (MQ=255)
              aaaaGCTTCATTTATCAAGAGTCCGTACAAAAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGATTAGCGTCCCTGGAAAg                                                                                          <  1:2306243/96‑1 (MQ=255)
                aaGCTTCATTTATCAAGAGTCCGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGATTAGCGTCCCTGGAAAGGTAACg                                                                                    <  2:3069529/100‑1 (MQ=255)
                   cTTCATTTATCAAGAGTCCGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGATTAGCGTCCCTGGAAAGGTAACGAAt                                                                                 <  2:2224332/100‑1 (MQ=255)
                   cTTCATTTATCAAGAGTCCGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGATTAGCGTCCCTGGAAAGGTAACGAAt                                                                                 <  2:146757/100‑1 (MQ=255)
                    ttCATTTATCAAGAGTCCGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGATTAGCGTCCCTGGAAAGGTAACGAAtt                                                                                <  1:1452391/100‑1 (MQ=255)
                     tCATTTATCAAGAGTCCGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGATTAGCGTCCCTGGAAAGGTAACGAATta                                                                               >  1:2381781/1‑100 (MQ=255)
                      cATTTATCAAGAGTCCGTACAACAAACAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGAGCAAAGATTAGCGTCCCTGGAAAGGTAACGAATtat                                                                              >  1:1301741/1‑100 (MQ=255)
                       aTTTATCAAGAGTCCGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGATTAGCGTCCCTGGAAAGGTAACGAATtata                                                                             >  2:2309071/1‑100 (MQ=255)
                       aTTTATCAAGAGTACGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGATTAGCGTCCCTGGAAAGGTAACGAATtata                                                                             <  1:2989308/100‑1 (MQ=255)
                        tttATCAAGAGGCCGTACAACAAAAAAAGAGACCATCGCTGTCCCGGAAACTTTCTTAAGGATCAAAGATTAGCGTCCCTGGAAAGGTAACGAATTATaa                                                                            <  1:1000749/100‑1 (MQ=255)
                        tttATCAAGAGGCCGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGATTAGCGTCCCTGGAAAGGTAACGAATTATaa                                                                            <  1:1001152/100‑1 (MQ=255)
                          tATCAAGAGTCCGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGATTAGCGCCCCTGGAAAGGTAACGAATTATaaaa                                                                          >  1:1551993/1‑100 (MQ=255)
                          tATCAAGAGTCCGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAACATTAGCGTCCCTGGAAAGGTAACGAATTATaaaa                                                                          >  1:1590220/1‑100 (MQ=255)
                          tATCAAGAGTCCGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTGAAGGATCAAAGATTAGCGTCCCTGGAAAGGTAACGAATTATaaaa                                                                          >  2:1290112/1‑100 (MQ=255)
                           aTCAAGAGTCCGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGATTAGCGTCCCTGGAAAGGTAACGAATTATaaaaa                                                                         >  1:1176585/1‑100 (MQ=255)
                           aTCAAGAGTCCGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGATTAGCGTCCCTGGAAAGGTAACGAATTATaaaaa                                                                         >  2:1493195/1‑100 (MQ=255)
                            tCAAGAGTCCGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGATTAGCGTCCCTGGAAAGGTAACGAATTATAAAAAg                                                                        >  2:1201719/1‑100 (MQ=255)
                            tCAAGAGTCCGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGATTAGCGTCCCTGGAAAGGTAACGAATTATAAAAAg                                                                        >  2:2845850/1‑100 (MQ=255)
                            tCAAGAGTCCGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGATTAGCGCCCCTGGAAAGGTAACGAATTATAAAAAg                                                                        >  2:2439565/1‑100 (MQ=255)
                             cAAGAGTCCGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGATTAGCGTCCCTGGAAAGGTAACGAATTATAAAAAgg                                                                       >  2:3774668/1‑100 (MQ=255)
                              aaGAGTCCGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGATTAGCGCCCCTGGAAAGGTAACGAATTATAAAAAGgc                                                                      <  2:3759632/100‑1 (MQ=255)
                                gagTCCGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGATTAGCGTCCCTGGAAAGGTAACGAATTATAAAAAGgcgc                                                                    <  1:1290108/100‑1 (MQ=255)
                                gagTCCGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGATTAGCGTCCCTGGAAAGGTAACGAATTATAAAAAGgcgc                                                                    <  1:3403623/100‑1 (MQ=255)
                                 agaccGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGATTAGCGCCCCTGGAAAGGTAACGAATTATAAAAAGGcgcg                                                                   <  1:1056853/97‑1 (MQ=255)
                                 agTCCGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGATTAGCGTCCCTGGAAAGGTAACGAATTATAAAAAGGcgcg                                                                   <  2:1854210/100‑1 (MQ=255)
                                 agTCCGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGATTAGCGCCCCTGGAAAGGTAACGAATTATAAAAAGGcgcg                                                                   <  2:258576/100‑1 (MQ=255)
                                     cGTATAACAAAAAAAGAGACCATCGCCGTCCCGGAAACTTTCTTAAGAATCAAAGATTAGCGCCCCTGGAAAGGTAACGAATTATAAAAAGGCGCGAATa                                                               <  1:3662800/100‑1 (MQ=255)
                                     cGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGATTAGCGTCCCTGGAAAGGTAACGAATTATAAAAAGGCGCGAATa                                                               >  2:744381/1‑100 (MQ=255)
                                     cGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGATTAGCGCCCCTGGAAAGGTAACGAATTATAAAAAGGCGCGAATa                                                               <  2:239069/100‑1 (MQ=255)
                                     cGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGATTAGCGCCCCTGGAAAGGTAACGAATTATAAAAAGGCGCGAATa                                                               <  2:238393/100‑1 (MQ=255)
                                     cGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGATTAGCGCCCCTGGAAAGGTAACGAATTATAAAAAGGCGCGAATa                                                               <  2:3628702/100‑1 (MQ=255)
                                      gTACAACAAAAAAAGAGACCATCGCGGTCCCGGCAACTTTCTTAAGGATCAAAGATTAGCGTCCCTGGAATGGTAACGAATTATAAAAAGGCGCGAATaa                                                              >  1:2571431/1‑100 (MQ=255)
                                       tACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGATTAGCGTCCCTGGAAAGGTAACGAATTATAAAAAGGCGCGAATAAc                                                             <  2:1919205/100‑1 (MQ=255)
                                       tACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGATTAGCGTCCCTGGAAAGGTAACGAATTATAAAAAGGCGCGAATAAc                                                             >  1:3119117/1‑100 (MQ=255)
                                        cgatctaaaaaaGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGATTAGCGTCCCTGGAAAGGTAACGAATTATAAAAAGGCGCGAATAACt                                                            <  1:451662/94‑1 (MQ=255)
                                         caacaaAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGATTAGCGCCCCTGGAAATGTAACGAATTATAAAAAGGCGCGAATAACtt                                                           <  1:750835/100‑1 (MQ=255)
                                          aacaaAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGATTAGCGTCCCTGGAAAGGTAACGAATTATAAAAAGGCGCGAATAACTTa                                                          <  2:1450931/100‑1 (MQ=255)
                                          aacaaAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGATTAGCGCCCCTGGAAAGGTAACGAATTATAAAAAGGCGCGAATAACTTa                                                          <  1:1511525/100‑1 (MQ=255)
                                          aacaaAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGATTAGCGCCCCTGGAAAGGTAACGAATTATAAAAAGGCGCGAATAACTTa                                                          <  1:2707018/100‑1 (MQ=255)
                                           acaaAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGATTAGCGCCCCTGGAAAGGTAACGAATTATAAAAAGGCGCGAATAACTTAg                                                         <  1:3174779/100‑1 (MQ=255)
                                           acaaAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGATTAGCGCCCCTGGAAAGGTAACGAATTATAAAAAGGCGCGAATAACTTAg                                                         <  2:1691921/100‑1 (MQ=255)
                                            caaAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGATTAGCGTCCCTGAAAAGGTAACGAATTATAAAAAGGCGCGAATAACTTAGc                                                        <  1:356483/100‑1 (MQ=255)
                                            caaAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGATTAGCGCCCCTGGAAAGGTAACGAATTATAAAAAGGCGCGAATAACTTAGc                                                        >  2:368946/1‑100 (MQ=255)
                                            caaAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGATTAGCGCCCCTGGAAAGGTAACGAATTATAAAAAGGCGCGAATAACTTAGc                                                        <  2:759923/100‑1 (MQ=255)
                                            caaAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGATTAGCGCCCCTGGAAAGGTAACGAATTATAAAAAGGCGCGAATAACTTAGc                                                        >  1:1156890/1‑100 (MQ=255)
                                             aaaaaaaGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGATTAGCGTCCCTGGAAAGGTAACGAATTATAAAAAGGCGCGAATAACTTAGCa                                                       <  2:44594/100‑1 (MQ=255)
                                              aaaaaaGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGATTAGCGTCCCTGGAAAGGTAACGAATTATAAAAAGGCGCGAATAACtagatcg                                                      >  2:451664/1‑94 (MQ=255)
                                               caaaaccgaCCATCGCGGTCCCCGGAACTTTCTTAGGCACCCAAGATTACTGGCCCTTGAAAGTGAACGAACTCTACAAAGGCGCGAATAACTTAGCAAt                                                     <  1:647884/93‑1 (MQ=255)
                                                aaaaGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGATTAGCGTCCCTGGAAAGGTAACGAATTATAAAAAGGCGCGAATAACTTAGCAATg                                                    <  1:2144252/100‑1 (MQ=255)
                                                aaaaGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGATTAGCGTCCCTGGAAAGGTAACGAATTATAAAAAGGCGCGAATAACTTAGCAATg                                                    >  1:2760003/1‑100 (MQ=255)
                                                aaaaGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGATTAGCGCCCCTGGAAAGGTAACGAATTATAAAAAGGCGCGAATAACTTAGCAATg                                                    <  1:926541/100‑1 (MQ=255)
                                                aaaaGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGATTAGCGCCCCTGGAAAGGTAACGAATTATAAAAAGGCGCGAATAACTTAGCAATg                                                    <  2:555547/100‑1 (MQ=255)
                                                     agaCCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGATTAGCGCCCCTGGAAAGGTAACGAATTATAAAAAGGCGCGAATAACTTAGCAATGTATTc                                               <  2:2422437/100‑1 (MQ=255)
                                                      gaCCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGATTAGCGTCCCTGGAAAGGTAACGAATTATAAAAAGGCGCGAATAACTTAGCAATGTATTCt                                              <  2:3171549/100‑1 (MQ=255)
                                                        ccATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGATTAGCGCCCCTGGAAAGGTAACGAATTATAAAAAGGCGCGAATAACTTAGCAATGTATTCTTa                                            >  2:921853/1‑100 (MQ=255)
                                                             gcgGTCCCGGAAACTTTCTTAAGGATCAAAGATTAGCGTCCCTGGAAAGGTAACGAATTATAAAAAGGCGCGAATAACTTAGCAATGTATTCTTATTTCa                                       >  2:1756859/1‑100 (MQ=255)
                                                             gcgGTCCCGGAAACTTTCTTAAGGATCAAAGATTAGCGTCCCTGGAAAGGTAACGAATTATAAAAAGGCGCGAATAACTTAGCAATGTATTCTTATTTCa                                       >  2:2741866/1‑100 (MQ=255)
                                                             gcgGTCCCGGAAACTTTCTTAAGGATCAAAGATTAGCGTCCCTGGAAAGGTAACGAATTATAAAAAGGCGCGAATAACTTAGCAATGTATTCTTATTTCa                                       >  1:2107955/1‑100 (MQ=255)
                                                             gcgGTCCCGGAAACTTTCTTAAGGATCAAAGATTAGCGTCCCTGGAAAGGTAACGAATTATAAAAAGGCGCGAATAACTTAGCAATGTATTCTTATTTCa                                       >  2:347586/1‑100 (MQ=255)
                                                                 tCCCGGAAACTTTCTTAAGGATCAAAGATTAGCGCCCCTGGAAAGGTAACGAATTATAAAAAGGCGCGAATAACTTAGCAATGTATTCTTATTTCAtttt                                   >  1:807761/1‑100 (MQ=255)
                                                                  cccGGAAACTTTCTTAAGGATCAAAGATTAGCGTCCCTGGAAAGGTAACGAATTATAAAAAGGCGCGAATAACTTAGCAATGTATTCTTATTTCAttttt                                  >  1:3879681/1‑100 (MQ=255)
                                                                   ccGGAAACTTTCTTAAGGATCAAAGATTAGCGCCCCTGGAAAGGTAACGAATTATAAAAAGGCGCGAATAACTTAGCAATGTATTCTTATTTCAtttttt                                 >  1:3961023/1‑100 (MQ=255)
                                                                   ccGGAAACTTTCTTAAGGATCAAAGATTAGCGCCCCTGGAAAGGTAACGAATTATAAAAAGGCGCGAATAACTTAGCAATGTATTCTTATTTCAtttttt                                 <  2:1545139/100‑1 (MQ=255)
                                                                    cGGAAACTTTCTTAAGGATCAAAGATTAGCGTCCCTGGAAAGGTAACGAATTATAAAAAGGCGCGAATAACTTAGCAATGTATTCTTATTTCATTTTTTg                                <  2:801843/100‑1 (MQ=255)
                                                                    cGGAAACTTTCTTAAGGATCAAAGATTAGCGTCCCTGGAAAGGTAACGAATTATAAAAAGGCGCGAATAACTTAGCAATGTATTCTTATTTCATTTTTTg                                >  1:2856924/1‑100 (MQ=255)
                                                                    cGGAAACTTTCTTAAGGATCAAAGATTAGCGTCCCTGGAAAGGTAACGAATTATAAAAAGGCGCGAATAACTTAGCAATGTATTCTTATTTCATTTTTTg                                >  2:308880/1‑100 (MQ=255)
                                                                    cGGAAACTTTCTTAAGGATCAAAGATTAGCGCCCCTGGAAAGGTAACGAATTATAAAAAGGCGCGAATAACTTAGCAATGTATTCTTATTTCATTTTTTg                                <  1:3027696/100‑1 (MQ=255)
                                                                       aaaCTTTCTTAAGGATCAAAGATTAGCGTCCCTGGAAAGGTAACGAATTATAAAAAGGCGCGAATAACTTAGCAATGTATTCTTATTTCATTTTTTGAAt                             <  2:3400818/100‑1 (MQ=255)
                                                                        aaCTTTCTTAAGGATCAAAGATTAGCGTCCCTGGAAAGGTAACGAATTATAAAAAGGCGCGAATAACTTAGCAATGTATTCTTATTTCATTTTTTGAATa                            <  1:3729475/100‑1 (MQ=255)
                                                                         aCTTTCTTAAGGATCAAAGATTAGCGCCCCTGGAAAGGTAACGAATTATAAAAAGGCGCGAATAACTTAGCAATGTATTCTTATTTCATTTTTTGAATaa                           <  2:2260385/100‑1 (MQ=255)
                                                                          cTTTCTTAAGGATCAAAGATTAGCGTCCCTGGAAAGGTAACGAATTATAAAAAGGCGCGAATAACTTAGCAATGTATTCTTATTTCATTTTTTGAATAAg                          >  1:992704/1‑100 (MQ=255)
                                                                          cTTTCTTAAGGATCAAAGATTAGCGTCCCTGGAAAGGTAACGAATTATAAAAAGGCGCGAATAACTTAGCAATGTATTCTTATTTCATTTTTTGAATAAg                          >  2:1173247/1‑100 (MQ=255)
                                                                          cTTTCTTAAGGATCAAAGATTAGCGTCCCTGGAAAGGTAACGAATTATAAAAAGGCGCGAATAACTTAGCAATGTATTCTTATTTCATTTTTTGAATAAg                          <  2:755867/100‑1 (MQ=255)
                                                                          cTTTCTTAAGGATCAAAGATTAGCGTCCCTGGAAAGGTAACGAATTATAAAAAGGCGCGAATAACTTAGCAATGTATTCTTATTTCATTTTTTGAATAAg                          >  2:1704004/1‑100 (MQ=255)
                                                                          cTTTCTTAAGGATCAAAGATTAGCGTCCCTGGAAACGTCACGAATTATAAAAAGGCGCCAATAACTTAGCACTGTATTCCTATTTCATTTTTTGAATaga                          >  1:3792300/1‑98 (MQ=255)
                                                                           tttCTTAAGGATCAAAGATTAGCGTCCCTGGAAAGGTAACGAATTATAAGAAGGCGCGAATAACTTAGCAATGTATTCTTATTTCATTTTTTGAATAAGc                         <  1:3704563/100‑1 (MQ=255)
                                                                           tttCTTAAGGATCAAAGATTAGCGTCCCTGGAAAGGTAACGAATTATAAGAAGGCGCGAATAACTTAGCAATGTATTCTTATTTCATTTTTTGAATAAGc                         >  2:3704579/1‑100 (MQ=255)
                                                                            ttCTTAAGGATCAAAGATTAGCGTCCCTGGAAAGGTAACGAATTATAAAAAGGCGCGAATAACTTAGCAATGTATTCTTATTTCATTTTTTGAATAAGCa                        <  2:1199281/100‑1 (MQ=255)
                                                                            ttCTTAAGGATCAAAGATTAGCGTCCCTGGAAAGGTAACGAATTATAAAAAGGCGCGAATAACTTAGCAATGTATTCTTATTTCATTTTTTGAATAAGCa                        <  2:2321841/100‑1 (MQ=255)
                                                                                tAAGGATCAAAGATTAGCGTCCCTGGAAAGGTAACGAATTATAAAAAGGCGCGAATAACTTAGCAATGTATTCTTATTTCATTTTTTGAATAAGCAtgtg                    <  2:4166771/100‑1 (MQ=255)
                                                                                tAAGGATCAAAGATTAGCGTCCCTGGAAAGGTAACGAATTATAAAAAGGCGCGAATAACTTAGCAATGTATTCTTATTTCATTTTTTGAATAAGCAtgtg                    <  2:4144754/100‑1 (MQ=255)
                                                                                  aGGATCAAAGATTAGCTCCCCTGGAAAGGTAACTAATTATAAAAAGGCGCTAATAACTTAGCTATGTATTCTTATTTCATTTTTTGAATAAGCATGTGGc                  <  1:1842670/100‑1 (MQ=255)
                                                                                   ggATCAAAGATTAGCGCCCCTGGAAAGGTAACGAATTATAAAAAGGCGCGAATAACTTAGCAATGTATTCTTATTTCATTTTTTGAATAAGCATGTGGCg                 >  1:3341727/1‑100 (MQ=255)
                                                                                    gATCAAAGATTAGCGCCCCTGGAAAGGTAACGAATTATAAAAAGGCGCGAATAACTTAGCAATGTATTCTTATTTCATTTTTTGAATAAGCATGTGGCGa                >  1:2870872/1‑100 (MQ=255)
                                                                                       cAAAGATTAGCGCCCCTGGAAAGGTAACGAATTATAAAAAGGCGCGAATAACTTAGCAATGTATTCTTATTTCATTTTTTGAATAAGCATGTGGCGaaaa             >  1:1249487/1‑100 (MQ=255)
                                                                                       cAAAGATTAGCGCCCCTGGAAAGGTAACGAATTATAAAAAGGCGCGAATAACTTAGCAATGTATTCTTATTTCATTTTTTGAATAAGCATGTGGCGaaaa             >  2:1581901/1‑100 (MQ=255)
                                                                                        aaaGATTAGCGCCCCTGGAAAGGTAACGAATTATAAAAAGGCGCGAATAACTTAGCAATGTATTCTTATTTCATTTTTTGAATAAGCATGTGGCGAAAAc            >  2:1909867/1‑100 (MQ=255)
                                                                                        aaaGATTAGCGCCCCTGGAAAGGTAACGAATTATAAAAAGGCGCGAATAACTTAGCAATGTATTCTTATTTCATTTTTTGAATAAGCATGTGGCGAAAAc            <  1:1202330/100‑1 (MQ=255)
                                                                                           gATTAGCGTCCCTGGAAAGGTAACGAATTATAAAAAGGCGCGAATAACTTAGCAATGTATTCTTATTTCATTTTTTGAATAAGCATGTGGCGAAAACAGa         <  2:2856935/100‑1 (MQ=255)
                                                                                                  gTCCCTGGAAAGGTAACGAATTATAAAAAGGCGCGAATAACTTAGCAATGTATTCTTATTTCATTTTTTGAATAAGCATGTGGCGAAAACAGATTtttat  >  2:3656452/1‑100 (MQ=255)
                                                                                                   |                                                                                                  
CAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACTTTCTTAAGGATCAAAGATTAGCGTCCCTGGAAAGGTAACGAATTATAAAAAGGCGCGAATAACTTAGCAATGTATTCTTATTTCATTTTTTGAATAAGCATGTGGCGAAAACAGATTTTTAT  >  NC_000913/1806059‑1806256

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