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 ← 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% 225.3 / 2.8 77V431A (GTC→GCC) ‡lgoTgalactonate:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base A (30/43);  new base G (2/2);  total (32/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.71e-01

AGATAGCGAGGAGAACTGGGAAACATTTTAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGCCCACCGACGCAGTCATGCGAGAAGCAACTGCGAC  >  NC_000913/4593282‑4593579
                                                                                                                                                     |                                                                                                                                                    
agttAGCTAGGAGAACTGGGAAACATTTTTATCGCCATGAAAATGCCTGATATGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCGGCGGTGa                                                                                                                                                      <  1:519467/147‑1 (MQ=255)
  aTAGCGAGGAGAACTGGGAAACATTTTAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACg                                                                                                                                                    <  1:535036/150‑1 (MQ=255)
  aTAGCGAGGAGAACTGGGAAACATTTTAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACg                                                                                                                                                    <  1:98920/150‑1 (MQ=255)
  aTAGCGAGGAGAACTGGGAAACATTTTAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGCCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACg                                                                                                                                                    <  2:1163960/150‑1 (MQ=255)
   tAGCGAGGAGAACTGGGAAACATTTTAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGc                                                                                                                                                   >  2:1249080/1‑150 (MQ=255)
       gaggagAACTGGGAAACATTTTAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAAcc                                                                                                                                               >  2:382550/1‑150 (MQ=255)
        aggagAACTGGGAAACATTTTAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCg                                                                                                                                              <  1:832079/150‑1 (MQ=255)
        aggagAACTGGGAAACATTTTAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCg                                                                                                                                              <  1:1582654/150‑1 (MQ=255)
              aCTGGGAAACATTTTAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGATGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAgatg                                                                                                                                        >  1:1446644/1‑150 (MQ=255)
              aCTGGGAAACATTTTAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAgatg                                                                                                                                        >  2:468272/1‑150 (MQ=255)
                tGGGAAACATTTTAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAATATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAgatgat                                                                                                                                      <  2:84093/150‑1 (MQ=255)
                tGGGAAACATTTTAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTTGCGTCACCGCAgcgggg                                                                                                                                      <  1:175068/150‑6 (MQ=255)
                tGGGAAACATTTTAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTTGCGTCACCGCAgcgggg                                                                                                                                      <  1:181968/150‑6 (MQ=255)
                   gAAACATTTTAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAg                                                                                                                                   >  2:1335095/1‑150 (MQ=255)
                   gAAACATTTTAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAg                                                                                                                                   >  1:414704/1‑150 (MQ=255)
                   gAAACATTTTAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAg                                                                                                                                   >  1:407415/1‑150 (MQ=255)
                    aaaCATTTTAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGt                                                                                                                                  >  1:603044/1‑150 (MQ=255)
                        aTTTTAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCa                                                                                                                              <  2:603044/150‑1 (MQ=255)
                         ttttAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCTGATGATTAGTGCCAg                                                                                                                             >  2:1318560/1‑150 (MQ=255)
                              aTCGCCATGAAAATGCCCGATAAGCAAAAAGTAAACGGGCATAAGGAGATTAATCTTTCCGTGGGAAGTTGATCGGATGACGAACCAGCAAGATGTACCCAAGTGCACCCGCTGCGGTGACGCAACCGCAGAAGATAAGTGACAGACCGa                                                                                                                        >  2:557281/1‑150 (MQ=255)
                                   cATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAg                                                                                                                   <  1:897306/150‑1 (MQ=255)
                                   cATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAg                                                                                                                   >  1:1737993/1‑150 (MQ=255)
                                    aTGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAgt                                                                                                                  <  1:1249080/150‑1 (MQ=255)
                                          aTGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAgtgggtgg                                                                                                            <  1:369408/150‑1 (MQ=255)
                                            gCCCGATAAGCAAAATGTTATCGTGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTa                                                                                                          <  2:1769471/150‑1 (MQ=255)
                                              ccGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATc                                                                                                        <  1:468272/150‑1 (MQ=255)
                                               cGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATca                                                                                                       <  2:441966/150‑1 (MQ=255)
                                                 aTAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATcaac                                                                                                     <  1:734087/150‑1 (MQ=255)
                                                  tAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATcaaca                                                                                                    <  1:160727/150‑1 (MQ=255)
                                                     gCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATa                                                                                                 >  2:931209/1‑150 (MQ=255)
                                                     gCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATa                                                                                                 <  1:1336740/150‑1 (MQ=255)
                                                      cAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATg                                                                                                                     <  2:1195321/129‑1 (MQ=255)
                                                      cAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATg                                                                                                                     >  1:1195321/1‑129 (MQ=255)
                                                       aaaaTGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATaaa                                                                                               >  2:793676/1‑150 (MQ=255)
                                                              tATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTa                                                                                        <  2:1045991/150‑1 (MQ=255)
                                                                tgccgCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAAt                                                                                      <  2:71042/146‑1 (MQ=255)
                                                                         aGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTTGTATGAACAATAAAACCAGTAATGATCGgcgc                                                                             <  2:634778/150‑1 (MQ=255)
                                                                          ggAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCa                                                                            <  2:1344540/150‑1 (MQ=255)
                                                                              aTTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAag                                                                        <  1:1541985/150‑1 (MQ=255)
                                                                                  aTCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGgc                                                                    >  2:664818/1‑150 (MQ=255)
                                                                                   tCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGgcc                                                                   >  1:20544/1‑149 (MQ=255)
                                                                                     tttACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGATATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCa                                                                 <  1:1335095/150‑1 (MQ=255)
                                                                                       tACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGa                                                               >  2:176718/1‑150 (MQ=255)
                                                                                       tACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCACAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGa                                                               >  2:496212/1‑150 (MQ=255)
                                                                                         cGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATg                                                             <  1:1874479/150‑1 (MQ=255)
                                                                                           tGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGaa                                                           <  2:414704/150‑1 (MQ=255)
                                                                                           tGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGaa                                                           <  2:407415/150‑1 (MQ=255)
                                                                                              gTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCt                                                        >  2:1379174/1‑150 (MQ=255)
                                                                                                cGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTgg                                                      <  1:15537/150‑1 (MQ=255)
                                                                                                     aacgccTGACGAACCAGGAAGATCTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTACAGCCAGACGCAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAg                                                 >  2:428420/6‑150 (MQ=255)
                                                                                                      tCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGt                                                >  1:1236678/1‑150 (MQ=255)
                                                                                                          cTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTg                                            <  1:931209/150‑1 (MQ=255)
                                                                                                           tGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTgg                                           <  2:1737993/150‑1 (MQ=255)
                                                                                                           tGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTgg                                           <  1:176718/150‑1 (MQ=255)
                                                                                                           tGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTgg                                           <  2:1933046/150‑1 (MQ=255)
                                                                                                           tGACGAACCAGGAAGATGTACGCCAGTACCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTgg                                           <  1:1845019/150‑1 (MQ=255)
                                                                                                            gatgAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGa                                          <  2:1446644/147‑1 (MQ=255)
                                                                                                            gACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAAcc                                                                                            >  2:858632/1‑100 (MQ=255)
                                                                                                            gACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAAcc                                                                                            <  1:858632/100‑1 (MQ=255)
                                                                                                            gACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGa                                          <  1:793676/150‑1 (MQ=255)
                                                                                                                  ccAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGc                                    >  1:565380/1‑150 (MQ=255)
                                                                                                                    aGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGccc                                  >  2:221182/1‑150 (MQ=255)
                                                                                                                    aGGAAGATGTACGACAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAAAAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGccc                                  >  2:1766804/1‑150 (MQ=255)
                                                                                                                     ggAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGCCCa                                 >  2:442661/1‑150 (MQ=255)
                                                                                                                            gTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGCCCACCGACGc                          >  1:811514/1‑150 (MQ=255)
                                                                                                                              aCGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGCCCACCGACGCAg                        <  2:1616193/150‑1 (MQ=255)
                                                                                                                              aCGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGCCCACCGACGCAg                        >  1:1438186/1‑150 (MQ=255)
                                                                                                                                 ccAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGCCCACCGACGCAGTCa                     <  1:256592/150‑1 (MQ=255)
                                                                                                                                   aGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGCCCACCGACGCAGTCATg                   <  2:1539724/150‑1 (MQ=255)
                                                                                                                                    gTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTGTGGATACTGCCCACCGACGCAGTCATGc                  <  2:1801818/150‑1 (MQ=255)
                                                                                                                                       cacccGCTGCGGTGACGCAACCTCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGCCGCAGATGAAGCTGGCAAAGTTCTGGATACTGCCCACCGACGCAGTCATGCgag               <  2:1460499/148‑1 (MQ=255)
                                                                                                                                        ccccGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGCCCACCGACGCAGTCATGCgaga              >  2:463201/1‑150 (MQ=255)
                                                                                                                                         tcatcTGCGGTTGCGTCACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGCCCACCGACGCAGTCATGCGAGaa             >  2:528370/5‑150 (MQ=255)
                                                                                                                                           cGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGCCCACCGACGCAGTCATGCGAGAAGc           <  1:382550/150‑1 (MQ=255)
                                                                                                                                           atcTGCGGTTGCGTCACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGCCCACCGACGCAGTCATGCGAGAAGc           >  2:109878/3‑150 (MQ=255)
                                                                                                                                            gCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGCCCACCGACGCAGTCATGCGAGAAGCa          <  2:1378675/150‑1 (MQ=255)
                                                                                                                                                 ggTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGCCCACCGACGCAGTCATGCGAGAAGCAACTGc     <  1:866563/150‑1 (MQ=255)
                                                                                                                                                    gACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGCCCACCGACGCAGTCATGCGAGAAGCAACTGCGAc  >  1:1697436/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AGATAGCGAGGAGAACTGGGAAACATTTTAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGCCCACCGACGCAGTCATGCGAGAAGCAACTGCGAC  >  NC_000913/4593282‑4593579

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