Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,593,431 A→G 5.2% V431A (GTC→GCC)  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,4310AG 5.2% 269.2 / 4.9 98V431A (GTC→GCC) lgoTputative L‑galactonate:H+ symporter
Reads supporting (aligned to +/- strand):  ref base A (38/52);  new base G (2/3);  total (43/55)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.20e-01

ATAGCGAGGAGAACTGGGAAACATTTTAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGCCCACCGACGCAGTCATGCGAGAAG  >  NC_000913/4593284‑4593569
                                                                                                                                                   |                                                                                                                                          
aTAGCGAGGAGAACTGGGAAACATTTTAATCGCCATGAAACTGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGAc                                                                                                                                           <  2:832117/149‑1 (MQ=255)
aTAGCGAGGAGAACTGGGAAACATTTTAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGAc                                                                                                                                           >  2:324230/1‑149 (MQ=255)
aTAGCGAGGAGAACTGGGAAACATTTTAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGAc                                                                                                                                           <  2:1235809/149‑1 (MQ=255)
aTAGCGAGGAGAACTGGGAAACATTTTAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGAc                                                                                                                                           <  1:645364/149‑1 (MQ=255)
  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                                                                                                                      <  2:212753/149‑1 (MQ=255)
                     cATTTTAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGc                                                                                                                      <  1:585079/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)
                                                                                                                                                   |                                                                                                                                          
ATAGCGAGGAGAACTGGGAAACATTTTAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGCCCACCGACGCAGTCATGCGAGAAG  >  NC_000913/4593284‑4593569

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: