Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,963,293:1 +C 5.3% coding (937/1485 nt) gpp ← guanosine pentaphosphatase/exopolyphosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,963,2931.C 5.3% 333.8 / 13.1 75coding (937/1485 nt)gppguanosine pentaphosphatase/exopolyphosphatase
Reads supporting (aligned to +/- strand):  ref base . (33/38);  new base C (2/2);  total (35/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.15e-01

GAAGCTGACAGGCGCTGATGAGCAAATCGCGGCTTATTGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAA  >  NC_000913/3963146‑3963437
                                                                                                                                                    |                                                                                                                                                
gAAGCTGACAGGCGCTGATGAGCAAATCGCGGCTTATTGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑c                                                                                                                                                 <  3:215206/149‑1 (MQ=255)
            cgcTGATGAGCAAATCGCGGCTTATTGCTTCAAGATGCCATTCGTTCTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCCGTAC‑GCTGCgt                                                                                                                                     <  1:202342/149‑2 (MQ=255)
             gcTGATGAGCAAATCGCGGCTTATTGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGaa                                                                                                                                    >  1:575870/1‑149 (MQ=255)
                  tgaGCAAATCGCGGCTTATTGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATcc                                                                                                                               <  2:445264/149‑1 (MQ=255)
                   gaGCAAATCGCGGCTTATTGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCt                                                                                                                              <  2:215290/149‑1 (MQ=255)
                    aGCAAATCGCGGCTTATTGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTg                                                                                                                             >  2:415201/1‑149 (MQ=255)
                      cAAATCGCGGCTTATTGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCt                                                                                                                           >  3:2167/1‑149 (MQ=255)
                              ggCTTATTGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGcc                                                                                                                   <  1:99632/149‑1 (MQ=255)
                               gCTTATTGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCCGTAC‑GCTGCgtaata                                                                                                                                 >  1:631765/1‑129 (MQ=255)
                               gCTTATTGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAGCCGTAC‑GCTGCgtaata                                                                                                                                 <  2:631765/134‑6 (MQ=255)
                                 ttATTGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGa                                                                                                                <  2:496429/149‑1 (MQ=255)
                                   aTTGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATg                                                                                                              >  1:11802/1‑149 (MQ=255)
                                     tGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCa                                                                                                            <  1:545833/149‑1 (MQ=255)
                                        ttCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCTGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACa                                                                                                         >  4:95875/1‑149 (MQ=255)
                                          cAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATg                                                                                                       >  4:621789/1‑149 (MQ=255)
                                           aaGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGc                                                                                                      <  2:11802/149‑1 (MQ=255)
                                             gATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCg                                                                                                    >  4:643808/1‑149 (MQ=255)
                                               tGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTa                                                                                                  <  2:575870/149‑1 (MQ=255)
                                                  cATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGAc                                                                                               <  3:645892/149‑1 (MQ=255)
                                                    ttCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCa                                                                                             >  2:620761/1‑149 (MQ=255)
                                                    ttCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCa                                                                                             >  2:674563/1‑149 (MQ=255)
                                                     tCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAg                                                                                            <  2:485624/149‑1 (MQ=255)
                                                     tCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAg                                                                                            <  3:57621/149‑1 (MQ=255)
                                                        ttttCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAAGATCGATCATAAAGCGGCGTTGGATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGcc                                                                                         >  2:578982/1‑149 (MQ=255)
                                                              aCCTGATCGAAAAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACg                                                                                   >  1:363636/1‑149 (MQ=255)
                                                                 tGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAg                                                                                            >  4:80492/1‑137 (MQ=255)
                                                                 tGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAg                                                                                            <  3:80492/137‑1 (MQ=255)
                                                                   aTCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCGGTAGACCAGGCCTTCACGCAgcg                                                                              >  4:178051/1‑149 (MQ=255)
                                                                   aTCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAgcg                                                                              >  1:396646/1‑149 (MQ=255)
                                                                      gaagaaGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCAc                                                                           >  1:608011/1‑149 (MQ=255)
                                                                      gaagaaGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCAc                                                                           >  4:114670/1‑149 (MQ=255)
                                                                       aagaagTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCAcc                                                                          >  2:360958/1‑149 (MQ=255)
                                                                       aagaagTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCAcc                                                                          <  3:271625/149‑1 (MQ=255)
                                                                                gCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAg                                                                 <  2:188930/149‑1 (MQ=255)
                                                                                gCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAg                                                                 >  2:624140/1‑149 (MQ=255)
                                                                                  cGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGcc                                                                                                                   <  2:474918/97‑1 (MQ=255)
                                                                                  cGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGcc                                                                                                                   >  1:474918/1‑97 (MQ=255)
                                                                                     aaCTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCa                                                            <  1:415201/149‑1 (MQ=255)
                                                                                       cTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATa                                                          <  3:253875/149‑1 (MQ=255)
                                                                                         ttGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATaca                                                        <  1:22374/149‑1 (MQ=255)
                                                                                             cTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACAt                                                                                                        <  1:159434/97‑1 (MQ=255)
                                                                                             cTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACAt                                                                                                        >  2:159434/1‑97 (MQ=255)
                                                                                              tACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGaa                                                   >  3:331138/1‑149 (MQ=255)
                                                                                              tACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGaa                                                   <  1:138087/149‑1 (MQ=255)
                                                                                              tACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGaa                                                   <  1:122460/149‑1 (MQ=255)
                                                                                                cgcgCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAata                                                 <  2:6974/149‑1 (MQ=255)
                                                                                                   gcTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTc                                              <  4:355212/149‑1 (MQ=255)
                                                                                                          cTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGaa                                                   >  4:262584/1‑137 (MQ=255)
                                                                                                          cTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGaa                                                   <  3:262584/137‑1 (MQ=255)
                                                                                                             aTCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTaa                                    >  1:354744/1‑149 (MQ=255)
                                                                                                             aTCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTaa                                    >  2:657878/1‑149 (MQ=255)
                                                                                                                aaTATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTaaaaa                                 <  2:396646/149‑1 (MQ=255)
                                                                                                                   atCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATgg                              <  4:331138/149‑1 (MQ=255)
                                                                                                                   atCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATgg                              <  3:621789/149‑1 (MQ=255)
                                                                                                                    tCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGc                             >  1:590586/1‑149 (MQ=255)
                                                                                                                    tCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGc                             <  2:590586/149‑1 (MQ=255)
                                                                                                                       aTCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGAt                          >  3:637127/1‑149 (MQ=255)
                                                                                                                         cATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCa                        <  4:2167/149‑1 (MQ=255)
                                                                                                                         cATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCa                        <  1:19044/149‑1 (MQ=255)
                                                                                                                          aTAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCGGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAg                       <  3:178051/149‑1 (MQ=255)
                                                                                                                             aaGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATa                                                          <  1:556162/111‑1 (MQ=255)
                                                                                                                             aaGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATa                                                          >  2:556162/1‑111 (MQ=255)
                                                                                                                              aGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATc                   <  1:656042/149‑1 (MQ=255)
                                                                                                                                  gcgTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAAAATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCCTGATCGCCa               <  2:608011/149‑1 (MQ=255)
                                                                                                                                     ttGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGAc            >  1:177490/1‑149 (MQ=255)
                                                                                                                                     ttGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGAc            >  4:615112/1‑149 (MQ=255)
                                                                                                                                     ttGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGAc            <  1:624140/149‑1 (MQ=255)
                                                                                                                                     ttGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATCCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGAc            >  3:127549/1‑149 (MQ=255)
                                                                                                                                         atatTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACt        <  1:657878/149‑1 (MQ=255)
                                                                                                                                          tatTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCCTGATCGCCatt             <  2:595756/143‑3 (MQ=255)
                                                                                                                                          tatTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCCTGATCGCCatt             >  1:595756/1‑141 (MQ=255)
                                                                                                                                          tatTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTc       >  3:541764/1‑149 (MQ=255)
                                                                                                                                           tattcGCAGCCGTAC‑GCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCg      >  4:414780/5‑149 (MQ=255)
                                                                                                                                           atTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCg      >  3:387559/1‑149 (MQ=255)
                                                                                                                                               cGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGaa  <  1:674563/149‑1 (MQ=255)
                                                                                                                                               cGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGaa  <  4:387559/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
GAAGCTGACAGGCGCTGATGAGCAAATCGCGGCTTATTGCTTCAAGATGCCATTCGTTTTCCACCTGATCGAAGAAGTTAGCCGCAACTTTGGCTACGCGCTGTGCCTGATCAATATCGATCATAAAGCGGCGTTGAATATTACGCAG‑CGTACGGCTGCGAATATCCTGCTCGACGGCCAGATGCAACATGCCGTAGACCAGGCCTTCACGCAGCGCACCGCCCGCCAGGGTCATACACTGAATATTCAGTTCGGTAAAAATGGCGATCAGGATCGCCAGACCACTCGGGAA  >  NC_000913/3963146‑3963437

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