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% 389.8 / NA 114A204T (GCC→ACC) aroH3‑deoxy‑D‑arabinoheptulosonate‑7‑phosphate synthase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (54/59);  minor base C (0/1);  total (54/60)
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 132 total aligned reads displayed.

ATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGT  >  NC_000913/1786978‑1787157
                                                                                          |                                                                                         
aTGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCacc                                                                                         <  2:2543598/93‑1 (MQ=255)
aTGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCacc                                                                                         >  1:1466363/1‑93 (MQ=255)
aTGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCacc                                                                                         >  1:346434/1‑93 (MQ=255)
   gCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCaccagc                                                                                      >  2:4475644/1‑93 (MQ=255)
    cTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCaccagcc                                                                                     >  2:410545/1‑93 (MQ=255)
    cTTCGGCACGCTCCTGTCCGGTAGTGTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCaccagcc                                                                                     <  1:2200706/93‑1 (MQ=255)
     ttCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCaccagcca                                                                                    >  2:1996898/1‑93 (MQ=255)
       cGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCata                                                                                  <  1:102333/93‑1 (MQ=255)
       cGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCata                                                                                  <  1:905686/93‑1 (MQ=255)
         gCACTCTCCTGTCCGTTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTAGCCGCGCAGCCCGCACCAGCCATATg                                                                                <  2:706859/93‑1 (MQ=255)
          ccctctCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGt                                                                               <  1:813488/91‑1 (MQ=255)
          cACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGt                                                                               <  1:2153668/93‑1 (MQ=255)
           aCTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGtt                                                                              >  1:458956/1‑93 (MQ=255)
           aCTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGtt                                                                              >  2:4370081/1‑93 (MQ=255)
            ctctCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTc                                                                             <  2:2023226/93‑1 (MQ=255)
            ctctCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTc                                                                             >  1:2157819/1‑93 (MQ=255)
            ctctCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTc                                                                             <  2:4623624/93‑1 (MQ=255)
             tctcCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTcc                                                                            >  2:4675678/1‑93 (MQ=255)
             tctcCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTcc                                                                            >  1:3538508/1‑93 (MQ=255)
             tctcCTGTCCGGTAGGTTTTAAAAATAGGACCGATGGCAATACGGGGAATGCTGTGGATGCTATTCGCGCAGTCTGTACCAGCCATATGTTcc                                                                            >  2:2949407/1‑93 (MQ=255)
               tcCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCctc                                                                          >  1:3926658/1‑93 (MQ=255)
               tcCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCctc                                                                          <  1:3706904/93‑1 (MQ=255)
                 cTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCtctc                                                                        >  2:247994/1‑93 (MQ=255)
                 cTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCtctc                                                                        >  1:721421/1‑93 (MQ=255)
                    tCCGGTAGGTTTTAAAAATGGTCCCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGcc                                                                     <  1:339559/93‑1 (MQ=255)
                      cGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATTCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAg                                                                   >  2:3584854/1‑93 (MQ=255)
                       ggTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGa                                                                  <  1:1655836/93‑1 (MQ=255)
                       ggTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGa                                                                  >  2:1907226/1‑93 (MQ=255)
                       ggTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGa                                                                  >  2:2853428/1‑93 (MQ=255)
                        gTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGAc                                                                 <  1:2862358/93‑1 (MQ=255)
                        gTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGAc                                                                 >  2:2099621/1‑93 (MQ=255)
                         tAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACa                                                                <  1:582725/93‑1 (MQ=255)
                          aGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACaa                                                               >  2:1993130/1‑93 (MQ=255)
                            gTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCTGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACaaaa                                                             >  1:791974/1‑93 (MQ=255)
                            gTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACaaaa                                                             >  1:4660308/1‑93 (MQ=255)
                            gTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACaaaa                                                             >  2:2123387/1‑93 (MQ=255)
                             ttttAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACaaaaa                                                            >  1:1308955/1‑93 (MQ=255)
                              tttAAAAATGGTCCCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAAt                                                           <  1:4425672/93‑1 (MQ=255)
                               ttAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATg                                                          >  2:1871951/1‑93 (MQ=255)
                                tAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATgg                                                         <  1:3025857/93‑1 (MQ=255)
                                tAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATgg                                                         <  1:2815216/93‑1 (MQ=255)
                                tAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATgg                                                         <  1:485549/93‑1 (MQ=255)
                                tAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATgg                                                         <  1:453888/93‑1 (MQ=255)
                                tAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATgg                                                         >  2:4332901/1‑93 (MQ=255)
                                  aaaaGGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTc                                                       <  2:3955921/93‑1 (MQ=255)
                                   aaaTGGTCCCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCa                                                      <  1:1194657/93‑1 (MQ=255)
                                    aaTGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAg                                                     >  2:1163044/1‑93 (MQ=255)
                                    aaTGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAg                                                     >  2:3916461/1‑93 (MQ=255)
                                     aTGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGa                                                    <  1:2492634/93‑1 (MQ=255)
                                      tGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGAt                                                   <  1:1730479/93‑1 (MQ=255)
                                      tGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGAt                                                   <  1:4440371/93‑1 (MQ=255)
                                       ggTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATg                                                  >  1:249318/1‑93 (MQ=255)
                                        gTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGa                                                 <  2:1132179/93‑1 (MQ=255)
                                          aCCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGagc                                               >  2:1222910/1‑91 (MQ=255)
                                           ccGATGGCAATACGCGGATGGCTGTGGATGCTATCCGCTCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGCCCa                                              <  2:1540918/93‑1 (MQ=255)
                                            cGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCAt                                             <  1:207245/93‑1 (MQ=255)
                                             gagggCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATc                                            <  2:4367179/90‑1 (MQ=255)
                                             gATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATc                                            <  2:4381473/93‑1 (MQ=255)
                                               tGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTa                                          >  2:3461896/1‑93 (MQ=255)
                                                  cAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCa                                       <  2:2848097/93‑1 (MQ=255)
                                                    aTACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGa                                     <  2:3294633/93‑1 (MQ=255)
                                                       cgcgGATTGCTGTGGATGCAATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCa                                  <  1:4119482/93‑1 (MQ=255)
                                                            aTTGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGc                             >  1:199597/1‑93 (MQ=255)
                                                              tGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCaa                           <  2:4481123/93‑1 (MQ=255)
                                                              tGCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCaa                           <  1:1843732/93‑1 (MQ=255)
                                                               gCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAAc                          <  2:3708219/93‑1 (MQ=255)
                                                               gCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAAc                          <  2:306465/93‑1 (MQ=255)
                                                               gCTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAAc                          >  1:1719357/1‑93 (MQ=255)
                                                                cTGTGGATGCTATCCGCGCAGCCCGCCCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAAcc                         <  1:3559147/93‑1 (MQ=255)
                                                                cTGTGGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAAcc                         >  1:3877761/1‑93 (MQ=255)
                                                                 tgtgGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAAccc                        >  2:753911/1‑93 (MQ=255)
                                                                 tgtgGATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAAccc                        >  2:854717/1‑93 (MQ=255)
                                                                    ggATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTa                     <  1:469733/93‑1 (MQ=255)
                                                                    ggATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTa                     <  1:2686770/93‑1 (MQ=255)
                                                                    ggATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTa                     <  2:41000/93‑1 (MQ=255)
                                                                     gATGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTAt                    >  1:3086860/1‑93 (MQ=255)
                                                                     gATGCTATCCGCGCAGCCCAGCCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTAt                    <  1:2311401/93‑1 (MQ=255)
                                                                      aTGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATg                   >  1:4303682/1‑93 (MQ=255)
                                                                       tGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATgg                  >  1:3384891/1‑93 (MQ=255)
                                                                       tGCTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATgg                  >  1:1522001/1‑93 (MQ=255)
                                                                         cTATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGcc                >  2:4019875/1‑93 (MQ=255)
                                                                          tATCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCca               >  1:3018076/1‑93 (MQ=255)
                                                                           aTCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCcac              >  1:4336898/1‑93 (MQ=255)
                                                                            tCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCcaca             >  2:2325997/1‑93 (MQ=255)
                                                                            tCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCcaca             <  1:4388358/93‑1 (MQ=255)
                                                                            tCCGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCcaca             >  1:745141/1‑93 (MQ=255)
                                                                             ccGCGCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACat            <  1:431544/93‑1 (MQ=255)
                                                                              cgcgCAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACatt           >  1:3655169/1‑93 (MQ=255)
                                                                               ccccAGCCCGCCCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACatta          <  1:3264203/90‑1 (MQ=255)
                                                                                cgcAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACattat         <  2:2392049/93‑1 (MQ=255)
                                                                                 gcAGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACattatt        >  1:3777962/1‑93 (MQ=255)
                                                                                  caccccGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCGACAttatta       <  2:1848912/90‑1 (MQ=255)
                                                                                  cAGCCCGCCCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACAttatta       <  1:3235559/93‑1 (MQ=255)
                                                                                  cAGCCCCCCCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACAttatta       <  2:1350253/93‑1 (MQ=255)
                                                                                  cAGCCCCCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACAttatta       <  1:1729511/93‑1 (MQ=255)
                                                                                   aGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATtattat      <  2:346434/93‑1 (MQ=255)
                                                                                   aGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATtattat      >  2:4625213/1‑93 (MQ=255)
                                                                                   aGCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATtattat      <  1:1178097/93‑1 (MQ=255)
                                                                                    gCCCGCACCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATg     <  2:4160330/93‑1 (MQ=255)
                                                                                       cgcaCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGt  <  2:3680220/93‑1 (MQ=255)
                                                                                          |                                                                                         
ATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGT  >  NC_000913/1786978‑1787157

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