Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,803,347 (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,803,3410T.97.1% 277.3 / 2.2 70coding (1994/2145 nt)nrdEribonucleoside‑diphosphate reductase 2, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base . (34/34);  total (34/36)
Fisher's exact test for biased strand distribution p-value = 4.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.19e-01

TTTACTACCCTGCCCCGTTTATGACTAACGAGAATCTGGCGCTGTATCAGGACGCTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACATCCGCCTGCGTCAGATGGCGCTGGAAGGCACTGAAATTGAAGGCTGC  >  NC_000913/2803206‑2803474
                                                                                                                                       |                                                                                                                                     
tttACTACCCTGCCCCGTTTATGACTAACGAGAATCTGGCGCTGTATCAGGACGCTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCtttt                                                                                                                                    >  2:200516/1‑139 (MQ=255)
     tacCCTGCCCCGTTTATGACTAACGAGAATCTGGCGCTGTATCAGGACGCTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTccc                                                                                                                               <  1:924196/139‑1 (MQ=255)
     tacCCTGCCCCGTTTATGACTAACGAGAATCTGGCGCTGTATCAGGACGCTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTccc                                                                                                                               <  1:200516/139‑1 (MQ=255)
        ccTGCCCCGTTTATGACTAACGAGAATCTGGCGCTGTATCAGGACGCTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGC‑TTTTTTCCCCGa                                                                                                                           <  2:1399895/139‑1 (MQ=255)
           gCCCCGTTTATGACTAACGAGAATCTGGCGCTGTATCAGGACGCTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGC‑TTTTTTCCCCGATAc                                                                                                                        <  2:1513956/139‑1 (MQ=255)
                 tttATGACTAACGAGAATCTGGCGCTGTATCAGGACGCTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGC‑TTTTTTCCCCGATACCGccac                                                                                                                  >  2:1019139/1‑139 (MQ=255)
                        cTAACGAGAATCTGGCGCTGTATCAGGACGCTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGC‑TTTTTTCCCCGATACCGCCACCACTcgc                                                                                                           <  1:710356/139‑1 (MQ=255)
                        cTAACGAGAATCTGGCGCTGTATCAGGACGCTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGC‑TTTTTTCCCCGATACCGCCACCACTcgc                                                                                                           <  2:361000/139‑1 (MQ=255)
                                   cTGGCGCTGTATCAGGACGCTTACGAAATTGGCGCAGAAAAGATCATCGAAACCTACGCGGAAGCGACTCGCCATGTCGATCAGGTGCTGTCGCTGACGC‑TTTTTTCCCCGATACCGCCACCACTCGCGATATcaacaa                                                                                                >  2:1117822/1‑139 (MQ=255)
                                       cgcTGTATCAGGACGCTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGC‑TTTCTTCCCCGATACCGCCACCACTCGCGATATCAACAAAgcg                                                                                            <  2:326472/139‑1 (MQ=255)
                                            tATCAGGACGCTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGC‑TTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGAt                                                                                       >  1:486195/1‑139 (MQ=255)
                                                   aCGCTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGC‑TTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGcc                                                                                >  1:1277490/1‑139 (MQ=255)
                                                   aCGCTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGC‑TTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGcc                                                                                <  2:628732/139‑1 (MQ=255)
                                                    cGCTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGC‑TTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCt                                                                               >  1:674220/1‑139 (MQ=255)
                                                      cTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGC‑TTTTTTCCCCGa                                                                                                                           <  2:1521983/93‑1 (MQ=255)
                                                      cTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGC‑TTTTTTCCCCGa                                                                                                                           >  1:1521983/1‑93 (MQ=255)
                                                      cTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGC‑TTTTTTCCCCGATa                                                                                                                         <  1:894633/95‑1 (MQ=255)
                                                      cTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGC‑TTTTTTCCCCGATACCg                                                                                                                      <  2:68425/98‑1 (MQ=255)
                                                      cTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGC‑TTTTTTCCCCGATACCg                                                                                                                      >  1:68425/1‑98 (MQ=255)
                                                      cTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGC‑TTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGcc                                                                                >  2:829063/1‑136 (MQ=255)
                                                      cTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGC‑TTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGcc                                                                                >  1:767211/1‑136 (MQ=255)
                                                      cTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGC‑TTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGcc                                                                                <  1:829063/136‑1 (MQ=255)
                                                      cTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGC‑TTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGcc                                                                                <  2:767211/136‑1 (MQ=255)
                                                      cTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGC‑TTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTgg                                                                             >  1:969923/1‑139 (MQ=255)
                                                      cTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGC‑TTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTgg                                                                             >  1:582147/1‑139 (MQ=255)
                                                      cTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGACTGTCGCTGACGC‑TTTTTTCCCCGATa                                                                                                                         >  2:894633/1‑95 (MQ=255)
                                                        tACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGC‑TTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGgcg                                                                           >  2:310074/1‑139 (MQ=255)
                                                         aCGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGC‑TTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGgcgc                                                                          >  1:684605/1‑139 (MQ=255)
                                                          cGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGC‑TTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCa                                                                         <  1:1019139/139‑1 (MQ=255)
                                                                   cgcAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGC‑TTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCa                                                                >  1:320144/1‑139 (MQ=255)
                                                                      aGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGC‑TTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCaaaa                                                             <  2:486195/139‑1 (MQ=255)
                                                                      aGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGC‑TTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCaaaa                                                             <  2:752726/139‑1 (MQ=255)
                                                                      aGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGC‑TTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCaaaa                                                             <  2:674220/139‑1 (MQ=255)
                                                                        aaaaGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGC‑TTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACg                                                           <  2:582147/139‑1 (MQ=255)
                                                                        aaaaGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGC‑TTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACg                                                           <  1:310074/139‑1 (MQ=255)
                                                                        aaaaGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGC‑TTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACg                                                           <  1:8924/139‑1 (MQ=255)
                                                                         aaaGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCATGGGCTGTCGCTGACGC‑TTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGc                                                          <  1:1117822/139‑1 (MQ=255)
                                                                                  cGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGC‑TTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACa                                                 >  1:295905/1‑139 (MQ=255)
                                                                                    acacCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCttttttcc                                                                                                                                <  1:1468008/59‑3 (MQ=255)
                                                                                    acacCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCttttttcc                                                                                                                                >  2:1468008/1‑57 (MQ=255)
                                                                                    acacCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGC‑TTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGcc                                                                                >  1:1121983/1‑106 (MQ=255)
                                                                                    acacCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGC‑TTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGcc                                                                                <  2:1121983/106‑1 (MQ=255)
                                                                                      acCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGC‑TTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCaa                                                                        <  2:120081/112‑1 (MQ=255)
                                                                                      acCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGC‑TTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCaa                                                                        >  1:120081/1‑112 (MQ=255)
                                                                                       ccTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGC‑TTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACg                                                           <  2:586042/124‑1 (MQ=255)
                                                                                       ccTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGC‑TTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACg                                                           >  1:586042/1‑124 (MQ=255)
                                                                                           cgcgGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGC‑TTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACATCCGCCTGc                                        >  1:1127882/1‑139 (MQ=255)
                                                                                           cgcgGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGC‑TTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACATCCGCCTGc                                        >  1:210251/1‑139 (MQ=255)
                                                                                                aaGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGC‑TTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAgcg                                                                                            <  2:1221866/82‑1 (MQ=255)
                                                                                                aaGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGC‑TTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAgcg                                                                                            >  1:1221866/1‑82 (MQ=255)
                                                                                                  gCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGC‑TTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACATCCGCCTGCGTCAGAt                                 >  2:504915/1‑139 (MQ=255)
                                                                                                      cTCGCCATGTCGATCAGGGGCTGTCGCTGACGC‑TTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACATCCGCCTGCGTCAGATGgcg                             <  1:504915/139‑1 (MQ=255)
                                                                                                           cATGTCGATCAGGGGCTGTCGCTGACGC‑TTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACATCCGCCTGCGTCAGATGGCGCTGGa                        <  2:1277490/139‑1 (MQ=255)
                                                                                                                 gATCAGGGGCTGTCGCTGACGC‑TTTTTTCCCCGATACCGCCACCACTc                                                                                                             >  1:32193/1‑48 (MQ=255)
                                                                                                                 gATCAGGGGCTGTCGCTGACGC‑TTTTTTCCCCGATACCGCCACCACTc                                                                                                             <  2:32193/48‑1 (MQ=255)
                                                                                                                   tCAGGGGCTGTCGCTGACGC‑TTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACATCCGCCTGCGTCAGa                                  <  2:77091/121‑1 (MQ=255)
                                                                                                                   tCAGGGGCTGTCGCTGACGC‑TTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACATCCGCCTGCGTCAGa                                  >  1:77091/1‑121 (MQ=255)
                                                                                                                   tCAGGGGCTGTCGCTGACGC‑TTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACATCCGCCTGCGTCAGATGGCGCTGGAAGGCACTg                <  2:210251/139‑1 (MQ=255)
                                                                                                                         gCTGTCGCTGACGC‑TTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACg                                                           >  2:1081275/1‑90 (MQ=255)
                                                                                                                         gCTGTCGCTGACGC‑TTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACg                                                           <  1:1081275/90‑1 (MQ=255)
                                                                                                                         gCTGTCGCTGACGC‑TTTTTTCCCCGATACCACCACCACTCGCGATATCAACAAAgcgc                                                                                           >  1:1218272/1‑58 (MQ=38)
                                                                                                                         gCTGTCGCTGACGC‑TTTTTTCCCCGATACCACCACCACTCGCGATATCAACAAAgcgc                                                                                           <  2:1218272/58‑1 (MQ=38)
                                                                                                                            gTCGCTGACGC‑TTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACATCCGCCTGCGTCAGATGGCGCTGGAAGGCACTGAAATTGAAg       >  1:156123/1‑139 (MQ=255)
                                                                                                                            gTCGCTGACGC‑TTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATT‑CATCCGCCTGCGTCAGATGGCGCTTGAAGGGCCTGAAAATGAAgg      >  1:557051/1‑139 (MQ=255)
                                                                                                                             tCGCTGACGC‑TTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGcc                                                                                >  1:1246450/1‑65 (MQ=255)
                                                                                                                             tCGCTGACGC‑TTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGcc                                                                                <  2:1246450/65‑1 (MQ=255)
                                                                                                                                 tGACGC‑TTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGcc                                                                                >  1:1227749/1‑61 (MQ=255)
                                                                                                                                 tGACGC‑TTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGcc                                                                                <  2:1227749/61‑1 (MQ=255)
                                                                                                                                 tGACGC‑TTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACATCCGCCt                                          <  1:1398096/99‑1 (MQ=255)
                                                                                                                                 tGACGC‑TTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACATCCGCCt                                          >  2:1398096/1‑99 (MQ=255)
                                                                                                                                 tGACGC‑TTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACATCCGCCTGCGTCAGATGGCGCTGGa                        >  2:1489440/1‑117 (MQ=255)
                                                                                                                                 tGACGC‑TTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACATCCGCCTGCGTCAGATGGCGCTGGa                        <  1:1489440/117‑1 (MQ=255)
                                                                                                                                 tGACGC‑TTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACATCCGCCTGCGTCAGATGGCGCTGGAAGGCACTGAAATTGAAGGCTGc  >  2:218085/1‑139 (MQ=255)
                                                                                                                                       |                                                                                                                                     
TTTACTACCCTGCCCCGTTTATGACTAACGAGAATCTGGCGCTGTATCAGGACGCTTACGAAATTGGCGCAGAAAAGATCATCGACACCTACGCGGAAGCGACTCGCCATGTCGATCAGGGGCTGTCGCTGACGCTTTTTTTCCCCGATACCGCCACCACTCGCGATATCAACAAAGCGCAGATTTACGCCTGGCGCAAGGGTATCAAAACGCTCTATTACATCCGCCTGCGTCAGATGGCGCTGGAAGGCACTGAAATTGAAGGCTGC  >  NC_000913/2803206‑2803474

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑