Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,801,368 (T)7→6 coding (2000/2145 nt) nrdE → ribonucleoside‑diphosphate reductase 2 alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,801,3620T.97.0% 1136.5 / 25.1 263coding (1994/2145 nt)nrdEribonucleoside‑diphosphate reductase 2 alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (4/4);  new base . (140/114);  total (144/119)
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.90e-01

Only 100 of 292 total aligned reads displayed.

GTATCAGGACGCTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACATC  >  NC_000913/2801270‑2801450
                                                                                            |                                                                                        
      ggACGCTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGCGCTGTCGCTGACGCttttttc                                                                                    >  1:7456606/1‑92 (MQ=255)
        aCGCTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTccc                                                                                  >  2:7811313/1‑93 (MQ=255)
          gCTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTccccg                                                                                <  2:1946110/93‑2 (MQ=255)
           cTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGa                                                                              <  1:962756/93‑1 (MQ=255)
           cTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGa                                                                              >  2:1899693/1‑93 (MQ=255)
            ttACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGAt                                                                             >  1:6761151/1‑93 (MQ=255)
            ttACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCCTTTTTCCCCGAt                                                                             >  2:7039169/1‑93 (MQ=255)
             tACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATa                                                                            <  1:6514453/93‑1 (MQ=255)
              aCGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATAc                                                                           <  1:7115152/93‑1 (MQ=255)
              aCGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATAc                                                                           <  2:4554851/93‑1 (MQ=255)
              aCGAAATGGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATAc                                                                           <  2:1661890/93‑1 (MQ=255)
                gAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCg                                                                         >  1:666835/1‑93 (MQ=255)
                 aaaTTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGc                                                                        >  2:1398715/1‑93 (MQ=255)
                 aaaTGGGCTCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGc                                                                        <  1:2443105/93‑1 (MQ=255)
                  aaTTGGCGCAGAAAAGATCATCGACACCTACGCGGCAGCGACTCGCCATGTCGATCAGGGGCTGTTGCTGACGCTCTTTTCCCCGATgccgcc                                                                       >  2:3468231/1‑93 (MQ=255)
                  aaTTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGcc                                                                       >  1:218130/1‑93 (MQ=255)
                    ttGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGccac                                                                     <  1:6717632/93‑1 (MQ=255)
                      ggCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGccacca                                                                   <  2:1737623/93‑1 (MQ=255)
                      ggCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGccacca                                                                   >  2:3257565/1‑93 (MQ=255)
                       gcgcAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCcaccac                                                                  <  2:3706167/93‑1 (MQ=255)
                       gcgcAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCcaccac                                                                  >  2:394055/1‑93 (MQ=255)
                         gcAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTc                                                                <  2:7225098/93‑1 (MQ=255)
                         gcAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTc                                                                <  2:604529/93‑1 (MQ=255)
                         gcAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTc                                                                <  2:5251361/93‑1 (MQ=255)
                         gcAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTc                                                                <  1:888547/93‑1 (MQ=255)
                           aGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTcgc                                                              >  1:7776568/1‑93 (MQ=255)
                             aaaaGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGa                                                            <  2:6561373/93‑1 (MQ=255)
                               aaGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGata                                                          <  2:3523626/93‑1 (MQ=255)
                               aaGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGata                                                          >  1:3446397/1‑93 (MQ=255)
                                aGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGatat                                                         >  1:568861/1‑93 (MQ=255)
                                 gATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATc                                                        <  1:895199/93‑1 (MQ=255)
                                 gATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATc                                                        >  1:2735940/1‑93 (MQ=255)
                                 gATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATc                                                        <  1:6385331/93‑1 (MQ=255)
                                  atcatcGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATca                                                       <  1:6847139/93‑1 (MQ=255)
                                   tcatcGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATcaa                                                      >  1:2662777/1‑93 (MQ=255)
                                    catcGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATcaac                                                     >  2:1117445/1‑93 (MQ=255)
                                    catcGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATcaac                                                     <  1:7302575/93‑1 (MQ=255)
                                         acacCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAgc                                                >  1:1705346/1‑93 (MQ=255)
                                         acacCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAgc                                                >  2:3334861/1‑93 (MQ=255)
                                          cgcgtACGCGGAAGCGTCCCGCCATGTCGATCAGGGGCGGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATACCAACAAAgcg                                               <  1:5107681/89‑1 (MQ=255)
                                           acCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAgcgc                                              <  2:5523086/93‑1 (MQ=255)
                                            ccTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCa                                             <  1:1034180/93‑1 (MQ=255)
                                            ccTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCa                                             >  1:2780861/1‑93 (MQ=255)
                                             cTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAg                                            >  2:5846553/1‑93 (MQ=255)
                                             cTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAg                                            >  1:4683328/1‑93 (MQ=255)
                                               aCGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGAt                                          <  2:4846746/93‑1 (MQ=255)
                                                 gcgGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGAttt                                        >  2:7599015/1‑93 (MQ=255)
                                                 gcgGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGAttt                                        >  1:409441/1‑93 (MQ=255)
                                                 gcgGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGAttt                                        <  1:3554335/93‑1 (MQ=255)
                                                 gcgGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGAttt                                        >  1:5999610/1‑93 (MQ=255)
                                                    gAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACg                                     >  2:4451100/1‑93 (MQ=255)
                                                    gAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACg                                     >  1:6126656/1‑93 (MQ=255)
                                                    gAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACg                                     >  1:3991376/1‑93 (MQ=255)
                                                    gAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACg                                     >  1:3713483/1‑93 (MQ=255)
                                                    gAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACg                                     >  2:2267694/1‑93 (MQ=255)
                                                     aaGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGc                                    >  1:2082417/1‑93 (MQ=255)
                                                        ggACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTg                                 <  2:5574615/92‑1 (MQ=255)
                                                        cGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTg                                 <  2:893512/93‑1 (MQ=255)
                                                         gACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTgg                                >  1:7545058/1‑93 (MQ=255)
                                                          gcTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGgc                               <  2:3078569/92‑1 (MQ=255)
                                                          aCTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGgc                               >  2:4173742/1‑93 (MQ=255)
                                                          aCTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGgc                               >  2:146792/1‑93 (MQ=255)
                                                            tCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGgcgc                             >  1:6281101/1‑93 (MQ=255)
                                                            tCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGgcgc                             >  1:2018112/1‑93 (MQ=255)
                                                              gcaaTGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACAACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCaa                           <  1:4922648/90‑1 (MQ=255)
                                                              gCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATGCCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCaa                           <  1:1361008/93‑1 (MQ=255)
                                                               ccATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCCACAAAGCGCAGATTTACGCCTGGCGCAAg                          <  2:2934850/93‑1 (MQ=255)
                                                               ccATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAg                          <  1:1267959/93‑1 (MQ=255)
                                                               ccATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAg                          <  1:5051987/93‑1 (MQ=255)
                                                               ccATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAg                          <  1:5516107/93‑1 (MQ=255)
                                                               ccATGTCGATCAGGGGCTGTCGCTGACGCTTTTTGCCCCGATACCGCCGCCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAg                          <  2:3898869/93‑1 (MQ=255)
                                                                cATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAgg                         <  1:3081966/93‑1 (MQ=255)
                                                                 aTGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAggg                        >  1:4394514/1‑93 (MQ=255)
                                                                 aTGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAggg                        <  1:5328095/93‑1 (MQ=255)
                                                                  tGTCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGt                       >  1:6993294/1‑93 (MQ=255)
                                                                    tCGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTAt                     <  2:1177893/93‑1 (MQ=255)
                                                                     cGATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATc                    <  1:1539419/93‑1 (MQ=255)
                                                                      gATCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCa                   >  2:6803455/1‑93 (MQ=255)
                                                                        tCAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCaaa                 >  2:6459138/1‑93 (MQ=255)
                                                                         cAGGGGCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCaaaa                <  1:1480221/93‑1 (MQ=255)
                                                                           ggggCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACg              >  2:1239185/1‑93 (MQ=255)
                                                                           ggggCTGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACg              >  2:2102331/1‑93 (MQ=255)
                                                                                tGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTa         <  2:7878053/93‑1 (MQ=255)
                                                                                tGTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTa         >  2:615839/1‑93 (MQ=255)
                                                                                 gTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTAt        <  1:4477572/93‑1 (MQ=255)
                                                                                 gTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTAt        <  2:5381831/93‑1 (MQ=255)
                                                                                 gTCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTAt        <  2:6393441/93‑1 (MQ=255)
                                                                                 gTCGCTGACGCTTTTTTCCCCGATACAGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTAt        >  2:26670/1‑93 (MQ=255)
                                                                                  tCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTAtt       >  1:822471/1‑93 (MQ=255)
                                                                                  tCGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTAtt       <  2:2422807/93‑1 (MQ=255)
                                                                                   cGCTGACGCTTTTTTCCCCGATCCCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTa      <  2:2888065/93‑1 (MQ=255)
                                                                                   cGCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTa      <  1:3403765/93‑1 (MQ=255)
                                                                                    gCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTAc     >  2:7128505/1‑93 (MQ=255)
                                                                                    gCTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTAc     >  2:3781865/1‑93 (MQ=255)
                                                                                     cTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACa    >  1:3601639/1‑93 (MQ=255)
                                                                                     cTGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACa    <  1:1019326/93‑1 (MQ=255)
                                                                                      tGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACAt   >  1:4038966/1‑93 (MQ=255)
                                                                                      tGACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACAt   >  1:22886/1‑93 (MQ=255)
                                                                                       gACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACATc  >  2:5037957/1‑93 (MQ=255)
                                                                                       gACGCTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACATc  <  1:3184462/93‑1 (MQ=255)
                                                                                            |                                                                                        
GTATCAGGACGCTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACATC  >  NC_000913/2801270‑2801450

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