Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,428,098 G→A G126G (GGC→GGT aroE ← dehydroshikimate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,428,0980GA99.5% 726.4 / ‑5.6 206G126G (GGC→GGTaroEdehydroshikimate reductase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (83/122);  total (83/123)
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.81e-01

Only 100 of 213 total aligned reads displayed.

TTCCGCGCGGGATACCGTCCGATTAGTTATTGTCACCGCACAGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCGCCGATAAGCAGAATACGTAAACCAGGGCGGATAAAAGACAGACGTTCCAGATCGCTTAACAAGCCTACACCATCGGTATTGTCACCCA  >  NC_000913/3428008‑3428186
                                                                                          |                                                                                        
ttCCGCGCGGGATACCGTCCGATTAGTTATTGTCACCGCACAGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCAcc                                                                                        >  1:3974546/1‑93 (MQ=255)
ttCCGCGCGGGATACCGTCCGATTAGTTATTGTCACCGCACAGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCAcc                                                                                        >  1:2294901/1‑93 (MQ=255)
  ccGCGCGGGATACCGTCCGATTAGTTATTGTCACCGCACAGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGa                                                                                      <  2:4614762/93‑1 (MQ=255)
   cgcgcgGGATACCGTCCGATTAGTTATTGTCACCGCACAGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGAt                                                                                     <  1:3031642/93‑1 (MQ=255)
     cgcgGGATACCGTCCGATTAGTTATTGTCACCGCACAGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATaa                                                                                   <  2:1177128/93‑1 (MQ=255)
     cgcgGGATACCGTCCGATTAGTTATTGTCACCGCACAGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATaa                                                                                   <  2:4480656/93‑1 (MQ=255)
     cgcgGGATACCGTCCGATTAGTTATTGTCACCGCACAGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATaa                                                                                   <  1:4328649/93‑1 (MQ=255)
      gcgGGATACCGTCCGATTAGTTATTGTCACCGCACAGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAg                                                                                  <  1:4665366/93‑1 (MQ=255)
       cgGGATACCGTCCGATTAGTTATTGTCACCGCACAGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGc                                                                                 >  1:1417348/1‑93 (MQ=255)
        gggATACCGTCCGATTAGTTATTGTCACCGCACAGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCa                                                                                <  1:4164780/93‑1 (MQ=255)
        gggATACCGTCCGATTAGTTATTGTCACCGCACAGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCa                                                                                <  1:1521026/93‑1 (MQ=255)
        gggATACCGTCCGATTAGTTATTGTCACCGCACAGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCa                                                                                <  2:4331666/93‑1 (MQ=255)
        gggATACCGTCCGATTAGTTATTGTCACCGCACAGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCa                                                                                <  1:2908693/93‑1 (MQ=255)
         ggATACCGTCCGATTAGTTATTGTCACCGCACAGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAg                                                                               >  1:3927010/1‑93 (MQ=255)
          gATACCGTCCGATTAGTTATTGTCACCGCACAGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGa                                                                              >  1:3557505/1‑93 (MQ=255)
           aTACCGTCCGATTAGTTATTGTCACCGCACAGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGaa                                                                             <  1:2694827/93‑1 (MQ=255)
            tACCGTCCGATTAGTTATTGTCACCGCACAGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAAt                                                                            >  2:4151405/1‑93 (MQ=255)
             aCCGTCCGATTAGTTATTGTCACCGCACAGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATa                                                                           <  1:2655722/93‑1 (MQ=255)
             aCCGTCCGATTAGTTATTGTCACCGCACAGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATa                                                                           <  2:2996961/93‑1 (MQ=255)
              ccGTCCGATTAGTTATTGTCACCGCACAGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATAc                                                                          >  2:640503/1‑93 (MQ=255)
              ccGTCCGATTAGTTATTGTCACCGCACAGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATAc                                                                          <  2:1692627/93‑1 (MQ=255)
               cGTCCGATTAGTTATTGTCACCGCACAGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACg                                                                         >  2:4176854/1‑93 (MQ=255)
                gTCCGATTAGTTATTGTCACCGCACAGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAACCAGAATACga                                                                        >  2:1096448/1‑92 (MQ=255)
                  ccGATTAGTTATTGTCACCGCACAGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTaa                                                                      <  2:187639/93‑1 (MQ=255)
                   cgcttagttagTGTCACCGTACAGTCGAGGTAAATGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTaaa                                                                     <  1:2073039/90‑1 (MQ=255)
                   cGATTAGTTATTGTCACCGCACAGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTaaa                                                                     <  2:52992/93‑1 (MQ=255)
                    gATTAGTTATTGTCACCGCACAGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAAc                                                                    >  1:2929903/1‑93 (MQ=255)
                    gATTAGTTATTGTCACCGCACAGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAAc                                                                    <  1:313893/93‑1 (MQ=255)
                    gATTAGTTATTGTCACCGCACAGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAAc                                                                    <  1:669525/93‑1 (MQ=255)
                     aTTAGTTATTGTCACCGCACAGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAAcc                                                                   >  2:3716264/1‑93 (MQ=255)
                      ttAGTTATTGTCACCGCACAGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCa                                                                  <  2:3667260/93‑1 (MQ=255)
                      ttAGTTATTGTCACCGCACAGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCa                                                                  >  1:2086546/1‑93 (MQ=255)
                       tAGTTATTGTCACCGCACAGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTCAACCAg                                                                 >  2:585074/1‑93 (MQ=255)
                       tAGTTATTGTCACCGCACAGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAg                                                                 <  2:344803/93‑1 (MQ=255)
                       tAGTTATTGTCACCGCACAGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAg                                                                 <  1:365597/93‑1 (MQ=255)
                        aGTTATTGTCACCGCACAGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAgg                                                                <  2:1991969/93‑1 (MQ=255)
                         gTTATTGTCACCGCACAGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAggg                                                               >  2:4198308/1‑93 (MQ=255)
                         gTTATTGTCACCGCACAGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAggg                                                               >  1:4220399/1‑93 (MQ=255)
                          ttATTGTCACCGCACAGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGc                                                              >  1:1193308/1‑93 (MQ=255)
                            aTTGTCACCGCACAGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGCTGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCgg                                                            >  2:2037002/1‑93 (MQ=255)
                              tGTCACCGCACAGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGAt                                                          >  1:1805383/1‑93 (MQ=255)
                                tCACCGCACAGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATaa                                                        <  2:3507589/93‑1 (MQ=255)
                                 cACCGCACAGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATaaa                                                       <  1:2261077/93‑1 (MQ=255)
                                  aCCGCACAGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATaaaa                                                      <  1:347858/93‑1 (MQ=255)
                                  aCCGCACAGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATaaaa                                                      <  2:2081916/93‑1 (MQ=255)
                                  aCCGCACAGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATaaaa                                                      <  1:748810/93‑1 (MQ=255)
                                     gCACAGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAagac                                                   >  2:2247135/1‑93 (MQ=255)
                                     gCACAGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAagac                                                   <  2:2043458/93‑1 (MQ=255)
                                       acaGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAagacag                                                 <  1:517248/93‑1 (MQ=255)
                                       acaGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAagacag                                                 >  2:150669/1‑93 (MQ=255)
                                        caGTCGAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAagacaga                                                <  1:2650529/93‑1 (MQ=255)
                                        caGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAagacaga                                                >  1:4037378/1‑93 (MQ=255)
                                        caGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAagacaga                                                <  2:3306965/93‑1 (MQ=255)
                                        caGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAagacaga                                                >  2:2960403/1‑93 (MQ=255)
                                        caGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAagacaga                                                <  2:2126678/93‑1 (MQ=255)
                                        caGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAagacaga                                                <  1:4645342/93‑1 (MQ=255)
                                         aGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAagacagac                                               <  1:2150399/93‑1 (MQ=255)
                                          gTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCCCCAGCACCGATAAGCAGAATACGTAAACAAGGGCGGATAAAAGACAGACg                                              >  2:1108987/1‑93 (MQ=255)
                                           tCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAAGACAGACGt                                             >  1:1888203/1‑93 (MQ=255)
                                           tCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAACCAGGGCGGATAAAAGACAGACGtt                                            <  1:376644/93‑1 (MQ=255)
                                            ccAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAAGACAGACGtt                                            <  1:4489220/93‑1 (MQ=255)
                                            ccAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAAGACAGACGtt                                            <  1:2467840/93‑1 (MQ=255)
                                               gggAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAAGACAGACGTTCCa                                         <  1:2626681/93‑1 (MQ=255)
                                               gggAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAAGACAGACGTTCCa                                         <  2:2696831/93‑1 (MQ=255)
                                                 gAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAAGACAGACGTTCCAGa                                       <  1:4273179/93‑1 (MQ=255)
                                                  aaaGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAAGACAGACGTTCCAGAt                                      >  1:1101509/1‑93 (MQ=255)
                                                  aaaGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAAGACAGACGTTCCAGAt                                      >  2:2559925/1‑93 (MQ=255)
                                                   aaGGAGTGGCAGAAGTACGCCGCGACATGCTCCAACAGCACCGATAAGAAGAATACGTAAACCAGGGCGGATAAAAGACAGACGTTCCAGATc                                     >  2:4407026/1‑93 (MQ=255)
                                                        gTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAAGACAGACGTTCCAGATCGCTTa                                <  1:1610049/93‑1 (MQ=255)
                                                        gTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAAGACAGACGTTCCAGATCGCTTa                                >  2:3023025/1‑93 (MQ=255)
                                                        gTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAAGACAGACGTTCCAGATCGCTTa                                >  1:4371974/1‑93 (MQ=255)
                                                        gTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGCATAAAAGACAGACGTTCCAGATCGCTTa                                >  1:2008273/1‑93 (MQ=255)
                                                         tGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAAGACAGACGTTCCAGATCGCTTaa                               <  1:122768/93‑1 (MQ=255)
                                                         tGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAAGACAGACGTTCCAGATCGCTTaa                               >  2:3210986/1‑93 (MQ=255)
                                                         tGGCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAAGACAGACGTTCCAGATCGCTTaa                               <  1:4593748/93‑1 (MQ=255)
                                                           gCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAAGACAGACGTTCCAGATCGCTTAACa                             >  1:2621748/1‑93 (MQ=255)
                                                           gCAGTAGTACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAAGACAGACGTTCCAGATCGCTTAACa                             >  2:4675115/1‑93 (MQ=255)
                                                             agtagtACGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAAGACAGACGTTCCAGATCGCTTAACAAg                           >  2:49340/1‑93 (MQ=255)
                                                                 gtaCGCCGCGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAAGACAGACGTTCCAGATCGCTTAACAAGCCTa                       <  1:3216651/93‑1 (MQ=255)
                                                                    cgccgcGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAAGACAGACGTTCCAGATCGCTTAACAAGCCTacac                    <  2:1048534/93‑1 (MQ=255)
                                                                      ccgcGAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAAGACAGACGTTCCAGATCGCTTAACAAGCCTACACCa                  <  1:2195164/93‑1 (MQ=255)
                                                                       cgcgAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAAGACAGACGTTCCAGATCGCTTAACAAGCCTACACCAt                 <  1:2343627/93‑1 (MQ=255)
                                                                       cgcgAGATGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAAGACAGACGTTCCAGATCGCTTAACAAGCCTACACCAt                 <  1:4516992/93‑1 (MQ=255)
                                                                          gagaTGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAAGACAGACGTTCCAGATCGCTTAACAAGCCTACACCATCgg              >  2:3555073/1‑93 (MQ=255)
                                                                          gagaTGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAAGACAGACGTTCCAGATCGCTTAACAAGCCTACACCATCgg              >  2:3989112/1‑93 (MQ=255)
                                                                          gagaTGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAAGACAGACGTTCCAGATCGCTTAACAAGCCTACACCATCgg              >  1:1299753/1‑93 (MQ=255)
                                                                            gaTGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAAGACAGACGTTCCAGATCGCTTAACAAGCCTACACCATCGGTa            >  2:4297451/1‑93 (MQ=255)
                                                                            gaTGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAAGACAGACGTTCCAGATCGCTTAACAAGCCTACACCATCGGTa            >  1:3870121/1‑93 (MQ=255)
                                                                             aTGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAAGACAGACGTTCCAGATCGCTTAACAAGCCTACACCATCGGTAt           <  2:205159/93‑1 (MQ=255)
                                                                             aTGCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAAGACAGACGTTCCAGATCGCTTAACAAGCCTACACCATCGGTAt           <  2:1528718/93‑1 (MQ=255)
                                                                               gCTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAAGACAGACGTTCCAGATCGCTTAACAAGCCTACACCATCGGTATTg         <  1:1858124/93‑1 (MQ=255)
                                                                                cTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAAGACAGACGTTCCAGATCGCTTAACAAGCCTACACCATCGGTATTGt        <  1:335833/93‑1 (MQ=255)
                                                                                cTCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAAGACAGACGTTCCAGATCGCTTAACAAGCCTACACCATCGGTATTGt        >  1:2482302/1‑93 (MQ=255)
                                                                                 tCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAAGACAGACGTTCCAGATCGCTTAACAAGCCTACACCATCGGTATTGTc       >  2:915015/1‑93 (MQ=255)
                                                                                 tCCACCAGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAAGACAGACGTTCCAGATCGCTTAACAAGCCTACACCATCGGTATTGTc       >  1:3295554/1‑93 (MQ=255)
                                                                                 cccaccaGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAAGACAGACGTTCCAGATCGCTTAACAAGCCTACACCATCGGTATTGTc       <  2:2818792/92‑1 (MQ=255)
                                                                                  ccaccaGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAAGACAGACGTTCCAGATCGCTTAACAAGCCTACACCATCGGTATTGTCa      <  1:2738725/93‑1 (MQ=255)
                                                                                  ccaccaGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAAGACAGACGTCCCAGATCGCTTAACAAGCCTACACCATCGGTATTGTCa      <  2:1287171/93‑1 (MQ=255)
                                                                                     ccaGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAAGACAGACGTTCCAGATCGCTTAACAAGCCTACACCATCGGTATTGTCAccc   <  2:2464674/93‑1 (MQ=255)
                                                                                      caGCACCGATAAGCAGAATACGTAAACCAGGGCGGATAAAAGACAGACGTTCCAGATCGCTTAACAAGCCTACACCATCGGTATTGTCACCca  <  2:4113340/93‑1 (MQ=255)
                                                                                          |                                                                                        
TTCCGCGCGGGATACCGTCCGATTAGTTATTGTCACCGCACAGTCCAGGGAAAGGAGTGGCAGTAGTACGCCGCGAGATGCTCCACCAGCGCCGATAAGCAGAATACGTAAACCAGGGCGGATAAAAGACAGACGTTCCAGATCGCTTAACAAGCCTACACCATCGGTATTGTCACCCA  >  NC_000913/3428008‑3428186

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