Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,232,033 G→T 75.7% A64S (GCG→TCG)  cdd → cytidine/deoxycytidine deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,232,0330GT75.7% 192.4 / 88.8 107A64S (GCG→TCG) cddcytidine/deoxycytidine deaminase
Reads supporting (aligned to +/- strand):  ref base G (13/13);  new base T (36/45);  total (49/58)
Fisher's exact test for biased strand distribution p-value = 6.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.91e-02

Only 100 of 107 total aligned reads displayed.

ATCTGCACTGGAACCTATTCTGGCAGACAAGTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGAACAAAGCGCGATCAGCCACGCCT  >  NC_000913/2231894‑2232177
                                                                                                                                           |                                                                                                                                                
aTCTGCACTGGAACCTATTCTGGCAGACAAGTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGc                                                                                                                                         >  4:606100/1‑149 (MQ=255)
       cTGGAACCTATTCTGGCAGACAAGTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTTCGCGTACGCCATTGTc                                                                                                                                  <  2:83852/149‑1 (MQ=255)
       cTGGAACCTATTCTGGCAGACAAGTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCACTGTc                                                                                                                                  <  3:606100/149‑1 (MQ=255)
          gAACCTATTCTGGCAGACAAGTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTTCGCGTACGCCATTGTCGaa                                                                                                                               >  1:876099/1‑149 (MQ=255)
            aCCTATTCTGGCAGACAAGTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTTCGCGTACGCCATTGTCGAAtt                                                                                                                             <  2:363334/149‑1 (MQ=255)
             ccTATTCTGGCAGACAAGTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTTCGCGTACGCCATTGTCGAAttt                                                                                                                            >  4:353516/1‑149 (MQ=255)
             ccTATTCTGGCAGACAAGTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAAttt                                                                                                                            >  1:827300/1‑149 (MQ=255)
              cTATTCTGGCAGACAAGTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTTCGCGTACGCCATTGTCGAAtttt                                                                                                                           >  4:414372/1‑149 (MQ=255)
              cTATTCTGGCAGACAAGTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTTCGCGTACGCCATTGTCGAAtttt                                                                                                                           <  3:432266/149‑1 (MQ=255)
                aTTCTGGCAGACAAGTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTTCGCGTACGCCATTGTCGAATTTTaa                                                                                                                         <  3:647111/149‑1 (MQ=255)
                aTTCTGGCAGACAAGTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTaa                                                                                                                         <  3:103912/149‑1 (MQ=255)
                aTTCTGGCAGACAAGTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTaa                                                                                                                         <  4:469855/149‑1 (MQ=255)
                aTTCTGGCAGACAAGTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCAGCTGGCGGCGGCCTGTTCGCGTTCGCCATTGTCGAATTTTaa                                                                                                                         <  3:731182/149‑1 (MQ=255)
                 ttCTGGCAGACAAGTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTTCGCGTACGCCATTGTCGAATTTTAAt                                                                                                                        <  3:353516/149‑1 (MQ=255)
                      gCAGACAAGTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTTCGCGTACGCCATTGTCGAATTTTAATGTTgg                                                                                                                   <  2:408552/149‑1 (MQ=255)
                            aaGTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGCCGGCCTGTTCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAAt                                                                                                             <  2:792777/149‑1 (MQ=255)
                              gTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTTCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTg                                                                                                           >  1:51418/1‑149 (MQ=255)
                               tACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGAAGCGCTGACATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCGATTGTCGAATTTTAATGTTGGCGCAATTgc                                                                                                          >  4:79462/1‑149 (MQ=255)
                                  ttCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTTCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTgcgcg                                                                                                       >  4:156762/1‑149 (MQ=255)
                                  ttCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTTCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTgcgcg                                                                                                       >  3:456011/1‑149 (MQ=255)
                                    ccccGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGcgcgcg                                                                                                     >  2:783127/1‑149 (MQ=255)
                                     cccGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTTCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCgg                                                                                                    <  4:636157/149‑1 (MQ=255)
                                     cccGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCgg                                                                                                    >  2:356932/1‑149 (MQ=255)
                                     accGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTTCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCCCGCgg                                                                                                    >  4:599884/2‑149 (MQ=255)
                                        gCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTTCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGgtgt                                                                                                 <  4:543244/149‑1 (MQ=255)
                                          tttGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGa                                                                                               <  1:455581/149‑1 (MQ=255)
                                           ttgttgACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAg                                                                                              <  1:783127/149‑1 (MQ=255)
                                              ttgACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTc                                                                                                                                  >  4:752355/1‑110 (MQ=255)
                                              ttgACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTc                                                                                                                                  <  3:752355/110‑1 (MQ=255)
                                                    ggggAGCAAGTCTCATCGCTGAAGAGCGCCACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTTCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTg                                                                                     <  2:581344/149‑1 (MQ=255)
                                                    ggggAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGTGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGCGAGCGGAACCTg                                                                                     <  4:781383/149‑1 (MQ=255)
                                                    ggggAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTTCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTg                                                                                     <  1:750951/149‑1 (MQ=255)
                                                    ggggAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTTCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTg                                                                                     <  3:336060/149‑1 (MQ=255)
                                                    ggggAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTTCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTg                                                                                     <  3:242316/149‑1 (MQ=255)
                                                      ggAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTTCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGt                                                                                   <  1:262878/149‑1 (MQ=255)
                                                      ggAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTTCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGt                                                                                   >  4:78954/1‑149 (MQ=255)
                                                        aGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTTCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTAt                                                                                 <  3:303688/149‑1 (MQ=255)
                                                          cAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTa                                                                                  <  2:537475/146‑1 (MQ=255)
                                                          cAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTa                                                                                  >  1:537475/1‑146 (MQ=255)
                                                             gTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTTCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTg                                                                                     <  1:415661/140‑1 (MQ=255)
                                                              tctcATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTTCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGt                                                                           <  2:682738/149‑1 (MQ=255)
                                                              tctcATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTTCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGt                                                                           <  4:456011/149‑1 (MQ=255)
                                                              tctcATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGt                                                                           <  1:672162/149‑1 (MQ=255)
                                                                tcATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTTCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGt                                                                           >  3:223834/1‑147 (MQ=255)
                                                                tcATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTTCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGt                                                                           <  4:640758/147‑1 (MQ=255)
                                                                tcATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTTCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGt                                                                           >  3:640758/1‑147 (MQ=255)
                                                                tcATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTTCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGt                                                                           <  4:223834/147‑1 (MQ=255)
                                                                tcATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTTCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGc                                                                         >  2:773094/1‑149 (MQ=255)
                                                                tcATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTTCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGc                                                                         >  4:133729/1‑149 (MQ=255)
                                                                tcATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGc                                                                         >  1:829884/1‑149 (MQ=255)
                                                                tcATCGCTGAAGAGCGCAACGAAGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTTCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGc                                                                         >  2:633984/1‑149 (MQ=255)
                                                                  aTCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTTCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCt                                                                                      >  1:357733/1‑134 (MQ=255)
                                                                  aTCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTTCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCt                                                                                      <  2:357733/134‑1 (MQ=255)
                                                                     gCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTTCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAata                                                                    >  3:283769/1‑149 (MQ=255)
                                                                      cTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTTCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAatat                                                                   >  4:833084/1‑149 (MQ=255)
                                                                       tGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTTCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATg                                                                  <  3:599884/149‑1 (MQ=255)
                                                                        gAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTTCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATgg                                                                 <  2:876099/149‑1 (MQ=255)
                                                                          agagCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTTCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGaa                                                               <  3:833084/149‑1 (MQ=255)
                                                                          agagCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTTCGCGTACGCCAATGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGaa                                                               <  2:79130/149‑1 (MQ=255)
                                                                              cgcAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTa                                                           <  1:58740/149‑1 (MQ=255)
                                                                               gcAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTTCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTAt                                                          <  4:510725/149‑1 (MQ=255)
                                                                                cAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTTCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTAtt                                                         <  4:283769/149‑1 (MQ=255)
                                                                                cAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTTCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAACTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTAtt                                                         <  3:78954/149‑1 (MQ=255)
                                                                                cAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTAtt                                                         <  3:824583/149‑1 (MQ=255)
                                                                                cAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTAtt                                                         <  2:343263/149‑1 (MQ=255)
                                                                                 aaCGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTTCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTg                                                        >  1:453998/1‑149 (MQ=255)
                                                                                 aaCGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTTCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTg                                                        >  2:347928/1‑149 (MQ=255)
                                                                                          ggACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTTCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACaa                                               >  3:385890/1‑149 (MQ=255)
                                                                                          ggACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTTCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACaa                                               >  1:473458/1‑149 (MQ=255)
                                                                                          ggACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTTCGACaa                                               >  3:817741/1‑149 (MQ=255)
                                                                                              gAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTTCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGgt                                                                                                   <  1:369450/93‑1 (MQ=255)
                                                                                              gAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTTCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGgt                                                                                                   >  2:369450/1‑93 (MQ=255)
                                                                                                  aCGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTTCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGt                                                                           <  1:521209/113‑1 (MQ=255)
                                                                                                  aCGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTTCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGt                                                                           >  2:521209/1‑113 (MQ=255)
                                                                                                         ggCATTCGCACTACTTCCGCTGGCGGCGGCCTGTTCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTc                                <  3:455618/149‑1 (MQ=255)
                                                                                                          gCATTCGCACTACTTCCGCTGGCGGCGGCCTGTTCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCa                               >  4:485613/1‑149 (MQ=255)
                                                                                                          gCATTCGCACTACTTCCGCTGGCGGCGGCCTGTTCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCa                               <  4:510953/149‑1 (MQ=255)
                                                                                                           cATTCGCACTACTTCCGCTGGCGGCGGCCTGTTCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCAt                              <  4:47668/149‑1 (MQ=255)
                                                                                                            aTTCGCACTACTTCCGCTGGCGGCGGCCTGTTCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCa                               >  4:557091/1‑147 (MQ=255)
                                                                                                            aTTCGCACTACTTCCGCTGGCGGCGGCCTGTTCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCa                               <  3:557091/147‑1 (MQ=255)
                                                                                                             ttCGCACTACTTCCGCTGGCGGCGGCCTGTTCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGc                            >  4:28311/1‑149 (MQ=255)
                                                                                                             ttCGCACTACTTCCGCTGGCGGCGGCCTGTTCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGc                            <  3:587739/149‑1 (MQ=255)
                                                                                                             ttCGCACTACTTCCGCTGGCGGCGGCCTGTTCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGc                            >  1:271033/1‑149 (MQ=255)
                                                                                                             ttCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGc                            >  4:151880/1‑149 (MQ=255)
                                                                                                                gCACTACTTCCGCTGGCGGCGGCCTGTTCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGa                         >  4:659390/1‑149 (MQ=255)
                                                                                                                  actactTCCGCTGGCGGCGGCCTGTTCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGAAc                       >  1:106150/1‑149 (MQ=255)
                                                                                                                  actactTCCGCTGGCGGCGGCCTGTTCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGAAc                       >  3:753919/1‑149 (MQ=255)
                                                                                                                      ctTCCGCTGGCGGCGGCCTGTTCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGAACAAAg                   <  1:128475/149‑1 (MQ=255)
                                                                                                                       ttCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGAACAAAgc                  >  1:89425/1‑149 (MQ=255)
                                                                                                                         ccGCTGGCGGCGGCCTGTTCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGAACAAAgcgc                <  2:51418/149‑1 (MQ=255)
                                                                                                                         ccGCTGGCGGCGGCCTGTTCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGAACAAAgcgc                <  3:801619/149‑1 (MQ=255)
                                                                                                                          cGCTGGCGGCGGCCTGTTCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCAt                              <  4:548895/134‑1 (MQ=255)
                                                                                                                          cGCTGGCGGCGGCCTGTTCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCAt                              >  3:548895/1‑134 (MQ=255)
                                                                                                                          cGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGAACAAAGcgcg               >  3:401751/1‑149 (MQ=255)
                                                                                                                            cTGGCGGCGGCCTGTTCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGAACAAAGCGCGAt             <  3:21407/149‑1 (MQ=255)
                                                                                                                              ggcggcggcCTGTTCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGAACAAAGCGCGATCa           <  2:674539/149‑1 (MQ=255)
                                                                                                                                cggcggcCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTTCGACAATGCAGCAAACCGTTCATGCCGAACAAAGCGCGATCAGc         <  4:817741/149‑1 (MQ=255)
                                                                                                                                     gcCTGTTCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGAACAAAGCGCGATCAGCCACgt    <  1:520732/149‑2 (MQ=255)
                                                                                                                                     gcCTGTTCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGAACAAAGCGCGATCAGCCACGc    >  3:150185/1‑149 (MQ=255)
                                                                                                                                       cTGTTCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGAACAAAGCGCGATCAGCCACGCCt  >  2:481474/1‑149 (MQ=255)
                                                                                                                                           |                                                                                                                                                
ATCTGCACTGGAACCTATTCTGGCAGACAAGTACTTCCCCGCTTTGTTGACCGGGGAGCAAGTCTCATCGCTGAAGAGCGCAACGGGGCTGGACGAAGACGCGCTGGCATTCGCACTACTTCCGCTGGCGGCGGCCTGTGCGCGTACGCCATTGTCGAATTTTAATGTTGGCGCAATTGCGCGCGGTGTGAGCGGAACCTGGTATTTCGGTGCCAATATGGAATTTATTGGTGCGACAATGCAGCAAACCGTTCATGCCGAACAAAGCGCGATCAGCCACGCCT  >  NC_000913/2231894‑2232177

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