Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,536,225 T→A 18.2% intergenic (+193/+390) yjhR → / ← nanS pseudogene, helicase family, putative frameshift suppressor/9‑O‑acetyl N‑acetylneuraminic acid esterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,536,2250TA18.2% 111.6 / 2.6 44intergenic (+193/+390)yjhR/nanSpseudogene, helicase family, putative frameshift suppressor/9‑O‑acetyl N‑acetylneuraminic acid esterase
Reads supporting (aligned to +/- strand):  ref base T (7/28);  new base A (6/2);  total (14/30)
Fisher's exact test for biased strand distribution p-value = 5.51e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATATATGGTAGTCTGGATCAACGTCAGCTTTAAATGTGATGAGCCGTTCATGATACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTGGCCCAAAACGTTCCAGTAAACTGTGTAGATGTTCGTAGATTTCCGAGCCATAATCGAAAAGGCGCGGATGTCTATTAAGGCATCCGGCTAATCTGTAAGCAATAACA  >  NC_000913/4535999‑4536450
                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                     
aaTATATGGTAGTCTGGATCAACGTCAGCTTTAAATGTGATGAGCCGTTCATGATACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGcc                                                                                                                                                                                                                           <  1:182233/239‑1 (MQ=255)
   atatGGTAGGCTGGATCAACGTCAGCTTTAAATGTGATGAGCCGTTCATGATACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTaa                                                                                                                                                                                                                        <  1:856422/239‑1 (MQ=255)
         tAGTCTGGATCAACGTCAGCTTTAAATGTGATGAGCCGTTCATGATACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGGTTTTTTGGATTAAAGAAAAATGAAAGCataatattatttttttctttttat                                                                                                                                                                                                                    >  2:186957/1‑223 (MQ=255)
            tCTGGATCAACGTCAGCTTTAAATGTGATGAGCCGTTCATGATACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGAGCTTTTGGGATGAAAAAAAAAGCCAACAATAACATTATAAATGCCTTATTTCttttt                                                                                                                                                                                                               >  1:255984/1‑237 (MQ=255)
            tCTGGATCAACGTCAGCTTTAAATGTGATGAGCCGTTCATGATACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGCTTTTTTTGATGAAAAAAAAAACAAGACATAAAAtaataaatttttcttttacatct                                                                                                                                                                                                               >  2:883889/1‑222 (MQ=255)
             cTGGATCAACGTCAGCTTTAAATGTGATGAGCCGTTCATGATACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGaa                                                                                                                                                                                                              <  2:402605/239‑1 (MQ=255)
              tGGATCAACGTCAGCTTTAAATGTGATGAGCCGTTCATGATACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGTTGGTTTTATTGATTAAAGAAAAAGAAAGCAAATATAAAAATAATGGGCAAtttttttttac                                                                                                                                                                                                              >  2:1089127/1‑236 (MQ=255)
                gATCAACGTCAGCTTTAAATGTGATGAGCCGTTCATGATACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGCTTTTTTTGTTGAAAAAAAAAGCAACCAATAATATAAAAATAGCCTAATttttttttacta                                                                                                                                                                                                           >  2:683162/1‑233 (MQ=255)
                   cAACGTCAGCTTTAAATGTGATGAGCCGTTCATGATACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGGTTTTTTTAATGTAAAAAAAATCAAGGAAAAAAATAAAAATAGCCTAATTTTTTTTTAcattat                                                                                                                                                                                                        >  2:76516/1‑239 (MQ=255)
                    aaCGTCAGCTTTAAATGTGATGAGCCGTTCATGATACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCAtt                                                                                                                                                                                                       <  1:683162/239‑1 (MQ=255)
                      cGTCAGCTTTAAATGTGATGAGCCGTTCATGATACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGaa                                                                                                                                                                                                              <  1:455772/230‑1 (MQ=255)
                       cgcAGCTTTAAATGTGATGAGCCGTTCATGATACAGGAGTAGTAAAGGGCATATTTTACCAATGGTAAGGCGGCTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATGAAg                                                                                                                                                                                                   <  2:1232949/237‑1 (MQ=255)
                         cAGCTTTAAATGTGATGAGCCGTTCATGATACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCa                                                                                                                                                                                                 <  1:76516/239‑1 (MQ=255)
                             tttAAATGTGATGAGCCGTTCATGATACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCatat                                                                                                                                                                                              <  1:836220/238‑1 (MQ=255)
                                aaaTGTGATGAGCCGTTCATGATACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGCTTTTTGTTATTTAAAAAAATATAAGAAATAATATAATAATTGGCTAATTTTTTTTTACATCATTaaacccgtttct                                                                                                                                                                                          >  2:455772/1‑229 (MQ=255)
                                         gaGCCGTTCATGATACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTAt                                                                                                                                                                                <  1:254868/239‑1 (MQ=255)
                                                       aCAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGt                                                                                                                                                                 <  1:1096828/239‑1 (MQ=255)
                                                         aGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGGTTTTTTGGATTAAAGAAAAATAAAAAAATAATATAAAAAAAGCCTAATTTTTTTTAACAACATTAAAAATATCATTATCATCTATTTTCCATAATcactc                                                                                                                                                                >  2:294508/1‑234 (MQ=255)
                                                              tagtaAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGGGCTTTTTGGAATGAAAAAAAAAACAAGCAATAAAATAATAAAAGCCTAAATTTTTTTTACAACATTAAAAATATCATAATATTCTATTTTCCAGAATCCCTCGACTta                                                                                                                                                          >  1:435029/1‑239 (MQ=255)
                                                                     gggCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGCCGGTTTTAGGAATAAAAGAAAATAGAAGCAATTATAAAATAATTACCCTCTTTCTTTTAACATCATTCAAACAATTATCCTAGTCTATTTTCTCGAATCCCTCCACTTCTAAttt                                                                                                                                                    >  1:1101114/1‑236 (MQ=255)
                                                                       gCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGTCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGc                                                                                                                                                   <  1:989394/237‑1 (MQ=255)
                                                                            ttttACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTGTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACg                                                                                                                                            <  1:172538/239‑1 (MQ=255)
                                                                             tttACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACgg                                                                                                                                           <  2:781996/239‑1 (MQ=255)
                                                                                     tGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGGTTTTTTGAATGAAAAAAAAAGCAAAAAATAATATAATAATAGCCTAATTTTTTTTAAAATAATTAAAAAAATTATCATACTTCTTTTTCTCGAATACGTCGACCTTTAATCACCACAACGGCATATCac                                                                                                                                   >  1:1085575/1‑238 (MQ=255)
                                                                                        tAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGa                                                                                                                                 <  1:294508/238‑1 (MQ=255)
                                                                                              gACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGtgc                                                                                                                           <  1:738449/238‑1 (MQ=255)
                                                                                                        tttATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCtg                                                                                                                <  2:67344/239‑1 (MQ=255)
                                                                                                                              gATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTGGCCCAAAACGTTCCAGTaaa                                                                                         <  1:622054/239‑1 (MQ=255)
                                                                                                                              gATTATAAATCTTATTGTTTACCTAATCATTTTCCTGTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTGGCCCAAAACGTTCCAGTaaa                                                                                         <  1:362747/239‑1 (MQ=255)
                                                                                                                                           aTTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTTTGAATTAAAAAAAATGCAAGCAATAAAATAAAAATGGGCTAATTTTTTTTTAAATCATTAAGAATATTATCAACCTCTATTTTCCCGAATCCCTCGACTTTTAATAGCCCAAACAGCATATCAAGCAGGGGCTTGGACCGCTGTGGGCCCAAAAAGTCCAGTaaaaaatgtaggtggt                                                                            >  2:794093/1‑226 (MQ=255)
                                                                                                                                                       aTCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTTTGTATGAAAAAAAAAACAAAGGATAAAATAAAAAATGCCTTATTTTTTTTTACATCATTTATCAAATTATCATCTTCTATTTTTTAGAATCCCTCCACTTATAAACTCCCAAAAGGCCTTTTAGGACGGGGCGTCGGACCCCTTGGGCCCAAACGCTCCCACAAACTGTTTAGACGTTCGTTGAttttc                                                                 >  2:435695/1‑236 (MQ=255)
                                                                                                                                                         aaTTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTGGCCCa                                                                                                        <  1:1079794/196‑1 (MQ=255)
                                                                                                                                                             ttCCTTTGAATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTGGCCCAAAACGTTCCAGTAAACTGTGTAGATGTTCGTAGATTTCCGAGCCAt                                                          <  1:794093/239‑1 (MQ=255)
                                                                                                                                                                    gCAAAAGGCTCCTGGGATTAAGGGCGTGGGTTTTTAGGAATAAAAAAAAATGAAAGAAATAATATAATAATGGCCTATTTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGATTTAAAATCGCCAGAACGGCATAGCAGGACGGTGCTGCGAACGCTGGGGCCAAAAACGTCCAAATAAACGGTGTAGATGTTCGGAGATTTCCGAGCCATAATCGaa                                                   <  2:980358/239‑1 (MQ=255)
                                                                                                                                                                           gcgcCTGGGATTAATGGCGGGTGTTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTGGCCCAAAACGTTCCAGTAAACTGTGTAGATGTTCGTAGATTTCCGAGCCATAATCGAAAAGgcgc                                            <  1:570398/239‑1 (MQ=255)
                                                                                                                                                                            cgcTTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATTCAAGCAATAATATAATAATGGCTTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAAAATCCGTCAACTTATAATCGCCAGAACGGCATATCAGGAGGGTGCTGCGAACGCTGTGGCCAAAAACGTTCCAGTAAACTGTGTAGATGTTCGTAGATTTCCGAGCCATAATCGAAAAGgcg                                             <  1:302346/237‑1 (MQ=255)
                                                                                                                                                                             gcCTGGAATTAATGGAGGTTGTTTTAAGAAATAAAAGAAAAGCCAAGAAATAAAATAATAAGGGCTTAATTTCTTTGAACATCATTAAGCATATCAACATCGTCATTTTGCCAGAATCCGTCGCCTTAAAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTGGCCAAAAACGTTCCAGTAAACTGTGTAGATGTTCGTAGTTTTCCGAGCCATAATCGAAAAGGCGCgg                                          <  2:1101114/239‑1 (MQ=255)
                                                                                                                                                                               gtGGGATAAATGGAGGGTGCTTTTAGGAATAAAAAAAAATGCAAGAAAAAATAAAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCGGGACGGTGCTGCGAACGCTGTGGCCAAAAACGTTCCAGTAAACTGTGTAGATGTTCGTAGATTTCCGAGCCATAATCGAAAAGGCGCGGAt                                        <  1:854941/238‑1 (MQ=255)
                                                                                                                                                                                tGGGATTAATGGCGGGTGCTTTTAGGTAGAAAATAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAAAATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGATTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCTAACGCTGTGGCCAAAAACGTTCCAGTAAACTGTGTAGATGTTCGTAGATTTCCGAGCCATAATCGAAAAGGCGCGGATg                                       <  1:240801/239‑1 (MQ=255)
                                                                                                                                                                                                            agaaaaaaaaaAGAAAGCAATAAAAAAAAAAGGGCCTAATTTCTTTAAACATCATAAAGCATACCATCATCGTCAATTTGCAAGAATCCGTCGACTTAAAATCGCAAGAACGGCATATCAGGAGGGTGCTGCGAACGCTGTGGCCAAAAACGTTCCAGAAAACTGTGTAGATTTTCGTAGATTTCCGAGCCATAATCGAAAAGGCGCGGATGTCTATTAAGGCATCCGGCTAATCTGTa            <  2:333381/236‑1 (MQ=255)
                                                                                                                                                                                                             tGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTGGCCCAAAACGTTCCAGTAAACTGTGTAGATGTTCGTAGATTTCCGAGCCATAATCGAAAAGGCGCGGATGTCTATTAAGGCATCCGGCTAATCTGTaa           >  2:203225/1‑239 (MQ=255)
                                                                                                                                                                                                              gtcagatgtgtataagaaacaaatataataATGGCCTAATTTCTTTGACCATCATTAAGCATATCATCAGCGTTTATTTGCCAGAATCCGGCGATTTATAATCGCAAGAACGTAATATGAGGAAGGTGCGGCAAACGCTGTGGCCAAAAACGTTCAAGAAAACTGGGAAGATGTTCGAAGATTTCCGAGCCATAATCGAAAAGGCGCGGATGTCTATTAAGGCATCCGGCTAATCTGTa            <  2:556863/218‑1 (MQ=255)
                                                                                                                                                                                                               aaagaaaagaCAAGCAATAATATAATAATGGCCTAATTTCTTTAAACATCATTAAGCAAACCATCATCGTCTTTTGGCCAGAATCGGTCGCTTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTGGCCCAAAACGTTCCAGTAAACTGTGTAGATGTTGGTAGATTTCCGAGCCATAATCGAAAAGGCGCGGATGTCTATTAAGGCATCCGGCTAATCTGTAAGc         <  1:110066/239‑1 (MQ=255)
                                                                                                                                                                                                               aaaaaaaaTGCAAGAAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCACCATCGTCTATTTGCAATAATCCGTCGACTTATAATCGCCAGAACGCCATAGCAGGACGGTGCTGGGAACGCTGTGGCCAAAAACGTTCCAGTAAACTGTGTAGATGTTCGTAGATTCCCGAGCCATAATCGAAAAGGCGCGGATGTCTATTAAGGCATCCGGCTAATCTGTAAGc         <  2:91214/239‑1 (MQ=255)
                                                                                                                                                                                                                aagaaaATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTGGCCCAAAACGTTCCAGTAAACTGTGTAGATGTTCGTAGATTTCCGAGCCATAATCGAAAAGGCGCGGATGTCTATTAAGGCATCCGGCTAATCTGTAAGCa        >  1:333381/1‑239 (MQ=255)
                                                                                                                                                                                                                       aTGCAAGCAATAATATAATAATGCCCTAATTTCTTTGAACATCATTAATCATATCATCATCGTCTATTTGCCTGAGTCCGTCGACTTATAACCGCCAGAACGGCATATCAGGTCGGTGCTGCGAACGCTGTGGCCCAAAACGTTCAAGTAAACTGTGTATATGTTCTTAGATTCCCGAGCCATAATCGAAAAGGCGCGGATGTCTATTAAGGCATCCGGCTAATCTGTAAGCAATaaca  <  1:439210/239‑1 (MQ=255)
                                                                                                                                                                                                                        aaagagaaagtaataaaataATGGCTTAATTTTTTTGAACATCATTAAGCATACAATCATGTTTATTTGGCAAGACTCCGGCGACTTATAACCGCCAGAACGGCATATCAGGAGGGTGCTGCGAACGCTGTGGCCCAAAACGTCCCAGTAACCTTTGTAGATGTTGGTAGATTCCCGAGCCATAATCGAAAAGGCGCGGATGTCTATTAAGGCATCCGGCTAATCTGTAAGCAATaaca  <  2:193256/229‑1 (MQ=255)
                                                                                                                                                                                                                                    atataataATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTGGCCCAAAACGTTCCAGTaaa                                                                                         >  2:79067/1‑140 (MQ=255)
                                                                                                                                                                                                                                    atataataATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTGGCCCAAAACGTTCCAGTaaa                                                                                         <  1:79067/140‑1 (MQ=255)
                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                     
AATATATGGTAGTCTGGATCAACGTCAGCTTTAAATGTGATGAGCCGTTCATGATACAGGAGTAGTAAAGGGCATATTTTACCAATGTTAAGGCGACTTTCCGTTTTATTGCGAGTCGCCTTTCCTGATTATAAATCTTATTGTTTACCTAATCATTTTCCTTTGCATAATGCGCCTGGGATTAATGGCGGGTGCTTTTAGGAATGAAAGAAAATGCAAGCAATAATATAATAATGGCCTAATTTCTTTGAACATCATTAAGCATATCATCATCGTCTATTTGCCAGAATCCGTCGACTTATAATCGCCAGAACGGCATATCAGGACGGTGCTGCGAACGCTGTGGCCCAAAACGTTCCAGTAAACTGTGTAGATGTTCGTAGATTTCCGAGCCATAATCGAAAAGGCGCGGATGTCTATTAAGGCATCCGGCTAATCTGTAAGCAATAACA  >  NC_000913/4535999‑4536450

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: