Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,808,078 A→G intergenic (+6/‑157) yniC → / → b1728 putative phosphatase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,808,0780AG99.4% 994.9 / ‑3.7 311intergenic (+6/‑157)yniC/b1728putative phosphatase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (165/144);  total (166/145)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.50e-01

Only 100 of 322 total aligned reads displayed.

AGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACC  >  NC_000913/1807988‑1808167
                                                                                          |                                                                                         
 gcgcAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCgggg                                                                                        <  1:309217/93‑1 (MQ=255)
  cgcAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGgggg                                                                                       <  2:2826653/93‑1 (MQ=255)
    cAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCa                                                                                     <  2:5293813/93‑1 (MQ=255)
      aaaTGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGt                                                                                   <  1:7375573/93‑1 (MQ=255)
        aTGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGa                                                                                 <  1:8743889/93‑1 (MQ=255)
         tGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGAt                                                                                >  2:26202/1‑93 (MQ=255)
         tGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGAt                                                                                <  1:5348129/93‑1 (MQ=255)
          gATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGAtc                                                                               >  1:500494/1‑93 (MQ=255)
          gATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGAtc                                                                               >  2:8821800/1‑93 (MQ=255)
          gATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGAtc                                                                               <  1:6548318/93‑1 (MQ=255)
          gATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGAtc                                                                               >  2:4774195/1‑93 (MQ=255)
          gATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGAtc                                                                               <  2:1765161/93‑1 (MQ=255)
             ccACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCg                                                                            >  1:7752501/1‑93 (MQ=255)
             ccACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCg                                                                            >  2:8566412/1‑93 (MQ=255)
              cACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGc                                                                           >  2:9032710/1‑93 (MQ=255)
              cACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGc                                                                           <  2:4187782/93‑1 (MQ=255)
               aCGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCt                                                                          >  2:10042298/1‑93 (MQ=255)
                cGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTg                                                                         >  2:9654349/1‑93 (MQ=255)
                cGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTg                                                                         >  1:1815178/1‑93 (MQ=255)
                 gTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGc                                                                        >  1:6652010/1‑93 (MQ=255)
                 gTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGc                                                                        <  1:4348213/93‑1 (MQ=255)
                  ttttGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGcc                                                                       >  2:8210822/1‑93 (MQ=255)
                   tttGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGccc                                                                      >  1:1602506/1‑93 (MQ=255)
                    ttGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGcccc                                                                     <  1:958637/93‑1 (MQ=255)
                    ttGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGcccc                                                                     >  1:1944932/1‑93 (MQ=255)
                    ttGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGcccc                                                                     <  2:3417303/93‑1 (MQ=255)
                     tGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCt                                                                    <  2:940236/93‑1 (MQ=255)
                     tGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCt                                                                    <  2:4490570/93‑1 (MQ=255)
                     tGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCt                                                                    <  2:2770844/93‑1 (MQ=255)
                       tATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTgg                                                                  >  1:1500580/1‑93 (MQ=255)
                        aTTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGt                                                                 <  2:3918531/93‑1 (MQ=255)
                         ttAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGtt                                                                >  1:5755816/1‑93 (MQ=255)
                           aGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCt                                                              <  2:169511/93‑1 (MQ=255)
                           aGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCt                                                              <  1:3883088/93‑1 (MQ=255)
                             cAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCttt                                                            >  1:3320313/1‑93 (MQ=255)
                              aGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTa                                                           >  2:7428752/1‑93 (MQ=255)
                               gACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTAt                                                          >  2:6988250/1‑93 (MQ=255)
                                   tCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGa                                                      <  1:9477968/93‑1 (MQ=255)
                                    cAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGaa                                                     <  2:5490708/93‑1 (MQ=255)
                                    cAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGaa                                                     >  1:5142552/1‑93 (MQ=255)
                                     aaaCTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGaaa                                                    <  1:3738312/93‑1 (MQ=255)
                                     aaaCTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGaaa                                                    <  2:3186789/93‑1 (MQ=255)
                                     aaaCTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGaaa                                                    <  2:8311194/93‑1 (MQ=255)
                                      aaCTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAAt                                                   <  1:9947028/93‑1 (MQ=255)
                                      aaCTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAAt                                                   >  2:7576253/1‑93 (MQ=255)
                                       aCTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAAtt                                                  <  1:9131284/93‑1 (MQ=255)
                                       aCTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAAtt                                                  <  2:4743312/93‑1 (MQ=255)
                                        cTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTg                                                 >  1:4633330/1‑93 (MQ=255)
                                         tttCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGc                                                >  2:6718339/1‑93 (MQ=255)
                                           tCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCAt                                              >  1:9674851/1‑93 (MQ=255)
                                           tCATCGCTGACAGAACGCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCAt                                              >  2:9555795/1‑93 (MQ=255)
                                            cATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCAtt                                             <  1:8323733/93‑1 (MQ=255)
                                             aTCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTc                                            <  2:2924209/93‑1 (MQ=255)
                                             aTCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTc                                            >  1:7060794/1‑93 (MQ=255)
                                              tCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCa                                           <  1:2272646/93‑1 (MQ=255)
                                               cGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCaa                                          <  2:7154258/93‑1 (MQ=255)
                                                gCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAAc                                         <  1:7940740/93‑1 (MQ=255)
                                                gCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAAc                                         <  1:2397505/93‑1 (MQ=255)
                                                gCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAAc                                         >  2:7999219/1‑93 (MQ=255)
                                                 cTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGCAATTGCATTCAACt                                        >  2:659791/1‑93 (MQ=255)
                                                 cTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACt                                        >  1:2672683/1‑93 (MQ=255)
                                                    aCAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGAc                                     <  2:1405214/93‑1 (MQ=255)
                                                       gAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGa                                  >  2:6248282/1‑93 (MQ=255)
                                                       gAACTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGa                                  <  2:130714/93‑1 (MQ=255)
                                                         aCTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGAtt                                <  2:3499542/93‑1 (MQ=255)
                                                          cTCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTa                               >  2:8448730/1‑93 (MQ=255)
                                                           tCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTaa                              <  1:7538253/93‑1 (MQ=255)
                                                           tCACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTaa                              >  1:3573304/1‑93 (MQ=255)
                                                            cACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAAt                             >  1:2027858/1‑93 (MQ=255)
                                                            cACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAAt                             <  2:2051396/93‑1 (MQ=255)
                                                            cACCGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAAt                             >  2:3235964/1‑93 (MQ=255)
                                                             cccGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATc                            <  2:1974850/92‑1 (MQ=255)
                                                              ccGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGGGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTACTCg                           >  1:1092881/1‑93 (MQ=255)
                                                               cGCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGt                          >  1:2897115/1‑93 (MQ=255)
                                                                gCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTc                         >  2:427239/1‑93 (MQ=255)
                                                                gCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTc                         >  1:9192041/1‑93 (MQ=255)
                                                                gCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTc                         >  2:5691095/1‑93 (MQ=255)
                                                                gCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTc                         >  1:4194632/1‑93 (MQ=255)
                                                                gCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTc                         <  1:5202625/93‑1 (MQ=255)
                                                                gCAAAAGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTc                         >  1:586170/1‑93 (MQ=255)
                                                                   aaaGACCTTCTTGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAAt                      >  1:8068772/1‑93 (MQ=255)
                                                                   aaaGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAAt                      <  2:3696625/93‑1 (MQ=255)
                                                                   aaaGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAAt                      >  2:6519859/1‑93 (MQ=255)
                                                                    aaGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAAtt                     <  1:634484/93‑1 (MQ=255)
                                                                    aaGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAAtt                     >  1:6804005/1‑93 (MQ=255)
                                                                     aGACCTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAAttt                    <  1:4300141/93‑1 (MQ=255)
                                                                        ccTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAag                 <  2:1082369/93‑1 (MQ=255)
                                                                         cTTCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAaga                <  1:7600930/93‑1 (MQ=255)
                                                                          ttCTCGGTTAATGACCGGGGGCAGTGATCTCGCTGCGCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAagag               <  2:4442068/93‑1 (MQ=255)
                                                                            ctcGGTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGaa             <  2:875677/93‑1 (MQ=255)
                                                                                gTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAagag         >  1:6251177/1‑93 (MQ=255)
                                                                                gTTAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAagag         <  2:6314876/93‑1 (MQ=255)
                                                                                 ttAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGt        <  2:9489235/93‑1 (MQ=255)
                                                                                 ttAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGt        >  2:5105025/1‑93 (MQ=255)
                                                                                  tAATGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGtt       >  1:4379238/1‑93 (MQ=255)
                                                                                    aTGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTac     >  1:7211803/1‑93 (MQ=255)
                                                                                    aTGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTac     >  2:2489786/1‑93 (MQ=255)
                                                                                     tGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTGTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGATaca    >  2:9163252/1‑93 (MQ=255)
                                                                                     tGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTaca    <  2:7805555/93‑1 (MQ=255)
                                                                                     tGACCGGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTaca    >  2:9809911/1‑93 (MQ=255)
                                                                                          |                                                                                         
AGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACC  >  NC_000913/1807988‑1808167

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