Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,141,904 C→G 5.0% V104V (GTG→GTC) ‡ 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,9040CG 5.0% 277.1 / 4.6 89V104V (GTG→GTC) ‡dcddeoxycytidine triphosphate deaminase; dCTP deaminase
Reads supporting (aligned to +/- strand):  ref base C (45/39);  new base G (2/2);  total (47/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.09e-01

CATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCAGCGTCACCGCCAGCGCCAGCTCTCCTGGGTGAAGATAGAACGCCTCGCCCTCGTCGAGAACGATCTCATCGCTCATCACGCGGTCAAGCGCGGCGCTCACTTCATCTTTGGGACCGCT  >  NC_000913/2141757‑2142053
                                                                                                                                                   |                                                                                                                                                     
cATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCccac                                                                                                                                                     >  2:1117227/1‑150 (MQ=255)
    cccGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGa                                                                                                                                                 >  1:1040183/1‑150 (MQ=255)
    cccGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGa                                                                                                                                                 <  1:1888632/150‑1 (MQ=255)
        gACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCCGCCAGCCCACCAGATCgg                                                                                                                                             >  1:1003472/1‑150 (MQ=255)
          cGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCCACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGTACGCCCGTCCAGCCAGCCCACCAGATCGGCt                                                                                                                                           >  2:300197/1‑150 (MQ=255)
          cGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCt                                                                                                                                           >  2:1705409/1‑150 (MQ=255)
          cGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGACCGTCCAGCCAGCCCACCAGATCGGCt                                                                                                                                           >  2:98943/1‑150 (MQ=255)
            cAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGCCGTGCCCCATCAGCCCCAGACGCGCCAGTGAGGACCGCCCGTCCAGCCAGCCCACCAGATCGGCggg                                                                                                                                         >  1:1448756/1‑150 (MQ=255)
            cAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTgg                                                                                                                                         >  1:567888/1‑150 (MQ=255)
            cAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTgg                                                                                                                                         <  1:753971/150‑1 (MQ=255)
              gcgcCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACAAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCCTCAGCGCCCGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGcc                                                                                                                                       >  2:1838967/1‑149 (MQ=255)
               cgcCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCGATCTGGTGGGCTGGctg                                                                                                                                      <  2:874970/150‑3 (MQ=255)
                 ccAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCg                                                                                                                                    <  1:388059/150‑1 (MQ=255)
                   aGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGCTGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAGCGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTc                                                                                                                                  >  1:268164/1‑150 (MQ=255)
                   aGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCAGTCCAGCCAGCCCACCAGATCGGCTGGCAGAGTc                                                                                                                                  >  1:1383275/1‑150 (MQ=255)
                   aGCGGCAGCTTACCGGAGTTGTAGAACTCAAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGAGAGGAACGCCCGTCCAGCCAGCCGATCTGGTGGGCTGGctggacg                                                                                                                                  <  2:1472640/150‑7 (MQ=255)
                    gCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCa                                                                                                                                 <  2:1836560/150‑1 (MQ=255)
                      ggCAGCTTACCGGAGTTTTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGTTCGATGCGGTGCGCGGTGACGTGCAGCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCAcc                                                                                                                               <  2:1233834/150‑1 (MQ=255)
                      ggCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTTGCAGCGTCAcc                                                                                                                               <  1:794542/150‑1 (MQ=255)
                       gCAGGTTACCGTAGTTGTAGAACTCCAGCACATTGCAACCAGACCAGCCCGGAGCGATGCGGTGCGCGGTGACGTGCACCATCAGCCTCAGACGCGCCAGTGAGGATGGCCCGTCCAGCCAGCCTATCAGATCTGCTGGGAGCGTCACCg                                                                                                                              <  2:283383/150‑1 (MQ=255)
                         aGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGCTCGGCTGGCAGCGTCACCgcc                                                                                                                            >  1:1327246/1‑148 (MQ=255)
                          gCTTACCGGAGTTGTAGAACTCCAGCACAATGCACCCAGCCCAGCCCGGAGCGATGCGGTGCGCGGGGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACt                                                                                                                           <  1:726983/150‑1 (MQ=255)
                          gCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACt                                                                                                                           >  2:780154/1‑150 (MQ=255)
                             tgctggAGTTCTACAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCa                                                                                                                        >  2:1677496/5‑150 (MQ=255)
                             tACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCa                                                                                                                        <  2:341303/150‑1 (MQ=255)
                             tACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCa                                                                                                                        >  1:863539/1‑150 (MQ=255)
                             tACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCa                                                                                                                        <  1:860122/150‑1 (MQ=255)
                               ccGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGCCTCCAGc                                                                                                                      >  1:741157/1‑150 (MQ=255)
                                 ggAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGACCGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCAGCGt                                                                                                                    >  2:140405/1‑150 (MQ=255)
                                      tGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCAGCGTCACcg                                                                                                               <  1:1544315/150‑1 (MQ=255)
                                      tGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCAGCGTCACcg                                                                                                               <  2:910366/150‑1 (MQ=255)
                                       gTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCCCCGACTCCAGCGTCACcgc                                                                                                              >  1:571000/1‑150 (MQ=255)
                                        tAGAACTCCAGCACAATGCAACAAGACCAGCCCGAATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCAGCGTCACcgcc                                                                                                             <  2:1834265/150‑1 (MQ=255)
                                         aGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCAGCGTCACcgccc                                                                                                            >  2:1255105/1‑149 (MQ=255)
                                          gAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGCCGTGCCCCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCAGCGTCCCcgccag                                                                                                           >  2:1023991/1‑150 (MQ=255)
                                            aCTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCAGCGTCACcgccagcg                                                                                                         >  1:1121546/1‑150 (MQ=255)
                                              tCCAGCACAATGCAACCAGACCAGCCCGGAGCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCAGCGTCACcgccagcgcc                                                                                                       <  1:858535/150‑1 (MQ=255)
                                               ccAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCAGCGTCACcgccagcgcca                                                                                                      >  1:1510859/1‑150 (MQ=255)
                                               ccAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGGGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCAGCGTCACcgccagcgcca                                                                                                      <  2:41672/150‑1 (MQ=255)
                                                cAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCAGCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCAGCGTCACcgccagcgccag                                                                                                     <  1:1330197/150‑1 (MQ=255)
                                                   cacaCTGCACCCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCAGCGTCACCGCCAGCGCCAGctc                                                                                                  <  1:1877824/150‑1 (MQ=255)
                                                    acaATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCAGCGTCACCGCCAGCGCCAGctct                                                                                                 <  2:805880/150‑1 (MQ=255)
                                                    acaATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCAGCGTCACCGCCAGCGCCAGctct                                                                                                 <  2:1148771/150‑1 (MQ=255)
                                                          cAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCAGCGTCACCGCCAGCGCCAGCTCTCCTggg                                                                                           <  1:804968/150‑1 (MQ=255)
                                                               agaccagCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCAGCGTCACCGCCAGCGCCAGCTCTCCTGGGTGAAg                                                                                      <  1:1524194/150‑1 (MQ=255)
                                                                gaccagCCCGGCTCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCAGCGTCACCGCCAGCGCCAGCTCTCCTGGGTGAAGa                                                                                     <  1:1886334/150‑1 (MQ=255)
                                                                gaccagCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCAGCGTCACCGCCAGCGCCAGCTCTCCTGGGTGAAGa                                                                                     <  2:784297/150‑1 (MQ=255)
                                                                  ccagCCCGGATCGATGCGGTGCGCGGTGACGTGCCCCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCCGCCAGCCCACCAGAGCGGCTGGCAGCGTCACCGACTCCAGCGGCACCGCCCGCGCCAGCTCGCCTGGGTGAAgcaa                                                                                   >  1:122700/1‑147 (MQ=255)
                                                                  ccagCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCAGCGTCACCGCCAGCGCCAGCTCTCCTGGGTGAAGATa                                                                                   >  1:1893252/1‑150 (MQ=255)
                                                                    agCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCAGCGTCACCGCCAGCGCCAGCTCTCCTGGGTGAAGATAGa                                                                                 >  2:1384434/1‑150 (MQ=255)
                                                                           aTCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCAGCGTCACCGCCAGCGCCAGCTCTCCTGGGTGAAGATAGAACGCCTc                                                                          >  2:563510/1‑150 (MQ=255)
                                                                             cGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCAGCGTCACCGCCAGCGCCAGCTCTCCTGGGTGAAGATAGAACGCCTCGc                                                                        <  2:84785/150‑1 (MQ=255)
                                                                             cGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCAGCGTCACCGCCAGCGCCAGCTCTCCTGGGTGAAGATAGAACGCCTCGc                                                                        <  2:244294/150‑1 (MQ=255)
                                                                              gATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCCCCGACTCCAGCGTCACCGCCAGCGCCAGCTCGCCTGGGGGAAGATAGAgcgccgcgcc                                                                       >  2:88972/1‑150 (MQ=255)
                                                                               aTGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCAGCGTCACCGCCAGCGCCAGCTCTCCTGGGTGAAGATAGAACGCCTCGccc                                                                      <  2:1851049/150‑1 (MQ=255)
                                                                                   ggTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCAGCGTCACCGCCAGCGCCAGCTCTCCTGGGTGAAGATAGAACGCCTCGCCCtcgt                                                                  >  1:1014559/1‑150 (MQ=255)
                                                                                       cgcgGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCAGCGTCACCGCCAGCGCCAGCTCTCCTGGGTGAAGATAGAACGCCTCGCCCTCGTCgag                                                              >  2:1471728/1‑150 (MQ=255)
                                                                                           gTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCTGCCCACCAGATCGGCTGGCCGCGTCCCCGCCTCCCGCGTCCCCGCCAGCGCCAGCTCCCCTGGGTGAAGCGAGCGCGCCTCGCCCTCGTCGAGCACg                                                          >  1:1161795/1‑150 (MQ=255)
                                                                                                 tGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCAGCGTCACCGCCAGCGCCAGCTCTCCTGGGTGAAGATAGAACGCCTCGCCCTCGTCGAGAACGATCTCa                                                    >  2:1340783/1‑150 (MQ=255)
                                                                                                  gCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCATCGTCACCGACTCCCGCGTCACCGACAGCGCCAGCTCTCCTGGGTGAAGATAGAACGCCTCGCCCTCGGCGAGACCGATCTCAt                                                   >  2:936414/1‑150 (MQ=255)
                                                                                                   cACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCAGCGTCACCGCCAGCGCCAGCTCTCCTGGGTGAAGATAGAACGCCTCGCCCTCGTCGAGAACGATCTCATc                                                  <  1:445482/150‑1 (MQ=255)
                                                                                                        tCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCAGCGTCACCGCCAGCGCCAGCTCTCCTGGGTGATGATAGAACGCCTCGCCCTCGTCGAGAACGCTCTCATCGCtca                                             >  1:790390/1‑150 (MQ=255)
                                                                                                        tCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCAGCGTCACCGCCAGCGCCAGCTCTCCTGGGTGAAGATAGAACGCCTCGCCCTCGTCGAGGACGATCTCATCGCtca                                             >  1:885441/1‑150 (MQ=255)
                                                                                                        tCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCAGCGTCACCGCCAGCGCCAGCTCTCCTGGGTGAAGATAGAACGCCTCGCCCTCGTCGAGAACGATCTCATCGCtca                                             <  1:648159/150‑1 (MQ=255)
                                                                                                          aGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCAGCGTCACCGCCAGCGCCAGCTCTCCTGGGTGAAGATAGAACGCCTCGCCCTCGTCGAGAACGATCTCATCGCtcatc                                           <  2:1510859/150‑1 (MQ=255)
                                                                                                          aGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCAGCGGCACCGCCAGCGCCAGCTCTCCTGGGTGAAGATAGAACGCCTCGCCCTCGTCGAGAACGATCTCATCGCtcatc                                           <  2:539577/150‑1 (MQ=255)
                                                                                                           gCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCAGCGTCACCGCCAGCGCCAGCTCTCCTGGGTGAAGATAGAACGCCTCGCCCTCGTCGAGAACGATCTCATCGCtcatca                                          <  1:1677496/150‑1 (MQ=255)
                                                                                                            ccccAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCAGCGTCACCGCCAGCGCCAGCTCTGCTGGGTGAAGATAGAACGCCTCGCCCTCGTCGAGAACGATCTCATCGctcctccc                                         >  2:1430066/1‑148 (MQ=255)
                                                                                                               cAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCAGCGTCACCGCCAGCGCCAGCTCTCCTGGGTGAAGATAGAACGCCTCGCCCTCGTCGAGAACGATCTCATCGCTCATCAcgcg                                      <  1:153411/150‑1 (MQ=255)
                                                                                                                   cgcgCCAGTGAGGAACGCGCGTCCAGCCAGCCCACCAGATCGGCTGGCTGCGTCACCGACTCCAGCTTCACCGCCATCTCCTGCTCTCCTGTGTGAAGATAGAACGACTCGCCCTGGTCGAGAACGATCTCTTCGCTCATCACGCGGTct                                  <  2:774316/150‑2 (MQ=255)
                                                                                                                   cgcgCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCAGCGTCACCGCCAGCGCCAGCTCTCCTGGGTGAAGATAGAACGCCTCGCCCTCGGCGAGAACGATCTCATCGCTCATCCCGCGGTCa                                  >  1:787802/1‑150 (MQ=255)
                                                                                                                     cgcCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCCCCGCCTCCAGCGTCACCGCCAGCGCCAGCTCTCCTGGGTGAAGATAGAACGCCTCGCCCTCGTCGAGAACGCTCTCATCGCTCATCACGCGGTCacg                                >  2:444644/1‑148 (MQ=255)
                                                                                                                     cgcCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCAGCGTCACCGCCAGCGCCAGCTCTCCTGGGTGAAGATAGAACGCCTCGCCCTCGTCGAGAACGATCTCATCGCTCATCACGCGGTCAAg                                <  2:60816/150‑1 (MQ=255)
                                                                                                                      gcCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCAGCGGCACCGCCAGCGCCAGCTCTCCTGGGTGAAGATAGAACGCCTCGCCCTCGTCGAGAACGATCTCATCGCTCATCACGCGGTCAAgc                               <  1:1150913/150‑1 (MQ=255)
                                                                                                                       ccAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCAGCGTCACCGCCAGCGCCAGCTCTCCTGGGTGAAGATAGAACGCCTCGCCCTCGTCGAGAACGATCTCATCGCTCATCACGCGGTCAAgcg                              >  2:154019/1‑150 (MQ=255)
                                                                                                                        cAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCCGCGTCACCGCCAGCGCCAGCTCTCCTGGGTGAAGATAGAACGCCTCGCCCTCGTCGAGAACGATCTCATCGCTCATCACGCGGTCAAgcgc                             >  2:1728189/1‑150 (MQ=255)
                                                                                                                              ggAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCAGCGTCACCGCCAGCGCCAGCTCTCCTGGGTGAAGATAGAACGCCTCGCCCTCGTCGAGCACGATCTCATCGCTCATCACGCGGTCAAGCGCGgcgcg                       >  1:66707/1‑149 (MQ=255)
                                                                                                                                aaCGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCCCCGACTCCCGCGTCAACGGCAGCGCCAGCTCTCCTGGGTGCAGATAGACCGCCTCGCCCTAGGCGCGAAGGATCTCAACGCTGATCACGCGGTCACGCGCGGCCCtgc                     >  2:1431827/1‑148 (MQ=255)
                                                                                                                                aaCGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCAGCGTCACCGCCAGCGCCAGCTCTCCTGGGTGAAGATAGAACGCCTCGCCCTCGTCGAGAACGATCTCATCGCTCATCACGCGGTCAAGCGCGGCGCTCa                     <  2:1040183/150‑1 (MQ=255)
                                                                                                                                   gCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCAGCGTCACCGCCAGCGCCAGCTCTCCTGGGTGAAGATAGAACGCCTCGCCCTCGTCGAGAACGATCTCATCGCTCATCACGCGGTCAAGCGCGGCGCTCACtt                  >  2:589971/1‑150 (MQ=255)
                                                                                                                                      cGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCAGCGTCACCGCCAGCGCCAGCTCTCCTGGGTGAAGATAGAACGCCTCGCCCTCGTCGAGAACGATCTCATCGCTCATCACGCGGTCAAGCGCGGCGCTCACTTCAt               >  2:1647061/1‑150 (MQ=255)
                                                                                                                                       gTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCAGCGTCACCGCCAGCGCCAGCTCTCCTGGGTGAAGATAGAACGCCTCGCCCTCGTCGAGAACGATCTCATCGCTCATCACGCGGTCAAGCGCGGCGCTCACTTCATc              <  2:1893252/150‑1 (MQ=255)
                                                                                                                                       gTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCAGCGTCACCGCCAGCGCCAGCTCTCCTGGGTGAAGATAGAACGCCTCGCCCTCGTCGAGAACGATCTCATCGCTCATCACGCGGTCAAGCGCGGCGCTCACTTCATc              <  2:1053997/150‑1 (MQ=255)
                                                                                                                                          ctgccagccGATCTGGTGGGCTGGCAGCGTCACCGACTCCAGCGTCACCGCCAGCGCCAGCTCTCCTGGGTGAAGATTGAACGCCTCGCCCTCGTCGAGAACGATCTCATCGCTCATCACGCGGTCAAGCGCGGCGCTCCCTTCATCttt           >  2:564782/3‑150 (MQ=255)
                                                                                                                                          ctgccagccGATCTGGTGGGCTGGCAGCGTCACCGACTCCAGCGTCACCGCCAGCGCCAGCTCTCCTGGGTGAAGATAGAACGCCTCGCCCTCGTCGAGAACGATCTCATCGCTCATCACGCGGTCAAGCGCGGCGCTCACTTCATCttt           >  2:1278984/3‑150 (MQ=255)
                                                                                                                                           agccagccCACCAGATCGGCTGGCAGCGTCACCGACTCCAGCGTCACCGCCAGCGCCAGCTCTCCTGGGTGAAGATAGAACGCCTCGCCCTCGACGAGAACGATATCATCGCTCATCACGCGGTCAACCGCGGCGCTCACTTCATCGTTg          >  2:618006/1‑150 (MQ=255)
                                                                                                                                            gccagccCACCAGATCGGCTGGCAGCGTCACCGACTCCAGCGTCACCGCCAGCGCCAGCTCTCCTGGGTGAAGATAGAACGCCTCGCCCTCGTCGAGCACGATCTCATCGCTCATCACGCGGTCAAGCGCGGCGCTCACTTCATCTTTgg         >  2:288954/1‑150 (MQ=255)
                                                                                                                                             ccagccCACCAGATCGGCTGGCAGCGTCACCGACTCCAGCGTCACCGCCAGCGCCAGCTCTCCTGGGTGAAGATAGAACGCCTCGCCCTCGTCGAGAACGATCTCATCGCTCATCACGCGGTCAAGCGCGGCGCTCACTTCATCTTTggg        <  1:98943/150‑1 (MQ=255)
                                                                                                                                              cagccCACCAGATCGGCTGGCAGCGTCACCGACTCCAGCGTCACCGCCAGCGCCAGCTCTCCTGGGTGAAGATAGAACGCCTCGCCCTCGTCGAGAACGATCTCATCGCTCATCACGCGGTCAAGCGCGGCGCTCACTTCATCTTTGGGa       <  2:443451/150‑1 (MQ=255)
                                                                                                                                              cagccCACCAGATCGGCTGGCAGCGTCACCGAATCCATCGTCACCGCCAGCCCCAGATCACCAGGGTGAAGATAGAACGCCTCGCCCACTTCGAGAACTATATCATAGAGCATCACTCGGTCAATCGCGGCGCTCACTTCATCATTGGGa       >  2:1567549/1‑150 (MQ=255)
                                                                                                                                                   caccaGATCGGCTGGCAGCGTCACCGACTCCAGCGTCACCGCCAGCGCCAGCTCTCCTGGGTGAAGATAGAACGCCTCGCCCTCGTCGAGAACGATCTCATCGCTCATCACGCGGGCAAGCGCGGCGCTCACTTCATCTTTGGGACCGcg  >  2:652211/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CATGCCCGGACGCAGCGCCAGCGGCAGCTTACCGGAGTTGTAGAACTCCAGCACAATGCAACCAGACCAGCCCGGATCGATGCGGTGCGCGGTGACGTGCACCATCAGCCCCAGACGCGCCAGTGAGGAACGCCCGTCCAGCCAGCCCACCAGATCGGCTGGCAGCGTCACCGACTCCAGCGTCACCGCCAGCGCCAGCTCTCCTGGGTGAAGATAGAACGCCTCGCCCTCGTCGAGAACGATCTCATCGCTCATCACGCGGTCAAGCGCGGCGCTCACTTCATCTTTGGGACCGCT  >  NC_000913/2141757‑2142053

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