Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 111,300 T→C L86P (CTG→CCG)  mutT → dGTP‑preferring nucleoside triphosphate pyrophosphohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913111,3000TC100.0% 316.2 / NA 83L86P (CTG→CCG) mutTdGTP‑preferring nucleoside triphosphate pyrophosphohydrolase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (38/45);  total (38/45)

GGCGGTAAAATTGAAATGGGTGAAACGCCGGAACAGGCGGTGGTGCGTGAACTTCAGGAAGAAGTCGGGATTACCCCCCAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCTGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAGGTCAGATAAGGCG  >  NC_000913/111152‑111446
                                                                                                                                                    |                                                                                                                                                  
ggcggTAAAATTGAAATGGGTGAAACGCCGGAACAGGCGGTGGTGCGTGAACTTCAGGAAGAAGTCGGGATTACCCCCCAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCgg                                                                                                                                                  <  1:397409/151‑1 (MQ=255)
          ttGAAATGGGTGAAACGCCGGAACAGGCGGTGGTGCGTGAACTTCAGGAAGAAGGCGGGATTACCCCCCAACATTTTTCGCTTTTTGAAAAACTGGGATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCt                                                                                                                                         >  2:418974/1‑150 (MQ=255)
          ttGAAATGGGTGAAACGCCGGAACAGGCGGTGGTGCGGGAACTTCAGGAAAAAGTCGGGATTACCCCCCAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAAAACTGTGGTTTTGGTCGGTCGAACGCt                                                                                                                                         >  1:1036659/1‑150 (MQ=255)
          ttGAAATGGGTGAAACGCCGGAACAGGCGGCCGAGCGTGAACTTCAGGGAGAAGGCGGGAGTACCCCCCAACATTTTTCCCAATTTGAAAAACTGGAATATGAATTCCCCGACAGGGATATAACCCTGCGGTTTTAGGCGGTCGaa                                                                                                                                             >  1:434102/1‑146 (MQ=255)
             aaaTGGGTGAAACGCCGGAACAGGCGGTGGTGCGTGGAATTCAGGAAGAAGTCGGGATTACCCCCCAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAAAACTGTGGTTTTGGCCGGGCGAACGCTggg                                                                                                                                      >  2:396410/1‑150 (MQ=255)
             aaaTGGGTGAAACGCCGGAACAGGCGGTGGTGCGTGAACTTCAGGAAGAAGTCGGGATTACCCCCCAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGa                                                                                                                                     <  1:341220/151‑1 (MQ=255)
                 gggTGAAACGCCGGAACAGGCGGTGGTGCGTGAACTTCAGGAAGAAGTCGGGATTACCCCCCAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAggg                                                                                                                                 <  1:1008824/151‑1 (MQ=255)
                  ggTGAAACGCCGGAACAGGCGGTGGTGCGTGAACTTCAGGAAGAAGTCGGGATTACCCCCCAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAgggg                                                                                                                                <  1:1027032/151‑1 (MQ=255)
                   gTGAAACGCCGGAACAGGCGGTGGTGCGTGAACTTCAGGAAGAAGTCGGGATTACCCCCCAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGa                                                                                                                               <  2:750876/151‑1 (MQ=255)
                    tGAAACGCCGGAACAGGCGGTGGTGCGGGAACTTTAGGGAGAAGGCGGGATTACCCCCCAAAATTTTTTGCTATTTGAAAAACTGGAATATGGATTCCCGGACAGGCATATAACACTGTGGTTTTGGGCGGTCGAAAGCTGGGAAgggggg                                                                                                                              >  2:690068/1‑149 (MQ=255)
                       aaCGCCGGAACAGGCGGTGGGGCGTGAACTTTAGGAAGAAGTCGGGATTACCCCCCCACCTTTTTCGCTCTTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAAAGCTGGGAAGGGGAg                                                                                                                              >  1:811182/1‑148 (MQ=255)
                        aCGCCGGAACAGGCGGTGGTGCGTGAACTTCAGGAAGAAGTCGGGATTACCCCCCAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGt                                                                                                                          <  1:220502/151‑1 (MQ=255)
                             ggAACAGGCGGTGGTGCGTGAACTTCAGGAAGAAGTCGGGATTACCCCCCAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGt                                                                                                                     <  1:421359/151‑1 (MQ=255)
                                      ggtggtGCGTGAACTTCAGGAAGAAGTCGGGATTACCCCCCAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAggg                                                                                                            >  2:1165338/1‑151 (MQ=255)
                                       gtggtgCGTGAACTTCAGGAAGAAGTCGGGATTACCCCCCAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGc                                                                                                           <  1:906165/151‑1 (MQ=255)
                                            gCGTGAACTTCAGGAAGAAGTCGGGATTACCCCCCAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAAccc                                                                                                      >  2:398092/1‑151 (MQ=255)
                                            gCGTGAACTTCAGGAAGAAGTCGGGATTACCCCCCAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAAcc                                                                                                       <  1:343637/150‑1 (MQ=255)
                                             cGTGAACTTCAGGAAGAAGTCGGGATTACCCCCCAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCg                                                                                                     >  2:628757/1‑151 (MQ=255)
                                              gTGAACTTCAGGAAGAAGTCGGGATTACCCCCCAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCgg                                                                                                    >  1:266753/1‑151 (MQ=255)
                                                gAACTTCAGGAAGAAGTCGGGATTACCCCCCAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGt                                                                                                   <  1:55947/150‑1 (MQ=255)
                                                  aCTTCAGGAAGAAGTCGGGATTACCCCCCAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTg                                                                                                  >  2:401444/1‑149 (MQ=255)
                                                     tCAGGAAGAAGTCGGGATTACCCCCCAATATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGt                                                                                               <  1:659184/149‑1 (MQ=255)
                                                     tCAGGAAGAAGTCGGGATTACCCCCCAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTg                                                                                              <  2:1144225/150‑1 (MQ=255)
                                                              agTCGGGATTACCCCCCAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTg                                                                                    >  2:215274/1‑151 (MQ=255)
                                                              agTCGGGATTACCCCCCAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTg                                                                                    >  2:635330/1‑151 (MQ=255)
                                                                 cGGGATTACCCCCCAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTggtc                                                                                 >  2:1113115/1‑151 (MQ=255)
                                                                 cGGGATTACCCCCCAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTggtc                                                                                 <  1:949870/151‑1 (MQ=255)
                                                                     aTTACCCCCCAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTggtcggtc                                                                             <  2:1049156/151‑1 (MQ=255)
                                                                       tACCCCCCAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCtt                                                                           >  1:607025/1‑151 (MQ=255)
                                                                        aCCCCCCAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCtt                                                                           >  2:459659/1‑150 (MQ=255)
                                                                        aCCCCCCAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCtt                                                                           <  1:698708/150‑1 (MQ=255)
                                                                        aCCCCCCAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCtt                                                                           <  1:820081/150‑1 (MQ=255)
                                                                        aCCCCCCAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCtt                                                                           >  1:879462/1‑150 (MQ=255)
                                                                         ccccccAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTaa                                                                         <  1:1108176/151‑1 (MQ=255)
                                                                              cAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCg                                                                    >  1:10344/1‑151 (MQ=255)
                                                                               aaCATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCg                                                                    >  2:4308/1‑150 (MQ=255)
                                                                                      ttCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGTAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGAttt                                                             <  2:174002/150‑1 (MQ=255)
                                                                                      ttCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGAtttt                                                            >  2:1139295/1‑151 (MQ=255)
                                                                                        cGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTcc                                                          >  1:861635/1‑151 (MQ=255)
                                                                                        cGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTc                                                           <  2:322188/150‑1 (MQ=255)
                                                                                         gCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCg                                                         >  1:119680/1‑151 (MQ=255)
                                                                                          cTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCgc                                                        >  2:1020063/1‑151 (MQ=255)
                                                                                            aTTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCgcca                                                      <  2:1171028/151‑1 (MQ=255)
                                                                                             tttGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCgcc                                                       >  1:1170547/1‑149 (MQ=255)
                                                                                              ttGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCgccag                                                     <  1:1173368/150‑1 (MQ=255)
                                                                                               tGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCgccagcc                                                   >  1:902020/1‑151 (MQ=255)
                                                                                               tGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCgccagc                                                    <  1:160511/150‑1 (MQ=255)
                                                                                               tGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCgccagc                                                    <  1:1024584/150‑1 (MQ=255)
                                                                                                gAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCgccagcca                                                  <  1:830151/151‑1 (MQ=255)
                                                                                                gAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCgccagcc                                                   >  1:958247/1‑150 (MQ=255)
                                                                                                 aaaaaCTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCgccagcca                                                  >  2:474423/1‑150 (MQ=255)
                                                                                                   aaaCTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATg                                               >  2:425575/1‑151 (MQ=255)
                                                                                                        ggAATATGAATTCCCGGACAGGCATATAACTCTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCg                                          <  1:659261/151‑1 (MQ=255)
                                                                                                        ggAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAAcc                                           <  1:612353/150‑1 (MQ=255)
                                                                                                           atatGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCgg                                         >  1:1207083/1‑149 (MQ=255)
                                                                                                           atatGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGt                                        >  1:1167458/1‑150 (MQ=255)
                                                                                                             atGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAAt                                     >  2:268922/1‑151 (MQ=255)
                                                                                                               gAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAAtt                                    <  2:1065842/150‑1 (MQ=255)
                                                                                                               gAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAAtt                                    <  2:807370/150‑1 (MQ=255)
                                                                                                               gAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTg                                   <  2:940933/151‑1 (MQ=255)
                                                                                                               gAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTg                                   <  2:240196/151‑1 (MQ=255)
                                                                                                                 atccccGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCg                                 <  2:438354/148‑1 (MQ=255)
                                                                                                                 aTTCCCGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGc                                  >  2:1155837/1‑150 (MQ=255)
                                                                                                                    cccGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGaa                               <  2:818579/150‑1 (MQ=255)
                                                                                                                    cccGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGaa                               <  1:1123054/150‑1 (MQ=255)
                                                                                                                    cccGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAg                              <  2:939046/151‑1 (MQ=255)
                                                                                                                     ccGGACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGa                                <  2:1178656/148‑1 (MQ=255)
                                                                                                                        gACAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCtt                           >  2:1166527/1‑150 (MQ=255)
                                                                                                                         aCAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTaa                         <  1:281155/151‑1 (MQ=255)
                                                                                                                         aCAGGCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTa                          >  2:471817/1‑150 (MQ=255)
                                                                                                                            ggCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACg                      <  1:1006753/151‑1 (MQ=255)
                                                                                                                             gCATATAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGt                     >  1:244227/1‑151 (MQ=255)
                                                                                                                               atatAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCt                   <  1:302700/151‑1 (MQ=255)
                                                                                                                               atatAACACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCt                   <  2:924672/151‑1 (MQ=255)
                                                                                                                                   aaCACTGTGGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAg               <  2:524385/151‑1 (MQ=255)
                                                                                                                                        tgtgGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAGGTc            >  1:1108164/1‑149 (MQ=255)
                                                                                                                                         gtgGTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAGGTCAg          <  2:562329/150‑1 (MQ=255)
                                                                                                                                         gtgGTTTTGGCCGGTCGAACGCTGGGAAGGGGAACCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAGGTc            >  1:26528/1‑148 (MQ=255)
                                                                                                                                           ggTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAGGTCAGAt        <  2:206481/150‑1 (MQ=255)
                                                                                                                                           ggTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAGGTCAGAt        >  1:23194/1‑150 (MQ=255)
                                                                                                                                           ggTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAGGTCAGAt        <  1:909150/150‑1 (MQ=255)
                                                                                                                                           ggTTTTGGCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAGGTCAGATa       <  1:124717/151‑1 (MQ=255)
                                                                                                                                                 ggCCGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAGGTCAGATAAGGCg  <  1:215283/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GGCGGTAAAATTGAAATGGGTGAAACGCCGGAACAGGCGGTGGTGCGTGAACTTCAGGAAGAAGTCGGGATTACCCCCCAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCTGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAGGTCAGATAAGGCG  >  NC_000913/111152‑111446

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