Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,141,832 A→T D128E (GAT→GAA)  dcd ← deoxycytidine triphosphate deaminase; dCTP deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,141,8320AT98.6% 254.2 / NA 71D128E (GAT→GAA) dcddeoxycytidine triphosphate deaminase; dCTP deaminase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base T (31/38);  total (33/38)
Fisher's exact test for biased strand distribution p-value = 4.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

ACGGCGGTTGTAAGGTCGCACCGCCGGGCCGGAAAGCGGCTCAAAGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCAGCGTCACCGCCAGCGCCAGCTCTCCTGGGTGAAGATAGAACGC  >  NC_000913/2141694‑2141978
                                                                                                                                          |                                                                                                                                                  
aCGGCGGTTGTAAGGTCGCACCGCCGGGCCGGAAAGCGGCTCAAAGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGTTCGATGCGGTgc                                                                                                                                        <  1:214733/151‑1 (MQ=255)
     ggTTGTAAGGTCGCACCGCCGGGCCGGAAAGCGGCTCAAAGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGTTCGATGCGGTGCGCgg                                                                                                                                    >  2:172756/1‑150 (MQ=255)
      gTTGTAAGGTCGCGCCGCCGGGCCGGAAAGCGGCTCAAAGCTCAGCGAACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGTTCGATGCGGTGCGCGGTg                                                                                                                                  <  1:662273/151‑1 (MQ=255)
      gTTGTAAGGTCGCACCGCCGGGCCGGAAAGCGGCTCAAAGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGTTCGATGCGGTGCGCGGTg                                                                                                                                  >  1:1335277/1‑151 (MQ=255)
       ttGTAAGGTCGCACCGCCGGGCCGGAAAGCGGCTCAAAGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGGTCGATGCGGTGCGCGGTg                                                                                                                                  >  1:898803/1‑150 (MQ=255)
        tGTAAGGTCGCACCGCCGGGCCGGAAAGCGGCTCAAAGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGTTCGATGCGGTGCGCGGTGa                                                                                                                                 <  2:1252715/150‑1 (MQ=255)
         gTAAGGTCGCACCGCCGGGCCGGAAAGCGGCTCAAAGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGTTCGATGCGGTGCGCGGTGAc                                                                                                                                <  2:202658/150‑1 (MQ=255)
         gTAAGGTCGCACCGCCGGGCCGGAAAGCGGCTCAAAGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGTTCGATGCGGTGCGCGGTGACg                                                                                                                               <  1:1185987/151‑1 (MQ=255)
          tAAGCTCGCACCGCCGGGCCCGAAAGCGGCTCAAAGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTCCCGGAGTTGTAGAACTCCAGCACAATGCAACCCGACCAGCCCGGTTCGATGCGGTGCGCGGTGACg                                                                                                                               <  1:1212819/150‑1 (MQ=255)
                cGCACCGCCGGGCCGGAAAGCGGCTCAAAGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGTTCGATGCGGTGCGCGGTGACGTGCACCa                                                                                                                        <  1:728663/151‑1 (MQ=255)
                  cACCGCCGGGCCGGAAAGCGGCTCAAAGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGTTCGATGCGGTGCGCGGTGACGTGCACCATc                                                                                                                      >  1:1406929/1‑151 (MQ=255)
                    ccgccgGGCCGGAAAGCGGCTCAAAGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGTTCGATGCGGTGCGCGGTGACGTGCACCATCAg                                                                                                                    <  1:181472/151‑1 (MQ=255)
                    ccgccgGGCCGGAAAGCGGCTCAAAGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGTTCGATGCGGTGCGCGGTGACGTGCACCATCAg                                                                                                                    <  2:205008/151‑1 (MQ=255)
                     cgccgGGCCGGAAAGCGGCTCAAAGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGTTCGATGCGGTGCGCGGTGACGTGCACCATCAGc                                                                                                                   >  2:219945/1‑151 (MQ=255)
                      gccgggccggAAAGCGGCTCAAAGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGTTCGATGCGGTGCGCGGAGACGTGCACCATCAGcc                                                                                                                  >  1:307690/1‑151 (MQ=255)
                        cgggccggAAAGCGGCTCAAAGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGTTCGATGCGGTGCGCGGTGACGTGCACCATCAGcccc                                                                                                                >  2:154014/1‑151 (MQ=255)
                               gAAAGCGGCTCAAAGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGTTCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGAcgcg                                                                                                         >  2:1094626/1‑151 (MQ=255)
                                  aGCGGCTCAAAGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGTTCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCa                                                                                                      <  1:4808/151‑1 (MQ=255)
                                   gCGGCTCAAAGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGTTCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGcc                                                                                                       >  2:44397/1‑149 (MQ=255)
                                    cGGCTCAAAGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGTTCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAg                                                                                                     >  1:1052983/1‑150 (MQ=255)
                                    cGGCTCAAAGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGTTCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGt                                                                                                    <  1:938062/151‑1 (MQ=255)
                                     ggCTCAAAGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGTGCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTg                                                                                                   >  1:917431/1‑151 (MQ=255)
                                      gCTCAAAGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGTTCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGa                                                                                                  <  1:1378129/151‑1 (MQ=255)
                                         cAAAGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGTTCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAgg                                                                                                <  2:322808/150‑1 (MQ=255)
                                         cAAAGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGTTCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAgg                                                                                                <  1:82685/150‑1 (MQ=255)
                                         cAAAGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGTTCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGa                                                                                               <  2:28032/151‑1 (MQ=255)
                                          aaaGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGTTCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGaa                                                                                              <  2:627057/151‑1 (MQ=255)
                                          aaaGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGTTCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGa                                                                                               >  2:146701/1‑150 (MQ=255)
                                          aaaGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACT‑CAGCACAATGCAACCAGACCAGCCCGGTTCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGa                                                                                               >  1:680118/1‑149 (MQ=255)
                                          aaaGCTCAGCGCACCAATTAACAGGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCATCCAGACCAGCCCGGTTCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGaa                                                                                              <  1:1231947/151‑1 (MQ=255)
                                           aaGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGTTCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAAc                                                                                             >  2:1210832/1‑151 (MQ=255)
                                           aaGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGTTCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAAc                                                                                             >  2:1009074/1‑151 (MQ=255)
                                              cTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGTTCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGcc                                                                                          >  1:594331/1‑151 (MQ=255)
                                                    gcACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGTTCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTcca                                                                                    <  2:505880/151‑1 (MQ=255)
                                                     cACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGTTCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTcc                                                                                     >  1:602483/1‑149 (MQ=255)
                                                      aCCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGTTCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTccag                                                                                   <  1:352283/150‑1 (MQ=255)
                                                        cAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGTTCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTccagcca                                                                                <  2:879619/151‑1 (MQ=255)
                                                        cAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGTTCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTccagcca                                                                                <  2:200574/151‑1 (MQ=255)
                                                                aTGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGTTCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCccacc                                                                         >  1:27719/1‑150 (MQ=255)
                                                                 tGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGTTCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCccacc                                                                         >  1:728652/1‑149 (MQ=255)
                                                                     cGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGTTCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCg                                                                   <  1:314173/151‑1 (MQ=255)
                                                                           cAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGTTCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGc                                                             >  2:171278/1‑151 (MQ=255)
                                                                            aGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGTTCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTggg                                                             >  1:1168833/1‑149 (MQ=255)
                                                                            aGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGTTCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCa                                                            <  2:39731/151‑1 (MQ=255)
                                                                              cgcCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGTTCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGc                                                          >  1:785520/1‑151 (MQ=255)
                                                                               gcCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGTTCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCg                                                         >  1:1271166/1‑151 (MQ=255)
                                                                                  aGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGTTCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGt                                                        >  1:1257593/1‑149 (MQ=255)
                                                                                        aGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGTTCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACt                                                <  2:1067922/151‑1 (MQ=255)
                                                                                          cTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGTTCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTcc                                              >  2:353237/1‑151 (MQ=255)
                                                                                             aCCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGTTCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCAg                                            >  1:1243899/1‑150 (MQ=255)
                                                                                              ccGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGTTCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCAGc                                           <  1:811457/150‑1 (MQ=255)
                                                                                                ggAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGTTCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCAGCGTc                                        <  2:29573/151‑1 (MQ=255)
                                                                                                 gAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGTTCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCAGCGTCa                                       <  2:1222590/151‑1 (MQ=255)
                                                                                                    ttGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGTTCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCAGCGTCAcc                                     >  2:59524/1‑150 (MQ=255)
                                                                                                     tGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGTTCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCAGCGTCACcg                                    <  1:544197/150‑1 (MQ=255)
                                                                                                        aGAACTCCAGCACAATGCAACCAGACCAGCCCGGTTCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCAGCGTCACcgcca                                 <  2:728775/150‑1 (MQ=255)
                                                                                                            cTCCAGCACAATGCAACCAGACCAGCCCGGTTCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCAGCGTCACcgccagcgcc                            >  2:505751/1‑151 (MQ=255)
                                                                                                            cTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGAGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCTTCCAGCCAGCCCACCAGAGCCGGTGGCCGCGTCAGCGACTCCTGCCTCACcgccagcg                              >  1:1065212/1‑149 (MQ=255)
                                                                                                              ccAGCACAATGCAACCAGACCAGCCCGGTTCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCGGCGTCACcgccagcgcca                           <  1:826428/150‑1 (MQ=255)
                                                                                                              ccAGCACAATGCAACCAGACCAGCCCGGTTCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCAGCGTCACcgccagcgcca                           <  2:1345907/150‑1 (MQ=255)
                                                                                                              ccAGCACAATGCAACCAGACCAGCCCGGTTCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCAGCGTCACcgccagcgcc                            >  2:310480/1‑149 (MQ=255)
                                                                                                                  cacaATGCAACCAGACCAGCCCGGTTCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCAGCGTCACCGCCAGCGCCAGctct                      >  2:1292206/1‑151 (MQ=255)
                                                                                                                  cacaATGCAACCAGACCAGCCCGGTTCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCAGCGTCACCGCCAGCGCCAGctct                      <  1:393644/151‑1 (MQ=255)
                                                                                                                       tGCAACCAGACCAGCCCGGTTCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCAGCGTCACCGCCAGCGCCAGCTCTCCTgg                 >  1:966954/1‑151 (MQ=255)
                                                                                                                              agaccagCCCGGTTCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCTAGCGTCACCGCCAGCGCCAGCTCTCCTGGGTGAAGa          <  1:4811/151‑1 (MQ=255)
                                                                                                                                 ccagCCCGGTTCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCCCCAGATCGGCTGGCAGCGTCACCGACTCCAGCGTCACCGCCAGCGCCAGCTCTCCTGGGTGAAGAt         <  1:857945/149‑1 (MQ=255)
                                                                                                                                  cagCCCGGTTCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCAGCGTCACCGCCAGCGCCAGCTCTCCTGGGTGAAGATAGa      <  1:172730/151‑1 (MQ=255)
                                                                                                                                     cccGGTTCGATGCGGTGCGCGGTGACGTGCCCCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCCCCAGATCGGCTGGCAGCGTCACCGACTCCAGCGTCACCGCCAGCGCCAGCTCTCCTGGGTGAAGATAGAAc    <  1:29179/150‑1 (MQ=255)
                                                                                                                                     cccGGTTCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCAGCGTCACCGCCAGCGCCAGCTCTCCTGGGTGAAGATAGAACg   <  2:352655/151‑1 (MQ=255)
                                                                                                                                     cccGGTCCGATGCGGTGCGCGGTCGCGCGCCCCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCCCCAGATCGGCTGGCAGCGTCTCAGACCCCAGCGTCGCCGCCAGCGCCAGCTCTCCTGGGTGAAGATAGAAc    <  1:46142/150‑1 (MQ=255)
                                                                                                                                      ccGGTTCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCAGCGTCACCGCCAGCGCCAGCTCTCCTGGGTGAAGATAGAACGc  >  2:224715/1‑151 (MQ=255)
                                                                                                                                          |                                                                                                                                                  
ACGGCGGTTGTAAGGTCGCACCGCCGGGCCGGAAAGCGGCTCAAAGCTCAGCGCACCAATTAACATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCAGCGTCACCGCCAGCGCCAGCTCTCCTGGGTGAAGATAGAACGC  >  NC_000913/2141694‑2141978

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