Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,364,998 C→A 15.4% intergenic (+19/+20) nudI → / ← ais nucleoside triphosphatase/putative LPS core heptose(II)‑phosphate phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,364,9980CA15.4% 459.2 / 60.1 162intergenic (+19/+20)nudI/aisnucleoside triphosphatase/putative LPS core heptose(II)‑phosphate phosphatase
Reads supporting (aligned to +/- strand):  ref base C (58/79);  new base A (10/15);  total (68/94)
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.71e-01

Only 100 of 165 total aligned reads displayed.

TGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTT  >  NC_000913/2364899‑2365097
                                                                                                   |                                                                                                   
   gTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTgg                                                                                                  >  1:2190553/1‑100 (MQ=255)
    tAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGc                                                                                                 <  1:257110/100‑1 (MQ=255)
     aaaaCCAGAAGATCGTGTGCATTATGATTTGAATCTCGCCACACGAAAAACGTTACGTTTGAAAGGTTTTCTGTAACAATGATGACAGCCATTACTGGc                                                                                                 <  2:1467934/99‑1 (MQ=255)
        aCCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTc                                                                                             <  1:3389409/100‑1 (MQ=255)
        aCCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTc                                                                                             >  2:2855826/1‑100 (MQ=255)
         ccTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCa                                                                                            <  2:1261439/100‑1 (MQ=255)
         ccTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCa                                                                                            >  1:29205/1‑100 (MQ=255)
         ccTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCAGTAATGGCTGTCa                                                                                            <  1:1755311/100‑1 (MQ=255)
          cTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCAGTAATGGCTGTCat                                                                                           <  1:2711490/100‑2 (MQ=255)
          cTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCAGTAATGGCTGTCat                                                                                           <  1:3618793/100‑2 (MQ=255)
          cTGAAGATCTGGTGCATATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCAGTAATGGCTGTCatc                                                                                          <  2:1518739/100‑3 (MQ=255)
           ttaaGATCTGGTGCATTATGATTTGTATGTCGCCACCCGAAAAACGTTACCTTTAACAGGTCTTCTGTAACAATGATGACATCCATTACTAGCGGTCAGa                                                                                          >  2:852906/3‑100 (MQ=255)
           tGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAg                                                                                           <  1:2998576/99‑1 (MQ=255)
           tGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCAGTAATGGCTGTCatc                                                                                          <  2:1507666/100‑3 (MQ=255)
           tGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCAGTAATGGCTGTCatc                                                                                          <  1:561834/100‑3 (MQ=255)
           tGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCAGTAATGGCTGTCatc                                                                                          <  1:1746412/100‑3 (MQ=255)
            gAAGATCTGGTGCATTATTATTTGAATGACAAAAAAAAAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCAGTAATGGCTGTCatc                                                                                          <  2:2471183/99‑3 (MQ=255)
            gAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGAt                                                                                         <  2:2944817/100‑1 (MQ=255)
             aaGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGAtt                                                                                        <  2:701813/100‑1 (MQ=255)
             aaGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGAtt                                                                                        <  1:1124676/100‑1 (MQ=255)
              aGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCAGTAATGGCTGTCatcatt                                                                                       <  2:2604834/100‑6 (MQ=39)
               gATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGa                                                                                      >  2:3503899/1‑100 (MQ=255)
               gATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGa                                                                                      >  1:1444020/1‑100 (MQ=255)
                 tCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGaaa                                                                                    <  2:2986352/100‑1 (MQ=255)
                   tggtggCTTATGATTCGAATTTCGCGTCCCGAATAACGTTACGTTTGAAAGGTCTTCTGTAACAAGGATGACAGCCATTACTGGCTGGCAGATTGaaa                                                                                    <  2:3266310/98‑1 (MQ=255)
                    ngTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAAtt                                                                                 >  2:2019407/2‑100 (MQ=255)
                    ggTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTg                                                                                       >  2:1226107/1‑94 (MQ=255)
                       gCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAg                                                                              <  2:2190539/100‑1 (MQ=255)
                        cATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGt                                                                             >  1:1150708/1‑100 (MQ=255)
                        cATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGt                                                                             <  1:1400431/100‑1 (MQ=255)
                        cATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGt                                                                             <  2:2940050/100‑1 (MQ=255)
                         aTTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTg                                                                            >  1:1720992/1‑100 (MQ=255)
                         aTTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTg                                                                            <  1:1048379/100‑1 (MQ=255)
                          ttATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTgg                                                                           <  2:2791329/100‑1 (MQ=255)
                            aTGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGtt                                                                         >  2:2160293/1‑100 (MQ=255)
                             tGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTa                                                                        <  1:1254551/100‑1 (MQ=255)
                             tGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTa                                                                        >  2:1570127/1‑100 (MQ=255)
                               aTTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAAc                                                                      >  1:3261509/1‑100 (MQ=255)
                               aTTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAAc                                                                      >  2:774195/1‑100 (MQ=255)
                                 ttGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGa                                                                    <  2:1869249/100‑1 (MQ=255)
                                  tGAATGTCGCCACCCGAAAAACGTTACGTTAGAAAGGTCTTATGTAACAATGATGACAGCCATTCCTGTCTGTCAGATTGAAAATTAAGTGGTTAAcaca                                                                   >  2:3525018/1‑97 (MQ=255)
                                   gAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAAt                                                                  <  1:1587379/100‑1 (MQ=255)
                                    aaTGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAAtt                                                                 >  2:560017/1‑100 (MQ=255)
                                    aaTGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAAtt                                                                 >  2:1825442/1‑100 (MQ=255)
                                      tGTCGTCACCCGAAATACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAATATTTAGTGTTTAACGAATTcc                                                               <  2:3382256/100‑1 (MQ=255)
                                        tCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTcccc                                                             >  1:2299539/1‑100 (MQ=255)
                                          gCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCAt                                                           <  1:109824/100‑1 (MQ=255)
                                           ccACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCAt                                                           <  2:1779154/99‑1 (MQ=255)
                                            cACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAAGAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATcc                                                         >  2:1946717/1‑100 (MQ=255)
                                             aCCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCa                                                        <  2:3225643/100‑1 (MQ=255)
                                              cccGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAg                                                       <  1:2964189/100‑1 (MQ=255)
                                              cccGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAg                                                       >  2:2679315/1‑100 (MQ=255)
                                                        gttacgttTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTg                                             <  2:1916290/100‑1 (MQ=255)
                                                         ttacgttTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACtagat                                            >  1:669201/1‑96 (MQ=255)
                                                         ttacgttTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGc                                            >  1:1004334/1‑100 (MQ=255)
                                                         ttacgttTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGc                                            <  2:1302276/100‑1 (MQ=255)
                                                          tacgttTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGcc                                           <  1:774200/100‑1 (MQ=255)
                                                          tacgttTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGcc                                           <  2:220341/100‑1 (MQ=255)
                                                           acgttTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCt                                          <  2:376416/100‑1 (MQ=255)
                                                           acgttTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCt                                          <  1:347593/100‑1 (MQ=255)
                                                            cgttTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCtt                                         <  1:2227275/100‑1 (MQ=255)
                                                             gttTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCttt                                        >  1:347282/1‑100 (MQ=255)
                                                               ttGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCttttt                                      <  1:2768559/100‑1 (MQ=255)
                                                                tGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTc                                     <  2:2310641/100‑1 (MQ=255)
                                                                   aaGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCcan                                  <  1:2019418/100‑2 (MQ=255)
                                                                     ggTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCAcat                                <  1:1257114/100‑1 (MQ=255)
                                                                     ggTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCAcat                                >  1:458060/1‑100 (MQ=255)
                                                                      gTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCAcatg                               <  1:1003440/100‑1 (MQ=255)
                                                                       tcttctGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCAcatgc                              <  2:2720799/100‑1 (MQ=255)
                                                                       tcttctGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGCTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTCCAcatgca                             <  2:589317/100‑1 (MQ=255)
                                                                        cttctGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCAcatgca                             <  2:528282/100‑1 (MQ=255)
                                                                        cttctGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCAcatgca                             <  1:355074/100‑1 (MQ=255)
                                                                          tctGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCAcatgcatg                           <  1:2911280/100‑1 (MQ=255)
                                                                            tGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGAc                         <  1:808450/100‑1 (MQ=255)
                                                                            tGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGAc                         >  1:3433294/1‑100 (MQ=255)
                                                                            tGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGAc                         >  1:1471973/1‑100 (MQ=255)
                                                                            tGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGAc                         <  1:1082319/100‑1 (MQ=255)
                                                                            tGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGAc                         >  1:2735316/1‑100 (MQ=255)
                                                                             gTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACt                        <  2:3253797/100‑1 (MQ=255)
                                                                              tAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTa                       <  1:2789233/100‑1 (MQ=255)
                                                                               aaCAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTaa                      <  1:1139257/100‑1 (MQ=255)
                                                                                 cAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAAc                    >  1:2055356/1‑100 (MQ=255)
                                                                                  caatctgaCAGCCAGTAATGGCTGTCAGATTGAAAATTTAGTGGTTAAAGAATTCTCCATCCAGATAAACTTTGCATTTTTCCACATGCATGCCTaaa                     >  2:1639560/6‑98 (MQ=37)
                                                                                    atctgaCAGCCAGTAATGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCAt                 >  1:37693/4‑100 (MQ=255)
                                                                                     gatgaCAGCCATTACTGGCTGTAAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATACATGACTCAACCATc                >  2:182923/1‑100 (MQ=255)
                                                                                      ctgaCAGCCAGTAATGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATcc               >  2:109093/2‑100 (MQ=255)
                                                                                      ctgaCAGCCAGTAATGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATcc               >  1:2241892/2‑100 (MQ=255)
                                                                                      ctgaCAGCCAGTAATGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATcc               >  1:195893/2‑100 (MQ=255)
                                                                                       tgaCAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATcc               <  2:2043394/99‑1 (MQ=255)
                                                                                       tgaCAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCa              >  1:2838026/1‑100 (MQ=255)
                                                                                       tgaCAGCCAGTAATGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCa              >  1:1683113/1‑100 (MQ=255)
                                                                                        gaCAGCCAGTAATGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAg             >  1:2622627/1‑100 (MQ=255)
                                                                                          cAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGAt           <  1:23534/100‑1 (MQ=255)
                                                                                           aGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATa          <  1:1657196/100‑1 (MQ=255)
                                                                                             ccATTACTGGCTGTTAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTCTCCAGATGCATGACTAAACCATCAAGATaa         <  2:2551362/99‑1 (MQ=255)
                                                                                             ccATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAAt        >  2:1901025/1‑100 (MQ=255)
                                                                                              cATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATc       >  2:2636180/1‑100 (MQ=255)
                                                                                                 tACTTGCTGTCAGATTGAATGTTTAGGTTCTCCTGAGTGCTCCATCCAGCTAAGCTTTGCCTTTTTCCACATGCATGACGAAACCATCCAGa            <  2:3587286/92‑1 (MQ=255)
                                                                                                 tACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAgg    <  2:1246093/100‑1 (MQ=255)
                                                                                                   cTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGtt  >  2:441907/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
TGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTT  >  NC_000913/2364899‑2365097

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: