Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,773,032 G→A G108D (GGC→GAC)  aroD → 3‑dehydroquinate dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,773,0320GA97.3% 674.9 / ‑0.6 187G108D (GGC→GAC) aroD3‑dehydroquinate dehydratase
Reads supporting (aligned to +/- strand):  ref base G (3/2);  new base A (89/93);  total (92/95)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 189 total aligned reads displayed.

TGTTTACCTTCCGCAGTGCCAAAGAAGGCGGCGAGCAGGCGATTTCCACCGAGGCTTATATTGCACTCAATCGTGCAGCCATCGACAGCGGCCTGGTTGATATGATCGATCTGGAGTTATTTACCGGTGATGATCAGGTTAAAGAAACCGTCGCCTACGCCCACGCGCATGATGTGAAAG  >  NC_000913/1772942‑1773121
                                                                                          |                                                                                         
 gTTTACCTTCCGCAGTGCCAAAGAAGGCGGCGAGCAGGCGATTTCCACCGAGGCTTATATTGCACTCAATCGTGCAGCCATCGACAGCGACCt                                                                                        >  2:1828127/1‑93 (MQ=255)
 gTTTACCTTCCGCAGTGCCAAAGAAGGCGGCGAGCAGGCGATTTCCACCGAGGCTTATATTGCACTCAATCGTGCAGCCATCGACAGCGACCt                                                                                        >  2:2712388/1‑93 (MQ=255)
 gTTTACCTTCCGCAGTGCCAAAGAAGGCGGCGAGCAGGCGATTTCCACCGAGGCTTATATTGCACTCAATCGTGCAGCCATCGACAGCGACCt                                                                                        >  2:4289972/1‑93 (MQ=255)
  tttACCTTCCGCAGTGCCAAAGAAGGCGGCGAGCAGGCGATTTCCACCGAGGCTTATATTGCACTCAATCGTGCGGCCATCGACAGCGACCTg                                                                                       >  2:215412/1‑93 (MQ=255)
  tttACCTTCCGCAGTGCCAAAGAAGGCGGCGAGCAGGCGATTTCCACCGAGGCTTATATTGCACTCAATCGTGCAGCCATCGACAGCGACCTg                                                                                       >  1:2722804/1‑93 (MQ=255)
    tACCTTCCGCAGTGCCAAAGAAGGCGGGGAGCAGGCGATTTCCACCGAGGCTTATATTGCACTCAATCGTGCAGCCATCGACAGCGACCTGGt                                                                                     <  1:3429/93‑1 (MQ=255)
    tACCTTCCGCAGTGCCAAAGAAGGCGGCGAGCAGGCGATTTCCACCGAGGCTTATATTGCACTCAATCGTGCAGCCATCGACAGCGACCTGGt                                                                                     <  2:2389995/93‑1 (MQ=255)
    tACCTTCCGCAGTGCCAAAGAAGGCGGCGAGCAGGCGATTTCCACCGAGGCTTATATTGCACTCAATCGTGCAGCCATCGACAGCGACCTGGt                                                                                     >  2:4654063/1‑93 (MQ=255)
     tccTTCCGCAGTGCCAAAGAAGGCGGCGAGCAGGCGATTTCCACCGAGGCTTATATTGCACTCAATCGTGCAGCCATCGACAGCGACCTGGtt                                                                                    <  1:3022001/92‑1 (MQ=255)
     tccTTCCGCAGTGCCAAAGAAGGCGGCGAGCAGGCGATTTCCACCGAGGCTTATATTGCACTCAATCGTGCAGCCATCGACAGCGACCTGGtt                                                                                    <  2:1650000/92‑1 (MQ=255)
      ccTTCCGCAGTGCCAAAGAAGGCGGCGAGCAGGCGATTTCCACCGAGGCTTATATTGCACTCAATCGTGCAGCCATCGACAGCGACCTGGTTg                                                                                   >  1:2169076/1‑93 (MQ=255)
       cTTCCGCAGTGCCAAAGAAGGCGGCGAGCAGGCGATTTCCACCGAGGCTTATATTGCACTCAATCGTGCAGCCATCGACAGCGACCTGGTTGa                                                                                  >  2:3427579/1‑93 (MQ=255)
         tCCGCAGTGCCAAAGAAGGCGGCGAGCAGGCGATTTCCACCGAGGCTTATATTGCACTCAATCGTGCAGCCATCGACAGCGACCTGGTTGata                                                                                <  1:1765177/93‑1 (MQ=255)
         tCCGCAGTGCCAAAGAAGGCGGCGAGCAGGCGATTTCCACCGAGGCTTATATTGCACTCAATCGTGCACGCATCGACAGCGACCTGGTTGata                                                                                >  1:3482400/1‑93 (MQ=255)
          ccGCAGTGCCAAAGAAGGCGGCGAGCAGGCGATTTCCACCGAGGCTTATATTGCACTCAATCGTGCAGCCATCGACAGCGGCCTGGTTGatat                                                                               <  2:1142161/93‑1 (MQ=255)
           cGCAGTGCCAAAGAAGGCGGCGAGCAGGCGATTTCCACCGAGGCTTATATTGCACTCAATCGTGCAGCCATCGACAGCGACCTGGTTGATATg                                                                              >  1:5166754/1‑93 (MQ=255)
            gCAGTGCCAAAGAAGGCGGCGAGCAGGCGATTTCCACCGAGGCTTATATTGCACTCAATCGTGCAGCCATCGACAGCGACCTGGTTGATATga                                                                             <  2:3304197/93‑1 (MQ=255)
            gCAGTGCCAAAGAAGGCGGCGAGCAGGCGATTTCCACCGAGGCTTATATTGCACTCAATCGTGCAGCCATCGACAGCGACCTGGTTGATATga                                                                             >  2:2712396/1‑93 (MQ=255)
             cAGTGCCAAAGAAGGCGGCGAGCAGGCGATTTCCACCGAGGCTTATATTGCACTCAATCGTGCAGCCATCGACAGCGACCTGGTTGATATgat                                                                            <  2:5658447/93‑1 (MQ=255)
             cAGTGCCAAAGAAGGCGGCGAGCAGGCGATTTCCACCGAGGCTTATATTGCACTCAATCGTGCAGCCATCGACAGCGACCTGGTTGATATgat                                                                            <  2:1826164/93‑1 (MQ=255)
             cAGTGCCAAAGAAGGCGGCGAGCAGGCGATTTCCACCGAGGCTTATATTGCACTCAATCGTGCAGCCATCGACAGCGACCTGGTTGATATgat                                                                            <  2:5085966/93‑1 (MQ=255)
              aGTGCCAAAGAAGGCGGCGAGCAGGCGATTTCCACCGAGGCTTATATTGCACTCAATCGTGCAGCCATCGACAGCGACCTGGTTGATATgatc                                                                           <  1:3214482/93‑1 (MQ=255)
                  ccAAAGAAGGCGGCGAGCAGGCGATTTCCACCGAGGCTTATATTGCACTCAATCGTGCAGCCATCGACAGCGACCTGGTTGATATgatcgatc                                                                       <  2:2472666/93‑1 (MQ=255)
                  ccAAAGAAGGCGGCGAGCAGGCGATTTCCACCGAGGCTTATATTGCACTCAATCGTGCAGCCATCGACAGCGACCTGGTTGATATgatcgatc                                                                       <  1:5597993/93‑1 (MQ=255)
                   cAAAGAAGGCGGCGAGCAGGCGATTTCCACCGAGGCTTATATTGCACTCAATCGTGCAGCCATCGACAGCGGCCTGGTTGATATGATCGATCt                                                                      >  1:5369374/1‑93 (MQ=255)
                   cAAAGAAGGCGGCGAGCAGGCGATTTCCACCGAGGCTTATATTGCACTCAATCGTGCAGCCATCGACAGCGACCTGGTTGATATGATCGATCt                                                                      >  2:2045408/1‑93 (MQ=255)
                   cAAAGAAGGCGGCGAGCAGGCGATTTCCACCGAGGCTTATATTGCACTCAATCGTGCAGCCATCGACAGCGACCTGGGTGATATGATCGATCt                                                                      >  1:4230420/1‑93 (MQ=255)
                    aaaGAAGGCGGCGAGCAGGCGATTTCCACCGAGGCTTATATTGCACTCAATCGTGCAGCCATCGACAGCGACCTGGTTGATATGATCGATCTg                                                                     >  2:3559722/1‑93 (MQ=255)
                     aagaagGCGGCGAGCAGGCGATTTCCACCGAGGCTTATATTGCACTCAATCGTGCAGCCATCGACAGCGACCTGGTTGATATGATCGATCTgg                                                                    >  1:2684296/1‑93 (MQ=255)
                      agaagGCGGCGAGCAGGCGATTTCCACCGAGGCTTATATTGCACTCAATCGTGCAGCCATCGACAGCGACCTGGTTGATATGATCGATCTGGa                                                                   <  2:5758609/93‑1 (MQ=255)
                       gaagGCGGCGAGCAGGCGATTTCCACCGAGGCTTATATTGCACTCAATCGTGCAGCCATCGACAGCGACCTGGTTGATATGATCGATCTGGAg                                                                  <  2:4638003/93‑1 (MQ=255)
                        aagGCGGCGAGCAGGCGATTTCCACCGAGGCTTATATTGCACTCAATCGTGCAGCCATCGACAGCGACCTGGTTGATATGATCGATCTGGAGt                                                                 <  2:4444789/93‑1 (MQ=255)
                              gcgAGCAGGCGATTTCCACCGAGGCTTATATTGCACTCAATCGTGCAGCCATCGACAGCGACCTGGTTGATATGATCGATCTGGAGTTATTTa                                                           >  1:5124670/1‑93 (MQ=255)
                              gcgAGCAGGCGATTTCCACCGAGGCTTATATTGCACTCAATCGTGCAGCCATCGACAGCGACCTGGTTGATATGATCGATCTGGAGTTATTTa                                                           <  2:5313755/93‑1 (MQ=255)
                                gAGCAGGCGATTTCCACCGAGGCTTATATTGCACTCAATCGTGCAGCCATCGACAGCGACCTGGTTGATATGATCGATCTGGAGTTATTTAcc                                                         <  1:5996804/93‑1 (MQ=255)
                                 aGCAGGCGATTTCCACCGAGGCTTATATTGCACTCAATCGTGCAGCCATCGACAGCGACCTGGTTGATATGATCGATCTGGAGTTATTTACCg                                                        >  2:567963/1‑93 (MQ=255)
                                  gCAGGCGATTTCCACCGAGGCTTATATTGCACTCAATCGTGCAGCCATCGACAGCGACCTGGTTGATATGATCGATCTGGAGTTATTTACCgg                                                       <  2:5883156/93‑1 (MQ=255)
                                    aGGCGATTTCCACCGAGGCTTATATTGCACTCAATCGTGCAGCCATCGACAGCGACCTGGTTGATATGATCGATCTGGAGTTATTTACCGGtg                                                     <  2:3801909/93‑1 (MQ=255)
                                     ggCGATTTCCACCGAGGCTTATATTGCACTCAATCGTGCAGCCATCGACAGCGACCTGGTTGATATGATCGATCTGGAGTTATTTACCGGtga                                                    >  2:2555115/1‑93 (MQ=255)
                                      gCGATTTCCACCGAGGCTTATATTGCACTCAATCGTGCAGCCATCGACAGCGACCTGGTTGATATGATCGATCTGGAGTTATTTACCGGtgat                                                   <  2:1348225/93‑1 (MQ=255)
                                      gCGATTTCCACCGAGGCTTATATTGCACTCAATCGTGCAGCCATCGACAGCGACCTGGTTGATATGATCGATCTGGAGTTATTTACCGGtgat                                                   >  2:1226172/1‑93 (MQ=255)
                                      gCGATTTCCACCGAGGCTTATATTGCACTCAATCGTGCAGCCATCGACAGCGACCTGGTTGATATGATCGATCTGGAGTTATTTACCGGtgat                                                   <  1:658767/93‑1 (MQ=255)
                                        gATTTCCACCGAGGCTTATATTGCACTCAATCGTGCAGCCATCGACAGCGACCTGGTTGATATGATCGATCTGGAGTTATTTACCGGtgatga                                                 <  2:4110399/93‑1 (MQ=255)
                                        gATTTCCACCGAGGCTTATATTGCACTCAATCGTGCAGCCATCGACAGCGACCTGGTTGATATGATCGATCTGGAGTTATTTACCGGtgatga                                                 <  2:4517191/93‑1 (MQ=255)
                                         aTTTCCACCGAGGCTTATATTGCACTCAATCGTGCAGCCATCGACAGCGACCTGGTTGATATGATCGATCTGGAGTTATTTACCGGTgatgat                                                <  1:1829377/93‑1 (MQ=255)
                                         aTTTCCACCGAGGCTTATATTGCACTCAATCGTGCAGCCATCGACAGCGACCTGGTTGATATGATCGATCTGGAGTTATTTACCGGTgatgat                                                <  2:6036483/93‑1 (MQ=255)
                                          tttCCACCGAGGCTTATATTGCACTCAATCGTGCAGCCATCGACAGCGACCTGGTTGATATGATCGATCTGGAGTTATTTACCGGTGATGATc                                               <  1:3296135/93‑1 (MQ=255)
                                          tttCCACCGAGGCTTATATTGCACTCAATCGTGCAGCCATCGACAGCGACCTGGTTGATATGATCGATCTGGAGTTATTTACCGGTGATGATc                                               <  1:3671978/93‑1 (MQ=255)
                                          tttCCACCGAGGCTTATATTGCACTCAATCGTGCAGCCATCGACAGCGACCTGGTTGATATGATCGATCTGGAGTTATTTACCGGTGATGATc                                               >  2:4818967/1‑93 (MQ=255)
                                          tttCCACCGAGGCTTATATTGCACTCAATCGTGCAGCCATCGACAGCGACCTGGTTGAGATGATCGATCTGGAGTTATTTACCGGTGATGATc                                               >  1:5151165/1‑93 (MQ=255)
                                           ttCCACCGAGGCTTATATTGCACTCAATCGTGCAGCCATCGACAGCGACCTGGTTGATATGATCGATCTGGAGTTATTTACCGGTGATGATCa                                              <  2:5505110/93‑1 (MQ=255)
                                              cACCGAGGCTTATATTGCACTCAATCGTGCAGCCATCGACAGCGACCTGGTTGATATGATCGATCTGGAGTTATTTACCGGTGATGATCAGGt                                           >  1:882916/1‑93 (MQ=255)
                                               aCCGAGGCTTATATTGCACTCAATCGTGCAGCCATCGACAGCGGCCTGGCTGATATGATCGATCTGGAGTTATTTACCGGTGATGATCAGGtt                                          >  1:1039898/1‑93 (MQ=255)
                                                  gAGGCTTATATTGCACTCAATCGTGCAGCCATCGACAGCGACCTGGTTGATATGATCGATCTGGAGTTATTTACCGGTGATGATCAGGTTaaa                                       >  1:1513340/1‑93 (MQ=255)
                                                    ggCTTATATTGCACTCAATCGTGCAGCCATCGACAGCGACCTGGTTGATATGATCGATCTGGAGTTATTTACCGGTGATGATCAGGTTAAAGa                                     >  1:4555581/1‑93 (MQ=255)
                                                    ggCTTATATTGCACTCAATCGTGCAGCCATCGACAGCGACCTGGTTGATATGATCGATCTGGAGTTATTTACCGGTGATGATCAGGTTAAAGa                                     <  1:1206544/93‑1 (MQ=255)
                                                    ggCTTATATTGCACTCAATCGTGCAGCCATCGACAGCGACCTGGTTGATATGATCGATCTGGAGTTATTTACCGGTGATGATCAGGTTAAAGa                                     >  1:2323215/1‑93 (MQ=255)
                                                     gCTTATATTGCACTCAATCGTGCAGCCATCGACAGCGACCTGGTTGATATGATCGATCTGGAGTTATTTACCGGTGATGATCAGGTTAAAGaa                                    >  2:4889773/1‑93 (MQ=255)
                                                      cTTATATTGCACTCAATCGTGCAGCCATCGACAGCGACCTGGTTGATATGATCGATCTGGAGTTATTTACCGGTGATGATCAGGTTAAAGaaa                                   <  1:4263424/93‑1 (MQ=255)
                                                      cTTATATTGCACTCAATCGTGCAGCCATCGACAGCGACCTGGTTGATATGATCGATCTGGAGTTATTTACCGGTGATGATCAGGTTAAAGaaa                                   <  2:3593977/93‑1 (MQ=255)
                                                      cTTATATTGCACTCAATCGTGCAGCCATCGACAGCGACCTGGTTGATATGATCGATCTGGAGTTATTTACCGGTGATGATCAGGTTAAAGaaa                                   >  2:2902808/1‑93 (MQ=255)
                                                       ttATATTGCACTCAATCGTGCAGCCATCGACAGCGACCTGGTTGATATGATCGATCTGGAGTTATTTACCGGTGATGATCAGGTTAAAGAAAc                                  <  1:4282254/93‑1 (MQ=255)
                                                       ttATATTGCACTCAATCGTGCAGCCATCGACAGCGACCTGGTTGATATGATCGATCTGGAGTTATTTACCGGTGATGATCAGGTTAAAGAAAc                                  <  1:3087146/93‑1 (MQ=255)
                                                        tataTTGCACTCAATCGTGCAGCCATCGACAGCGACCTGGTTGATATGATCGATCTGGAGTTATTTACCGGTGATGATCAGGTTAAAGAAAcc                                 >  1:471471/1‑93 (MQ=255)
                                                        tataTTGCACTCAATCGTGCAGCCATCGACAGCGACCTGGTTGATATGATCGATCTGGAGTTATTTACCGGTGATGATCAGGTTAAAGAAAcc                                 >  2:5520665/1‑93 (MQ=255)
                                                         atatTGCACTCAATCGTGCAGCCATCGACAGCGACCTGGTTGATATGATCGATCTGGAGTTATTTACCGGTGATGATCAGGTTAAAGAAACCg                                >  2:3786524/1‑93 (MQ=255)
                                                         atatTGCACTCAATCGTGCAGCCATCGACAGCGACCTGGTTGATATGATCGATCTGGAGTTATTTACCGGTGATGATCAGGTTAAAGAAACCg                                <  1:5696485/93‑1 (MQ=255)
                                                         atatTGCACTCAATCGTGCAGCCATCGACAGCGACCTGGTTGATATGATCGATCTGGAGTTATTTACCGGTGATGATCAGGTTAAAGAAACCg                                >  1:2944336/1‑93 (MQ=255)
                                                          tatTGCACTCAATCGTGCAGCCATCGACAGCGACCTGGTTGATATGATCGATCTGGAGTTATTTACCGGTGATGATCAGGTTAAAGAAACCGt                               >  1:2821276/1‑93 (MQ=255)
                                                          tatTGCACTCAATCGTGCAGCCATCGACAGCGACCTGGTTGATATGATCGATCTGGAGTTATTTACCGGTGATGATCAGGTTAAAGAAACCGt                               <  2:3472860/93‑1 (MQ=255)
                                                          tatTGCACTCAATCGTGCAGCCATCGACAGCGACCTGGTTGATATGATCGATCTGGAGTTATTTACCGGTGATGATCAGGTTAAAGAAACCGt                               <  2:3875745/93‑1 (MQ=255)
                                                           atTGCACTCAATCGTGCAGCCATCGACAGCGACCTGGTTGATATGATCGATCTGGAGTTATTTACCGGTGATGATCAGGTTAAAGAAACCGTc                              >  1:4434849/1‑93 (MQ=255)
                                                           atTGCACTCAATCGTGCAGCCATCGACAGCGACCTGGTTGATATGATCGATCTGGAGTTATTTACCGGTGATGATCAGGTTAAAGAAACCGTc                              <  1:6121726/93‑1 (MQ=255)
                                                            ttGCACTCAATCGTGCAGCCATCGACAGCGACCTGGTTGATATGATCGATCTGGAGTTATTTACCGGTGATGATCAGGTTAAAGAAACCGTCg                             >  2:2452204/1‑93 (MQ=255)
                                                            ttGCACTCAATCGTGCAGCCATCGACAGCGACCTGGTTGATATGATCGATCTGGAGTTATTTACCGGTGATGATCAGGTTAAAGAAACCGTCg                             >  2:411499/1‑93 (MQ=255)
                                                            ttGCACTCAATCGTGCAGCCATCGACAGCGACCTGGTTGATATGATCGATCTGGAGTTATTTACCGGTGATGATCAGGATAAAGAAACCGTCg                             >  2:5445211/1‑93 (MQ=255)
                                                               cACTCAATCGTGCAGCCATCGACAGCGACCTGGTTGATATGATCGATCTGGAGTTATTTACCGGTGATGATCAGGTTAAAGAAACCGTCGCCt                          >  2:1690533/1‑93 (MQ=255)
                                                                ccTCAATCGTGCAGCCATCGACAGCGACCTGGTTGATATGATCGATCTGGAGTTATTTACCGGTGATGATCAGGTTAAAGAAACCGTCGCCTa                         <  1:4864704/92‑1 (MQ=255)
                                                                aCTCAATCGTGCAGCCATCGACAGCGACCTGGTTGATATGATCGATCTGGAGTTATTTACCGGTGATGATCAGGTTAAAGAAACCGTCGCCTa                         >  1:25093/1‑93 (MQ=255)
                                                                aCTCAATCGTGCAGCCATCGACAGCGACCTGGTTGATATGATCGATCTGGAGTTATTTACCGGTGATGATCAGGTTAAAGAAACCGTCGCCTa                         <  1:3908618/93‑1 (MQ=255)
                                                                   cAATCGTGCAGCCATCGACAGCGACCTGGTTGATATGATCGATCTGGAGTTATTTACCGGTGATGATCAGGTTAAAGAAACCGTCGCCTACGc                      >  2:969955/1‑93 (MQ=255)
                                                                    aaTCGTGCAGCCATCGACAGCGACCTGGTTGATATGATCGATCTGGAGTTATTTACCGGTGATGATCAGGTTAAAGAAACCGTCGCCTACGcc                     <  2:911426/93‑1 (MQ=255)
                                                                      tCGTGCAGCCATCGACAGCGACCTGGTTGATATGATCGATCTGGAGTTATTTACCGGTGATGATCAGGTTAAAGAAACCGTCGCCTACGCCCa                   <  1:1026363/93‑1 (MQ=255)
                                                                      tCGTGCAGCCATCGACAGCGACCTGGTTGATATGATCGATCTGGAGTTATTTACCGGTGATGATCAGGTTAAAGAAACCGTCGCCTACGCCCa                   <  1:1892354/93‑1 (MQ=255)
                                                                      tCGTGCAGCCATCGACAGCGACCTGGTTGATATGATCGATCTGGAGTTATTTACCGGTGATGATCAGGTTAAAGAAACCGTCGCCTAAGcccc                   >  2:685373/1‑92 (MQ=255)
                                                                       cGTGCAGCCATCGACAGCGACCTGGTTGATATGATCGATCTGGAGTTATTTACCGGTGATGATCAGGTTAAAGAAACCGTCGCCTACGCCCAc                  <  2:1269123/93‑1 (MQ=255)
                                                                        gTGCAGCCATCGACAGCGACCTGGTTGATATGATCGATCTGGAGTTATTTACCGGTGATGATCAGGTTAAAGAAACCGTCGCCTACGCCCAcg                 >  1:4403658/1‑93 (MQ=255)
                                                                         tGCAGCCATCGACAGCGACCTGGTTGATATGATCGATCTGGAGTTATTTACCGGTGATGATCAGGTTAAAGAAACCGTCGCCTACGCCCAcgc                >  1:5331190/1‑93 (MQ=255)
                                                                          gCAGCCATCGACAGCGACCTGGTTGATATGATCGATCTGGAGTTATTTACCGGTGATGATCAGGTTAAAGAAACCGTCGCCTACGCCCAcgcg               >  1:654716/1‑93 (MQ=255)
                                                                            aGCCATCGACAGCGACCTGGTTGATATGATCGATCTGGAGTTATTTACCGGTGATGATCAGGTTAAAGAAACCGTCGCCTACGCCCACGCGCa             <  2:5047500/93‑1 (MQ=255)
                                                                             gCCATCGACAGCGACCTGGTTGATATGATCGATCTGGAGTTATTTACCGGTGATGATCAGGTTAAAGAAACCGTCGCCTACGCCCACGCGCat            >  1:338051/1‑93 (MQ=255)
                                                                                aTCGACAGCGACCTGGTTGATATGATCGATCTGGAGTTATTTACCGGTGATGATCAGGTTAAAGAAACCGTCGCCTACGCCCACGCGCatgat         >  2:3006130/1‑93 (MQ=255)
                                                                                aTCGACAGCGACCTGGTTGATATGATCGATCTGGAGTTATTTACCGGTGATGATCAGGTTAAAGAAACCGTCGCCTACGCCCACGCGCatgat         >  1:1652315/1‑93 (MQ=255)
                                                                                aTCGACAGCGACCTGGTTGATATGATCGATCTGGAGTTATTTACCGGTGATGATCAGGTTAAAGAAACCGTCGCCTACGCCCACGCGCatgat         <  2:2234344/93‑1 (MQ=255)
                                                                                aTCGACAGCGACCTGGTTGATATGATCGATCTGGAGTTATTTACCGGTGATGATCAGGTTAAAGAAACCGTCGCCTACGCCCACGCGCatgat         >  2:4000481/1‑93 (MQ=255)
                                                                                 tCGACAGCGACCTGGTTGATATGATCGATCTGGAGTTATTTACCGGTGATGATCAGGTTAAAGAAACCGTCGCCTACGCCCACGCGCatgatg        >  1:228315/1‑93 (MQ=255)
                                                                                    aCAGCGACCTGGTTGATATGATCGATCTGGAGTTATTTACCGGTGATGATCAGGTTAAAGAAACCGTCGCCTACGCCCACGCGCATGATGTGa     <  1:5672073/93‑1 (MQ=255)
                                                                                      aGCGACCTGGTTGATATGATCGATCTGGAGTTATTTACCGGTGATGATCAGGTTAAAGAAACCGTCGCCTACGCCCACGCGCATGATGTGaaa   <  1:3619217/93‑1 (MQ=255)
                                                                                      aGCGACCTGGTTGATATGATCGATCTGGAGTTATTTACCGGTGATGATCAGGTTAAAGAAACCGTCGCCTACGCCCACGCGCATGATGTGaaa   <  2:5404269/93‑1 (MQ=255)
                                                                                       gCGACCTGTTTGATATGATCGATCTGGAGTTATTTACCGGTGATGATCAGGTTAAAGAAACCGTCGCCTACGCCCACGCGCATGATGTGAAag  <  1:4656846/93‑1 (MQ=255)
                                                                                          |                                                                                         
TGTTTACCTTCCGCAGTGCCAAAGAAGGCGGCGAGCAGGCGATTTCCACCGAGGCTTATATTGCACTCAATCGTGCAGCCATCGACAGCGGCCTGGTTGATATGATCGATCTGGAGTTATTTACCGGTGATGATCAGGTTAAAGAAACCGTCGCCTACGCCCACGCGCATGATGTGAAAG  >  NC_000913/1772942‑1773121

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