Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,797,391 A→C 5.4% L185R (CTG→CGG)  pheT ← phenylalanine tRNA synthetase, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,797,3910AC 5.4% 525.4 / 17.9 149L185R (CTG→CGG) pheTphenylalanine tRNA synthetase, beta subunit
Reads supporting (aligned to +/- strand):  ref base A (77/64);  new base C (2/6);  total (79/70)
Fisher's exact test for biased strand distribution p-value = 1.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.96e-01

Only 100 of 149 total aligned reads displayed.

ATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCAGGTATTCACGGATATCGGTGCCAATCGGCGCATCCGCAGGCAG  >  NC_000913/1797244‑1797536
                                                                                                                                                   |                                                                                                                                                 
aTCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAg                                                                                                                                                  >  2:717765/1‑149 (MQ=255)
   cacaGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTg                                                                                                                                               <  1:127153/149‑1 (MQ=255)
     caGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCTGCCGCTGGt                                                                                                                                             <  4:355321/149‑1 (MQ=255)
      aGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCTGCCGCTGGtt                                                                                                                                            <  2:1242565/149‑1 (MQ=255)
       gCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTc                                                                                                                                           >  1:731704/1‑149 (MQ=255)
          gAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCTGCCGCTGGTTCaac                                                                                                                                        <  4:390716/149‑3 (MQ=255)
          gAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGc                                                                                                                                        <  4:1241046/149‑1 (MQ=255)
           aggtGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAGCGCGGGCAGGCTCCGGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGAGGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCa                                                                                                                                       <  2:512612/146‑1 (MQ=255)
             ttGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACg                                                                                                                                     >  3:1165708/1‑149 (MQ=255)
               ggCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGc                                                                                                                                   >  1:821072/1‑149 (MQ=255)
               ggCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGc                                                                                                                                   >  1:340090/1‑149 (MQ=255)
                gcgcTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCa                                                                                                                                  >  2:519871/1‑149 (MQ=255)
                  gcTTTAACGTTAATGCCTTTTACCACACGGTCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAAc                                                                                                                                <  1:1201773/149‑1 (MQ=255)
                  gcTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAAc                                                                                                                                <  1:519871/149‑1 (MQ=255)
                    tttAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGt                                                                                                                              >  3:551122/1‑149 (MQ=255)
                    tttAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGt                                                                                                                              >  4:1166623/1‑149 (MQ=255)
                         cGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCAcgc                                                                                                                         >  1:70496/1‑149 (MQ=255)
                          gTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCAcgcg                                                                                                                        >  2:818858/1‑149 (MQ=255)
                          gTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCAcgcg                                                                                                                        >  1:418666/1‑149 (MQ=255)
                          gTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCAcgcg                                                                                                                        >  1:1167594/1‑149 (MQ=255)
                           ttAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACgcgc                                                                                                                       >  4:263587/1‑149 (MQ=255)
                           ttAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACgcgc                                                                                                                       >  2:3510/1‑149 (MQ=255)
                              aTGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAac                                                                                                                    >  1:1376690/1‑149 (MQ=255)
                               tGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAaca                                                                                                                   >  1:464741/1‑149 (MQ=255)
                                  cTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCa                                                                                                                >  2:1104956/1‑149 (MQ=255)
                                   ttttACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCaa                                                                                                               <  1:983813/149‑1 (MQ=255)
                                   ttttACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCaa                                                                                                               >  1:992499/1‑149 (MQ=255)
                                   ttttACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCaa                                                                                                               >  3:112954/1‑149 (MQ=255)
                                    tttACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAat                                                                                                              <  1:226269/149‑1 (MQ=255)
                                      tACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAatga                                                                                                            <  1:1462788/149‑1 (MQ=255)
                                      tACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAatga                                                                                                            <  2:841368/149‑1 (MQ=255)
                                       aCCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAatgat                                                                                                           <  4:358978/149‑1 (MQ=255)
                                       aCCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAatgat                                                                                                           >  3:348360/1‑149 (MQ=255)
                                        ccACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAatgatg                                                                                                          <  1:1104956/149‑1 (MQ=255)
                                        ccACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAatgatg                                                                                                          <  2:418666/149‑1 (MQ=255)
                                         gacacGGCCAAGATAACGCGGTCAGGCTTCCGGCGCTTAGACTGTAATCGGCAGAGTGTCGCCGATGGTCGCACCAACCGGAACGATTTAAGGTTTAACCAGCGGCAGCTGGTTCAGAACGGCCACGTCACGAGCAACACCAATGATGc                                                                                                         <  2:229733/148‑1 (MQ=255)
                                         cacaCGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGc                                                                                                         >  1:295708/1‑149 (MQ=255)
                                          acacGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGcc                                                                                                        >  3:132662/1‑149 (MQ=255)
                                          acacGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGcc                                                                                                        <  2:1480101/149‑1 (MQ=255)
                                          acacGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACGGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGCTCAGCACGGCAACGTCACGCGCAACACCAATGATGcc                                                                                                        >  3:951170/1‑149 (MQ=255)
                                             cGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTaa                                                                                                     <  4:80256/149‑1 (MQ=255)
                                              ggCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAg                                                                                                    <  1:676511/149‑1 (MQ=255)
                                                ccAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCa                                                                                                  <  1:295656/149‑1 (MQ=255)
                                                  aaGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTaa                                                                                                     <  3:920865/144‑1 (MQ=255)
                                                  aaGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTaa                                                                                                     >  4:920865/1‑144 (MQ=255)
                                                  aaGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGt                                                                                                >  1:134567/1‑149 (MQ=255)
                                                   aGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTc                                                                                               <  4:968395/149‑1 (MQ=255)
                                                    gATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCg                                                                                              >  4:1407058/1‑149 (MQ=255)
                                                        aCGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCAc                                                                                          >  2:1467256/1‑149 (MQ=255)
                                                        aCGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCAc                                                                                          >  2:927193/1‑149 (MQ=255)
                                                        aCGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCAc                                                                                          >  4:510045/1‑149 (MQ=255)
                                                            gggCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGtt                                                                                      <  3:50903/149‑1 (MQ=255)
                                                             ggCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGttt                                                                                     >  4:602637/1‑149 (MQ=255)
                                                              gCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTg                                                                                    <  1:1002147/149‑1 (MQ=255)
                                                                  gCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGt                                                                                >  2:1405640/1‑149 (MQ=255)
                                                                     tCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCAc                                                                             <  4:410506/149‑1 (MQ=255)
                                                                                gACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGa                                                                  >  3:486737/1‑149 (MQ=255)
                                                                                 aCTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGAt                                                                 <  2:946837/149‑1 (MQ=255)
                                                                                   tGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAAc                                                                                                                                >  2:515729/1‑84 (MQ=255)
                                                                                   tGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAAc                                                                                                                                <  1:515729/84‑1 (MQ=255)
                                                                                   tGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATgg                                                               >  2:708381/1‑149 (MQ=255)
                                                                                     tAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGgtg                                                             <  2:1149256/149‑1 (MQ=255)
                                                                                      aaTCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGgtgt                                                            >  3:1385085/1‑149 (MQ=255)
                                                                                       aTCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGtgtt                                                           >  1:1122353/1‑149 (MQ=255)
                                                                                       aTCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGtgtt                                                           <  1:249044/149‑1 (MQ=255)
                                                                                        tCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGtgttg                                                          >  1:754163/1‑149 (MQ=255)
                                                                                         cGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGtgttgt                                                         <  3:767257/149‑1 (MQ=255)
                                                                                         cGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGtgttgt                                                         >  1:829725/1‑149 (MQ=255)
                                                                                           gCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGtca                                                       >  1:1265067/1‑149 (MQ=255)
                                                                                                 tgtCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGt                                                 >  2:1384724/1‑149 (MQ=255)
                                                                                                   tcgtcgATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGttt                                               >  3:264944/1‑149 (MQ=255)
                                                                                                     gtcgATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCGACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCa                                             >  4:338822/1‑149 (MQ=255)
                                                                                                     gtcgATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCa                                             >  3:663163/1‑149 (MQ=255)
                                                                                                     gtcgATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCa                                             <  1:175733/149‑1 (MQ=255)
                                                                                                      tcgATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCAg                                            <  4:132662/149‑1 (MQ=255)
                                                                                                      tcgATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCAg                                            >  3:896087/1‑149 (MQ=255)
                                                                                                           ggTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCAGGTAtt                                       >  1:802617/1‑149 (MQ=255)
                                                                                                                cACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGAt                                                                 <  2:772269/118‑1 (MQ=255)
                                                                                                                 aCCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGCTGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCAGGTATTCACggt                                 >  1:1307761/1‑148 (MQ=255)
                                                                                                                   cAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCAGGTATTCACGGata                               >  3:488575/1‑149 (MQ=255)
                                                                                                                    aaCCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCAGGTATTCACGGatat                              >  4:143128/1‑149 (MQ=255)
                                                                                                                     aCCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCAGGTATTCACGGATATc                             >  1:461114/1‑149 (MQ=255)
                                                                                                                      ccGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCAGGTATTCACGGATATCg                            <  3:842965/149‑1 (MQ=255)
                                                                                                                       cGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCAGGTATTCACGGATATCgg                           <  2:1394918/149‑1 (MQ=255)
                                                                                                                        ggAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCAGGTATTCACGGATATCGGt                          >  3:370338/1‑149 (MQ=255)
                                                                                                                         gAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCAGGTATTCACGGATATCGGTg                         <  4:1023699/149‑1 (MQ=255)
                                                                                                                                ttCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCAGGTATTCACGGATATCGGTGCCAATCg                  <  1:811469/149‑1 (MQ=255)
                                                                                                                                  ccGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCAGGTATTCACGGATATCGGTGCCAATCGgc                <  3:263587/149‑1 (MQ=255)
                                                                                                                                  ccGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCAGGTATTCACGGATATCGGTGCCAATCGgc                <  2:137628/149‑1 (MQ=255)
                                                                                                                                     gTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCAGGTATTCACGGATATCGGTGCCAATCGGCGCa             >  2:1100356/1‑149 (MQ=255)
                                                                                                                                       tGAACCAGCTGCCGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCAGGTATTCACGGATATCGGTGCCAATCGGCGCATc           >  2:830235/1‑149 (MQ=255)
                                                                                                                                       tGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCAGGTATTCACGGATATCGGTGCCAATCGGCGCATc           <  2:821072/149‑1 (MQ=255)
                                                                                                                                       tGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCAGGTATTCACGGATATCGGTGCCAATCGGCGCATc           <  2:75925/149‑1 (MQ=255)
                                                                                                                                         aaCCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCAGGTATTCACGGATATCGGTGCCAATCGGCGCATCCg         >  3:697383/1‑149 (MQ=255)
                                                                                                                                           ccAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGCTGGTGTTGTCATCAAGTTTCAGGTATTCACGGATATCGGTGCCAATCGGCGCATCCgca       >  1:1360125/1‑149 (MQ=255)
                                                                                                                                           ccAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCAGGTATTCACGGATATCGGTGCCAATCGGCGCATCCgca       >  3:79836/1‑149 (MQ=255)
                                                                                                                                           ccAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCAGGTATTCACGGATATCGGTGCCAATCGGCGCATCCgca       >  2:1126445/1‑149 (MQ=255)
                                                                                                                                            cAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCAGGTATTCACGGATATCGGTGCCAATCGGCGCATCCgcag      >  2:20006/1‑149 (MQ=255)
                                                                                                                                            cAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCAGGTATTCACGGATATCGGTGCCAATCGGCGCATCCgcag      <  2:177108/149‑1 (MQ=255)
                                                                                                                                              gCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCAGGTATTCACGGATATCGGTGCCAATCGGCGCATCCgcaggc    >  4:468948/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
ATCCACAGCGGAGTTGGCGCTTTAACGTTAATGCCTTTTACCACACGGCCAAGATAACGCGGGCAGGCTTCCGGCGCTTCGACTGTAATCGGCAGCGTGTCGTCGATGGTCGCACCAACCGGAACGATTTCCGGTTGAACCAGCGGCAGCTGGTTCAGCACGGCAACGTCACGCGCAACACCAATGATGCCTAAGCAGTCGGCACGGTTTGGCGTCACGCTGATTTCGATGGTGTTGTCATCAAGTTTCAGGTATTCACGGATATCGGTGCCAATCGGCGCATCCGCAGGCAG  >  NC_000913/1797244‑1797536

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