Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,776,126 C→G 7.7% S185S (TCC→TCG) ‡ mtlR → mannitol operon repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,776,1260CG 7.7% 873.0 / 45.8 247S185S (TCC→TCG) ‡mtlRmannitol operon repressor
Reads supporting (aligned to +/- strand):  ref base C (115/113);  new base G (9/10);  total (124/123)
Fisher's exact test for biased strand distribution p-value = 8.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.22e-01

Only 100 of 261 total aligned reads displayed.

GGCGTTACCACCGCCGCCACAGTTTGAACCAGCAGACTCCAGTTTGTATGCAATGCAAATTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCG  >  NC_000913/3776006‑3776273
                                                                                                                        |                                                                                                                                                   
ggcgTTACCACCGCCGCCACAGTTTGAACCAGCAGACTCCAGTTTGTATGCAATGCAAATTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTcc                                                                                                                                                     >  3:386098/1‑121 (MQ=255)
       ccaccGCCGCCACAGTTTGAACCAGCAGACTCCAGTTTGTATGCAATGCAAATTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCaa                                                                                                                                                   <  2:612449/116‑1 (MQ=255)
         accGCCGCCACAGTTTGAACCAGCAGACTCCAGTTTGTATGCAATGCAAATTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTTGGAAATCAg                                                                                                                                             >  4:817518/1‑120 (MQ=255)
         accGCCGCCACAGTTTGAACCAGCAGACTCCAGTTTGTATGCAATGCAAATTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCa                                                                                                                                              >  2:604712/1‑119 (MQ=255)
         accGCCGCCACAGTTTGAACCAGCAGACTCCAGTTTGTATGCAATGCAAATTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCa                                                                                                                                              <  1:604712/119‑1 (MQ=255)
         accGCCGCCACAGTTTGAACCAGCAGACTCCAGTTTGTATGCAATGCAAATTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTAAGCTGATTTTGGAAATCAg                                                                                                                                             >  4:906732/1‑120 (MQ=255)
             ccgccACAGTTTGAACCAGCAGACTCCAGTTTGTATGCAATGCAAATTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGt                                                                                                             >  4:545786/1‑148 (MQ=255)
                 cacaGTTTGAACCAGCAGACTCCAGTTTGTATGCAATGCAAATTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCaaa                                                                                                                         <  2:804717/132‑1 (MQ=255)
                   caGTTTGAACCAGCAGACTCCAGTTTGTATGCAATGCAAATTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCaaa                                                                                                                         <  4:44983/130‑1 (MQ=255)
                   caGTTTGAACCAGCAGACTCCAGTTTGTATGCAATGCAAATTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTa                                                                                                                     >  1:496316/1‑134 (MQ=255)
                         gAACCAGCAGACTCCAGTTTGTATGCAATGCAAATTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGcc                                                                                                 <  3:184110/148‑1 (MQ=255)
                          aaCCAGCAGACTCCAGTTTGTATGCAATGCAAATTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAATTCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTg                                                                                                          >  2:68175/1‑138 (MQ=255)
                           aCCAGCAGACTCCAGTTTGTATGCAATGCAAATTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTaaaaaa                                                                                                                                    >  4:346347/1‑111 (MQ=255)
                              agcagACTCCAGTTTGTATGCAATGCAAATTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCtt                                                                                                                              >  2:207897/1‑114 (MQ=255)
                                   aCTCCAGTTTGTATGCAATGCAAATTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGCCCCTCTACCTg                                                                                       <  1:245269/148‑1 (MQ=255)
                                      ccAGTTTGTATGCAATGCAAATTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTTGGAAATCAGCtc                                                                                                                                          >  3:922023/1‑93 (MQ=255)
                                      ccAGTTTGTATGCAATGCAAATTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTaaaaaa                                                                                                                                    >  2:494279/1‑100 (MQ=255)
                                      ccAGTTTGTATGCAATGCAAATTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGt                                                                                                             >  3:82846/1‑123 (MQ=255)
                                      ccAGTTTGTATGCAATGCAAATTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTg                                                                                                          >  4:113792/1‑126 (MQ=255)
                                      ccAGTTTGTATGCAATGCAAATTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTg                                                                                                          <  3:113792/126‑1 (MQ=255)
                                           ttGTATGCAATGCAAATTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGcc                                                                                                 <  1:153139/130‑1 (MQ=255)
                                              tATGCAATGCAAATTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCa                                                                                                                                              <  1:683244/82‑1 (MQ=255)
                                               atgcaatgcaAATTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTTGGAAATCAGCt                                                                                                                                           >  1:556197/1‑84 (MQ=255)
                                               atgcaatgcaAATTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCaaa                                                                                                                         >  1:655646/1‑102 (MQ=255)
                                               atgcaatgcaAATTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCaaa                                                                                                                         >  1:283351/1‑102 (MQ=255)
                                               atgcaatgcaAATTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCaaa                                                                                                                         <  2:655646/102‑1 (MQ=255)
                                                 gcaatgcaaATTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATaaa                                                                         >  2:1093902/1‑148 (MQ=255)
                                                  caatgcaaATTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCaaa                                                                                                                         <  1:120589/99‑1 (MQ=255)
                                                    atgcaaATTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCtt                                                                                                                              >  3:847978/1‑92 (MQ=255)
                                                     tgcaaATTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTTGGAAATCAGCtc                                                                                                                                          >  4:1027037/1‑78 (MQ=38)
                                                     tgcaaATTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTaaaaaa                                                                                                                                    >  1:1026399/1‑85 (MQ=255)
                                                      gcaaATTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTg                                                                                                   <  2:516086/117‑1 (MQ=255)
                                                       caaATTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGCCCCTCTACCTGCt                                                                                     <  4:1077564/130‑1 (MQ=255)
                                                        aaaTTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCcgacg                                                                  <  4:830063/148‑1 (MQ=255)
                                                         aaTTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGCCCctc                                                                                             >  1:615955/1‑120 (MQ=255)
                                                         aaTTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGCCCct                                                                                              >  4:1045550/1‑119 (MQ=255)
                                                           ttCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGata                                                               >  2:152758/1‑148 (MQ=255)
                                                           ttCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTACACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGCCCCTCTACCTGCTTCgg                                                                                 >  4:955608/1‑130 (MQ=255)
                                                            tCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCa                                                                                                                 <  4:292188/97‑1 (MQ=255)
                                                            tCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCa                                                                                                                 >  3:292188/1‑97 (MQ=255)
                                                             cAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGataat                                                             <  1:1093902/148‑1 (MQ=255)
                                                                   tATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCa                                                       <  1:984518/148‑1 (MQ=255)
                                                                    aTCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTaaaaaa                                                                                                                                    <  3:49853/70‑1 (MQ=255)
                                                                    aTCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCtt                                                                                                                              >  3:1025919/1‑76 (MQ=255)
                                                                    aTCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCaaa                                                                                                                         <  3:1114835/81‑1 (MQ=255)
                                                                     tCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTc                                                         <  1:1017075/144‑1 (MQ=255)
                                                                      cAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTa                                                                                                                     >  3:103073/1‑83 (MQ=255)
                                                                        aCAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTg                                                                                                          <  1:358025/92‑1 (MQ=255)
                                                                        aCAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCg                                                                              <  2:739900/120‑1 (MQ=255)
                                                                          aGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGa                                                                             >  2:1078478/1‑119 (MQ=255)
                                                                          aGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTg                                                <  3:636852/148‑1 (MQ=255)
                                                                             cTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGc                                                                                                                  >  2:433303/1‑79 (MQ=255)
                                                                               gtgCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGat                                           <  3:544487/148‑1 (MQ=255)
                                                                                  cGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGCCCCTCTAc                                                                                          <  2:314065/98‑1 (MQ=255)
                                                                                  cGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGCCCCTCTAc                                                                                          >  1:314065/1‑98 (MQ=255)
                                                                                   gATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGCCCCTCTACCTGCTTCGGc                                                                                <  3:969791/107‑1 (MQ=255)
                                                                                     tCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTg                                                                                                          <  3:193940/79‑1 (MQ=255)
                                                                                     tCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTg                                                                                                          >  4:193940/1‑79 (MQ=255)
                                                                                      caacaaTGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACgg                                    >  1:436212/1‑148 (MQ=255)
                                                                                        acaaTGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCt                                                                                                                               >  4:896812/1‑55 (MQ=255)
                                                                                        acaaTGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCt                                                                                                                               <  3:708931/55‑1 (MQ=255)
                                                                                        acaaTGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCt                                                                                                                               >  4:708931/1‑55 (MQ=255)
                                                                                        acaaTGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTa                                                                                                      <  3:899375/80‑1 (MQ=255)
                                                                                        acaaTGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGCCCCTCTACCTg                                                                                       >  2:338339/1‑95 (MQ=255)
                                                                                        acaaTGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGCCCCTCTACCTGCtt                                                                                    >  2:567617/1‑98 (MQ=255)
                                                                                        acaaTGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCCTAAAAAAAGCCttt                                                                                                                             >  2:119920/1‑57 (MQ=255)
                                                                                         caaTGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTg                                                                                                          >  4:799518/1‑75 (MQ=255)
                                                                                         caaTGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTg                                                                                                          <  3:799518/75‑1 (MQ=255)
                                                                                         caaTGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGt                                               <  1:881299/134‑1 (MQ=255)
                                                                                         caaTGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGt                                               <  3:337635/134‑1 (MQ=255)
                                                                                         caaTGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGt                                               >  2:881299/1‑134 (MQ=255)
                                                                                           aTGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCaaa                                                                                                                         >  3:739878/1‑58 (MQ=255)
                                                                                           aTGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCt                                                                                                        >  1:849874/1‑75 (MQ=255)
                                                                                           aTGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCt                                                                                                        <  2:849874/75‑1 (MQ=255)
                                                                                           aTGGTCCTTTCACTGACTGAGCTGAGTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGGAAGCAACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGc                               >  2:264295/1‑148 (MQ=255)
                                                                                            tGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTAc                                                                                                     >  1:587481/1‑77 (MQ=255)
                                                                                               tCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGCCCCTCTACCTGc                                                                                      >  4:501866/1‑89 (MQ=255)
                                                                                                   ttCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAg                                                                                                                   >  4:485060/1‑56 (MQ=255)
                                                                                                   ttCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTg                                                                                                          >  3:451567/1‑65 (MQ=255)
                                                                                                   ttCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCcg                                                                     <  3:667032/102‑1 (MQ=255)
                                                                                                   ttCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCg                       >  1:830317/1‑148 (MQ=255)
                                                                                                       ctgactgaGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATaa                                                                          >  4:160023/1‑93 (MQ=255)
                                                                                                         gactgaGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGCCCCTCTACCTGCTTCgg                                                                                 <  1:432286/84‑1 (MQ=255)
                                                                                                           ctgaGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGCCCCTCTACCTGCTTCGGc                                                                                <  4:331510/83‑1 (MQ=255)
                                                                                                             gaGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACgg                                    >  4:234300/1‑125 (MQ=255)
                                                                                                             gaGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACgg                                    <  3:234300/125‑1 (MQ=255)
                                                                                                                 tGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGCCCCTCTACCTGCtt                                                                                    >  1:769134/1‑73 (MQ=255)
                                                                                                                 tGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCcggacgg         >  2:404028/1‑148 (MQ=255)
                                                                                                                  gATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGata                                                               >  4:209575/1‑93 (MQ=255)
                                                                                                                  gATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGata                                                               <  3:209575/93‑1 (MQ=255)
                                                                                                                     ttCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGAtata                                        <  2:507155/113‑1 (MQ=255)
                                                                                                                     ttCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGAtata                                        >  1:507155/1‑113 (MQ=255)
                                                                                                                      tCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCCGGACGGAATcg    <  4:732284/148‑1 (MQ=255)
                                                                                                                       ccAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAg                                                      >  1:777539/1‑97 (MQ=255)
                                                                                                                       ccAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAg                                                      <  2:777539/97‑1 (MQ=255)
                                                                                                                       ccAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGa                                                     >  2:93577/1‑98 (MQ=255)
                                                                                                                       ccAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGa                                                     <  1:93577/98‑1 (MQ=255)
                                                                                                                       ccAAAACCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCcgacg                                                                  <  3:1077172/85‑1 (MQ=255)
                                                                                                                        cAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACt                                                          <  4:656919/92‑1 (MQ=255)
                                                                                                                        cAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCg                       <  1:108554/127‑1 (MQ=255)
                                                                                                                        |                                                                                                                                                   
GGCGTTACCACCGCCGCCACAGTTTGAACCAGCAGACTCCAGTTTGTATGCAATGCAAATTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCCGGACGGAATCGCG  >  NC_000913/3776006‑3776273

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