Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,787,068 G→A A204T (GCC→ACC)  aroH → 3‑deoxy‑D‑arabinoheptulosonate‑7‑phosphate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,787,0680GA99.1% 379.1 / NA 117A204T (GCC→ACC) aroH3‑deoxy‑D‑arabinoheptulosonate‑7‑phosphate synthase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  major base A (58/57);  minor base C (0/1);  total (59/58)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 137 total aligned reads displayed.

ATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGT  >  NC_000913/1786978‑1787157
                                                                                          |                                                                                         
aTGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGGTGTGGATGCTATCCGCGCAGCCCGCacc                                                                                         >  2:1712997/1‑93 (MQ=255)
   gCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCaccagc                                                                                      <  1:3466768/93‑1 (MQ=255)
    cTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCaccagcc                                                                                     >  1:2819696/1‑93 (MQ=255)
     ttCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCaccagcca                                                                                    >  1:278132/1‑93 (MQ=255)
     ttCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTCTGGATGCTATCCGCGCAGCCCGCaccagcca                                                                                    >  1:496293/1‑93 (MQ=255)
      tCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCat                                                                                   >  1:131837/1‑93 (MQ=255)
       cGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCata                                                                                  <  1:2066438/93‑1 (MQ=255)
       cGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCata                                                                                  <  2:3461403/93‑1 (MQ=255)
       cGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCata                                                                                  <  2:1873360/93‑1 (MQ=255)
        ggCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCacat                                                                                 >  1:3072947/1‑93 (MQ=255)
         gCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATg                                                                                >  1:2793962/1‑93 (MQ=255)
          cACGCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGt                                                                               <  2:2535894/93‑1 (MQ=255)
               tcCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCctc                                                                          >  1:1043959/1‑93 (MQ=255)
                ccTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTACCCGCGCAGCCCGCACCAGCCATATGTTCctct                                                                         <  2:3507265/93‑1 (MQ=255)
                 cTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCtctc                                                                        >  2:3942606/1‑93 (MQ=255)
                  tGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCg                                                                       >  1:4514384/1‑93 (MQ=255)
                  tGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCg                                                                       >  2:2230002/1‑93 (MQ=255)
                  tGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTATGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCg                                                                       >  2:581831/1‑93 (MQ=255)
                   gTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGc                                                                      >  1:1251099/1‑93 (MQ=255)
                    tCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGcc                                                                     >  1:2478321/1‑93 (MQ=255)
                       ggTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGGCAGa                                                                  >  2:2953219/1‑93 (MQ=255)
                         tAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACa                                                                <  1:1255885/93‑1 (MQ=255)
                           ggTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACaaa                                                              <  1:2391518/93‑1 (MQ=255)
                             ttttAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACaaaaa                                                            <  1:3265408/93‑1 (MQ=255)
                             ttttAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACaaaaa                                                            >  2:2105887/1‑93 (MQ=255)
                             ttttAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACaaaaa                                                            >  2:388202/1‑93 (MQ=255)
                               ttAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATg                                                          >  2:4346492/1‑93 (MQ=255)
                                 aaaaaTGGTACCGGTGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGt                                                        <  2:2886662/93‑1 (MQ=255)
                                 aaaaaTGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGt                                                        <  2:3496738/93‑1 (MQ=255)
                                  aaaaTGGTCCCGATGGCAATACGCGGATTGCTGTGGATTCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTc                                                       <  1:2643003/93‑1 (MQ=255)
                                   aaaTGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCa                                                      >  1:2911307/1‑93 (MQ=255)
                                     agggTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGa                                                    <  2:1709535/91‑1 (MQ=255)
                                     aTGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGa                                                    <  1:1835327/93‑1 (MQ=255)
                                      tGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGAt                                                   >  1:232660/1‑93 (MQ=255)
                                        gTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGa                                                 >  1:4499816/1‑93 (MQ=255)
                                          aCCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGAcc                                               >  2:4399886/1‑93 (MQ=255)
                                          aCCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGAcc                                               <  1:240290/93‑1 (MQ=255)
                                          aCCGATGGAAATACGCGGGTTGCTGTGGATGCTAGCCGCGCACCCCCCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGCTGAcc                                               <  1:1641451/93‑1 (MQ=255)
                                           ccGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCa                                              <  1:4500326/93‑1 (MQ=255)
                                           ccGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCACACAAAAATGGTCAGATGACCa                                              >  1:3743367/1‑93 (MQ=255)
                                             gATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATc                                            >  2:1361786/1‑93 (MQ=255)
                                             gATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATc                                            <  2:1321463/93‑1 (MQ=255)
                                             gATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATc                                            >  1:516759/1‑93 (MQ=255)
                                             cgtGGCAATGCGCAGATTGCTGTGGATGCTATCCGCACAGCCCGCCCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATc                                            <  2:677824/91‑1 (MQ=255)
                                              aTGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCt                                           <  2:2961563/93‑1 (MQ=255)
                                               tGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTa                                          <  1:2040362/93‑1 (MQ=255)
                                                ggCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTAt                                         <  1:1256351/93‑1 (MQ=255)
                                                 gCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATc                                        <  1:4647661/93‑1 (MQ=255)
                                                  cAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCa                                       <  2:253307/93‑1 (MQ=255)
                                                   aaTACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAg                                      >  1:4368737/1‑93 (MQ=255)
                                                   aaTACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAg                                      >  1:18512/1‑93 (MQ=255)
                                                   aaTACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAg                                      >  1:2695369/1‑93 (MQ=255)
                                                    aTACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGa                                     >  2:834060/1‑93 (MQ=255)
                                                    aTACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGa                                     >  1:4376249/1‑93 (MQ=255)
                                                    aTACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGa                                     >  1:1274559/1‑93 (MQ=255)
                                                     tACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGAc                                    <  1:434289/93‑1 (MQ=255)
                                                       cgcgGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCa                                  <  1:2583949/93‑1 (MQ=255)
                                                        gcgGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAg                                 <  1:558281/93‑1 (MQ=255)
                                                        gcgGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAg                                 <  2:4367209/93‑1 (MQ=255)
                                                         cgGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGc                                <  1:708444/93‑1 (MQ=255)
                                                         cgGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGc                                <  2:2679246/93‑1 (MQ=255)
                                                          ggATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCg                               >  2:1453063/1‑93 (MQ=255)
                                                          ggATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCg                               <  2:2854423/93‑1 (MQ=255)
                                                           gATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCgg                              >  1:689266/1‑93 (MQ=255)
                                                           gATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCgg                              <  1:549988/93‑1 (MQ=255)
                                                             ttGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCa                            >  1:334914/1‑93 (MQ=255)
                                                              tGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCaa                           <  1:238285/93‑1 (MQ=255)
                                                               gCTGTGGATGCTATCCGCGCAGCCCGCAACAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAAc                          >  1:2256424/1‑93 (MQ=255)
                                                                cTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAAcc                         <  1:124285/93‑1 (MQ=255)
                                                                  gtgGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCg                       <  1:3856939/93‑1 (MQ=255)
                                                                   tgGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGt                      >  2:25578/1‑93 (MQ=255)
                                                                    ggATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTa                     <  2:2045154/93‑1 (MQ=255)
                                                                     gATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTAt                    <  2:426941/93‑1 (MQ=255)
                                                                       tGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATgg                  <  2:2326499/93‑1 (MQ=255)
                                                                       tGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATgg                  >  2:1390405/1‑93 (MQ=255)
                                                                        gCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGc                 >  2:607010/1‑93 (MQ=255)
                                                                         cTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGcc                <  1:1744323/93‑1 (MQ=255)
                                                                          tATCCGCGCAGCCCCCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCca               <  2:3377320/93‑1 (MQ=255)
                                                                           aTCTGCGCGGCCCATCCCAGCTATATGTTCCTCTCGCCAGACAAAAAGGGTCAGATTACCATCTAGCAGACCAGCGGCAACCCGTATGGCcac              <  2:1302967/93‑1 (MQ=255)
                                                                            tCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCcaca             >  2:1759323/1‑93 (MQ=255)
                                                                             ccGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAAACCGTATGGCCACat            <  1:1684330/93‑1 (MQ=255)
                                                                              cgcgCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACatt           <  1:4544298/93‑1 (MQ=255)
                                                                                cgcAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACattat         <  2:883683/93‑1 (MQ=255)
                                                                                 gcAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACattatt        >  2:2157365/1‑93 (MQ=255)
                                                                                 gcAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCTACCTGTATGGCCACattatt        >  1:4070992/1‑93 (MQ=255)
                                                                                 gcAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACattatt        <  1:1917825/93‑1 (MQ=255)
                                                                                 gcAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACattatt        >  1:3437833/1‑93 (MQ=255)
                                                                                 gcAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACattatt        >  1:1042342/1‑93 (MQ=255)
                                                                                 gcAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACattatt        >  2:3420434/1‑93 (MQ=255)
                                                                                  cAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACAttatta       <  2:4541398/93‑1 (MQ=255)
                                                                                   aGCCCGCCCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCAtattattat      <  2:3172936/93‑1 (MQ=255)
                                                                                   aGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATtattat      >  1:3882076/1‑93 (MQ=255)
                                                                                    gCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCGACCCGTATGGCCACATTATTATg     >  1:993142/1‑93 (MQ=255)
                                                                                    gCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATg     >  2:1561539/1‑93 (MQ=255)
                                                                                    gCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATg     >  1:4165062/1‑93 (MQ=255)
                                                                                    gCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATg     >  1:4165063/1‑93 (MQ=255)
                                                                                     cccGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGc    <  2:4388499/93‑1 (MQ=255)
                                                                                      ccGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCg   <  2:1867537/93‑1 (MQ=255)
                                                                                       cgccCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGt  <  2:1611950/93‑1 (MQ=255)
                                                                                       cgcaCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGt  <  1:3302075/93‑1 (MQ=255)
                                                                                          |                                                                                         
ATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGT  >  NC_000913/1786978‑1787157

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