Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,287,420 G→C 5.7% intergenic (‑92/+118) tyrV ← / ← tyrT tRNA‑Tyr/tRNA‑Tyr

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,287,4200GC 5.7% 248.5 / 3.3 70intergenic (‑92/+118)tyrV/tyrTtRNA‑Tyr/tRNA‑Tyr
Reads supporting (aligned to +/- strand):  ref base G (35/31);  new base C (2/2);  total (37/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCG‑AAGTTACCACATCGCTGTGGTGAATTATGGTGGTGGGGGAAGGATTCGAACCTTCGAA  >  NC_000913/1287278‑1287568
                                                                                                                                              |                                                                                                                                                     
gATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCAGTTGCCTGTGGCaa                                                                                                                                                 >  2:665257/1‑149 (MQ=255)
 aTGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAAc                                                                                                                                                <  4:512105/149‑1 (MQ=255)
 aTGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTCGAGCTGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAAc                                                                                                                                                >  2:472881/1‑149 (MQ=255)
       gCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTtctc                                                                                                                                          >  4:160741/1‑149 (MQ=255)
       gCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCACAGGCAACgg                                                                                                                                              <  2:314589/145‑2 (MQ=255)
       gCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCACAGGCAACgg                                                                                                                                              >  1:314589/1‑144 (MQ=255)
           aTTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGc                                                                                                                                      <  2:451073/149‑1 (MQ=255)
            tttACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCt                                                                                                                                     >  4:315864/1‑149 (MQ=255)
                   tCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCtc                                                                                                                              >  1:339569/1‑149 (MQ=255)
                   tCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCtc                                                                                                                              >  2:416491/1‑149 (MQ=255)
                          ccTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCtcgaatcga                                                                                                                       >  4:361506/1‑149 (MQ=255)
                           cTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCtcgaatcgaa                                                                                                                      <  1:537289/149‑1 (MQ=255)
                            tttGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCtcgaa                                                                                                                           >  3:342576/1‑143 (MQ=255)
                            tttGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAAc                                                                                                                     <  3:315864/149‑1 (MQ=255)
                            tttGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAAc                                                                                                                     >  3:680828/1‑149 (MQ=255)
                            tttGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGCTCTCTCGCTTTCGCtcgaa                                                                                                                           <  4:342576/143‑1 (MQ=255)
                                  cGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGt                                                                                                               <  3:361506/149‑1 (MQ=255)
                                       gggAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAg                                                                                                          <  3:160741/149‑1 (MQ=255)
                                       gggAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAg                                                                                                          <  4:66979/149‑1 (MQ=255)
                                        ggAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAgg                                                                                                         >  2:439298/1‑149 (MQ=255)
                                             cccACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTtctc                                                                                                    >  1:14996/1‑149 (MQ=255)
                                                          aTTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGAt                                                                                       <  2:14996/149‑1 (MQ=255)
                                                           ttCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATg                                                                                      >  3:681058/1‑149 (MQ=255)
                                                           ttCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGAGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATg                                                                                      >  3:571094/1‑149 (MQ=255)
                                                            tCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGa                                                                                     <  1:439298/149‑1 (MQ=255)
                                                                 ttttGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCa                                                                                >  2:84143/1‑149 (MQ=255)
                                                                        gTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTc                                                                         >  1:463767/1‑149 (MQ=255)
                                                                        gTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTc                                                                         <  4:680828/149‑1 (MQ=255)
                                                                          aaTGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTcac                                                                       >  1:481073/1‑149 (MQ=255)
                                                                          aaTGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTcac                                                                       >  2:466239/1‑149 (MQ=255)
                                                                               ttGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAAtct                                                                  >  1:450263/1‑149 (MQ=255)
                                                                                    aTGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACttt                                                                          >  3:633539/1‑136 (MQ=255)
                                                                                    aTGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACttt                                                                          <  4:633539/136‑1 (MQ=255)
                                                                                      ggtggtggGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCG‑aa                                                          >  4:101755/1‑149 (MQ=255)
                                                                                         ggtggGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCG‑AAGtt                                                       <  4:681058/149‑1 (MQ=255)
                                                                                         ggtggGAGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCG‑AAGtt                                                       <  4:571094/149‑1 (MQ=255)
                                                                                             ggggAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCG‑AAGTTACca                                                   <  1:665257/149‑1 (MQ=255)
                                                                                               ggAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCG‑AAGTTACcaca                                                 >  4:503775/1‑149 (MQ=255)
                                                                                                    gATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCG‑AAGTTACCACATCGCt                                            >  3:232422/1‑149 (MQ=255)
                                                                                                    gATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCG‑AAGTTACCACATCGCt                                            <  3:503775/149‑1 (MQ=255)
                                                                                                       tattCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCG‑AAGTTACCACATCGCtgtg                                         <  1:416491/149‑1 (MQ=255)
                                                                                                         ttCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCG‑AAGTTACCACATCGCTgtggt                                       >  1:291376/1‑149 (MQ=255)
                                                                                                            gtcgCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCG‑AAGTTACCACATCGCTGTGGTGaa                                    >  2:629146/1‑149 (MQ=255)
                                                                                                            gtcgCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCG‑AAGTTACCACATCGCTGTGGTGaa                                    <  2:339569/149‑1 (MQ=255)
                                                                                                               gcttcgctCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTca                                                                        <  1:131359/111‑1 (MQ=255)
                                                                                                               gcttcgctCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTca                                                                        >  2:131359/1‑111 (MQ=255)
                                                                                                                cttcgctCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCG‑AAGTTACCACATCGCTGTGGTGAATTAt                                <  2:291376/149‑1 (MQ=255)
                                                                                                                 ttcgctCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCG‑AAGTTACCACATCGCTGTGGTGAATTAtg                               <  2:463767/149‑1 (MQ=255)
                                                                                                                   cgctCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCG‑AAGTTACCACATCGCTGTGGTGAATTAtggt                             <  1:84143/149‑1 (MQ=255)
                                                                                                                       cctcACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCG‑AAGTTACCACATCGCTGTGGTGAATTAtggtggtg                         <  2:481073/149‑1 (MQ=255)
                                                                                                                        ctcACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCG‑AAGTTACCACATCGCTGTGGTGAATTAtggtggtgg                        <  2:450263/149‑1 (MQ=255)
                                                                                                                         tcACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGATCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCG‑AAGTTACCACATCGCTGTGGTGAATTATGGTGGTggg                       >  1:181854/1‑149 (MQ=255)
                                                                                                                          cACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCG‑AAGTTACCACATCGCTGTGGTGAATTATGGTGGTgggg                      >  1:429576/1‑149 (MQ=255)
                                                                                                                           aCCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCG‑AAGTTACCACATCGCt                                            <  2:289795/126‑1 (MQ=255)
                                                                                                                           aCCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCG‑AAGTTACCACATCGCt                                            >  1:289795/1‑126 (MQ=255)
                                                                                                                             ccTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCG‑AAGTTACCACATCGCTGTGGTGAATTAtggt                             <  1:305009/139‑1 (MQ=255)
                                                                                                                             ccTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCG‑AAGTTACCACATCGCTGTGGTGAATTAtggt                             >  2:305009/1‑139 (MQ=255)
                                                                                                                             ccTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCG‑AAGTTACCACATCGCTGTGGTGAATTATGGTGGTGGGGGaa                   >  3:326722/1‑149 (MQ=255)
                                                                                                                             ccTTCGGGACGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCG‑AAGTTACCACATCGCTGTGGTGAATTATGGTGGTGGGGGaa                   >  3:636653/1‑149 (MQ=255)
                                                                                                                              cTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCG‑AAGTTACCACATCGCTGTGGTGAATTATGGTGGTGGGGGAAg                  >  3:671747/1‑149 (MQ=255)
                                                                                                                                 gagaacGTTGCCACAGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCG‑AAGTTACCACATCGCTGTGGTGAATTATGGTGGTGGGGGaa                   <  4:600505/140‑1 (MQ=255)
                                                                                                                                 gagaacGTTGCCACAGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCG‑AAGTTACCACATCGCTGTGGCGAATTATGGTGGTGGGGGaa                   >  3:600505/6‑145 (MQ=255)
                                                                                                                                 cGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCG‑AAGTTACCACATCGCTGTGGTGAATTATGGTGGTGGGGGAAGGAt               <  2:429576/149‑1 (MQ=255)
                                                                                                                                   ggCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCATAATTCACCACAGCGATgtggt                                       <  1:312655/124‑1 (MQ=255)
                                                                                                                                     acGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCG‑AAGTTACCACATCGCTGTGGTGAATTATGGTGGTGGGGGAAGGATTCGa           >  4:242893/2‑149 (MQ=255)
                                                                                                                                      cGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCG‑AAGTTACCACATCGCTGTGGTGAATTATGGTGGTGGGGGAAGGATTCGaa          <  4:671747/149‑1 (MQ=255)
                                                                                                                                        ttGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCG‑AAGTTACCACATCGCTGTGGTGAATTATGGTGGTGGGGGAAGGATTCGAAcc        <  3:101755/149‑1 (MQ=255)
                                                                                                                                           ccTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCG‑AAGTTACCACATCGCTGTGGTGAATTATGGTGGTGGGGGAAGGATTCGAACCTTc     <  1:466239/149‑1 (MQ=255)
                                                                                                                                            cTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCG‑AAGTTACCACATCGCTGTGGTGAATTATGGTGGTGGGGGAAGGATTCGAACCTTCg    >  2:144272/1‑149 (MQ=255)
                                                                                                                                            cTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGATCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCG‑AAGTTACCACATCGCTGTGGTGAATTAGGGTGGTGGGGGAAGGATTCGAACCTTCg    >  2:642880/1‑149 (MQ=255)
                                                                                                                                              gtgGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCG‑AAGTTACCACATCGCTGTGGTGAATTATGGTGGTGGGGGAAGGATTCGAACCTTCGaa  <  4:636653/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                     
GATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCGTTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCG‑AAGTTACCACATCGCTGTGGTGAATTATGGTGGTGGGGGAAGGATTCGAACCTTCGAA  >  NC_000913/1287278‑1287568

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