Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,844,629 A→G E222E (GAA→GAG adeD → cryptic adenine deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,844,6290AG100.0% 451.4 / NA 114E222E (GAA→GAGadeDcryptic adenine deaminase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (48/66);  total (48/66)

Only 100 of 114 total aligned reads displayed.

CCTGGCGGGACCATCCGCAGGTTACCGGCCTTGCAGAAATGATGGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAACGCGCTGGCACCGTTGAT  >  NC_000913/3844481‑3844775
                                                                                                                                                    |                                                                                                                                                  
   ggCGGGACCATCCGCAGGTTACCGGCCTTGCAGAAATGATGGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAGCtt                                                                                                                                                 <  4:426555/149‑1 (MQ=255)
   ggCGGGACCATCCGCAGGTTACCGGCCTTGCAGAAATGATGGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAGCtt                                                                                                                                                 <  4:275927/149‑1 (MQ=255)
      gggACCATCCGCAGGTTACCGGCCTTGCAGAAATGATGGACTACCCTGGCGTAATTAGCGGGCAGCATGCGCTGATCGATAAACTGGATGCATTTCGCCACCTGACGCTGGCCGGTCACTGCGCGGGTTTGGGTGGTAAAGAGCtt                                                                                                                                                 >  2:456912/1‑146 (MQ=255)
         aCCATCCGCAGGTTACCGGCCTTGCAGAAATGATGGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAGCTTAACGCCt                                                                                                                                          >  2:280116/1‑150 (MQ=255)
            aTCCGCAGGTTACCGGCCTTGCAGAAATGATGGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAGCTTAACGCCTatat                                                                                                                                      <  2:1074133/151‑1 (MQ=255)
              ccGCAGGTTACCGGCCTTGCAGAAATGATGGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAGCTTAACGCCTATAtt                                                                                                                                     >  4:34528/1‑150 (MQ=255)
              ccGCAGGTTACCGGCCTTGCAGAAATGATGGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAGCTTAACGCCTATAtt                                                                                                                                     >  2:187078/1‑150 (MQ=255)
              ccGCAGGTTACCGGCCTTGCAGAAATGATGGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAGCTTAACGCCTATAtt                                                                                                                                     <  1:1065768/150‑1 (MQ=255)
              ccGCAGGTTACCGGCCTTGCAGAAATGATGGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAGCTTAACGCCTATATTa                                                                                                                                    <  1:1032123/151‑1 (MQ=255)
               cGCAGGTTACCGGCCTTGCAGAAATGATGGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAGCTTAACGCCTATATTAc                                                                                                                                   >  2:206418/1‑151 (MQ=255)
                  aGGTTACCGGCCTTGCAGAAATGATGGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAGCTTAACGCCTATATTACTGc                                                                                                                                <  3:218876/151‑1 (MQ=255)
                      tACCGGCCTTGCAGAAATGATGGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAGCTTAACGCCTATATTACTGCGGGt                                                                                                                            >  1:431240/1‑151 (MQ=255)
                        ccGGCCTTGCAGAAATGATGGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAGCTTAACGCCTATATTACTGCGGGTAt                                                                                                                          <  1:12212/151‑1 (MQ=255)
                          ggCCTTGCAGAAATGATGGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAGCTTAACGCCTATATTACTGCGGGTATTg                                                                                                                        <  2:976459/151‑1 (MQ=255)
                            ccTTGCAGAAATGATGGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAATAGCTTAACGCCTATATTACTGCGGGTATTGaa                                                                                                                      <  2:1060728/151‑1 (MQ=255)
                            ccTTGCAGAAATGATGGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAGCTTAACGCCTATATTACTGCGGGTATTGaa                                                                                                                      <  1:414957/151‑1 (MQ=255)
                            ccTTGCAGAAATGATGGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAGCTTAACGCCTATATTACTGCGGGTATTGaa                                                                                                                      <  2:9835/151‑1 (MQ=255)
                               tGCAGAAATGATGGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAGCTTAACGCCTATATTACTGCGGGTATTGaaaa                                                                                                                    <  1:1063429/150‑1 (MQ=255)
                                 cAGAAATGATGGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAGCTTAACGCCTATATTACTGCGGGTATTGAAAACTg                                                                                                                 <  4:257648/151‑1 (MQ=255)
                                   gAAATGATGGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAGCTTAACGCCTATATTACTGCGGGTATTGAAAACTGcc                                                                                                               <  4:425184/151‑1 (MQ=255)
                                    aaaTGATGGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAGCTTAACGCCTATATTACTGCGGGTATTGAAAACTGcc                                                                                                               >  1:1009323/1‑150 (MQ=255)
                                    aaaTGATGGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAGCTTAACGCCTATATTACTGCGGGTATTGAAAACTGcc                                                                                                               >  2:1265995/1‑150 (MQ=255)
                                    aaaTGATGGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAGCTTAACGCCTATATTACTGCGGGTATTGAAAACTGcc                                                                                                               >  4:124632/1‑150 (MQ=255)
                                    aaaTGATGGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACAGTCACTGCCCGGGTTTGGGTGGTAAAGAGCTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCa                                                                                                              <  3:245075/151‑1 (MQ=255)
                                     aaTGATGGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAGCTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCAc                                                                                                             <  1:69481/151‑1 (MQ=255)
                                     aaTGATGGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAGCTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCAc                                                                                                             <  2:1090619/151‑1 (MQ=255)
                                     aaTGATGGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCGGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAGCTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCAc                                                                                                             <  3:19674/151‑1 (MQ=255)
                                      atgatgGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAGCTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACg                                                                                                            <  3:358634/151‑1 (MQ=255)
                                        gatgGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAGCTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGa                                                                                                           >  3:71105/1‑150 (MQ=255)
                                            gACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAGCTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGtt                                                                                                      <  1:307740/151‑1 (MQ=255)
                                              cTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAGCTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTa                                                                                                     <  2:606941/150‑1 (MQ=255)
                                                aCCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAGCTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTAt                                                                                                    <  4:237547/149‑1 (MQ=255)
                                                  ccTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAGCTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGc                                                                                                >  2:827675/1‑151 (MQ=255)
                                                     ggCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAGCTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTgg                                                                                             <  1:5929/151‑1 (MQ=255)
                                                      gCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAGCTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTgg                                                                                             >  3:406103/1‑150 (MQ=255)
                                                      gCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAGCTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGga                                                                                            <  1:351100/151‑1 (MQ=255)
                                                       cGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAGCTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGgaa                                                                                           <  2:427399/151‑1 (MQ=255)
                                                       cGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAGCTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGgaa                                                                                           <  4:272606/151‑1 (MQ=255)
                                                          aaTTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAGCTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGgaaga                                                                                         <  1:1133475/150‑1 (MQ=255)
                                                           aTTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAGCTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGaagaag                                                                                       <  4:225611/151‑1 (MQ=255)
                                                             tAGCGGGCAGAATGCGCTGCTCGCTAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAGCTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAgg                                                                                      >  1:344477/1‑150 (MQ=255)
                                                             tAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAGCTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAgg                                                                                      >  3:62455/1‑150 (MQ=255)
                                                              aGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAGCTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGAc                                                                                    <  4:297006/151‑1 (MQ=255)
                                                                 gggCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAGCTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGcc                                                                                 <  3:286903/151‑1 (MQ=255)
                                                                   gCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAGCTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCgg                                                                               >  1:49349/1‑151 (MQ=255)
                                                                    cAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAGCTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGa                                                                              <  1:92938/151‑1 (MQ=255)
                                                                     aGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAGCTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGa                                                                              >  2:953578/1‑150 (MQ=255)
                                                                       aaTGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAGCTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAAt                                                                           <  1:223309/151‑1 (MQ=255)
                                                                          gcgcTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAGCTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTAc                                                                        >  4:121662/1‑151 (MQ=255)
                                                                          gcgcTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAGCTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTAc                                                                        <  2:425727/151‑1 (MQ=255)
                                                                          gcgcTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAGCTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTAc                                                                        <  1:558399/151‑1 (MQ=255)
                                                                           cgcTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAGCTTACCGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTAcg                                                                       <  2:745447/151‑2 (MQ=255)
                                                                           cgcTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAGCTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTAc                                                                        >  4:438950/1‑150 (MQ=255)
                                                                            gctgctCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAGCTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACa                                                                       >  2:1213803/1‑150 (MQ=255)
                                                                              tgctCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAGCTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACt                                                                    >  1:408116/1‑151 (MQ=255)
                                                                               gctCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAGCTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTc                                                                   <  2:70491/151‑1 (MQ=255)
                                                                                 tCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAGCTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCgg                                                                 >  1:158035/1‑151 (MQ=255)
                                                                                  cGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAGCTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCgg                                                                 <  1:533872/150‑1 (MQ=255)
                                                                                   gATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAGCTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGc                                                                >  1:862489/1‑150 (MQ=255)
                                                                                   gATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAGCTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCa                                                               <  2:60522/151‑1 (MQ=255)
                                                                                    aTAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAGCTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCAt                                                              <  2:611515/151‑1 (MQ=255)
                                                                                     tAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAGCTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATg                                                             >  4:154816/1‑151 (MQ=255)
                                                                                      aaaCTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAGCTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATg                                                             <  2:799470/150‑1 (MQ=255)
                                                                                            gATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAGCTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGtt                                                        <  2:1087957/149‑1 (MQ=255)
                                                                                             aTGCATTTTGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAGCTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTtgat                                                     >  1:282957/1‑151 (MQ=255)
                                                                                             aTGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAGCTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTtga                                                      <  1:180447/150‑1 (MQ=255)
                                                                                              tGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAGCTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTtgatg                                                    >  1:1213584/1‑151 (MQ=255)
                                                                                                  tttCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAGCTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATc                                                 >  2:1045554/1‑150 (MQ=255)
                                                                                                    tCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAGCTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCcgc                                              >  1:364331/1‑151 (MQ=255)
                                                                                                     cGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAGCTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCcgcg                                             <  1:527028/151‑1 (MQ=255)
                                                                                                       ccACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAGCTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGaa                                           <  4:140504/151‑1 (MQ=255)
                                                                                                       ccACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAGCTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGa                                            >  1:159427/1‑150 (MQ=255)
                                                                                                        cACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAGCTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAg                                          <  4:118471/151‑1 (MQ=255)
                                                                                                         aCCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAGCTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAgg                                         >  2:741959/1‑151 (MQ=255)
                                                                                                         aCCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAGCTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAgg                                         <  3:389711/151‑1 (MQ=255)
                                                                                                             gACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAGCTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGt                                       <  2:1163279/149‑1 (MQ=255)
                                                                                                             gACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAGCTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTcc                                     <  1:201449/151‑1 (MQ=255)
                                                                                                             gACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAGCTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTcc                                     >  2:159250/1‑151 (MQ=255)
                                                                                                             gACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAGCTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTc                                      <  1:239531/150‑1 (MQ=255)
                                                                                                              aCGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAGCTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTcc                                     >  1:995272/1‑150 (MQ=255)
                                                                                                                    gACGGTCACTGCCCGGGTTTGGGTGGTAAAGAGCTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGccc                              <  2:1261204/151‑1 (MQ=255)
                                                                                                                       ggTCACTGCCCGGGTTTGGGTGGTAAAGAGCTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAGCTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGc                            <  2:514279/150‑1 (MQ=255)
                                                                                                                       ggTCACTGCCCGGGTTTGGGTGGTAAAGAGCTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGc                            >  3:210698/1‑150 (MQ=255)
                                                                                                                            cTGCCCGGGTTTGGGTGGTAAAGAGCTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAAtct                       <  2:306405/150‑1 (MQ=255)
                                                                                                                                ccGGGTTTGGGTGGTAAAGAGCTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAAcg                  <  1:300039/151‑1 (MQ=255)
                                                                                                                                 tgggTTTGGGTGGTAAAGAGCTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAAcgc                 >  2:858312/2‑151 (MQ=255)
                                                                                                                                   ggTTTGGGTGGTAAAGAGCTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAACgcgc               >  4:61300/1‑151 (MQ=255)
                                                                                                                                   ggTTTGGGTGGTAAAGAGCTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAACgcgc               >  2:959100/1‑151 (MQ=255)
                                                                                                                                    gTTTGGGTGGTAAAGAGCTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAACGCGCt              <  2:859900/151‑1 (MQ=255)
                                                                                                                                      ttGGGTGGTAAAGAGCTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAACGCGCTgg            >  3:317927/1‑151 (MQ=255)
                                                                                                                                      ttGGGTGGTAAAGAGCTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAACGCGCTgg            >  1:994267/1‑151 (MQ=255)
                                                                                                                                      ttGGGTGGTAAAGAGCTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACACGGCATGTCGTTGATGATCCGCGAAGGGTCCGCGGCCCGCAATCTCAACGCGCTg             <  3:378551/150‑1 (MQ=255)
                                                                                                                                       tGGGTGGTAAAGAGCTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAACGCGCTGGc           >  1:944959/1‑151 (MQ=255)
                                                                                                                                          gtggtAAAGAGCTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAACGCGCTGGCAcc        <  1:182865/151‑1 (MQ=255)
                                                                                                                                             gtAAAGAGCTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAACGCGCTGGCACCGtt     <  1:1012170/151‑1 (MQ=255)
                                                                                                                                              tAAAGAGCTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAACGCGCTGGCACCGTTg    >  4:299800/1‑151 (MQ=255)
                                                                                                                                               aaaGAGCTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAACGCGCTGGCACCGTTg    <  2:492479/150‑1 (MQ=255)
                                                                                                                                               aaaGAGCTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAACGCGCTGGCACCGTTg    >  4:320717/1‑150 (MQ=255)
                                                                                                                                               aaaGAGCTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAACGCGCTGGCACCGTTGa   <  2:1030865/151‑1 (MQ=255)
                                                                                                                                                aaGAGCTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAACGCGCTGGCACCGTTGAt  <  4:338172/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CCTGGCGGGACCATCCGCAGGTTACCGGCCTTGCAGAAATGATGGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAACGCGCTGGCACCGTTGAT  >  NC_000913/3844481‑3844775

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