Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,593,435 A→C 5.1% C430G (TGC→GGC)  lgoT ← galactonate:H(+) symporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,593,4350AC 5.1% 248.3 / 4.7 78C430G (TGC→GGC) lgoTgalactonate:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base A (30/44);  new base C (2/2);  total (32/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.38e-01

GAGGAGAACTGGGAAACATTTTAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGCCCACCGACGCAGTCATGCGAGAAGCAACTGCGACGTGG  >  NC_000913/4593289‑4593583
                                                                                                                                                  |                                                                                                                                                    
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:181968/150‑6 (MQ=255)
         tGGGAAACATTTTAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTTGCGTCACCGCAgcgggg                                                                                                                                          <  1:175068/150‑6 (MQ=255)
            gAAACATTTTAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAg                                                                                                                                       >  1:414704/1‑150 (MQ=255)
            gAAACATTTTAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAg                                                                                                                                       >  1:407415/1‑150 (MQ=255)
            gAAACATTTTAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAg                                                                                                                                       >  2:1335095/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:1737993/1‑150 (MQ=255)
                            cATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAg                                                                                                                       <  1:897306/150‑1 (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                                                                                                     <  1:1336740/150‑1 (MQ=255)
                                              gCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATa                                                                                                     >  2:931209/1‑150 (MQ=255)
                                               cAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATg                                                                                                                         >  1:1195321/1‑129 (MQ=255)
                                               cAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATg                                                                                                                         <  2:1195321/129‑1 (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:1933046/150‑1 (MQ=255)
                                                                                                    tGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTgg                                               <  1:176718/150‑1 (MQ=255)
                                                                                                    tGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTgg                                               <  2:1737993/150‑1 (MQ=255)
                                                                                                    tGACGAACCAGGAAGATGTACGCCAGTACCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTgg                                               <  1:1845019/150‑1 (MQ=255)
                                                                                                     gatgAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGa                                              <  2:1446644/147‑1 (MQ=255)
                                                                                                     gACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAAcc                                                                                                <  1:858632/100‑1 (MQ=255)
                                                                                                     gACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAAcc                                                                                                >  2:858632/1‑100 (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                            >  1:1438186/1‑150 (MQ=255)
                                                                                                                       aCGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGCCCACCGACGCAg                            <  2:1616193/150‑1 (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)
                                                                                                                                               cGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGCCCACCGACGCAGTCATGCGAGAAGCAACTGCGACGt    >  1:475010/1‑150 (MQ=255)
                                                                                                                                               cGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGCCCACCGACGCAGTCATGCGAGAAGCAACTGCGACGt    <  2:455094/150‑1 (MQ=255)
                                                                                                                                                gCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGCCCACCGACGCAGTCATGCGAGAAGCAACTGCGACGTg   <  1:442661/150‑1 (MQ=255)
                                                                                                                                                gCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGCCCACCGACGCAGTCATGCGAGAAGCAACTGCGACGTg   <  1:1464666/150‑1 (MQ=255)
                                                                                                                                                 cAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGCCCACCGACGCAGTCATGCGAGAAGCAACTGCGACGTgg  <  1:1379174/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GAGGAGAACTGGGAAACATTTTAATCGCCATGAAAATGCCCGATAAGCAAAATGTTATCGGGCATAAGGAGATTAATCTTTACGTGGGTCGTTGATCGGCTGACGAACCAGGAAGATGTACGCCAGTGCCCCCGCTGCGGTGACGCAACCGCAGATGATTAGTGCCAGACGGAATGAGTGGGTGGTATCAACAATAAAACCAGTAATGATCGGCGCAAAAGAGGCGCAGATGAAGCTGGCAAAGTTCTGGATACTGCCCACCGACGCAGTCATGCGAGAAGCAACTGCGACGTGG  >  NC_000913/4593289‑4593583

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