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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,593,4340CT 5.4% 258.7 / 4.6 95C430Y (TGC→TAC) lgoTputative L‑galactonate:H+ symporter
Reads supporting (aligned to +/- strand):  ref base C (40/48);  new base T (2/3);  total (44/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.48e-01

AGCGAGGAGAACTGGGAAACATTTTAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGCCCACCGACGCAGTCATGCGAGAAGCAACTGCGACGT  >  NC_000913/4593286‑4593581
                                                                                                                                                    |                                                                                                                                                   
atcgcggagAACTTGGTATCATTTTAATCGCCATGAAAATGCCCGATAAGCAAAATGTTAGCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGc                                                                                                                                                     <  2:1489097/144‑1 (MQ=255)
  cGAGGAGAACTGGGAAACATTTTAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCaa                                                                                                                                                   >  1:706979/1‑149 (MQ=255)
   gaggagAACTGGGAAACATTTTAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGCCATAAGGAGATAAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAAc                                                                                                                                                  >  2:603019/1‑149 (MQ=255)
    aggagAACTGGGAAACATTTTAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAAcc                                                                                                                                                 >  2:2038/1‑149 (MQ=255)
    aggagAACTGGGAAACATTTTAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCAGAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAAcc                                                                                                                                                 >  1:1262291/1‑149 (MQ=255)
      gagaACTGGGAAACATTTTAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGc                                                                                                                                               >  1:110274/1‑149 (MQ=255)
       agaACTGTGAAACATTTTAATCGCCATGATAATGCCCGATAAGCAAAATGTTAGCGGGCATAAGGAGATTAATCTTTACGTGGGTCTTTGATCGGCTGACGAACCAGGAAGATGAACGCCAGTGCCCCCGCTGCTGTGACGCAACCGCa                                                                                                                                              <  2:1427154/149‑1 (MQ=255)
        gaACTGGGAAACATTTTAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAg                                                                                                                                             <  1:343668/149‑1 (MQ=255)
         aaCTGGGAAACATTTTAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAga                                                                                                                                            <  2:1274905/149‑1 (MQ=255)
         aaCTGGGAAACATTTTAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGGGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTTGCGTCACCGCAgc                                                                                                                                            <  1:612845/149‑2 (MQ=255)
          aCTGGGAAACATTTTAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAgat                                                                                                                                           <  2:191370/149‑1 (MQ=255)
           cTGGGAAACATTTTAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTTGCGTCACCGCAgcgg                                                                                                                                          <  1:1191560/149‑4 (MQ=255)
           cTGGGAAACATTTTAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGCCGTTGATCGGCTGCCGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTTGCGTCACCGCAgcgg                                                                                                                                          <  2:243230/149‑4 (MQ=255)
            tGGGAAACATTTTAATCGCCATGAAAATGCCCGATAAACAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAgatga                                                                                                                                         <  1:1102594/149‑1 (MQ=255)
             gggAAACAATTTAATAGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGAAGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAgatgat                                                                                                                                        <  2:1418756/149‑1 (MQ=255)
              ggAAACATTTTAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCTCTGCGGTGACGCAACCGCAGATGAtt                                                                                                                                       <  2:273865/149‑1 (MQ=255)
              ggAAACATTTTAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGCCGCAACCGCAGATGAtt                                                                                                                                       >  2:636460/1‑149 (MQ=255)
              ggAAACATTTTAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGAtt                                                                                                                                       <  2:71177/149‑1 (MQ=255)
              ggAAACATTTTAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGAtt                                                                                                                                       >  2:198134/1‑149 (MQ=255)
               gAAACATTTTAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCCGGCATAAGGAGATAAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGCTGACGCAACCGCAGATGATTa                                                                                                                                      >  1:965619/1‑149 (MQ=255)
                aaaCATTTTAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGGCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGGGACGCAACCGCAGATGATTAg                                                                                                                                     <  1:809776/149‑1 (MQ=255)
                 aaCATTTTAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGt                                                                                                                                    <  1:1126336/149‑1 (MQ=255)
                   cATTTTAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGc                                                                                                                                  <  1:585079/149‑1 (MQ=255)
                   cATTTTAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGc                                                                                                                                  <  2:212753/149‑1 (MQ=255)
                    aTTTTAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGcc                                                                                                                                 >  2:131976/1‑149 (MQ=255)
                      tttAATCGCCATGAAAATGCCCGATATGCTAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGAAACCGCCGATGATTTGTGcctg                                                                                                                               >  1:1056047/1‑147 (MQ=255)
                      tttAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGCCCCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCTCCCGCAGATGATTAGTGCCAg                                                                                                                               >  1:348829/1‑149 (MQ=255)
                      tttAATCGCCATGAAAATGCCCGAGAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGAGGGACGTTGATCGGCTGACGAACAAGGAAGATGTCCGCCAGTGCCCCCGCTGCGGTGACGCAACCCCAGATGATTAGTGCCAg                                                                                                                               >  2:872510/1‑149 (MQ=255)
                           tCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGa                                                                                                                          <  1:131976/149‑1 (MQ=255)
                             gCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAAt                                                                                                                        <  1:509613/149‑1 (MQ=255)
                             gCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAAt                                                                                                                        <  2:519401/149‑1 (MQ=255)
                              ccATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCACCCGCAGATGATTAGTGCCAGACGGAATg                                                                                                                       >  2:716359/1‑149 (MQ=255)
                              ccATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATg                                                                                                                       >  2:1128223/1‑149 (MQ=255)
                               cATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGa                                                                                                                      <  2:392267/149‑1 (MQ=255)
                                     aaTGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAgtgggt                                                                                                                <  2:1290114/149‑1 (MQ=255)
                                     aaTGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAgtgggt                                                                                                                <  1:62833/149‑1 (MQ=255)
                                      aTGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTCCGTGGGGCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAgtgggtg                                                                                                               <  1:144088/149‑1 (MQ=255)
                                      aTGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAgtgggtg                                                                                                               <  2:23253/149‑1 (MQ=255)
                                       tGCGCGATATGCAAAATGTTATGGGGCATAAGGAGATTATTCTTTACGGTGGGCTTTTATCGGCTTACGAAGCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGAAGCAACCGCAGATGATTAGTGCCAGACGGAATGAgtgggttg                                                                                                              <  1:219233/149‑3 (MQ=255)
                                        gCCCGATAATCAAAATGTTATCGGGAATAAGGAGATTAATCTTTACGTGGGGCGTTGATCGGCTCACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGggtggt                                                                                                             <  2:269765/149‑1 (MQ=255)
                                         cccGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTa                                                                                                            >  1:543284/1‑149 (MQ=255)
                                         cccGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACAAGGAAGATGTACGCCAGTGCCCCCGCTGCGCTGACGCAACCGCAGATGATTAGTGCCAGACGAAATGAGTGGGTGGTa                                                                                                            >  1:847772/1‑149 (MQ=255)
                                           cGATAAGCAAAATGTTATCGGGCATAAGGAGATAAATCTTTACGTGGGTCGTTGATCGGCTGAGGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATc                                                                                                          <  2:1252778/149‑1 (MQ=255)
                                            gATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTCCGCCAGTGCCCCCGCTGCGGTGACGCAGCCGCAGATGATTAGTGCCAGAGGGAATGAGTGGGTGGTATca                                                                                                         >  2:674290/1‑149 (MQ=255)
                                             aTAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATcaa                                                                                                        <  2:495828/149‑1 (MQ=255)
                                               aaGCAAATTGTTATCGGGTATAAGGAGATTTATCTTTACGTGGGTAGTTGATCGGGTGAAGAACCAGGTAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGTTGATTAGTGCCAGACGGAATGAGTGGGTGGTTTcaaca                                                                                                      <  1:651569/149‑1 (MQ=255)
                                                gataaaaTGTTCTCGTTCTTAAGGAGATTAATCTTTTCGGGGGTCGTTGATCTGCTTACGAACCATGAAGATGTACTACATTGCCCACTCTGCGGTGACGCAACCGCAGATTATTAGTTCCAGACGGAATGATTTTGGGGTATctacaa                                                                                                     <  2:252947/146‑1 (MQ=255)
                                                 gCAAAATGTTATCGGGCATAAGGAGTTTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAAt                                                                                                    <  1:964224/149‑1 (MQ=255)
                                                     aaTGTTATCGGGCATAAGGAGATTAATGTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATaaaa                                                                                                <  2:346747/149‑1 (MQ=255)
                                                      aTGTTATCGGGCATAAGGAGATTAATCTTTACTTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCAAGTGCACCCGCTTCGGTGACGCAACCGCAGATTATTAGTTCCAGACGGAATGAGTGGGTTGTATcaacaaaagatc                                                                                               >  1:252947/1‑143 (MQ=255)
                                                            tCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGTAGGTATCAACAATAAAACCAGTaa                                                                                         >  1:703092/1‑149 (MQ=255)
                                                             tgggAATAAGGAGCTTAATCTTTCCGTGGGGCGTTGATCGGCTGACGAACCAGGCAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAAt                                                                                        <  2:1352342/148‑1 (MQ=255)
                                                             cGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAAt                                                                                        <  1:538203/149‑1 (MQ=255)
                                                              gggCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATg                                                                                       <  1:716359/149‑1 (MQ=255)
                                                                    aaGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGgc                                                                                 >  1:200597/1‑149 (MQ=255)
                                                                      ggAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGGAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGgcgc                                                                               >  1:20135/1‑149 (MQ=255)
                                                                      ggAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGgcgc                                                                               >  2:109058/1‑149 (MQ=255)
                                                                       gagaTTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCa                                                                              >  2:621647/1‑149 (MQ=255)
                                                                          aTTAATCTTTACGTGGGTCGTTGCTCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGCCGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCaaaa                                                                           >  2:325740/1‑149 (MQ=255)
                                                                            tactCTTTACGGGGGTCCTTGATCGGCTTACGAACCAGGATGATGTACGCCAGTGCCCCCGCTGCGGTTACGCAACAGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAaga                                                                         <  1:1457463/146‑1 (MQ=255)
                                                                            tAATCTCTCCGTGGGTCGTTGATCGGCTGACGCACCAGGACGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAaga                                                                         <  2:1436126/149‑1 (MQ=255)
                                                                             aaTCTTTACGGTGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTTGTATCAACAATAAAACCAGTAATGATCGGCGCAAAagag                                                                        <  1:100819/149‑1 (MQ=255)
                                                                                   tACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGCCGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAg                                                                  >  2:1176128/1‑149 (MQ=255)
                                                                                    aCGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCTGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGa                                                                 >  2:1118088/1‑149 (MQ=255)
                                                                                        gggTCGTGTATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCGCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGaa                                                             <  2:1334515/149‑1 (MQ=255)
                                                                                           tCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCt                                                          >  1:1532273/1‑149 (MQ=255)
                                                                                             gTTGAGCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCGGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTgg                                                        <  1:721188/149‑1 (MQ=255)
                                                                                               tGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCa                                                      >  1:1463973/1‑149 (MQ=255)
                                                                                                    ggCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACAAGAAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGtt                                                 <  2:27131/149‑1 (MQ=255)
                                                                                                    ggCTGACGAACCAGGAAGATGTACGCCAGTGCCACCGCTGCGGTGACGCAACCGCAGAGGATTAGTTCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATTAAGCTGTCAAAGtt                                                 >  1:343088/1‑149 (MQ=255)
                                                                                                       tGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTg                                              >  1:142512/1‑149 (MQ=255)
                                                                                                        gACGAACCAGGAAGATGTACGCCAGGGCCTCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTgg                                             <  1:426219/149‑1 (MQ=255)
                                                                                                           gAACCAGGAAGCTGAACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGTATGAGTTGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATa                                          <  2:1119007/149‑1 (MQ=255)
                                                                                                           gAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGAt                                                                                     >  2:672524/1‑106 (MQ=255)
                                                                                                           gAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGAt                                                                                     <  1:672524/106‑1 (MQ=255)
                                                                                                           gAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGGGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATa                                          <  2:921405/149‑1 (MQ=255)
                                                                                                             aCCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGCCGCAACCGCAGAGGAGTAGTGCCAGACGGAATGACTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGACGCTGGCAAAGTTCTGGATACt                                        >  1:79768/1‑149 (MQ=255)
                                                                                                             aCCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACt                                        >  2:763855/1‑149 (MQ=255)
                                                                                                              ccAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTg                                       >  2:1209633/1‑149 (MQ=255)
                                                                                                              ccAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACAGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTg                                       <  1:324230/149‑1 (MQ=255)
                                                                                                              ccAGGAAGATGTACGCAAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTg                                       >  1:1464684/1‑149 (MQ=255)
                                                                                                               cAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGAGCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGGTCTGGATACTGc                                      <  1:1352020/149‑1 (MQ=255)
                                                                                                                  gAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGCCCa                                   >  2:1218002/1‑149 (MQ=255)
                                                                                                                   aaGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGCCCAc                                  >  2:1375539/1‑149 (MQ=255)
                                                                                                                         tACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGCATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGCCCACCGACGc                            >  2:690796/1‑149 (MQ=255)
                                                                                                                              cAGTACCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGCCCACCGACGCAGTCa                       >  2:943615/1‑149 (MQ=255)
                                                                                                                                gTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGCCCACCGACGCAGTCATg                     <  2:110274/149‑1 (MQ=255)
                                                                                                                                  gCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGCCCACCGACGCAGTCATGCg                   <  1:1097789/149‑1 (MQ=255)
                                                                                                                                   cccccGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGCCCACCGACGCAGTCATGCga                  <  2:66901/149‑1 (MQ=255)
                                                                                                                                     tcatcTGCGGTTGCGTCACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGCTCGGCGCTAAAGAGGCGCAGATGAAGCTGGCAATGTTCTGGATACTGCCCACCGACGCCGTCATGCgaga                >  1:1293639/5‑149 (MQ=255)
                                                                                                                                      ccGCTTCGGGGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGCCCACCGACGCAGTCATGCGAGaa               <  2:25011/149‑1 (MQ=255)
                                                                                                                                      ccGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTCGATACTGCCCACCGACGCAGTCATGCGAGaa               >  1:343003/1‑149 (MQ=255)
                                                                                                                                      ccGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAAAGAGTGGGTGGTATCAACAATAAAACCAGTAATGAGCGGCGCAAAAGAGGCGCAGATGAAGCTTGCAAAGGTCTGGATACTGCCCACCGACGCAGTCATGCGAGaa               <  2:1262291/149‑1 (MQ=255)
                                                                                                                                       atcTGCGGTTGCGTCACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGCCCACCGACGCAGTCATGCGAGAAg              >  2:427341/3‑149 (MQ=255)
                                                                                                                                                   gCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGCCCACCGACGCAGTCATGCGAGAAGCAACTGCTACGt  >  2:72567/1‑149 (MQ=255)
                                                                                                                                                   gCAACCGCAGATGATTAGTGCCAGACGCAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGAAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGCCCACCGACGCAGGCATGCGAGAAGCAACTGCGACGt  >  2:1438884/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
AGCGAGGAGAACTGGGAAACATTTTAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGCCCACCGACGCAGTCATGCGAGAAGCAACTGCGACGT  >  NC_000913/4593286‑4593581

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