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.5% 682.9 / NA 197A204T (GCC→ACC) aroH3‑deoxy‑D‑arabinoheptulosonate‑7‑phosphate synthase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (102/93);  total (102/95)
Fisher's exact test for biased strand distribution p-value = 4.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 225 total aligned reads displayed.

ATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGT  >  NC_000913/1786978‑1787157
                                                                                          |                                                                                         
aTGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCacc                                                                                         <  2:6402318/93‑1 (MQ=255)
     ttCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCGGCGCAGCCCGCaccagcca                                                                                    <  1:6543459/93‑1 (MQ=255)
     ttCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCaccagcca                                                                                    >  2:6002855/1‑93 (MQ=255)
     ttCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCaccagcca                                                                                    >  1:1947834/1‑93 (MQ=255)
       cGGCACTCTCCTGTCCGTTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCata                                                                                  <  2:3261694/93‑1 (MQ=255)
       cGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCata                                                                                  <  2:100560/93‑1 (MQ=255)
       cGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCata                                                                                  <  2:27627/93‑1 (MQ=255)
       cGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCata                                                                                  <  2:6049916/93‑1 (MQ=255)
        ggCGCTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCatat                                                                                 <  2:5330256/93‑1 (MQ=255)
         gCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATg                                                                                <  1:3179658/93‑1 (MQ=255)
         gCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATg                                                                                <  2:2579019/93‑1 (MQ=255)
          cACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGt                                                                               <  1:3606281/93‑1 (MQ=255)
           aCTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGtt                                                                              >  2:2437068/1‑93 (MQ=255)
           aCTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGtt                                                                              >  1:5832392/1‑93 (MQ=255)
             tctcCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTcc                                                                            <  1:2630955/93‑1 (MQ=255)
             tctcCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTcc                                                                            <  2:3200073/93‑1 (MQ=255)
                ccTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCCCCAGCCATATGTTCctct                                                                         <  1:5059206/93‑1 (MQ=255)
                   gTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGc                                                                      >  2:6307508/1‑93 (MQ=255)
                    tCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGcc                                                                     >  2:3366272/1‑93 (MQ=255)
                    tCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGcc                                                                     >  2:4534023/1‑93 (MQ=255)
                     ccGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCa                                                                    >  1:5283956/1‑93 (MQ=255)
                       ggTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGa                                                                  >  1:6510223/1‑93 (MQ=255)
                        gTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGAc                                                                 >  2:4277427/1‑93 (MQ=255)
                        gTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGAc                                                                 >  2:4672765/1‑93 (MQ=255)
                         tAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACa                                                                <  1:1471069/93‑1 (MQ=255)
                         tAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACa                                                                >  2:2092064/1‑93 (MQ=255)
                          aGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACaa                                                               >  1:4943440/1‑93 (MQ=255)
                          aGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGCTCCTCTCGCCAGACaa                                                               >  1:2599754/1‑93 (MQ=255)
                           ggTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACaaa                                                              <  2:5535824/93‑1 (MQ=255)
                            gTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACaaaa                                                             >  2:3005438/1‑93 (MQ=255)
                            gTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACaaaa                                                             <  2:3910199/93‑1 (MQ=255)
                            gTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACaaaa                                                             >  2:6523933/1‑93 (MQ=255)
                            gTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACaaaa                                                             <  1:4417153/93‑1 (MQ=255)
                            gTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACaaaa                                                             >  2:3466170/1‑93 (MQ=255)
                             ttttAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACaaaaa                                                            <  1:818211/93‑1 (MQ=255)
                             ttttAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACaaaaa                                                            >  2:147671/1‑93 (MQ=255)
                             ttttAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACaaaaa                                                            <  1:2531934/93‑1 (MQ=255)
                              tttAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAAt                                                           <  2:443031/93‑1 (MQ=255)
                               ttAAAAATGGTACCGATGGCAATACTCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTTGCCAGACAAAAATg                                                          >  1:6731142/1‑93 (MQ=255)
                               ttAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATg                                                          >  1:4825861/1‑93 (MQ=255)
                                 aaaaaTGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGt                                                        >  1:2449515/1‑93 (MQ=255)
                                  aaaaTGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTc                                                       >  2:6632061/1‑93 (MQ=255)
                                   aaaTGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATTGTCa                                                      <  2:4029679/93‑1 (MQ=255)
                                   aaaTGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCa                                                      >  2:5830905/1‑93 (MQ=255)
                                   aaaTGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCa                                                      >  2:5098943/1‑93 (MQ=255)
                                   aaaTGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCa                                                      >  1:1519447/1‑93 (MQ=255)
                                     aTGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGa                                                    >  2:1778487/1‑93 (MQ=255)
                                     aTGGTACCGATGGCAATACGCGGATTGCGGTGGATGCTACCCGGGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGa                                                    <  2:4963107/93‑1 (MQ=255)
                                       ggTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATg                                                  >  2:2732926/1‑93 (MQ=255)
                                       ggTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATg                                                  >  1:5707626/1‑93 (MQ=255)
                                        gTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGa                                                 <  2:1336726/93‑1 (MQ=255)
                                            cGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATAGTCAGATGACCAt                                             >  1:1074360/1‑93 (MQ=255)
                                             gATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATc                                            <  1:3525285/93‑1 (MQ=255)
                                              aTGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCCCCAGACAAAAATGGTCAGAGGACCATCt                                           >  1:6381293/1‑93 (MQ=255)
                                               tGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTa                                          >  1:6163838/1‑93 (MQ=255)
                                                 gCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGTTGACCATCTATc                                        <  2:3591421/93‑1 (MQ=255)
                                                   aaTACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAg                                      >  2:841168/1‑93 (MQ=255)
                                                    agaCGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGa                                     >  1:4203593/3‑93 (MQ=255)
                                                    aTACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGa                                     >  1:3065379/1‑93 (MQ=255)
                                                     tACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGAc                                    >  1:935607/1‑93 (MQ=255)
                                                     tACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGAc                                    >  1:1675990/1‑93 (MQ=255)
                                                      aCGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGAcc                                   >  2:2881525/1‑93 (MQ=255)
                                                       cgcgGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCa                                  >  1:3686106/1‑93 (MQ=255)
                                                       cgcgGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCa                                  <  1:2217711/93‑1 (MQ=255)
                                                       cgcgGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCa                                  <  1:4920905/93‑1 (MQ=255)
                                                        gcgGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAg                                 <  2:6218866/93‑1 (MQ=255)
                                                        ccgGATTGCTGTGGATGCTATCCGCCCAGCCCGCCCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAg                                 <  1:5592615/92‑1 (MQ=255)
                                                         cgGATTGCTGTGGATGCTACCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGACGACCATCTATCAGACCAGc                                <  1:3278968/93‑1 (MQ=255)
                                                          ggATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCg                               <  2:3121578/93‑1 (MQ=255)
                                                            aTTGCTGTGGATCCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGc                             <  1:6267038/93‑1 (MQ=255)
                                                                cTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAAcc                         <  1:5441752/93‑1 (MQ=255)
                                                                 tgtgGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAAccc                        >  2:4742832/1‑93 (MQ=255)
                                                                 tgtgGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAAccc                        >  2:1138268/1‑93 (MQ=255)
                                                                     gATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTAt                    <  2:6046027/93‑1 (MQ=255)
                                                                     gATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTAt                    <  1:1135118/93‑1 (MQ=255)
                                                                      aTGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATg                   <  2:3980657/93‑1 (MQ=255)
                                                                        gCTACCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGc                 <  2:861644/93‑1 (MQ=255)
                                                                         cTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGcc                <  1:202725/93‑1 (MQ=255)
                                                                           acccGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCcac              <  1:3438300/91‑1 (MQ=255)
                                                                            tCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCcaca             >  1:467873/1‑93 (MQ=255)
                                                                             ccGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACat            <  2:407645/93‑1 (MQ=255)
                                                                             ccGCGCACCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACat            <  2:2833034/93‑1 (MQ=255)
                                                                               gcgcAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACatta          >  1:6637132/1‑93 (MQ=255)
                                                                               gcgcAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACatta          >  2:6267018/1‑93 (MQ=255)
                                                                               gcgcAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACatta          >  1:2869219/1‑93 (MQ=255)
                                                                                cgcAGCCCTCACCAGCCATATTTTCCTTTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGGAACCAGTATGGCCACattat         <  1:116273/93‑1 (MQ=255)
                                                                                cgcAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACattat         >  2:230900/1‑93 (MQ=255)
                                                                                cgcAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACattat         <  2:5013054/93‑1 (MQ=255)
                                                                                 gcAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACattatt        <  1:1236539/93‑1 (MQ=255)
                                                                                 ccAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACattatt        <  2:2489488/92‑1 (MQ=255)
                                                                                  cAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACAttatta       <  2:4572921/93‑1 (MQ=255)
                                                                                  cAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACAttatta       <  1:3421925/93‑1 (MQ=255)
                                                                                   aGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATtattat      >  2:4883545/1‑93 (MQ=255)
                                                                                   aGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATtattat      <  1:5372023/93‑1 (MQ=255)
                                                                                   aGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATtattat      >  2:1933549/1‑93 (MQ=255)
                                                                                     cccGCGCCAGTCATAAGTTCCTCTCGCCAGCCAAAAATGGTCAGATGACCATCTATAAGACCAGCGGAAACCAGTATGGCCACATTATTATGc    <  2:2612694/93‑1 (MQ=255)
                                                                                     cccGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGc    >  1:4290907/1‑93 (MQ=255)
                                                                                      ccGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCg   >  1:2101253/1‑93 (MQ=255)
                                                                                      ccGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCg   <  1:1355497/93‑1 (MQ=255)
                                                                                      ccGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCg   <  2:3792003/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: