Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 111,305 G→T E88* (GAA→TAA)  mutT → dGTP‑preferring nucleoside triphosphate pyrophosphohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913111,3050GT100.0% 250.4 / NA 63E88* (GAA→TAA) mutTdGTP‑preferring nucleoside triphosphate pyrophosphohydrolase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (32/31);  total (32/31)

AAAATTGAAATGGGTGAAACGCCGGAACAGGCGGTGGTGCGTGAACTTCAGGAAGAAGTCGGGATTACCCCCCAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCTGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAGGTCAGATAAGGC  >  NC_000913/111158‑111445
                                                                                                                                                   |                                                                                                                                            
aaaaTTGAAATGGGTGAAACGCCGGAACAGGCGGTGGTGCGTGAACTTCAGGAAGAAGTCGGGATTACCCCCCAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCTGGTCTaa                                                                                                                                            <  2:75635/150‑1 (MQ=255)
    ttGAAATGGGTGAAACGCCGGAACAGGCGGTGGTGCGTGAACTTCAGGAAGAAGTCGGGATTACCCCCCAACATTTTTCGCTATTTGGAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCTGGTCTAAAGCt                                                                                                                                        >  1:267277/1‑150 (MQ=255)
           gggTGAAACGCCGGAACAGGCGGTGGTGCGTGAACTTCAGGAAGAAGTCGGGATTACCCCCCAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCTGGTCTAACGCTGGGAAgg                                                                                                                                 <  1:27441/150‑1 (MQ=255)
             gTGAAACGCCGGAACAGGCGGTGGTGCGTGAACTTCAGGAAGAAGTCGGGATTACCCCCCAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCTGGTCTAACGCTGGGAAGGGGa                                                                                                                              <  2:302270/151‑1 (MQ=255)
                    gCCGGAACAGGCGGTGGTGCGTGAACTTCAGGAAGAAGTCGGGATTACCCCCCAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCTGGTCTAACGCTGGGAAGGGGAGCCGTgg                                                                                                                       >  2:149903/1‑151 (MQ=255)
                      cGGAACAGGCGGTGGTGCGTGAACTTCAGGAAGAAGGCGGGGTTACCCCCCAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCTGGGCTAACGCTGGGAAGGGGAACCGAgggg                                                                                                                     >  2:825404/1‑151 (MQ=255)
                       ggAACAGGCGGTGGTGCGTGAACTTCAGGAAGAAGTCGGGATTACCCCCCAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCTGGTCTAACGCTGGGAAGGGGAGCCGTgggg                                                                                                                     >  2:289313/1‑150 (MQ=255)
                         aaCAGGCGGTGGTGCGTGAACTTCAGGAAGAAGTCGGGATTACCCCCCAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCTGGTCTAACGCTGGGAAGGGGAGCCGTGGGGTaa                                                                                                                  <  2:820198/151‑1 (MQ=255)
                           cAGGCGGTGGTGCGTGAACTTCAGGAAGAAGTCGGGATTACCCCCCAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCTGGTCTAACGCTGGGAAGGGGAGCCGTGGGGTAaag                                                                                                                <  2:579529/151‑1 (MQ=255)
                                ggtggtGCGTGAACTTCAGGAAGAAGTCGGGATTACCCCCCAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCTGGTCTAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAggg                                                                                                           <  2:176343/151‑1 (MQ=255)
                                 gtggtgCGTGAACTTCAGGAAGAAGTCGGGATTACCCCCCAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCTGGTCTAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGc                                                                                                          >  1:581166/1‑151 (MQ=255)
                                   ggtgCGTGAACTTCAGGAAGAAGTCGGGATTACCCCCCAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCTGGTCTAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCa                                                                                                         <  1:343749/150‑1 (MQ=255)
                                         tGAACTTCAGGAAGAAGTCGGGATTACCCCCCAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCTGGTCTAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCgg                                                                                                   >  1:480024/1‑150 (MQ=255)
                                          gAACTTCAGGAAGAAGTCGGGATTACCCCCCAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCTGGTCTAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGt                                                                                                  <  1:461940/150‑1 (MQ=255)
                                          gAACTTCAGGAAGAAGTCGGGATTACCCCCCAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCTGGTCTAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTg                                                                                                 <  1:721319/151‑1 (MQ=255)
                                          gAACTTCAGGAAGAAGTCGGGATTACCCCCCAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCTGGTCTAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTg                                                                                                 <  1:708698/151‑1 (MQ=255)
                                           aaCTTCAGGAAGAAGTCGGGATTACCCCCCAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCTGGTCTAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTg                                                                                                 >  2:801206/1‑150 (MQ=255)
                                               tCAGGAAGAAGTCGGGATTACCCCCCAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCTGGTCTAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTg                                                                                             <  2:136074/150‑1 (MQ=255)
                                                  ggAAGAAGTCGGGATTACCCCCCAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCTGGTCTAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATg                                                                                         <  1:711495/151‑1 (MQ=255)
                                                      gaagTCGGGATTACCCCCCAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCTGGTCTAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGc                                                                                     >  2:433084/1‑151 (MQ=255)
                                                                ttACCCCCCAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCTGGTCTAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCt                                                                           >  2:243657/1‑151 (MQ=255)
                                                                 tACCCCCCAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCTGGTCTAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCt                                                                           <  1:408417/150‑1 (MQ=255)
                                                                  aCCCCCCAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCTGGTCTAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCtt                                                                          <  1:239990/150‑1 (MQ=255)
                                                                  aCCCCCCAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCTGGTCTAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCtt                                                                          >  2:827883/1‑150 (MQ=255)
                                                                    cccccAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCTGGTCTAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTa                                                                         >  1:569093/1‑149 (MQ=255)
                                                                    cccccAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCTGGTCTAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAAt                                                                       >  1:214463/1‑151 (MQ=255)
                                                                      cccAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCTGGTCTAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGc                                                                     <  1:414962/151‑1 (MQ=255)
                                                                         aaCATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCTGGTCTAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCg                                                                   >  2:564963/1‑150 (MQ=255)
                                                                          aCATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCTGGTCTAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCgat                                                                 >  2:812120/1‑151 (MQ=255)
                                                                          aCATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCTGGTCTAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCgat                                                                 <  2:70208/151‑1 (MQ=255)
                                                                            aTTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCTGGTCTAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCgatg                                                                >  2:721730/1‑150 (MQ=255)
                                                                               tttCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCTGGTCTAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGAttt                                                            >  2:747802/1‑151 (MQ=255)
                                                                                ttCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCTGGTCTAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGAtttt                                                           >  1:257798/1‑151 (MQ=255)
                                                                                ttCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATACAACACTGTGGTTTTGGCTGGTCTAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGAtttt                                                           >  1:387408/1‑151 (MQ=255)
                                                                                  cGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCTGTTCTAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAAGGCAACCCGGTGAGTGGATGCCGCTGGTCGGTCTTAATGCAGATGAtttt                                                           >  1:130315/1‑149 (MQ=255)
                                                                                  cGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCTGGTCTAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTcc                                                         <  1:903935/151‑1 (MQ=255)
                                                                                  cGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCTGGTCTAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTc                                                          <  1:140263/150‑1 (MQ=255)
                                                                                      aTTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCTGGTCTAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCgcca                                                     <  2:444897/151‑1 (MQ=255)
                                                                                      aTTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCTGGTCTAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCgcc                                                      >  2:84403/1‑150 (MQ=255)
                                                                                           aaaaaCTGGAATATGAATTCCCGGAAAGGCATATAACACTGTGGTTTTGGCTGGTCTAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCgccagcca                                                 >  1:615927/1‑150 (MQ=255)
                                                                                               aCTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCTGGTCTAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGa                                             >  1:590083/1‑150 (MQ=255)
                                                                                                  ggAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCTGGTCTAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGTCAATGAAcc                                          <  1:300914/150‑1 (MQ=255)
                                                                                                  ggAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCTGGTCTAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCg                                         <  1:921218/151‑1 (MQ=255)
                                                                                                       atGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCTGGTCTAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAAt                                    >  1:919194/1‑151 (MQ=255)
                                                                                                         gAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCTGGTCTAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTg                                  <  2:872633/151‑1 (MQ=255)
                                                                                                          aaTTCCCGGACAGGCATATAACACTGTGGTTTTGGCTGGTCTAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGc                                 >  2:480329/1‑151 (MQ=255)
                                                                                                              cccGGACAGGCATATAACACTGTGGTTTTGGCTGGTCTAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGa                               <  2:629102/149‑1 (MQ=255)
                                                                                                                cGGACAGGCATATAACACTGTGGTTTTGGCTGGTCTAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCt                           >  2:822296/1‑151 (MQ=255)
                                                                                                                   aCAGGCATATAACACTGTGGTTTTGGCTGGTCTAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTa                         <  1:485740/150‑1 (MQ=255)
                                                                                                                    cAGGCATATAACACTGTGGTTTTGGCTGGTCTAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGGTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTa                         <  2:250543/149‑1 (MQ=255)
                                                                                                                     aGGCATATAACACTGTGGTTTTGGCTGGTCTAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAAc                      <  2:825701/151‑1 (MQ=255)
                                                                                                                        cATATAACACTGTGGTTTTGGCTGGTCTAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTc                   >  2:130187/1‑151 (MQ=255)
                                                                                                                         atatAACACTGTGGTTTTGGCTGGTCTAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCt                  >  1:837955/1‑151 (MQ=255)
                                                                                                                          tataACACTGTGGTTTTGGCTGGTCTAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTg                 >  1:731211/1‑151 (MQ=255)
                                                                                                                                acTGTGGTTTTGGCTGGTCTAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAGGTc           >  2:217691/1‑151 (MQ=255)
                                                                                                                                   gtgGTTTTGGCTGGTCTAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAGGTc           >  2:263139/1‑148 (MQ=255)
                                                                                                                                   gtgGTTTTGGCTGGTCTAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAGGTCa          <  2:205/149‑1 (MQ=255)
                                                                                                                                   gtgGTTTTGGCTGGTCTAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAGGTCAGa        <  2:766152/151‑1 (MQ=255)
                                                                                                                                   gtgGTTTTGGCTGGTCTAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAGGTCAGa        <  2:753751/151‑1 (MQ=255)
                                                                                                                                     ggTTTTGGCTGGTCTAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAGGTCAGAt       >  2:51770/1‑150 (MQ=255)
                                                                                                                                        tttGGCTGGTCTAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAGGTCAGATAAg    <  1:448024/150‑1 (MQ=255)
                                                                                                                                          tGGCTGGTCTAACGCTGGGAAGTGGAGCCGTGGGGTAAAGAATGGCTACCCGGTGAGTGGGTGTCGCTGGTCGGGCTTACTGCAGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAGGTCAGATaa     >  1:764889/1‑147 (MQ=255)
                                                                                                                                          tGGCTGGTCTAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAGGTCAGATAAGGc  <  1:234411/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                            
AAAATTGAAATGGGTGAAACGCCGGAACAGGCGGTGGTGCGTGAACTTCAGGAAGAAGTCGGGATTACCCCCCAACATTTTTCGCTATTTGAAAAACTGGAATATGAATTCCCGGACAGGCATATAACACTGTGGTTTTGGCTGGTCGAACGCTGGGAAGGGGAGCCGTGGGGTAAAGAAGGGCAACCCGGTGAGTGGATGTCGCTGGTCGGTCTTAATGCCGATGATTTTCCGCCAGCCAATGAACCGGTAATTGCGAAGCTTAAACGTCTGTAGGTCAGATAAGGC  >  NC_000913/111158‑111445

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