Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 694,941 G→A L5L (CTC→CTT miaB ← isopentenyl‑adenosine A37 tRNA methylthiolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913694,9410GA91.1% 109.8 / 2.8 45L5L (CTC→CTTmiaBisopentenyl‑adenosine A37 tRNA methylthiolase
Reads supporting (aligned to +/- strand):  ref base G (2/2);  new base A (19/22);  total (21/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.81e-01

CGATGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGGAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCATAGTGTAATGCTTTGCTGCCGTTGTGACCAGTATGAGCGTTATCAGCCCTTAGGGGTAAAAATCCTGTAAACTTAAAGCATATTGCTAACAGGATGATTGACCATGACAAATCAACCAACGGAAATTGCCATTGTCGGCGGAGGAATGGTCGGCGGCGCACTGGCGCTGGGGCTGGCACAGCACGGA  >  NC_000913/694698‑695184
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                   
cGATGCCACGCAGCCACCGACTCCTGTAATCAGGTCTGGATTCTTCTCTGTTAACAGTCTCCAGCGACCCAACTGATAGAAGACTTTTTCCTGAGCCTTCTCGCGGAGTGAGCAGGTGTACAGCAGCAGCACACGCGCTTCTTCTGCCACGTCGGTCAGTTGATAGCCTTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGttttt                                                                                                                                                                                                                                              <  1:172244/251‑1 (MQ=255)
       aCGCAGCCACCGACGCCGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCAt                                                                                                                                                                                                                                       >  1:302263/1‑251 (MQ=255)
           aGCCACCGACGCCGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTGTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCGCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTTAGTTGATAGCCGTGGTTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTATTATGGAGTTTTTTGGTCATCGAc                                                                                                                                                                                                                                   <  2:22025/251‑1 (MQ=255)
           aGCCACCGACGCCGATAATCAGGCCTGGATTCTTCTCTTTTAACCGTTTCCAGCGACCCAACTGATGGACGACTTTTTCCTGAGCCTTCTCGCGGCTTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGAAGATCGGCCATATTCGATGCCTCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCTTCgcc                                                                                                                                                                                                                                   >  2:360104/1‑249 (MQ=255)
              cACCGACGCCGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCTACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGGAGTTTTTTGGTCATCGACTTg                                                                                                                                                                                                                                >  1:63501/1‑251 (MQ=255)
                  gACGCCGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGctct                                                                                                                                                                                                                            <  1:134096/251‑1 (MQ=255)
                     gCCGATAATCGGGTCTGGTTTCTTCGCTTTTAACATTGTCCAGCGCTCCAATGGCTGTAATACGTTTTCCTGATCCTTCCCGCGTTTTGATCAGGTGTTCAGCATCAGCCCATCTGCTTCTTCCGCCACGTCGTTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGTCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGGTTTAATATGAATTTTTTTTGTCATCGACTTCCTCTTGc                                                                                                                                                                                                                         <  2:247316/251‑1 (MQ=255)
                     cccGATAATCAGTTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGc                                                                                                                                                                                                                         <  1:199565/250‑1 (MQ=255)
                      ccGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTAATCGACTTGCTCTTGCg                                                                                                                                                                                                                        >  1:372482/1‑251 (MQ=255)
                             tCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCTCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATATGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAAtagt                                                                                                                                                                                                                 >  2:131638/1‑251 (MQ=255)
                                  tCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGcc                                                                                                                                                                                                             >  1:216713/1‑250 (MQ=255)
                                   cTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGCCTTg                                                                                                                                                                                                                                >  2:125619/1‑230 (MQ=255)
                                   cTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTg                                                                                                                                                                                                                                <  1:125619/230‑1 (MQ=255)
                                       aTTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAAt                                                                                                                                                                                                       <  1:80274/251‑1 (MQ=255)
                                                 tttAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGgcgc                                                                                                                                                                                              >  1:216948/1‑250 (MQ=255)
                                                    aaCAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGGAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCAt                                                                                                                                                                                            >  1:112591/1‑249 (MQ=255)
                                                      caggtcctcgCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGTATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATTAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCATAgtg                                                                                                                                                                                        <  1:356709/242‑1 (MQ=255)
                                                             cAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAAtag                                                                                                                                                                                                                  >  1:201290/1‑218 (MQ=255)
                                                                      aaCTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTGCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTGTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCATAGTGTAATGCTTTGCTGCCg                                                                                                                                                                        <  1:360104/251‑1 (MQ=255)
                                                                         ggATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGTCGCATAGTGTAATGCTTTGCTGCCGTtg                                                                                                                                                                     <  1:287839/250‑1 (MQ=255)
                                                                                        tCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCATCTCATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCAGCGACTTTCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCATAGTGTAATGCTTTGCTGCCGTTGTGACCAGTATGAGc                                                                                                                                                       <  1:23651/250‑1 (MQ=255)
                                                                                         ccTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCATAGTGTAATGCTTTGCTGCCGTTGTGACCAGTATGAGCGt                                                                                                                                                     >  1:156452/1‑251 (MQ=255)
                                                                                            gAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGTGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCAGCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCATAGTGTAATGCTTTGCTGCCGTTGTGACCAGTATGAGCGTTa                                                                                                                                                   <  1:367629/250‑1 (MQ=255)
                                                                                                    ccgcggtgttaGCAGGTGTTCAGCAGCAGCCCAGCCGCTTATTCCGCTACTTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTGTAATATGGAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCATAGTGTAATGCTTTGCTGCCGTTGTGACCAGTATGAGCGTTATCAGCCCtt                                                                                                                                          <  2:63501/241‑1 (MQ=255)
                                                                                                     cgcgGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCATAGTGTAATGCTTTGCTGCCGTTGTGACCAGTATGAGCGTTATCAGCCCtt                                                                                                                                          <  1:152399/250‑1 (MQ=255)
                                                                                                       cgGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCATAGTGTAATGCTTTGCTGCCGTTGTGACCAGTATGAGCGTTATCAGCCCTTAgg                                                                                                                                       <  1:170184/251‑1 (MQ=255)
                                                                                                              ggCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCATAGTGTAATGCTTTGCTGCCGTTGTGCCCAGTATGAGCGTTATCAGCCCTTAGGGGTaaaa                                                                                                                                <  1:189774/250‑1 (MQ=255)
                                                                                                              aGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCATAGTGTAATGCTTTGCTGCCGTTGTGACCAGTATGAGCGTTATCAGCCCTTAGGGGTaaa                                                                                                                                 >  1:200551/1‑250 (MQ=255)
                                                                                                                 aGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCATAGTGTAATGCTTTGCTGCCGTTGTGACCAGTATGAGCGTTATCAGCCCTTAGGGGTAAAAATc                                                                                                                             >  1:303489/1‑251 (MQ=255)
                                                                                                                         agcaacagcAAATACGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGCGGCATCCAGCAGATCGGCCATCTTCGATGAATCGAACTCGTACATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTACACAGGAATGCAGGGCGCATAGTGTAATGCTTTGCTGCCGTTGTGACCAGTATGAGCGTTCACAGCCCTTAGGGGTAAAAAACCTTTaaa                                                                                                                      >  2:56033/1‑250 (MQ=255)
                                                                                                                                   caTCCGCCTCTCCCGCCCCGTCGGTCATGTGATAGCAGTGGGGGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACCTGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCGTGCGAAATAGTAGCCAGGAATGCAGGGCGCATAGTGTAATGCTTTGCTGCCGTTGTGACCAGTATGAGCGTTATCAGCCCTTAGGGGTAAAAATCCTGTAAACTTAAAGCata                                                                                                           <  2:200551/251‑1 (MQ=255)
                                                                                                                                         cttcttCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTAGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCATAGTGTAATGCTATGCTGCCGTTGTGACCAGTATGAGCGTTATCAGCCCTTAGGGGTAAAAATCCTGTAAACTTAAAGCATATTGCt                                                                                                      >  1:379015/1‑250 (MQ=255)
                                                                                                                                           ttttcCGCCACGTCGGTCAGTTGATAGTCGTTGGTGGCCTCCAGCAGAGCGGCCATCTTCGATTAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCATAGTGTAATGCTTTGCTGCCGTTGTGACCAGTATGAGCGTTATCAGCCCTTAGGGGTAAAAATCCTGTAAACTTAAAGCATATTGCt                                                                                                      <  2:331334/246‑1 (MQ=255)
                                                                                                                                               ccGCCACGTCGGGCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCTGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCATAGTGTAATGCTTTGCTGCCGTTGTGACCAGTATGAGCGTTATCAGCCCTTAGGGGTAAAAATCCTGTAAACTTAAAGCATATTGCTAACAg                                                                                                 <  2:192233/249‑1 (MQ=255)
                                                                                                                                                      gTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCATAGTGTAATGCTTTGCTGCCGTTGTGACCAGTATGAGCGTTATCa                                                                                                                                                <  1:173540/195‑1 (MQ=255)
                                                                                                                                                      gTCGGTCAGTTGATAGCCGTGGGTGGCATACAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCATAGTGTAATGCTTTGCTGCCGTTGTGACCAGTATGAGCGTTATCa                                                                                                                                                >  2:173540/1‑195 (MQ=255)
                                                                                                                                                                  aTAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCATAGTGTAATGCTTTGCTGCCGTTGTGACCAGTATGAGCGTTATCAGCCCTTAGGGGTAAAAATCCTGTAAACTTAAAGCATATTGCTACCAGGATGATTGACCATGACAAATc                                                                            >  1:134739/1‑251 (MQ=255)
                                                                                                                                                                             ggtcATCCAGGAGATCGGCCATCTGCTATGAATCGTACTCGTTCATCTTACAGCCCCAGGTGTTAATATGAAGTTTTTTGTGCAGCGACTTGCTCTTGCGAAATAGTAGCCATTAATGCAGGGCGCATAGTGTAATGCTTTGCTGCCGTTGAGACCAGTATTAGCGTTATCAGCCCTTAGGGGTAAAAATCCTGTAAACTTAAAGAATATTTCTAACATGATGATTGACCAAGACAAATCAACCAACGGaa                                                                 <  2:216713/248‑1 (MQ=255)
                                                                                                                                                                                     agcagAGCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCATAGTGTAATGCTTTGCTGCCGTTGTGACCAGTATGAGCGTTATCAGCCCTTAGGGGTAAAAATCCTGTAAACTTAAAGCATATTGCTAACAGGATGATTGACCATGACAAATCAACCAACGGAAATTGCCAt                                                         <  1:28666/251‑1 (MQ=255)
                                                                                                                                                                                                     tCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCATAGTGTAATGCTTTGCTGCCGTTGTGACCAGTATGAGCGTTATCAGCCCTTAGGGGTAAAAATCCTGTAAACTTAAAGCATATTGCTAACAGGATGATTGACCATGACAAATCAACCAACGGAAATTGCCATTGTCGGCGGAGGaa                                           <  1:275472/249‑1 (MQ=255)
                                                                                                                                                                                                                       tCATCTGACAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCATAGTGTAATGCTTTGCTGCCGTTGTGACCAGTATGAGCGTTATCAGCCCTTAGGGGTAAAAATCCTGTAAACTTAAAGCATATTGCTAACAGGATGATTGACCATGACAAATCAACCAACGGAAATTGCCATTGTCGGCGGAGGAATGGTCGGCGGCGCACTGgcg                       >  1:38057/1‑251 (MQ=255)
                                                                                                                                                                                                                              aCAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCATAGTGTAATGCTTTGCTGCCGTTGTGACCAGTATGAGCGTTATCAGCCCTTAGGGGTAAAAATCCTGTAAACTTAAAGCATATTGCTAACAGGATGATTGACCATGACAAATCAACCAACGGAAATTGCCATTGTCGGCggag                                              <  2:174843/221‑1 (MQ=255)
                                                                                                                                                                                                                              aCAGCCCCAGGTTTTAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCATAGTGTAATGCTTTGCTGCCGTTGTGACCAGTATGAGCGTTATCAGCCCTTAGGGGTAAAAATCCTGTAAACTTAAAGCATATTGCTAACAGGATGATTGACCATGACAAATCAACCAACGGAAATTGCCATTGTCGGCggag                                              >  1:174843/1‑221 (MQ=255)
                                                                                                                                                                                                                                  tcccAGGTTTTAATATGAAGTTTTTTGGTCATCGTCTTGCTCTTGCGAAATCGTAGCCAGGATTTCAGGGCGCATAGTGTAATGCTTTGCTGCCGTTGTGACCAGTATGAGCGTTATCAGCCCTTAGGGGTAAAACTCCTGTAAACTTAAAGCATATTGCTAACAGGATGATTGACAATGCCAAATCACCCAACGGAAATTGCAATTGTCGGCGTAAGAATGGTCGGAGGCGCACTGGCGCTGGGGCTGGc            >  2:100593/2‑251 (MQ=255)
                                                                                                                                                                                                                                            tAATATGAAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCATAGTGTAATGCTTTGCTGCCGTTGTGACCAGTATGAGCGTTATCAGCCCTTAGGGGTAAAAATCCTGTAAACTTAAAGCATATTGCTAACAGGATGATTGACCATGACAAATCAACCAACGGAAATTGCCATTGTCGGCGGAGGAATGGTAGGCGGCGCACTGGCGCTGGGGCTGGCACAGCACggg  >  1:175900/1‑250 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                   
CGATGCCACGCAGCCACCGACGCCGATAATCAGGTCTGGATTCTTCTCTTTTAACAGTTTCCAGCGACCCAACTGATGGAAGACTTTTTCCTGAGCCTTCTCGCGGATTGAGCAGGTGTTCAGCAGCAGCACATCCGCTTCTTCCGCCACGTCGGTCAGTTGATAGCCGTGGGTGGCATCCAGCAGATCGGCCATCTTCGATGAATCGTACTCGTTCATCTGACAGCCCCAGGTTTTAATATGGAGTTTTTTGGTCATCGACTTGCTCTTGCGAAATAGTAGCCAGGAATGCAGGGCGCATAGTGTAATGCTTTGCTGCCGTTGTGACCAGTATGAGCGTTATCAGCCCTTAGGGGTAAAAATCCTGTAAACTTAAAGCATATTGCTAACAGGATGATTGACCATGACAAATCAACCAACGGAAATTGCCATTGTCGGCGGAGGAATGGTCGGCGGCGCACTGGCGCTGGGGCTGGCACAGCACGGA  >  NC_000913/694698‑695184

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