Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,776,125 C→T 7.6% S185F (TCC→TTC) ‡ mtlR → mannitol operon repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,776,1250CT 7.6% 797.7 / 37.8 247S185F (TCC→TTC) ‡mtlRmannitol operon repressor
Reads supporting (aligned to +/- strand):  ref base C (115/113);  new base T (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 = 7.63e-01

Only 100 of 255 total aligned reads displayed.

GGCGTTACCACCGCCGCCACAGTTTGAACCAGCAGACTCCAGTTTGTATGCAATGCAAATTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCCGGACGGAATCG  >  NC_000913/3776006‑3776271
                                                                                                                       |                                                                                                                                                   
       ccaccGCCGCCACAGTTTGAACCAGCAGACTCCAGTTTGTATGCAATGCAAATTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCaa                                                                                                                                                  >  1:612449/1‑116 (MQ=255)
         accGCCGCCACAGTTTGAACCAGCAGACTCCAGTTTGTATGCAATGCAAATTCAGCGCTATCAACAGGCTGTGCTATCAACAATGGTCCTTTCACTGACTAAGCTGATTTTGGAAATCAg                                                                                                                                            <  3:906732/120‑1 (MQ=255)
         accGCCGCCACAGTTTGAACCAGCAGACTCCAGTTTGTATGCAATGCAAATTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTTGGAAATCAGCtc                                                                                                                                         >  2:475639/1‑122 (MQ=255)
         accGCCGCCACAGTTTGAACCAGCAGACTCCAGTTTGTATGCAATGCAAATTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCa                                                                                                                                             <  1:604712/119‑1 (MQ=255)
         accGCCGCCACAGTTTGAACCAGCAGACTCCAGTTTGTATGCAATGCAAATTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAg                                                                                                                                            >  3:263640/1‑120 (MQ=255)
         accGCCGCCACAGTTTGAACCAGCAGACTCCAGTTTGTATGCAATGCAAATTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTAAGCTGATTTTGGAAATCAg                                                                                                                                            >  4:906732/1‑120 (MQ=255)
             ccgccACAGTTTGAACCAGCAGACTCCAGTTTGTATGCAATGCAAATTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGt                                                                                                           >  4:545786/1‑148 (MQ=255)
               gccACAGTTTGAACCAGCAGACTCCAGTTTGTATGCAATGCAAATTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCttt                                                                                                                            <  2:609025/130‑1 (MQ=255)
               gccACAGTTTGAACCAGCAGACTCCAGTTTGTATGCAATGCAAATTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCt                                                                                                         <  3:972014/148‑1 (MQ=255)
                 cacaGTTTGAACCAGCAGACTCCAGTTTGTATGCAATGCAAATTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCaaa                                                                                                                        >  1:804717/1‑132 (MQ=255)
                   caGTTTGAACCAGCAGACTCCAGTTTGTATGCAATGCAAATTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTa                                                                                                                    >  1:496316/1‑134 (MQ=255)
                       ttGAACCAGCAGACTCCAGTTTGTATGCAATGCAAATTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTg                                                                                                 >  4:184110/1‑148 (MQ=255)
                         gAACCAGCAGACTCCAGTTTGTATGCAATGCAAATTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGcc                                                                                               <  3:184110/148‑1 (MQ=255)
                          aaCCAGCAGACTCCAGTTTGTATGCAATGCAAATTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAATTCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTg                                                                                                        <  1:68175/138‑1 (MQ=255)
                           aCCAGCAGACTCCAGTTTGTATGCAATGCAAATTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTaaaaaa                                                                                                                                   >  4:346347/1‑111 (MQ=255)
                            ccAGCAGACTCCAGTTTGTATGCAATGCAAATTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTTGGAAATCAg                                                                                                                                            <  1:1070978/101‑1 (MQ=255)
                              agcagACTCCAGTTTGTATGCAATGCAAATTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCtt                                                                                                                             <  1:207897/114‑1 (MQ=255)
                              agcagACTCCAGTTTGTATGCAATGCAAATTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCtt                                                                                                                             >  2:207897/1‑114 (MQ=255)
                                      ccAGTTTGTATGCAATGCAAATTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTTGGAAATCAGCtc                                                                                                                                         >  3:922023/1‑93 (MQ=255)
                                      ccAGTTTGTATGCAATGCAAATTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTg                                                                                                        <  2:977070/126‑1 (MQ=255)
                                      ccAGTTTGTATGCAATGCAAATTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTg                                                                                                        >  4:113792/1‑126 (MQ=255)
                                      ccAGTTTGTATGCAATGCAAATTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTg                                                                                                        >  1:977070/1‑126 (MQ=255)
                                        aGTTTGTATGCAATGCAAATTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTaaaaaa                                                                                                                                   >  2:91385/1‑98 (MQ=255)
                                         gTTTGTATGCAATGCAAATTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTTGGAAATCAg                                                                                                                                            <  2:671117/88‑1 (MQ=255)
                                          tttGTATGCAATGCAAATTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTa                                                                                                                    >  2:130841/1‑111 (MQ=255)
                                               atgcaatgcaAATTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTTGGAAATCAGCt                                                                                                                                          >  1:556197/1‑84 (MQ=255)
                                               atgcaatgcaAATTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTTGGAAATCAGCt                                                                                                                                          >  3:329849/1‑84 (MQ=255)
                                               atgcaatgcaAATTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCaaa                                                                                                                        >  1:655646/1‑102 (MQ=255)
                                                 gcaatgcaaATTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATaaa                                                                       >  1:846750/1‑148 (MQ=255)
                                                  caatgcaaATTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAg                                                                      >  3:830063/1‑148 (MQ=255)
                                                   aatgcaaATTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCg                                                                            >  2:648875/1‑141 (MQ=255)
                                                     tgcaaATTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTTGGAAATCAGCtc                                                                                                                                         <  3:1027037/79‑2 (MQ=38)
                                                     tgcaaATTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTTGGAAATCAGCtc                                                                                                                                         >  4:1027037/1‑78 (MQ=38)
                                                     tgcaaATTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTaaaaaa                                                                                                                                   <  2:1026399/85‑1 (MQ=255)
                                                     tgcaaATTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGccc                                                                                              <  4:401361/121‑1 (MQ=255)
                                                     tgcaaATTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGccc                                                                                              >  3:401361/1‑121 (MQ=255)
                                                      gcaaATTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCagacgt                                                                                                           <  2:306354/108‑6 (MQ=255)
                                                      gcaaATTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCagacgt                                                                                                           >  1:306354/1‑103 (MQ=255)
                                                       caaATTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAATGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGCCCCTCTACCTGCt                                                                                   >  3:1077564/1‑130 (MQ=255)
                                                       caaATTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAGACGTTGCTTAc                                                                                                   >  2:697638/1‑114 (MQ=255)
                                                       caaATTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTc                                                                                                          <  1:856436/107‑1 (MQ=255)
                                                       caaATTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTc                                                                                                          >  2:856436/1‑107 (MQ=255)
                                                       caaATTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGCCCCTCTACCTGCt                                                                                   <  4:1077564/130‑1 (MQ=255)
                                                       caaATTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATa                                                                         >  3:72166/1‑140 (MQ=255)
                                                       caaATTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCcgac                                                                 >  4:1053966/1‑148 (MQ=255)
                                                        aaaTTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCcgacg                                                                <  4:830063/148‑1 (MQ=255)
                                                         aaTTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGCCCct                                                                                            <  3:1045550/119‑1 (MQ=255)
                                                           ttCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTACACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGCCCCTCTACCTGCTTCgg                                                                               >  4:955608/1‑130 (MQ=255)
                                                            tCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCa                                                                                                                <  4:292188/97‑1 (MQ=255)
                                                             cAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGataat                                                           <  1:1093902/148‑1 (MQ=255)
                                                                    aTCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTaaaaaa                                                                                                                                   >  2:833412/1‑70 (MQ=255)
                                                                    aTCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTaaaaaa                                                                                                                                   >  4:49853/1‑70 (MQ=255)
                                                                    aTCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTaaaaaa                                                                                                                                   <  3:49853/70‑1 (MQ=255)
                                                                     tCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTc                                                       <  1:1017075/144‑1 (MQ=255)
                                                                      cAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGCCCCTCTACCTGc                                                                                    >  1:37101/1‑114 (MQ=255)
                                                                      cAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGCCCCTCTACCTGc                                                                                    <  2:37101/114‑1 (MQ=255)
                                                                        aCAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGTTTAAAAAAAGCCTTTCaaa                                                                                                                        <  2:59343/77‑1 (MQ=255)
                                                                        aCAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCg                                                                            >  1:739900/1‑120 (MQ=255)
                                                                        aCAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCg                                                                            >  3:808959/1‑120 (MQ=255)
                                                                         cAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCttt                                                                                                                            >  2:794604/1‑72 (MQ=255)
                                                                         cAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGCCCctc                                                                                           >  3:142489/1‑104 (MQ=255)
                                                                         cAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGtt                                               >  3:85626/1‑148 (MQ=255)
                                                                          aGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTg                                              <  3:636852/148‑1 (MQ=255)
                                                                             cTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGc                                                                                                                 >  2:433303/1‑79 (MQ=255)
                                                                              tgtgCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAg                                                                                                                  <  2:1087967/77‑1 (MQ=255)
                                                                              tgtgCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGa                                          <  1:152758/148‑1 (MQ=255)
                                                                               gtgCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGat                                         <  3:544487/148‑1 (MQ=255)
                                                                                     tCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTg                                                                                                        >  4:193940/1‑79 (MQ=255)
                                                                                      caacaaTGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACgg                                  >  1:436212/1‑148 (MQ=255)
                                                                                        acaaTGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCttt                                                                                                                            <  1:119920/57‑1 (MQ=255)
                                                                                        acaaTGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCt                                                                                                                              >  4:896812/1‑55 (MQ=255)
                                                                                        acaaTGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCt                                                                                                                              >  4:708931/1‑55 (MQ=255)
                                                                                        acaaTGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGCCCCTCTACCTg                                                                                     >  2:338339/1‑95 (MQ=255)
                                                                                        acaaTGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGCCCCTCTACCTg                                                                                     <  1:338339/95‑1 (MQ=255)
                                                                                        acaaTGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGCCCCTCTACCTGCtt                                                                                  <  1:567617/98‑1 (MQ=255)
                                                                                        acaaTGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCCTAAAAAAAGCCttt                                                                                                                            >  2:119920/1‑57 (MQ=255)
                                                                                         caaTGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTg                                                                                                        >  4:799518/1‑75 (MQ=255)
                                                                                         caaTGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGt                                             <  3:337635/134‑1 (MQ=255)
                                                                                          aaTGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACg                              <  1:9176/148‑1 (MQ=255)
                                                                                           aTGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAATCCTTTCaaa                                                                                                                        <  4:739878/58‑1 (MQ=255)
                                                                                           aTGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCaaa                                                                                                                        >  3:739878/1‑58 (MQ=255)
                                                                                           aTGGTCCTTTCACTGACTGAGCTGAGTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGGAAGCAACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGc                             >  2:264295/1‑148 (MQ=255)
                                                                                               tCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTg                                                                                                        >  2:558896/1‑69 (MQ=255)
                                                                                                  tttCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCt                                                                                                      >  2:497751/1‑68 (MQ=255)
                                                                                                   ttCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAg                                                                                                                  >  4:485060/1‑56 (MQ=255)
                                                                                                   ttCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTg                                                                                                        >  3:451567/1‑65 (MQ=255)
                                                                                                   ttCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCcg                                                                   >  4:667032/1‑102 (MQ=255)
                                                                                                   ttCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCg                     <  1:264295/148‑1 (MQ=255)
                                                                                                    tCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCcg                                                                   >  2:510756/1‑101 (MQ=255)
                                                                                                       ctgactgaGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATc                      <  3:203779/143‑1 (MQ=255)
                                                                                                           ctgaGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGCCCCTCTACCTGCTTCGGc                                                                              <  4:331510/83‑1 (MQ=255)
                                                                                                             gaGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACgg                                  >  4:234300/1‑125 (MQ=255)
                                                                                                             gaGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCcgg           <  4:818729/148‑1 (MQ=255)
                                                                                                               gCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCcggac         >  3:670411/1‑148 (MQ=255)
                                                                                                                cTGATTTTGGAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGata                                                             <  3:112923/95‑1 (MQ=255)
                                                                                                                 tGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGCCCCTCTACCTGCtt                                                                                  <  2:769134/73‑1 (MQ=255)
                                                                                                                  gATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGata                                                             <  3:209575/93‑1 (MQ=255)
                                                                                                                  gATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGata                                                             >  4:209575/1‑93 (MQ=255)
                                                                                                                     ttCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGAtata                                      >  1:507155/1‑113 (MQ=255)
                                                                                                                       ccAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAg                                                    >  1:777539/1‑97 (MQ=255)
                                                                                                                       |                                                                                                                                                   
GGCGTTACCACCGCCGCCACAGTTTGAACCAGCAGACTCCAGTTTGTATGCAATGCAAATTCAGCGCTATCAACAGGCTGTGCGATCAACAATGGTCCTTTCACTGACTGAGCTGATTTCCAAAATCAGCTTAAAAAAAGCCTTTCAAAAGTAAGCAACGTCTGCTTACTGCCCCTCTACCTGCTTCGGCCGATAAAGCCGACGATAATACTCCAGACGTTGTAGATATAACGGCACGCTTTCATCGGGTATTCCGGACGGAATCG  >  NC_000913/3776006‑3776271

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: