Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,231,994 G→T A51S (GCG→TCG)  cdd → cytidine/deoxycytidine deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,231,9940GT99.0% 338.7 / ‑3.1 97A51S (GCG→TCG) cddcytidine/deoxycytidine deaminase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base T (51/45);  total (52/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.28e-01

TTCAAACCGCTTTTGCCCAACTTGCGGATAACTTGCAATCTGCACTGGAACCTATTCTGGCAGACAAGTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGA  >  NC_000913/2231857‑2232129
                                                                                                                                         |                                                                                                                                       
ttCAAACCGCTTTTGCCCAACTTGCGGATAACTTGCAATCTGCACTGGAACCTATTCTGGCAGACAAGTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGActcg                                                                                                                                       >  1:7518/1‑140 (MQ=255)
  cAAACCGCTTTTGCCCAACTTGCGGATAACTTGCAATCTGCACTGGGACCTATTCTGGCAGACAAGTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACTCGCt                                                                                                                                     <  1:1660429/140‑1 (MQ=255)
  cAAACCGCTTTTGCCCAACTTGCGGATAACTTGCAATCTGCACTGGAACCTATTCTGGCAGACAAGTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACTCGCt                                                                                                                                     <  2:706240/140‑1 (MQ=255)
         cTTTTGCCCAACTTGCGGATAACTTGCAATCTGCACTGGAACCTATTCTGGCAGACAAGTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACTCGCTGGCATTc                                                                                                                              >  1:8391/1‑140 (MQ=255)
         cTTTTGCCCAACTTGCGGATAACTTGCAATCTGCACTGGAACCTATTCTGGCAGACAAGTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACTCGCTGGCATTc                                                                                                                              >  2:1687672/1‑140 (MQ=255)
            ttGCCCAACTTGCGGATAACTTGCAATCTGCACTGGAACCTATTCTGGCAGACAAGTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACTCGCTGGCATTCGCa                                                                                                                           <  2:850138/140‑1 (MQ=255)
                   aCTTGCGGATAACTTGCAATCTGCACTGGAACCTATTCTGGCAGACAAGTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACTCGCTGGCATTCGCACTACTTc                                                                                                                    >  2:350176/1‑140 (MQ=255)
                    cTTGCGGATAACTTGCAATCTGCACTGGAACCTATTCTGGCAGACAAGTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACTCGCTGGCATTCGCACTACTTcc                                                                                                                   >  1:74305/1‑140 (MQ=255)
                      tGCGGATAACTTGCAATCTGCACTGGAACCTATTCTGGCAGACAAGTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACTCGCTGGCATTCGCACTACTTCCGc                                                                                                                 <  1:772989/140‑1 (MQ=255)
                      tGCGGATAACTTGCAATCTGCACTGGAACCTATTCTGGCAGACAAGTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACTCGCTGGCATTCGCACTACTTCCGc                                                                                                                 <  2:1020537/140‑1 (MQ=255)
                      tGCGGATAACTTGCAATCTGCACTGGAACCTATTCTGGCAGACAAGTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACTCGCTGGCATTCGCACTACTTCCGc                                                                                                                 <  2:114514/140‑1 (MQ=255)
                        cGGATAACTTGCAATCTGCACTGGAACCTATTCTGGCAGACAAGTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACTCGCTGGCATTCGCACTACTTCCGCTg                                                                                                               >  2:797720/1‑140 (MQ=255)
                            tAACTTGCAATCTGCACTGGAACCTATTCTGGCAGACAAGTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACTCGCTGGCATTCGCACTACTTCCGCTggcgg                                                                                                           <  2:1518944/140‑1 (MQ=255)
                             aaCTTGCAATCTGCACTGGAACCTATTCTGGCAGACAAGTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACTCGCTGGCATTCGCACTACTTCCGCTggcggc                                                                                                          <  1:70008/140‑1 (MQ=255)
                              aCTTGCAATCTGCACTGGAACCTATTCTGGCAGACAAGTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACTCGCTGGCATTCGCACTACTTCCGCTggcggcg                                                                                                         >  2:579507/1‑140 (MQ=255)
                              aCTTGCAATCTGCACTGGAACCTATTCTGGCAGACAAGTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACTCGCTGGCATTCGCACTACTTCCGCTggcggcg                                                                                                         <  1:797720/140‑1 (MQ=255)
                               cTTGCAATCTGCACTGGAACCTATTCTGGCAGACAAGTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACTCGCTGGCATTCGCACTACTTCCGCTggcggcgg                                                                                                        >  1:1301162/1‑140 (MQ=255)
                                    aaTCTGCACTGGAACCTATTCTGGCAGACAAGTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACTCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCtgt                                                                                                   >  2:1140973/1‑140 (MQ=255)
                                        tGCACTGGAACCTATTCTGGCAGACAAGTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACTCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTgcgc                                                                                               >  1:431356/1‑140 (MQ=255)
                                        tGCACTGGAACCTATTCTGGCAGACAAGTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACTCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTgcgc                                                                                               >  1:1301563/1‑140 (MQ=255)
                                        tGCACTGGAACCTATTCTGGCAGACAAGTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACTCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTgcgc                                                                                               >  1:451335/1‑140 (MQ=255)
                                          cACTGGAACCTATTCTGGCAGACAAGTACTTCCACGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACTCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGt                                                                                             >  2:1277891/1‑140 (MQ=255)
                                              ggAACCTATTCTGGCAGACAAGTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACTCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGc                                                                                         <  1:217609/140‑1 (MQ=255)
                                              ggAACCTATTCTGGCAGACAAGTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACTCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGc                                                                                         <  1:1497946/140‑1 (MQ=255)
                                                  ccTATTCTGGCAGACAAGTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACTCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCAtt                                                                                     >  1:1181244/1‑140 (MQ=255)
                                                    tATTCTGGCAGACAAGTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACTCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGt                                                                                   >  2:1278589/1‑140 (MQ=255)
                                                    tATTCTGGCAGACAAGTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACTCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGt                                                                                   >  1:1415221/1‑140 (MQ=255)
                                                    tATTCTGGCAGACAAGTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACTCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGt                                                                                   >  1:730725/1‑140 (MQ=255)
                                                     aTTCTGGCAGACATGTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACTCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTc                                                                                  <  2:1755358/140‑1 (MQ=255)
                                                     aTTCTGGCAGACAAGTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACTCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTc                                                                                  <  2:1301162/140‑1 (MQ=255)
                                                     aTTCTGGCAGACAAGTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACTCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTc                                                                                  <  2:8391/140‑1 (MQ=255)
                                                      ttCTGGCAGACAAGTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACTCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCg                                                                                 <  2:624049/140‑1 (MQ=255)
                                                      ttCTGGCAGACAAGTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACTCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCg                                                                                 <  2:431356/140‑1 (MQ=255)
                                                          ggCAGACAAGTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACTCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAAtt                                                                             >  1:1243881/1‑140 (MQ=255)
                                                           gCAGACAAGTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACTCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCg                                                                                 <  1:246888/135‑1 (MQ=255)
                                                           gCAGACAAGTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACTCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCg                                                                                 >  2:246888/1‑135 (MQ=255)
                                                           gCAGACAAGTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACTCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAAttt                                                                            >  1:227007/1‑140 (MQ=255)
                                                           gCAGACAAGTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACTCGCTGGCATTCGCACTACTTCCGCTGGCGGCGACCTGTGCGCGTACGCCATTGTCGAAttt                                                                            >  1:1959632/1‑140 (MQ=255)
                                                            cAGACAAGTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACTCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAAtttt                                                                           >  1:566551/1‑140 (MQ=255)
                                                              gACAAGTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACTCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTaa                                                                         >  2:760689/1‑140 (MQ=255)
                                                               aCAAGTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACTCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAAt                                                                        <  2:517407/140‑1 (MQ=255)
                                                               aCAAGTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACTCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAAt                                                                        >  1:1854366/1‑140 (MQ=255)
                                                                cAAGTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACTCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATg                                                                       >  1:1082739/1‑140 (MQ=255)
                                                                       ttCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACTCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGgcgc                                                                <  1:514113/140‑1 (MQ=255)
                                                                         ccccGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACTCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCaa                                                              >  1:2031831/1‑140 (MQ=255)
                                                                           ccGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACTCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAAtt                                                            >  1:679148/1‑140 (MQ=255)
                                                                             gCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACTCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTgc                                                          >  2:528662/1‑140 (MQ=255)
                                                                             gCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACTCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTgc                                                          >  2:1077626/1‑140 (MQ=255)
                                                                              cTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACTCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTgcg                                                         >  1:861976/1‑140 (MQ=255)
                                                                                ttgttgACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACTCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTgcgcg                                                       <  1:1053829/140‑1 (MQ=255)
                                                                                ttgttgACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACTCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTgcgcg                                                       <  1:981246/140‑1 (MQ=255)
                                                                                ttgttgACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACTCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTgcgcg                                                       <  2:1181244/140‑1 (MQ=255)
                                                                                ttgttgACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACTCGCTGGCATTCGCACTACTTCCGCTGGCGACGGCCTGTGCGCGTACACCATTGTCGAATTTTAATGTTGGCGCAATTgcgcg                                                       <  2:7518/140‑1 (MQ=255)
                                                                                 tgttgACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACTCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTgcgcgc                                                      <  2:1305162/140‑1 (MQ=255)
                                                                                    tgACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACTCGCTGGCATTCGCACTACTTCCGCTGGCGGCGACCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGgt                                                   <  2:1959632/140‑1 (MQ=255)
                                                                                     gACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACTCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGgtg                                                  <  2:1082739/140‑1 (MQ=255)
                                                                                        cGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACTCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGa                                               >  1:1582055/1‑140 (MQ=255)
                                                                                         ggggAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACTCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAg                                              >  1:1940110/1‑140 (MQ=255)
                                                                                         ggggAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACTCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAg                                              <  1:1278589/140‑1 (MQ=255)
                                                                                           ggAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACTCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCg                                            >  2:1447394/1‑140 (MQ=255)
                                                                                           ggAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACTCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCg                                            >  1:1432325/1‑140 (MQ=255)
                                                                                                aaGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACTCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAAcc                                       <  2:1940110/140‑1 (MQ=255)
                                                                                                 aGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACTCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCt                                      >  1:968125/1‑140 (MQ=255)
                                                                                                       aTCGCTGAAGAGCGCAACGGGGCTGGACGAAGACTCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTAtt                                >  2:1205824/1‑140 (MQ=255)
                                                                                                           cTGAAGAGCGCAACGGGGCTGGACGAAGACTCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCgg                            <  1:855186/140‑1 (MQ=255)
                                                                                                            tGAAGAGCGCAACGGGGCTGGACGAAGACTCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGt                           >  2:1587856/1‑140 (MQ=255)
                                                                                                             gAAGAGCGCAACGGGGCTGGACGAAGACTCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTg                          <  2:1301563/140‑1 (MQ=255)
                                                                                                             gAAGAGCGCAACGGGGCTGGACGAAGACTCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTg                          >  2:979764/1‑140 (MQ=255)
                                                                                                              aaGAGCGCAACGGGGCTGGACGAAGACTCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGc                         <  2:1415221/140‑1 (MQ=255)
                                                                                                              aaGAGCGCAACGGGGCTGGACGAAGACTCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGc                         <  1:1277891/140‑1 (MQ=255)
                                                                                                              aaGAGCGCAACGGGGCTGGACGAAGACTCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGc                         <  2:730725/140‑1 (MQ=255)
                                                                                                              aaGAGCGCAACGGGGCTGAACGAAGACTCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGc                         >  1:824237/1‑140 (MQ=255)
                                                                                                                   cgcAACGGGGCTGGACGAAGACTCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTg                          >  2:110470/1‑134 (MQ=255)
                                                                                                                   cgcAACGGGGCTGGACGAAGACTCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTg                          <  1:110470/134‑1 (MQ=255)
                                                                                                                   cgcAACGGGGCTGGACGAAGACTCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAata                    >  1:936894/1‑140 (MQ=255)
                                                                                                                   cgcAACGGGGCTGGACGAAGACTCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAata                    >  1:601824/1‑140 (MQ=255)
                                                                                                                   cgcAACGGGGCTGGACGAAGACTCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAata                    >  2:532240/1‑140 (MQ=255)
                                                                                                                   cgcAACGGGGCTGGACGAAGACTCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAata                    >  2:640286/1‑140 (MQ=255)
                                                                                                                    gcAACGGGGCTGGACGAAGACTCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAatat                   <  1:760689/140‑1 (MQ=255)
                                                                                                                     cAACGGGGCTGGACGAAGACTCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATg                  <  1:1420605/140‑1 (MQ=255)
                                                                                                                     cAACGGGACTGGACGAAGACTCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATg                  >  2:1651512/1‑140 (MQ=255)
                                                                                                                       aCGGGGCTGGACGAAGACTCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGa                <  1:96160/140‑1 (MQ=255)
                                                                                                                       aCGGGGCTGGACGAAGACTCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGa                <  2:2031831/140‑1 (MQ=255)
                                                                                                                        cGGGGCTGGACGAAGACTCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGaa               >  1:861292/1‑140 (MQ=255)
                                                                                                                        cGGGGCTGGACGAAGACTCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGaa               >  1:721548/1‑140 (MQ=255)
                                                                                                                           ggCTGGACGAAGACTCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAAttt            >  1:1788539/1‑140 (MQ=255)
                                                                                                                             cTGGACGAAGACTCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTAt          >  2:696816/1‑140 (MQ=255)
                                                                                                                             cTGGACGAAGACTCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTAt          <  2:74305/140‑1 (MQ=255)
                                                                                                                             cTGGACGAAGACTCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTAt          <  1:1587856/140‑1 (MQ=255)
                                                                                                                             cTGGACGAAGACTCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTAt          >  2:1959911/1‑140 (MQ=255)
                                                                                                                                gACGAAGACTCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTgg       <  1:1649748/140‑1 (MQ=255)
                                                                                                                                gACGAAGACTCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTgg       >  2:952474/1‑140 (MQ=255)
                                                                                                                                  cGAAGACTCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTg     <  1:532240/140‑1 (MQ=255)
                                                                                                                                  cGAAGACTCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTg     <  2:1678170/140‑1 (MQ=255)
                                                                                                                                     aGACTCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGa  >  2:713917/1‑140 (MQ=255)
                                                                                                                                     aGACTCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGa  <  1:1447394/140‑1 (MQ=255)
                                                                                                                                     aGACTCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGa  <  2:227007/140‑1 (MQ=255)
                                                                                                                                     aGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGa  >  1:689464/1‑140 (MQ=255)
                                                                                                                                         |                                                                                                                                       
TTCAAACCGCTTTTGCCCAACTTGCGGATAACTTGCAATCTGCACTGGAACCTATTCTGGCAGACAAGTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGA  >  NC_000913/2231857‑2232129

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: