Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,972,288 G→A T354M (ACG→ATG)  cheA ← CheA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,972,2880GA98.5% 468.9 / ‑3.8 135T354M (ACG→ATG) cheACheA
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (67/66);  total (68/67)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.56e-01

Only 100 of 141 total aligned reads displayed.

ATCGAGGCTATTGCGTACCAGGTGGGTCAGCGGGTCGATAATGCGTTCTATCAGGCTTTTGTCGAGTTCAGTAGAACTGCCCACCAGCGTCAGTTCTACCTGCTTGCCGAGTTTTCCCGCCAGATCACGCACCAGCCGGGGATAGCGACTAAAAACATATTCCATCGGCATCATGC  >  NC_000913/1972200‑1972375
                                                                                        |                                                                                       
aTCGAGGCTATTGCGTACCAGGTGGGTCAGCGGGTCGATAATGCGTTCTATCAGGCTTTTGTCGAGTTCAGTAGAACTGCCCACCAGCATCAg                                                                                     <  2:3004098/93‑1 (MQ=255)
aTCGAGGCTATTGCGTACCAGGTGGGTCAGCGGGTCGATAATGCGTTCTATCAGGCTTTTGTCGAGTTCAGTAGAACTGCCCACCAGCATCAg                                                                                     <  2:279570/93‑1 (MQ=255)
  cGAGGCTATTGCGTACCAGGTGGGTCAGCGGGTCGATAATGCGTTCTATCAGGCTTTTGTCGAGTTCAGTAGAACTGCCCACCAGCATCAGtt                                                                                   >  2:3091792/1‑93 (MQ=255)
   gAGGCTATTGCGTACCAGGTGGGTCAGCGGGTCGATAATGCGTTCTATCAGGCTTTTGTCGAGTTCAGTAGAACTGCCCACCAGCATCAGTTc                                                                                  <  2:4285735/93‑1 (MQ=255)
   gAGGCTATTGCGTACCAGGTGGGTCAGCGGGTCGATAATGCGTTCTATCAGGCTTTTGTCGAGTTCAGTAGAACTGCCCACCAGCATCAGTTc                                                                                  <  1:4419519/93‑1 (MQ=255)
     ggCTATTGCGTACCAGGTGGGTCAGCGGGTCGATAATGCGTTCTATCAGGCTTTTGTCGAGTTCAGTAGAACTGCCCACCAGCATCAGTTCTa                                                                                >  2:479805/1‑93 (MQ=255)
      gCTATTGCGTACCAGGTGGGTCAGCGGGTCGATAATGCGTTCTATCAGGCTTTTGTCGAGTTCAGTAGAACTGCCCACCAGCATCAGTTCTAc                                                                               >  1:373024/1‑93 (MQ=255)
      gCTATTGCGTACCAGGTGGGTCAGCGGGTCGATAATGCGTTCTATCAGGCTTTTGTCGAGTTCAGTAGAACTGCCCACCAGCATCAGTTCTAc                                                                               <  2:530754/93‑1 (MQ=255)
        tATTGCGTACCAGGTGGGTCAGCGGGTCGATAATGCGTTCTATCAGGCTTTTGTCGAGTTCAGTAGAACTGCCCACCAGCATCAGTTCTACCt                                                                             <  1:221276/93‑1 (MQ=255)
            gCGTACCAGGTGGGTCAGCGGGTCGATAATGCGTTCTATCAGGCTTTTGTCGAGTTCAGTAGAACTGCCCACCAGCATCAGTTCTACCTGCtt                                                                         <  2:3625575/93‑1 (MQ=255)
            gCGTACCAGGTGGGTCAGCGGGTCGATAATGCGTTCTATCAGGCTTTTGTCGAGTTCAGTAGAACTGCCCACCAGCATCAGTTCTACCTGCtt                                                                         >  2:2798321/1‑93 (MQ=255)
             cGTACCAGGTGGGTCAGCGGGTCGATAATGCGTTCTATCAGGCTTTTGTCGAGTTCAGTAGAACTGCCCACCAGCATCAGTTCTACCTGCTTg                                                                        >  2:3465716/1‑93 (MQ=255)
              gTACCAGGTGGGTCAGCGGGTCGATAATGCGTTCTATCAGGCTTTTGTCGAGTTCAGTAGAACTGCCCACCAGCATCAGTTCTACCTGCTTGc                                                                       >  2:2638144/1‑93 (MQ=255)
              gTACCAGGTGGGTCAGCGGGTCGATAATGCGTTCTATCAGGCTTTTGTCGAGTTCAGTAGAACTGCCCACCAGCATCAGTTCTACCTGCTTGc                                                                       >  2:1357798/1‑93 (MQ=255)
               tcccAGGTTGTTCAGCGGGTCGATAATGCGTTCTATCAGGCTTTTGTCGAGTTCAGTAGAACTGCCCCCCAGCATCAGTTCTACCTGCTTGcc                                                                      <  1:23269/91‑1 (MQ=255)
               tACCAGGTGGGTCAGCGGGTCGATAATGCGTTCTATCAGGCTTTTGTCGAGTTCAGTAGAACTGCCCACCAGCATCAGTTCTACCTGCTTGcc                                                                      >  1:2189024/1‑93 (MQ=255)
               tACCAGGTGGGTCAGCGGGTCGATAATGCGTTCTATCAGGCTTTTGTCGAGTTCAGTAGAACTGCCCACCAGCATCAGTTCTACCTGCTTGcc                                                                      <  2:1340296/93‑1 (MQ=255)
                  cAGGTGGGTCAGCGGGTCGATAATGCGTTCTATCAGGCTTTTGTCGAGTTCAGTAGAACTGCCCACCAGCATCAGGTCTATCTGCTTGCCGAg                                                                   >  2:3592100/1‑93 (MQ=255)
                    ggTGGGTCAGCGGGTCGGTCATTCGTTCTATCAGGGTTTTGTCAAGTTCACTAGTACTGCCCACCAGCATCAGTTCTACCGGCTTGGCGAGtt                                                                 >  2:3605671/1‑93 (MQ=255)
                     gTGGGTCAGCGGGTCGATAATGCGTTCTATCAGGCTTTTGTCGAGTTCAGTAGAACTGCCCACCAGCATCAGTTCTACCTGCTTGCCGAGttt                                                                >  2:395800/1‑93 (MQ=255)
                          tCAGCGGGTCGATAATGCGTTCTATCAGACTTTTGTCGAGTTCAGTAGAACTGCCCACCAGCATCAGTTCTACCTGCTTGCCGAGTTTTCCCg                                                           <  2:3026753/93‑1 (MQ=255)
                             gCGGGTCGATAATGCGTTCTATCAGGCTTTTGTCGAGTTCAGTAGAACTGCCCACCAGCATCAGTTCTACCTGCTTGCCGAGTTTTCCCGCCa                                                        <  2:1249258/93‑1 (MQ=255)
                             gCGGGTCGATAATGCGTTCTATCAGGCTTTTGTCGAGTTCAGTAGAACTGCCCACCAGCATCAGTTCTACCTGCTTGCCGAGTTTTCCCGCCa                                                        <  1:3970475/93‑1 (MQ=255)
                             gCGGGTCGATAATGCGTTCTATCAGGCTTTTGTCGAGTTCAGTAGAACTGCCCACCAGCATCAGTTCTACCTGCTTGCCGAGTTTTCCCGCCa                                                        >  2:2010051/1‑93 (MQ=255)
                               gggTCGATAATGCGTTCTATCAGGCTTTTGTCGAGTTCAGTAGAACTGCCCACCAGCATCAGTTCTACCTGCTTGCCGAGTTTTCCCGCCAGa                                                      >  1:1658262/1‑93 (MQ=255)
                                 gTCGATAATGCGTTCTATCAGGCTTTTGTCGAGTTCAGTAGAACTGCCCACCAGCATCAGTTCTACCTGCTTGCCGAGTTTTCCCGCCAGATc                                                    <  2:2149082/93‑1 (MQ=255)
                                 gTCGATAATGCGTTCTATCAGGCTTTTGTCGAGTTCAGTAGAACTGCCCACCAGCATCAGTTCTACCTGCTTGCCGAGTTTTCCCGCCAGATc                                                    <  1:491381/93‑1 (MQ=255)
                                 gTCGATAATGCGTTCTATCAGGCTTTTGTCGAGTTCAGTAGAACTGCCCACCAGCATCAGTTCTACCTGCTTGCCGAGTTTTCCCGCCAGATc                                                    >  2:25651/1‑93 (MQ=255)
                                 gTCGATAATGCGTTCTATCAGGCTTTTGTCGAGTTCAGTAGAACTGCCCACCAGCATCAGTTCTACCTGCTTGCCGAGTTTTCCCGCCAGATc                                                    >  1:3622716/1‑93 (MQ=255)
                                  tCGATAATGCGTTCTATCAGGCTTTTGTCGAGTTCAGTAGAACTGCCCACCAGCATCAGTTCTACCTGCTTGCCGAGTTTTCCCGCCAGATCa                                                   <  1:1243219/93‑1 (MQ=255)
                                    gATAATGCGTTCTATCAGGCTTTTGTCGAGTTCAGTAGAACTGCCCACCAGCATCAGTTCTACCTGCTTGCCGAGTTTTCCCGCCAGATCACg                                                 >  1:2034370/1‑93 (MQ=255)
                                    gATAATGCGTTCTATCAGGCTTTTGTCGAGTTCAGTAGAACTGCCCACCAGCATCAGTTCTACCTGCTTGCCGAGTTTTCCCACCAGATCACg                                                 >  1:1677289/1‑93 (MQ=255)
                                     aTAATGCGTTCTATCAGGCTTTTTTCGAGTTCAGTAGAACTGCCCACCAGCATCAGTTCTACCTGCTTGCCGAGTTTTCCCGCCAGATCACGc                                                >  2:3782989/1‑93 (MQ=255)
                                     aTAATGCGTTCTATCAGGCTTTTGTCGAGTTCAGTAGAACTGCCCACCAGCATCAGTTCTACCTGCTTGCCGAGTTTTCCCGCCAGATCACGc                                                >  1:204853/1‑93 (MQ=255)
                                     aTAATGCGTTCTATCAGGCTTTTGTCGAGTTCAGTAGAACTGCCCACCAGCATCAGTTCTACCTGCTTGCCGAGTTTTCCCGCCAGATCACGc                                                <  1:160152/93‑1 (MQ=255)
                                      tAATGCGTTCTATCAGGCTTTTGTCGAGTTCAGTAGAACTGCCCACCAGCATCAGTTCTACCTGCTTGCCGAGTTTTCCCGCCAGATCACGCa                                               >  1:1947806/1‑93 (MQ=255)
                                      tAATGCGTTCTATCAGGCTTTTGTCGAGTTCAGTAGAACTGCCCACCAGCATCAGTTCTACCTGCTTGCCGAGTTTTCCCGCCAGATCACGCa                                               >  2:504886/1‑93 (MQ=255)
                                      tAATGCGTTCTATCAGGCTTTTGTCGAGTTCAGTAGAACTGCCCACCAGCATCAGTTCTACCTGCTTGCCGAGTTTTCCCGCCAGATCACGCa                                               >  2:2371985/1‑93 (MQ=255)
                                       aaTGCGTTCTATCAGGCTTTTGTCGAGTTCAGTAGAACTGCCCACCAGCATCAGTTCTACCTGCTTGCCGAGTTTTCCCGCCAGATCACGCAc                                              >  1:2255369/1‑93 (MQ=255)
                                       aaTGCGTTCTATCAAGCTTTTGTCGAGTTCAGCAGTACTGCTCACCAGCATCATGTTTACCTGGTTGCCGGCGCTTTCCGCCAGAACACGCAc                                              >  2:3954186/1‑93 (MQ=255)
                                        aTGCGTTCTATCAGGCTTTTGTCGAGTTCAGTAGAACTGCCCACCAGCATCAGTTCTACCTGCTTGCCGAGTTTTCCCGCCAGATCACGCAcc                                             <  2:3184615/93‑1 (MQ=255)
                                        aTGCGTTCTATCAGGCTTTTGTCGAGTTCAGTAGAACTGCCCACCAGCATCAGTTCTACCTGCTTGCCGAGTTTTCCCGCCAGATCACGCAcc                                             <  1:3975606/93‑1 (MQ=255)
                                        aTGCGTTCTATCAGGCTTTTGTCGAGTTCAGTAGAACTGCCCACCAGCATCAGTTCTACCTGCTTGCCGAGTTTTCCCGCCAGATCACGCAcc                                             >  2:4401183/1‑93 (MQ=255)
                                        aTGCGTTCTATCAGGCTTTTGTCGAGTTCAGTAGAACTGCCCACCAGCATCAGTTCTACCTGCTTGCCGAGTTTTCCCGCCAGATCACGCAcc                                             >  1:1220956/1‑93 (MQ=255)
                                        aTGCGTTCTATCAGGCTTTTGTCGAGTTCAGTAGAACTGCCCACCAGCATCAGTTCTACCTGCTTGCCGAGTTTTCCCGCCAGATCACGCAcc                                             >  1:1941068/1‑93 (MQ=255)
                                        aTGCGTTCTATCAGGCTTTTGTCGAGTTCAGTAGAACTGCCCACCAGCATCAGTTCTACCTGCTTGCCGAGTTTTCCCGCCAGATCACGCAcc                                             <  1:4307148/93‑1 (MQ=255)
                                         tGCGTTCTATCAGGCTTTTGTCGAGTTCAGTAGAACTGCCCACCAGCATCAGTTCTACCTGCTTGCCGAGTTTTCCCGCCAGATCACGCACCa                                            >  1:4094950/1‑93 (MQ=255)
                                          gCGTTCTATCAGGCTTTTGTCGAGTTCAGTAGAACTGCCCACCAGCATCAGTTCTACCTGCTTGCCGAGTTTTCCCGCCAGATCACGCACCAg                                           <  1:3731132/93‑1 (MQ=255)
                                           cGTTCTATCAGGCTTTTGTCGAGTTCAGTAGAACTGCCCACCAGCATCAGTTCTACCTGCTTGCCGAGTTTTCCCGCCAGATCACGCACCAGc                                          >  1:2993724/1‑93 (MQ=255)
                                            gTTCTATCAGGCTTTTGTCGAGTTCAGTAGAACTGCCCACCAGCATCAGTTCTACCTGCTTGCCGAGTTTTCCCGCCAGATCACGCACCAGcc                                         >  2:2194497/1‑93 (MQ=255)
                                             ttCTATCAGGCTTTTGTCGAGTTCAGTAGAACTGCCCACCAGCATCAGTTCTACCTGCTTGCCGAGTTTTCCCGCCAGATCACGCACCAGCCg                                        <  1:4582733/93‑1 (MQ=255)
                                              tCTATCAGGCTTTTGTCGAGTTCAGTAGAACTGCCCACCAGCATCAGTTCTACCTGCTTGCCGAGTTTTCCCGCCAGATCACGCACCAGCCgg                                       <  2:1191922/93‑1 (MQ=255)
                                                tATCAGGCTTTTGTCGAGTTCAGTAGAACTGCCCACCAGCATCAGTTCTACCTGCTTGCCGAGTTTTTCCGCCAGATGACGCACCAGCCggga                                     >  1:257130/1‑92 (MQ=255)
                                                tATCAGGCTTTTGTCGAGTTCAGTAGAACTGCCCACCAGCATCAGTTCTACCTGCTTGCCGAGTTTTCCCGCCAGATCACGCACCAGCCgggg                                     <  1:3392913/93‑1 (MQ=255)
                                                 aTCAGGCTTTTGTCGAGTTCAGTAGAACTGCCCACCAGCATCAGTTCTACCTGCTTGCCGAGTTTTCCCGCCAGATCACGCACCAGCCGGGGa                                    >  2:386776/1‑93 (MQ=255)
                                                  tCAGGCTTTTGTCGAGTTCAGTAGAACTGCCCACCAGCATCAGTTCTACCTGCTTGCCGAGTTTTCCCGCCAGATCACGCACCAGCCGGGGAt                                   >  2:238809/1‑93 (MQ=255)
                                                   cAGGCTTTTGTCGAGTTCAGTAGAACTGCCCACCAGCATCAGTTCTACCTGCTTGCCGAGTTTTCCCGCCAGATCACGCACCAGCCGGGGATa                                  >  1:1332089/1‑93 (MQ=255)
                                                   cAGGCTTTTGTCGAGTTCAGTAGAACTGCCCACCAGCATCAGTTCTACCTGCTTGCCGAGTTTTCCCGCCAGATCACGCACCAGCCGGGGATa                                  <  1:3106304/93‑1 (MQ=255)
                                                   cAGGCTTTTGTCGAGTTCAGTAGAACTGCCCACCAGCATCAGTTCTACCTGCTTGCCGACTTTTCCCGCCAGATCACGCACCAGCCGGGGATa                                  <  1:3161804/93‑1 (MQ=255)
                                                    aGGCTTTTGTCGAGTTCAGTAGAACTGCCCACCAGCATCAGTTCTACCTGCTTGCCGAGTTTTCCCGCCAGATCACGCACCAGCCGGGGATAg                                 <  2:4120479/93‑1 (MQ=255)
                                                     ggCTTTTGTCGAGTTCAGTAGAACTGCCCACCAGCATCAGTTCTACCTGCTTGCCGAGTTTTCCCGCCAGATCACGCACCAGCCGGGGATAGc                                >  1:651454/1‑93 (MQ=255)
                                                     ggCTTTTGTCGAGTTCAGTAGAACTGCCCACCAGCATCAGTTCTACCTGCTTGCCGAGTTTTCCCGCCAGATCACGCACCAGCCGGGGATAGc                                <  1:2631871/93‑1 (MQ=255)
                                                     ggCTTTTGTCGAGTTCAGTAGAACTGCCCACCAGCATCAGTTCTACCTGCTTGCCGAGTTTTCCCGCCAGATCACGCACCAGCCGGGGATAGc                                <  2:2189688/93‑1 (MQ=255)
                                                        ttttGTCGAGTTCAGTAGAACTGCCCACCAGCATCAGTTCTACCTGCTTGCCGAGTTTTCCCGCCAGATCACGCACTAGCCGGGGATAGCGAc                             >  1:1322829/1‑93 (MQ=255)
                                                        ttttGTCGAGTTCAGTAGAACTGCCCACCAGCATCAGTTCTACCTGCTTGCCGAGTTTTCCCGCCAGATCACGCACCAGCCGGGGATAGCGAc                             <  1:2312107/93‑1 (MQ=255)
                                                        ttttGTCGAGTTCAGTAGAACTGCCCACCAGCATCAGTTCTACCTGCTTGCCGAGTTTTCCCGCCAGATCACGCACCAGCCGGGGATAGCGAc                             <  1:3113002/93‑1 (MQ=255)
                                                         tgtGTCGAGTTCAGTAGAACTGCCCACCAGCATCAGTTCTACCTGCTTGCCGAGTTTTCCCGCCAGATCACGCACCAGCCGGGGATAGCGACt                            <  1:2906022/91‑1 (MQ=255)
                                                           tGTCGAGTTCAGTAGAACTGCCCACCAGCATCAGTTCTACCTGCTTGCCGAGTTTTCCCGCCAGATCACGCACCAGCCGGGGATAGCGACTaa                          <  1:4715679/93‑1 (MQ=255)
                                                               gAGTTCAGTAGAACTGCCCACCAGCATCAGTTCTACCTGCTTGCCGAGTTTTCCCGCCAGATCACGCACCAGCCGGGGATAGCGACTAAAAAc                      >  1:1472698/1‑93 (MQ=255)
                                                               gAGTTCAGTAGAACTGCCCACCAGCATCAGTTCTACCTGCTTGCCGAGTTTTCCCGCCAGATCACGCACCAGCCGGGGATAGCGACTAAAAAc                      >  2:2547139/1‑93 (MQ=255)
                                                               gAGTTCAGTAGAACTGCCCACCAGCATCAGTTCTACCTGCTTGCCGAGTTTTCCCGCCAGATCACGCACCAGCCGGGGATAGCGACTAAAAAc                      >  1:2837983/1‑93 (MQ=255)
                                                                aGTTCAGTAGAACTGCCCACCAGCATCAGTTCTACCTGCTTGCCGAGTTTTCCCGCCAGATCACGCACCAGCCGGGGATAGCGACTAAAAACa                     <  2:957441/93‑1 (MQ=255)
                                                                 gTTCAGTAGAACTGCCCACCAGCATCAGTTCTACCTGCTTGCCGAGTTTTCCCGCCAGATCACGCACCAGCCGGGGATAGCGACTAAAAACat                    <  1:1120881/93‑1 (MQ=255)
                                                                 gTTCAGTAGAACTGCCCACCAGCATCAGTTCTACCTGCTTGCCGAGTTTTCCCGCCAGATCACGCACCAGCCGGGGATAGCGACTAAAAACat                    >  2:1008212/1‑93 (MQ=255)
                                                                 gTTCAGTAGAACTGCCCACCAGCATCAGTTCTACCTGCTTGCCGAGTTTTCCCGCCAGATCACGCACCAGCCGGGGATAGCGACTAAAAACat                    >  1:1543420/1‑93 (MQ=255)
                                                                 gTTCAGTAGAACTGCCCACCAGCATCAGTTCTACCTGCTTGCCGAGTTTTCCCGCCAGATCACGCACCAGCCGGGGATAGCGACTAAAAACat                    <  2:4418357/93‑1 (MQ=255)
                                                                 gTTCAGTAGAACTGCCCACCAGCATCAGTTCTACCTGCTTGCCGAGTTTTCCCGCCAGATCACGCACCAGCCGGGGATAGCGACTAAAAACat                    >  2:4341332/1‑93 (MQ=255)
                                                                 gTTCAGTAGAACTGCCCACCAGCATCAGTTCTACCTGCTTGCCGAGTTTTCCCGCCAGATCACGCACCAGCCGGGGATAGCGACTAAAAACat                    <  2:4303082/93‑1 (MQ=255)
                                                                 gTTCAGTAGAACTGCCCACCAACATCAGTTCTACCTGCTTGCCGAGTTTTCCCGCCAGATCACGCACCAGCCGGGGATAGCGACTAAAAACat                    >  1:4159834/1‑93 (MQ=255)
                                                                   tCAGTAGAACTGCCCACCAGCATCAGTTCTACCTGCTTGCCGAGTTTTCCCGCCAGATCACGCACCAGCCGGGGATAGCGACTAAAAACatat                  <  1:4635420/93‑1 (MQ=255)
                                                                   tCAGTAGAACTGCCCACCAGCATCAGTTCTACCTGCTTGCCGAGTTTTCCCGCCAGATCACGCACCAGCCGGGGATAGCGACTAAAAACatat                  <  2:3997067/93‑1 (MQ=255)
                                                                    cAGTAGAGCTGCCCACCAGCATCAGTTCTACCTGCTTGCCGAGTTTTCCCGCCAGATCACGCACCAGCCGGGGATAGCGACTAAAAACATAtt                 <  2:131548/93‑1 (MQ=255)
                                                                    cAGTAGAACTGCCCCCCAGCATCAGTTCTACCTGCTTGCCGAGTTTTCCCGCCAGATCACGCACCAGCCGGGGATAGCGACTAAAAACATAtt                 <  1:4509783/93‑1 (MQ=255)
                                                                        aGAACTGCCCACCAGCATCAGTTCTACCTGCTTGCCGAGTTTTCCCGCCAGATCACGCACCAGCCGGGGATAGCGACTAAAAACATATTCCAt             >  2:4375852/1‑93 (MQ=255)
                                                                        aGAACTGCCCACCAGCATCAGTTCTACCTGCTTGCCGAGTTTTCCCGCCAGATCACGCACCAGCCGGGGATAGCGACTAAAAACATATTCCAt             <  1:745768/93‑1 (MQ=255)
                                                                        aGAACTGCCCACCAGCATCAGTTCTACCTGCTTGCCGAGTTTTCCCGCCAGATCACGCACCAGCCGGGGTAGCGACTAAAAACATATTCCATc            >  2:1483596/1‑93 (MQ=255)
                                                                         gAACTGCCCACCAGCATCAGTTCTACCTGCTTGCCGAGTTTTCCCGCCAGATCACGCACCAGCCGGGGATAGCGACTAAAAACATATTCCATc            <  2:716609/93‑1 (MQ=255)
                                                                          aaCTGCCCACCAGCATCAGTTCTCCCTGCTTGCCGAGTTTTCCCGCAAGATCACGCACCAGCCGGGGATAGCGACTAAAAACATATTCCATCg           <  1:1868347/93‑1 (MQ=255)
                                                                           aCTGCCCACCAGCATCAGTTCTACCTGCTTGCCGAGTTTTCCCGCCAGATCACGCACCAGCCGGGGATAGCGACTAAAAACATATTCCATCgg          >  2:4235315/1‑93 (MQ=255)
                                                                           aCTGCCCACCAGCATCAGTTCTACCTGCTTGCCGAGTTTTCCCGCCAGATCACGCACCAGCCGGGGATAGCGACTAAAAACATATTCCATCgg          <  2:942144/93‑1 (MQ=255)
                                                                            cTGCCCACCAGCATCAGTTCTACCTGCTTGCCGAGTTTTCCCGCCAGATCACGCACCAGCCGGGGATAGCGACTAAAAACATATTCCATCgtc         >  1:4586801/1‑91 (MQ=255)
                                                                               cccACCAGCATCAGTTCTACCTGCTTGCCGAGTTTTCCCGCCAGATCACGCACCAGCCGGGGATAGCGACTAAAAACATATTCCATCGGcatc      <  1:2844319/93‑1 (MQ=255)
                                                                                  cccaGCATCAGTTCTACCTGCTTGCCGAGTTTTCCCGCCAGATCACGCACCAGCCGGGGATAGCGACTAAAAACATATTCCATCGGCATCATg   <  1:3368934/92‑1 (MQ=255)
                                                                                  accaGCATCAGTTCTACCTGCTTGCCGAGTTTTCCCGCCAGATCACGCACCAGCCGGGGATAGCGACTAAAAACATATTCCATCGGCATCATg   >  2:773985/1‑93 (MQ=255)
                                                                                  accaGCATCAGTTCTACCTGCTTGCCGAGTTTTCCCGCCAGATCACGCACCAGCCGGGGATAGCGACTAAAAACATATTCCATCGGCATCATg   >  1:4386995/1‑93 (MQ=255)
                                                                                   ccaGCATCAGTTCTACCTGCTTGCCGAGTTTTCCCGCCAGATCACGCACCAGCCGGGGATAGCGACTAAAAACATATTCCATCGGCATCATGc  >  1:587274/1‑93 (MQ=255)
                                                                                   ccaGCATCAGTTCTACCTGCTTGCCGAGTTTTCCCGCCAGATCACGCACCAGCCGGGGATAGCGACTAAAAACATATTCCATCGGCATCATGc  <  2:4525125/93‑1 (MQ=255)
                                                                                   ccaGCATCAGTTCTACCTGCTTGCCGAGTTTTCCCGCCAGATCACGCACCAGCCGGGGATAGCGACTAAAAACATATTCCATCGGCATCATGc  >  1:296378/1‑93 (MQ=255)
                                                                                   ccaGCATCAGTTCTACCTGCTTGCCGAGTTTTCCCGCCAGATCACGCACCAGCCGGGGATAGCGACTAAAAACATATTCCATCGGCATCATGc  >  1:587273/1‑93 (MQ=255)
                                                                                   ccaGCATCAGTTCTACCTGCTTGCCGAGTTTTCCCGCCAGATCACGCACCAGCCGGGGATAGCGACTAAAAACATATTCCATCGGCATCATGc  <  1:4392865/93‑1 (MQ=255)
                                                                                        |                                                                                       
ATCGAGGCTATTGCGTACCAGGTGGGTCAGCGGGTCGATAATGCGTTCTATCAGGCTTTTGTCGAGTTCAGTAGAACTGCCCACCAGCGTCAGTTCTACCTGCTTGCCGAGTTTTCCCGCCAGATCACGCACCAGCCGGGGATAGCGACTAAAAACATATTCCATCGGCATCATGC  >  NC_000913/1972200‑1972375

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