Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,593,434 C→T 5.1% C430Y (TGC→TAC) ‡ lgoT ← galactonate:H(+) symporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,593,4340CT 5.1% 318.7 / 4.0 98C430Y (TGC→TAC) ‡lgoTgalactonate:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base C (44/49);  new base T (2/3);  total (46/52)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.31e-01

TAGCGAGGAGAACTGGGAAACATTTTAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGCCCACCGACGCAGTCATGCGAGAAGCAACTGCGACGTGG  >  NC_000913/4593285‑4593583
                                                                                                                                                     |                                                                                                                                                     
tAGCGAGGAGAACTGGGAAACATTTTAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGGTTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGc                                                                                                                                                       <  1:159615/150‑1 (MQ=255)
  gCGAGGAGAACTGGGAAACATTTTAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCaa                                                                                                                                                     >  1:629345/1‑150 (MQ=255)
    gaggagAACTGGGAAACATTTTAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGGCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAAcc                                                                                                                                                   <  1:159017/150‑1 (MQ=255)
      ggagAACTGGGAAACATTTTAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGc                                                                                                                                                 >  1:595897/1‑150 (MQ=255)
       gagaACTGGGAAACATTTTAATCGCCATGAAAATGCCCGATAATCAAAATTTTATCGGTCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCATTTCCCCCGCTGCGGTGACGCAACCGCa                                                                                                                                                <  2:901094/150‑1 (MQ=255)
       gagaACTGGGAAACATTTTAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCa                                                                                                                                                <  1:1174223/150‑1 (MQ=255)
       gagaACTGGGAAACATTTTAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCa                                                                                                                                                >  1:922830/1‑150 (MQ=255)
       gagaACTGGGAAACATTTTAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTGGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCa                                                                                                                                                <  1:641549/150‑1 (MQ=255)
        agaACTGGGAAACATTTTAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAg                                                                                                                                               >  1:268237/1‑150 (MQ=255)
         gaACTGGGAAACATTTTAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTTGCGTCACCGCAgc                                                                                                                                              <  2:896733/150‑2 (MQ=255)
         gaACTGGGAAACATTTTAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAga                                                                                                                                              >  1:45107/1‑150 (MQ=255)
           aCTGGGAAACATTTTAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAgatg                                                                                                                                            <  2:1304459/150‑1 (MQ=255)
           aCTGGGAAACATTTTAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAgatg                                                                                                                                            <  1:1428462/150‑1 (MQ=255)
           aCTGGGAAACATTTTAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAgatg                                                                                                                                            >  1:646713/1‑150 (MQ=255)
            cTGGGAAACATTTTAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTTGCGTCACCGCAgcgg                                                                                                                                            >  1:630563/1‑146 (MQ=255)
            cTGGGAAACATTTTAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTTGCGTCACCGCAgcgg                                                                                                                                            <  2:630563/149‑4 (MQ=255)
            cTGGGAAACATTTTAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAgatga                                                                                                                                           >  1:442122/1‑150 (MQ=255)
             tGGGAAACATTTTAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTTGCGTCACCGCAgcgggg                                                                                                                                          <  1:560956/150‑6 (MQ=255)
             tGGGAAACATTTTAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAgatgat                                                                                                                                          <  1:373863/150‑1 (MQ=255)
              gggAAACATTTTAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGAtt                                                                                                                                         <  1:681719/150‑1 (MQ=255)
              gggAAACATTTTAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGAtt                                                                                                                                         >  1:1376408/1‑150 (MQ=255)
               ggAAACATTTTAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTa                                                                                                                                        <  1:1061115/150‑1 (MQ=255)
                    cATTTTAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGcc                                                                                                                                   <  1:1036422/150‑1 (MQ=255)
                         tAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACg                                                                                                                              >  1:783981/1‑150 (MQ=255)
                          aaTCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACgg                                                                                                                             >  1:905621/1‑150 (MQ=255)
                          aaTCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACgg                                                                                                                             >  2:1506456/1‑150 (MQ=255)
                             cGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAAt                                                                                                                          >  2:1017936/1‑150 (MQ=255)
                             cGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAAt                                                                                                                          >  1:221004/1‑150 (MQ=255)
                                   gAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAgtgg                                                                                                                    <  2:426549/150‑1 (MQ=255)
                                      aaTGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAgtgggtg                                                                                                                 <  2:642810/150‑1 (MQ=255)
                                       aTGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAgtgggtgg                                                                                                                >  1:378510/1‑150 (MQ=255)
                                         gCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCACCCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTa                                                                                                              >  1:403228/1‑150 (MQ=255)
                                         gCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTa                                                                                                              >  1:757546/1‑150 (MQ=255)
                                          cccGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTAt                                                                                                             <  2:48169/150‑1 (MQ=255)
                                           ccGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATc                                                                                                            <  2:757546/150‑1 (MQ=255)
                                            cgttAAGCAAAATTTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGTTCGGCTGACGAACCAGGAAGATGTTCGCCAGTTCCCGCGCGGGGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTCTca                                                                                                           <  2:1522310/147‑1 (MQ=255)
                                            cGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATca                                                                                                           <  2:480333/150‑1 (MQ=255)
                                                  gCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGGCGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATa                                                                                                     <  2:629791/150‑1 (MQ=255)
                                                  gCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATa                                                                                                     <  2:758708/150‑1 (MQ=255)
                                                  gCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATa                                                                                                     >  1:250194/1‑150 (MQ=255)
                                                  ccAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATa                                                                                                     <  2:1319000/149‑1 (MQ=255)
                                                    aaaaTGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGGAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATaaa                                                                                                   >  2:98148/1‑150 (MQ=255)
                                                     aaaTGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATaaaa                                                                                                  <  1:374293/150‑1 (MQ=255)
                                                     aaaTGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATaaaa                                                                                                  <  1:368793/150‑1 (MQ=255)
                                                      aaTGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGGAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAAc                                                                                                 <  1:98148/150‑1 (MQ=255)
                                                       aTGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATAGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAAcc                                                                                                <  2:646713/150‑1 (MQ=255)
                                                        tGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCa                                                                                               >  1:17169/1‑150 (MQ=255)
                                                         gTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAg                                                                                              >  1:405921/1‑150 (MQ=255)
                                                            aTCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTaa                                                                                           >  2:475294/1‑150 (MQ=255)
                                                             tCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAAt                                                                                          <  2:929304/150‑1 (MQ=255)
                                                                 gCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATc                                                                                      >  2:1320369/1‑150 (MQ=255)
                                                                   aTAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCgg                                                                                    <  2:1523760/150‑1 (MQ=255)
                                                                   aTAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCgg                                                                                    <  2:45107/150‑1 (MQ=255)
                                                                   aTAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCgg                                                                                    >  1:1523760/1‑150 (MQ=255)
                                                                    tAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGATGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGgc                                                                                   >  2:870726/1‑150 (MQ=255)
                                                                     aaGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGgcg                                                                                  <  2:1475790/150‑1 (MQ=255)
                                                                          gaTTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCaaaa                                                                             >  1:350209/1‑150 (MQ=255)
                                                                            ttAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAaga                                                                           >  1:275252/1‑150 (MQ=255)
                                                                                 cTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGgcgc                                                                      <  2:449778/150‑1 (MQ=255)
                                                                                  tttACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCa                                                                     >  1:1237670/1‑150 (MQ=255)
                                                                                     aCGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGAt                                                                  <  1:1174005/150‑1 (MQ=255)
                                                                                           gTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCt                                                            <  1:697008/150‑1 (MQ=255)
                                                                                             cGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTgg                                                          <  2:1416000/150‑1 (MQ=255)
                                                                                             cGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTgg                                                          <  2:17169/150‑1 (MQ=255)
                                                                                              gTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGc                                                         <  2:848549/150‑1 (MQ=255)
                                                                                                 gATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCaaa                                                      >  2:993205/1‑150 (MQ=255)
                                                                                                 gATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCaaa                                                      >  2:62004/1‑150 (MQ=255)
                                                                                                 gATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCaaa                                                      >  1:394264/1‑150 (MQ=255)
                                                                                                    cGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCACGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGCTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGtt                                                   >  2:1391395/1‑150 (MQ=255)
                                                                                                       cTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTg                                                <  1:62004/150‑1 (MQ=255)
                                                                                                        tGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTgg                                               <  2:868504/150‑1 (MQ=255)
                                                                                                        tGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTgg                                               <  2:905621/150‑1 (MQ=255)
                                                                                                         gACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCaa                                                       >  2:703518/1‑141 (MQ=255)
                                                                                                         gACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCaa                                                       <  1:703518/141‑1 (MQ=255)
                                                                                                         gACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGa                                              <  2:783981/150‑1 (MQ=255)
                                                                                                              aCCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTg                                         >  2:1455025/1‑150 (MQ=255)
                                                                                                                cAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGcc                                       >  1:395445/1‑150 (MQ=255)
                                                                                                                 aGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAGAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGccc                                      <  1:1017936/150‑1 (MQ=255)
                                                                                                                 aGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGccc                                      <  1:1506456/150‑1 (MQ=255)
                                                                                                                  ggAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGTAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGCCCa                                     <  2:466616/150‑1 (MQ=255)
                                                                                                                   gAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGCCCAc                                    <  2:1250093/150‑1 (MQ=255)
                                                                                                                   gAAGATGTACGCCAGTGCCCCCGCTGCGGGGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGCCCAc                                    <  2:1255593/150‑1 (MQ=255)
                                                                                                                       aTGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGCCCACCGAc                                >  2:1333758/1‑150 (MQ=255)
                                                                                                                         gTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGCCCACCGACGc                              >  2:487672/1‑150 (MQ=255)
                                                                                                                          tACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAACCCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGCCCACCGACGCa                             >  1:943923/1‑150 (MQ=255)
                                                                                                                          tACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGCCCACCGACGCa                             >  2:342016/1‑150 (MQ=255)
                                                                                                                          tACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGCCCACCGACGCa                             <  2:595897/150‑1 (MQ=255)
                                                                                                                              ccAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGCCCACCGACGCAGTCa                         >  1:760418/1‑150 (MQ=255)
                                                                                                                              ccAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGCCCACCGACGCAGTCa                         <  2:1237443/150‑1 (MQ=255)
                                                                                                                                   gCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGCCCACCGACGCAGTCATGCga                    <  1:993205/150‑1 (MQ=255)
                                                                                                                                    cccccGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGCCCACCGACGCAGTCATGCgag                   >  2:14267/1‑150 (MQ=255)
                                                                                                                                      cccGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGCCCACCGACGCAGTCATGCGAGaa                 >  1:693049/1‑150 (MQ=255)
                                                                                                                                        cGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGCCCACCGACGCAGTCATGCGAGAAGc               <  2:442122/150‑1 (MQ=255)
                                                                                                                                        atcTGCGGTTGCGTCACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGCCCACCGACGCAGTCATGCGAGAAGc               >  2:980990/3‑150 (MQ=255)
                                                                                                                                             cGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGCCCACCGACGCAGTCATGCGAGAAGCAACTg          <  2:221004/150‑1 (MQ=255)
                                                                                                                                              ggTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGCCCACCGACGCAGTCATGCGAGAAGCAACTGc         >  2:1265833/1‑150 (MQ=255)
                                                                                                                                              ggTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGCCCACCGACGCAGTCATGCGAGAAGCAACTGc         >  1:91948/1‑150 (MQ=255)
                                                                                                                                                  aCGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGCCCACCGACGCAGTCATGCGAGAAGCAACTGCGACg     <  2:281508/150‑1 (MQ=255)
                                                                                                                                                   cGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGCCCACCGACGCAGTCATGCGAGAAGCAACTGCGACGt    >  1:120987/1‑150 (MQ=255)
                                                                                                                                                     cAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGCCCACCGACGCAGTCATGCGAGAAGCGACTGCGACGTgg  <  1:497277/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TAGCGAGGAGAACTGGGAAACATTTTAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGCCCACCGACGCAGTCATGCGAGAAGCAACTGCGACGTGG  >  NC_000913/4593285‑4593583

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