Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 219,011 G→A 5.2% P195L (CCG→CTG)  tsaA ← tRNA‑Thr(GGU) m(6)t(6)A37 methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913219,0110GA 5.2% 258.5 / 3.4 77P195L (CCG→CTG) tsaAtRNA‑Thr(GGU) m(6)t(6)A37 methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base G (40/33);  new base A (2/2);  total (42/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

AAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAGCTGCTTTTCGACCTCTGCGGTAAAACTCACCGCCATCTCTGCAGCTGGCGCACTTTGCGCATAGCTGGCAC  >  NC_000913/218868‑219149
                                                                                                                                               |                                                                                                                                          
aaaaGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTcc                                                                                                                                       >  2:1392906/1‑149 (MQ=255)
aaaaGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTcc                                                                                                                                       >  4:6845/1‑149 (MQ=255)
  aaGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTg                                                                                                                                     <  3:6845/149‑1 (MQ=255)
           gCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTAc                                                                                                                            >  2:1096639/1‑149 (MQ=255)
             gTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACtt                                                                                                                          >  2:1065814/1‑149 (MQ=255)
                     aaCGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTa                                                                                                                             >  1:149331/1‑138 (MQ=255)
                     aaCGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTa                                                                                                                             <  2:149331/138‑1 (MQ=255)
                             tCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGt                                                                                                          <  1:231575/149‑1 (MQ=255)
                                    cAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGc                                                                                                   <  1:1392906/149‑1 (MQ=255)
                                         aCTTCAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATa                                                                                              >  1:236180/1‑149 (MQ=255)
                                               aaaCCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTa                                                                                                                             >  4:563220/1‑112 (MQ=255)
                                               aaaCCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTa                                                                                                                             <  3:563220/112‑1 (MQ=255)
                                               aaaCCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCtt                                                                                        <  1:1065814/149‑1 (MQ=255)
                                               aaaCCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCtt                                                                                        <  2:564935/149‑1 (MQ=255)
                                                       gTCGGTGACGCGCCAGCGAACGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAg                                                                                                                               <  3:169338/102‑1 (MQ=255)
                                                       gTCGGTGACGCGCCAGCGAACGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAg                                                                                                                               >  4:169338/1‑102 (MQ=255)
                                                                    cAGCGAACGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAGCTGCtt                                                                   >  4:220225/1‑149 (MQ=255)
                                                                       cGAACGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAGCTGCTTTTc                                                                <  3:220225/149‑1 (MQ=255)
                                                                       cGAACGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAGCTGCTTTTc                                                                <  2:236180/149‑1 (MQ=255)
                                                                            gTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAAc                                                                                                               >  1:291864/1‑97 (MQ=255)
                                                                            gTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAAc                                                                                                               <  2:291864/97‑1 (MQ=255)
                                                                            gTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAGCTGCTTTTc                                                                <  4:893770/144‑1 (MQ=255)
                                                                            gTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAGCTGCTTTTc                                                                >  3:893770/1‑144 (MQ=255)
                                                                             ttAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACtt                                                                                                                          >  2:1112682/1‑85 (MQ=255)
                                                                             ttAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACtt                                                                                                                          <  1:1112682/85‑1 (MQ=255)
                                                                              tAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTtctc                                                                                     <  4:215230/121‑1 (MQ=255)
                                                                              tAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTtctc                                                                                     >  3:215230/1‑121 (MQ=255)
                                                                                 aaTCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATaaa                                                                                                                >  3:251002/1‑91 (MQ=255)
                                                                                 aaTCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATaaa                                                                                                                <  4:251002/91‑1 (MQ=255)
                                                                                 aaTCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGc                                                                                                            >  1:919540/1‑95 (MQ=255)
                                                                                 aaTCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGc                                                                                                            <  2:919540/95‑1 (MQ=255)
                                                                                     aTGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAGCTGCTTTTCGACCTCTGCGGTaa                                                  >  4:1326372/1‑149 (MQ=255)
                                                                                      tGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGc                                                                                                   <  4:1091685/99‑1 (MQ=255)
                                                                                      tGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGc                                                                                                   >  3:1091685/1‑99 (MQ=255)
                                                                                          gCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAGCTGCtt                                                                   >  2:1350067/1‑127 (MQ=255)
                                                                                          gCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAGCTGCtt                                                                   <  1:1350067/127‑1 (MQ=255)
                                                                                             aGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAGCTGCTTTTCGACCTCTGCGGTAAAACTCa                                            >  4:142812/1‑147 (MQ=255)
                                                                                             aGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAGCTGCTTTTCGACCTCTGCGGTAAAACTCa                                            <  3:142812/147‑1 (MQ=255)
                                                                                                  gCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAGCTGCTTTTCGACCTCTGCGGTAAAACTCACCGCCAt                                     <  3:1326372/149‑1 (MQ=255)
                                                                                                    aTAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTtc                                                                                       <  2:507387/97‑1 (MQ=255)
                                                                                                    aTAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTtc                                                                                       >  1:507387/1‑97 (MQ=255)
                                                                                                    aTAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAgctgc                                                                     <  2:733449/115‑1 (MQ=255)
                                                                                                    aTAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAgctgc                                                                     >  1:733449/1‑115 (MQ=255)
                                                                                                    aTAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAGCTGCTTTTCGACCTCTGCGGTAAAACTCACCGCCAtct                                   >  1:1239274/1‑149 (MQ=255)
                                                                                                      aGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAGCTGCTTTTCGACCTCTGCGGTAAAACTCACCGCCATctct                                 >  3:1279437/1‑149 (MQ=255)
                                                                                                        gTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAGCTGCTTTTCGACCTCTGCGGTAAAACTCACCGCCAtct                                   >  4:1469801/1‑145 (MQ=255)
                                                                                                        gTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAGCTGCTTTTCGACCTCTGCGGTAAAACTCACCGCCAtct                                   <  3:1469801/145‑1 (MQ=255)
                                                                                                         ttttGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAGCTGCTTTTCGACCTCTGCGGTAAAACTCACCGCCATCTCTGc                               >  1:1028916/1‑148 (MQ=255)
                                                                                                         ttttGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAGCTGCTTTTCGACCTCTGCGGTAAAACTCACCGCCATCTCTGc                               <  2:1028916/148‑1 (MQ=255)
                                                                                                               ccGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAGCTGCTTTTCGACCTCTGCGGTAAAACTCACCGCCATCTCTGCAGCTGgc                        >  4:453323/1‑149 (MQ=255)
                                                                                                                 gTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAGCTGCTTTTCGACCTCTGCGGTaaaa                                                <  1:1401147/123‑1 (MQ=255)
                                                                                                                 gTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAGCTGCTTTTCGACCTCTGCGGTaaaa                                                >  2:1401147/1‑123 (MQ=255)
                                                                                                                 gTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAGCTGCTTTTCGACCTCTGCGGTAAAACTCACCGCCATCTCTGCAGCTGgcgc                      <  4:1279437/149‑1 (MQ=255)
                                                                                                                      cTCACCTTTACGATAGGCCGGGCGCGGGTCCTACGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAGCTGCTTTTCGACCTCTGCGGTAAAACTCACCGCCATCTCTGCAGCTGGCGCACttt                 >  2:826541/1‑149 (MQ=255)
                                                                                                                       tCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGc                                                                                 <  2:1419977/84‑1 (MQ=255)
                                                                                                                       tCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGc                                                                                 >  1:1419977/1‑84 (MQ=255)
                                                                                                                       tCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAGCTGCTTTTCGACCTCTGCGGTAAAACTCACCGCCATCTCTGCAGCTGGCGCACTTTg                >  4:610853/1‑149 (MQ=255)
                                                                                                                        cACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAGCTGCTTTTCGACCTCTGCGGTAAAACTCACCGCCATCTCTGCAGCTGGCGCACTTTgc               <  3:453323/149‑1 (MQ=255)
                                                                                                                           cTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAGCTGCTTTTCGACCTCTGCGGTAAAACTCACCGCCATCTCTGCAGCTGGCGCACTTTGCGCa            >  3:519811/1‑149 (MQ=255)
                                                                                                                           cTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAGCTGCTTTTCGACCTCTGCGGTAAAACTCACCGCCATCTCTGCAGCTGGCGCACTTTGCGCa            >  1:940699/1‑149 (MQ=255)
                                                                                                                            tttACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAGCTGCTTTTCGACCTCTGCGGTAAAACTCACCGCCATCTCTGCAGCTGGCGCACTTTGCGCAt           >  2:517409/1‑149 (MQ=255)
                                                                                                                              tACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAGCTGCTTTTCGACCTCTGCGGTAAAACTCAc                                           <  3:1167778/115‑1 (MQ=255)
                                                                                                                              tACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAGCTGCTTTTCGACCTCTGCGGTAAAACTCAc                                           >  4:1167778/1‑115 (MQ=255)
                                                                                                                               aCGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAGCTGCTTTTCGACCTCTGCGGTAAAACTCACCGCCATCTCTGCAGCTGGCGCACTTTGCGCATAGc        >  2:251579/1‑149 (MQ=255)
                                                                                                                                 gATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAg                                                                         >  1:939349/1‑82 (MQ=255)
                                                                                                                                 gATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAg                                                                         <  2:939349/82‑1 (MQ=255)
                                                                                                                                 gATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAGCTGCTTTTc                                                                >  3:622587/1‑91 (MQ=255)
                                                                                                                                 gATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAGCTGCTTTTc                                                                <  4:622587/91‑1 (MQ=255)
                                                                                                                                 gATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAGCTGCTTTTCGACCTCTGCGGTAACACTCACCGCCATCTCTGCAGCTGGCGCACTTTGCGCATAGctg      >  3:1411551/1‑149 (MQ=255)
                                                                                                                                  aTAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAGCTGCTTTTCGACCTCTGCGGTAAAACTCACCGCCATCTCTGCAGCTGGCGCACTTTGCGCATAGctgg     >  3:271832/1‑149 (MQ=255)
                                                                                                                                     ggCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAGCTGCTTTTCGACCTCTGCGGTAAAACTCACCGCCATCTCTGCAGCTGGCGCACTTTGCGCATAGctggcac  <  1:251579/149‑1 (MQ=255)
                                                                                                                                        ctggCGCAGGACCCGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAGCTGCTTTTCGAc                                                             <  4:1288311/85‑1 (MQ=255)
                                                                                                                                        ctggCGCAGGACCCGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAGCTGCTTTTCGAc                                                             >  3:1288311/3‑87 (MQ=255)
                                                                                                                                        ctggCGCAGGACCCGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAGCTGCTTTTCGACCTCTGCGGTaaa                                                 <  1:389825/97‑1 (MQ=255)
                                                                                                                                        ctggCGCAGGACCCGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAGCTGCTTTTCGACCTCTGCGGTaaa                                                 >  2:389825/3‑99 (MQ=255)
                                                                                                                                          ggCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAGCTGCTTTTCGACCTCTGCGGTAAAACTCACCGCCAtctc                                  >  3:448990/1‑112 (MQ=255)
                                                                                                                                          ggCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAGCTGCTTTTCGACCTCTGCGGTAAAACTCACCGCCAtctc                                  <  4:448990/112‑1 (MQ=255)
                                                                                                                                             gcgGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAGCTGCTTTTCGACCTCTGCGGTAAAACTCACCGCCATCTCTGCAGCTGgc                        <  3:1067941/119‑1 (MQ=255)
                                                                                                                                             gcgGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAGCTGCTTTTCGACCTCTGCGGTAAAACTCACCGCCATCTCTGCAGCTGgc                        >  4:1067941/1‑119 (MQ=255)
                                                                                                                                               |                                                                                                                                          
AAAAGAGAAGGGCGTGAATTTAACGCGGTTCCAGCGCAAAGACTTCAAAACCTGCGTCGGTGACGCGCCAGCGAACGTTAAAATCATGTAGCCAGACGGCATAGGTTTTGCCCGTTTCCTCACCTTTACGATAGGCCGGGCGCGGGTCCTGCGCCAGTACTTCGCGGATAAACAGCGTTAACTGCGGATAACGCTTCTCCAGCGTCAAAAGCTGCTTTTCGACCTCTGCGGTAAAACTCACCGCCATCTCTGCAGCTGGCGCACTTTGCGCATAGCTGGCAC  >  NC_000913/218868‑219149

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: ‑