Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 256,913 T→A 5.5% A129A (GCT→GCA frsA → fermentation‑respiration switch protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913256,9130TA 5.5% 350.7 / 7.4 109A129A (GCT→GCAfrsAfermentation‑respiration switch protein
Reads supporting (aligned to +/- strand):  ref base T (60/43);  new base A (2/4);  total (62/47)
Fisher's exact test for biased strand distribution p-value = 3.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.61e-01

Only 100 of 110 total aligned reads displayed.

TACCGACGATGATTTATACGATACGGTGATTGGCTACCGTGGCGGCAACTGGATTTATGAGTGGGCCACCCAGGCGATGGTGTGGCAACAAAAAGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTGTCAAACCGCGCCTATGAAGAGGCCGCTCAGCGTCTACCGGGCACGATGCGGCAGATGGAGTTTACCGTACCCGGCGGTGC  >  NC_000913/256766‑257062
                                                                                                                                                   |                                                                                                                                                     
tACCGACGATGATTTATACGATACGGTGATTGGCTACCGTGGCGGCAACTGGATTTATGAGTGGGCCACCCAGGCGATGGTGTGGCAACAAAAAGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTAc                                                                                                                                                     >  2:1154975/1‑150 (MQ=255)
  ccGACGATGATTTATACGATACGGTGATTGGCTACCGTGGCGGCAACTGGATTTATGAGTGGGCCACCCAGGCGATGGTGTGGCAACAAAAAGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGt                                                                                                                                                   <  1:1529885/150‑1 (MQ=255)
    gacgaTGATTTATACGATACGGTGATTGGCTACCGTGGCGGCAACTGGATTTATGAGTGGGCCACCCAGGCGATGGTGTGGCAACAAAAAGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTg                                                                                                                                                 <  2:1468389/150‑1 (MQ=255)
    gacgaTGATTTATACGATACGGTGATTGGCTACCGTGGCGGCAACTGGATTTATGAGTGGGCCACCCAGGCGATGGTGTGGCAACAAAAAGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTg                                                                                                                                                 >  1:930844/1‑150 (MQ=255)
      cgaTGATTTATACGATACGGTGATTGGCTACCGTGGCGGCAACTGGATTTATGAGTGGGCCACCCAGGCGATGGTGTGGCAACAAAAAGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTa                                                                                                                                               <  2:1451767/150‑1 (MQ=255)
       gatgatTTATACGATACGGTGATTGGCTACCGTGGCGGCAACTGGATTTATGAGTGGGCCACCCAGGCGATGGTGTGGCAACAAAAAGCCTGTGCGGAAGACGATCCGCAATTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTac                                                                                                                                              <  1:661210/150‑1 (MQ=255)
       gatgatTTATACGATACGGTGATTGGCTACCGTGGCGGCAACTGGATTTATGAGTGGGCCACCCAGGCGATGGTGTGGCAACAAAAAGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTac                                                                                                                                              >  1:591932/1‑150 (MQ=255)
       gatgatTTATACGATACGGTGATTGGCTACCGTGGCGGCAACTGGATTTATGAGTGGGCCACCCAGGCGATGGTGTGGCAACAAAAAGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTac                                                                                                                                              >  1:1488765/1‑150 (MQ=255)
       gatgatTTATACGATACGGTGATTGGCTACCGTGGCGGCAACTGGATTTATGAGGGGGCCACCCAGGCGATGGGGGGGCAACAAAAAGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTac                                                                                                                                              <  2:520589/150‑1 (MQ=255)
             ttATACGATACGGTGATTGGCTACCGTGGCGGCAACTGGATTTATGAGTGGGCCACCCAGGCGATGGTGTGGCAACAAAAAGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACAtt                                                                                                                                        >  2:491470/1‑150 (MQ=255)
               atacgatacgGTGATTGGCTACCGTGGCGGCAACTGGATTTATGAGTGGGCCACCCAGGCGATGGTGTGGCAACAAAAAGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTgc                                                                                                                                      >  1:361968/1‑150 (MQ=255)
               atacgatacgGTGATTGGCTACCGTGGCGGCAACTGGATTTATGAGTGGGCCACCCAGGCGATGGTGTGGCAACAAAAAGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTgc                                                                                                                                      >  1:622642/1‑150 (MQ=255)
                   gatacgGTGATTGGCTACCGTGGCGGCAACTGGATTTATGAGTGGGCCACCCAGGCGATGGTGTGGCAACAAAAAGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTgccgcc                                                                                                                                  <  2:1530558/150‑1 (MQ=255)
                    atacgGTGATTGGCTACCGTGGCGGCAACTGGATTTATGAGTGGGCCACCCAGGCGATGGTGTGGCAACAAATAGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTGGTCACTGGCTGCATGCGGCAATGTTGTACAACGTAgccgcat                                                                                                                                 <  2:69806/150‑3 (MQ=255)
                    atacgGTGATTGGCTACCGTGGCGGCAACTGGATTTATGAGTGGGCCACCCAGGCGATGGTGTGGCAACAAAAAGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCAATGTTGTACAACGTAgccgcat                                                                                                                                 <  1:1065869/150‑3 (MQ=255)
                     tacgGTGATTGGCTACCGTGGCGGCAACTGGATTTATGAGTGGGCCACCCAGGCGATGGTGTGGCAACAAAATGCCTGTGCGGAAGACGATCCGCAACTCTGTGGTCGTCACTGGCTGCATGCGGGTACGTTGTACAACATTGCCGCCTa                                                                                                                                <  2:1451961/150‑1 (MQ=255)
                      acgGTGATTGGCTACCGTGGCGGCAACTGGATTTATGAGTGGGCCACCCAGGCGATGGTGTGGCAACAAAAAGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCAATGTTGTACAACGTAgccgcatgc                                                                                                                               <  1:20488/150‑5 (MQ=255)
                       cgGTGATTGGCTACCGTGGCGGCAACTGGATTTATGAGTGGGCCACCCAGGCGATGGTGTGGCAACAAAAAGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATc                                                                                                                              >  1:1367578/1‑150 (MQ=255)
                       cgGTGATTGGCTACCGTGGCGGCAACTGGATTTATGAGTGGGCCACCCAGGCGATGGTGTGGCAACAAAAAGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATc                                                                                                                              <  2:1629431/150‑1 (MQ=255)
                          tGATTGGCTACCGTGGCGGCAACTGGATTTATGAGTGGGCCACCCAGGCGATGGTGTGGCAACAAAAAGCCTGTGCGGAAGACGATCCGCACCTCAGGGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCtg                                                                                                                           >  2:312748/1‑149 (MQ=255)
                          tGATTGGCTACCGTGGCGGCAACTGGATTTATGAGTGGGCCACCCAGGCGATGGTGTGGCAACAAAAAGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTc                                                                                                                           >  2:1531358/1‑150 (MQ=255)
                            aTTGGCTACCGTGGCGGCAACTGGATTTATGAGTGGGCCACCCAGGCGATGGTGTGGCAACAAAAAGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCAt                                                                                                                         >  1:1292241/1‑150 (MQ=255)
                               ggCTACCGTGGCGGCAACTGGATTTATGAGTGGGCCACCCAGGCGATGGTGTGGCAACAAAAAGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTg                                                                                                                      >  2:1668266/1‑150 (MQ=255)
                               ggCTACCGTGGCGGCAACTGGATTTATGAGTGGGCCACCCAGGCGATGGTGTGGCAACAAAAAGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTg                                                                                                                      >  1:803799/1‑150 (MQ=255)
                               ggCTACCGTGGCGGCAACTGGATTTATGAGTGGGCCACACAGGCGATGGTGTGGCAACAAAAAGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTg                                                                                                                      >  2:1225912/1‑150 (MQ=255)
                                    ccGTGGCGGCAACTGGATTTATGAGTGGGCCACCCAGGCGATGGTGTGGCAACAAAAAGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAgg                                                                                                                 >  1:858592/1‑150 (MQ=255)
                                       tGGCGGCAACTGGATTTATGAGTGGGCCACCCAGGCGATGGTGTGGCAACAAAAAGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGgaga                                                                                                              <  2:1542192/150‑1 (MQ=255)
                                        ggcggcAACTGGATTTATGAGTGGGCCACCCAGGCGATGGTGTGGCAACAAAAAGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTGACTGGCTGCATTCGGCTACGTTGTACAACATTGCCGCGTATCCTCATCTGATAGGAGAt                                                                                                             <  2:825432/150‑1 (MQ=255)
                                         gcggcAACTGGATTTATGAGTGGGCCACCCAGGCGATGGTGTGGCAACAAAAAGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATg                                                                                                            >  2:420878/1‑150 (MQ=255)
                                           ggcAACTGGATTTATGAGTGGGCCACCCAGGCGATGGTGTGGCAACAAAAAGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGAc                                                                                                          >  1:343028/1‑150 (MQ=255)
                                            gcAACTGGATTTATGAGTGGGCCACCCAGGCGATGGTGTGGCAACAAAAAGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGAcc                                                                                                         >  2:189833/1‑150 (MQ=255)
                                             cAACTGGATTTATGAGTGGGCCACCCAGGCGATGGTGTGGCAACAAAAAGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCt                                                                                                        >  1:34950/1‑150 (MQ=255)
                                               aCTGGATTTATGAGTGGGCCACCCAGGCGATGGTGTGGCAACAAAAAGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTgg                                                                                                      >  1:533548/1‑150 (MQ=255)
                                                cTGGATTTATGAGTGGGCCACCCAGGCGATGGTGTGGCAACAAAAAGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGc                                                                                                     <  1:835156/150‑1 (MQ=255)
                                                cTGGATTTATGAGTGGGCCACCCAGGCGATGGTGTGGCAACAAAAAGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGc                                                                                                     >  1:800287/1‑150 (MQ=255)
                                                cTGGATTTATGAGTGGGCCACCCAGGCGATGGTGTGGCAACAAAAAGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGc                                                                                                     <  1:914526/150‑1 (MQ=255)
                                                  ggATTTATGAGTGGGCCACCCAGGCGATGGTGTGGCAACAAAAAGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCg                                                                                                   >  1:1536617/1‑150 (MQ=255)
                                                   gATTTATGAGTGGGCCACGCAGGCGATGGTGCGGCAACAAAAAGCCTGTGCGGAAGACGATGCGCAACTCAGTGTTCGTCACTGGCTGGATGCGGCTACGTTGTTCAACATTGACGCCTATCCTGATCTGAAAGGAGATGACATGGCCGa                                                                                                  <  2:856444/150‑1 (MQ=255)
                                                   gATTTATGAGTGGGCCACCCAGGCGATGGTGTGGCAACAAAAAGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGa                                                                                                  >  2:793390/1‑150 (MQ=255)
                                                     tttATGAGTGGGCCACCCAGGCGATGGTGTGGCAACAAAAAGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGc                                                                                                >  1:174851/1‑150 (MQ=255)
                                                          gAGTGGGCCACCCAGGCGATGGTGTGGCAACAAAAAGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTGCGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAgcg                                                                                           >  1:1142450/1‑150 (MQ=255)
                                                            gTGGGCCACCCAGGCGATGGTGTTGCAACAAAAAGCCTGGGCGGAAGACGATCCGCAACTCAGTGGGCGTCACTGGCTGCAGGCGGCTACGTCGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCa                                                                                         <  1:1567845/150‑1 (MQ=255)
                                                            gTGGGCCACCCAGGCGATGGTGTGGCAACAAAAAGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCa                                                                                         >  2:241662/1‑150 (MQ=255)
                                                            gTGGGCCACCCAGGCGATGGTGTGGCAACAAAAAGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCa                                                                                         >  2:276477/1‑150 (MQ=255)
                                                               ggCCACCCAGGCGATGGTGTGGCAACAAAAAGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGc                                                                                      >  1:160510/1‑150 (MQ=255)
                                                               ggCCACCCAGGCGATGGTGTGGCAACAAAAAGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGc                                                                                      >  2:56392/1‑150 (MQ=255)
                                                                gCCACCCAGGCGATGGTGTGGCAACAAAAAGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCt                                                                                     >  2:437790/1‑150 (MQ=255)
                                                                gCCACCCAGGCGATGGTGTGGCAACAAAAAGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCt                                                                                     <  2:858592/150‑1 (MQ=255)
                                                                     ccAGGCGATGGTGTGGCAACAAAAAGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTGTc                                                                                <  1:1488193/150‑1 (MQ=255)
                                                                          cGATGGTGTGGCAACAAAAAGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTGTCAAAcc                                                                           <  1:1225912/150‑1 (MQ=255)
                                                                              ggTGTGGCAACAAAAAGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTGTCAAACCgcgc                                                                       >  1:1218818/1‑150 (MQ=255)
                                                                               gtgtGGCAACAAAAAGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTGTCAAACCGCGcc                                                                      <  2:1536617/150‑1 (MQ=255)
                                                                                   ggCAACAAAAAGCCTGTGCGGAAGCCGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTGTCAAACCGCGCCTATg                                                                  >  1:1434278/1‑150 (MQ=255)
                                                                                   ggCAACAAAAAGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTGTCAAACCGCGCCTATg                                                                  <  2:1596309/150‑1 (MQ=255)
                                                                                    gCAACAAAAAGCCTGTGCGGAAGACGCTCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTGTCAAACCGCGCCTATGa                                                                 >  2:582942/1‑150 (MQ=255)
                                                                                    gCAACAAAAAGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTGTCAAACCGCGCCTATGa                                                                 <  1:442203/150‑1 (MQ=255)
                                                                                      aacaaAAAGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTGTCAAACCGCGCCTATGAag                                                               <  2:693928/150‑1 (MQ=255)
                                                                                        caaAAAGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAg                                                                                             >  2:1646934/1‑118 (MQ=255)
                                                                                        caaAAAGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAg                                                                                             <  1:1646934/118‑1 (MQ=255)
                                                                                         aaaaaGCCTGTGCGGAAGACGATCCGCAACTCAGTCCTCGTCACTGGCTGCATGCCGCTACGTTGTACAACATTGCCGCCTATCCTCAACAGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTGTCAAACCGCGCCTATGAAGAgg                                                            >  2:17621/1‑150 (MQ=255)
                                                                                           aaaGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTGTCAAACCGCGCCTATGAAGAGGcc                                                          <  1:96132/150‑1 (MQ=255)
                                                                                           aaaGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTGTCAAACCGCGCCTATGAAGAGGcc                                                          >  1:1528952/1‑150 (MQ=255)
                                                                                           aaaGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTGTCAAACCGCGCCTATGAAGAGGcc                                                          <  1:1531358/150‑1 (MQ=255)
                                                                                                 tgtgCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCGGCATGCGGCTACGTTGTACAACATTGCCGCCTATGCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTGTCAAACCGCGCCTATGAAGAGGCCGCTCAg                                                    <  2:869294/150‑1 (MQ=255)
                                                                                                  gtgCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTGTCAAACCGCGCCTATGAAGAGGCCGCTCAGc                                                   >  1:1067413/1‑150 (MQ=255)
                                                                                                   tgCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTGTCAAACCGCGCCTATGAAGAGGCCGCTCAGCg                                                  >  1:661344/1‑150 (MQ=255)
                                                                                                     cGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGACTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTGTCAAACCGCGCCTATGAAGAGGCCGCTCAGCGTc                                                >  1:1161/1‑150 (MQ=255)
                                                                                                      ggAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTGTCAAACCGCGCCTATGAAGAGGCCGCTCAGCGTCt                                               >  2:1662587/1‑150 (MQ=255)
                                                                                                       gAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTGTCAAACCGCGCCTATGAAGAGGCCGCTCAGCGTCTa                                              >  1:1530444/1‑150 (MQ=255)
                                                                                                         aGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTGTCAAACCGCGCCTATGAAGAGGCCGCTCAGCGTCTAcc                                            >  2:680038/1‑150 (MQ=255)
                                                                                                         aGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTGTCAAACCGCGCCTATGAAGAGGCCGCTCAGCGTCTAcc                                            <  1:1471055/150‑1 (MQ=255)
                                                                                                             gATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTGTCAAACCGCGCCTATGAAGAGGCCGCTCAGCGTCTACCGGGc                                        >  2:109456/1‑150 (MQ=255)
                                                                                                             gATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTGTCAAACCGCGCCTATGAAGAGGCCGCTCAGCGTCTACCGGGc                                        <  1:965792/150‑1 (MQ=255)
                                                                                                                ccGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTGTCAAACCGCGCCTATGAAGAGGCCGCTCAGCGTCTACCGGGCACg                                     >  1:1454981/1‑150 (MQ=255)
                                                                                                                 cGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTGTCAAACCGCGCCTATGAAGAGGCCGCTCAGCGTCTACCGGGCACGa                                    <  2:361968/150‑1 (MQ=255)
                                                                                                                      cTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTGTCAAACCGCGCCTATGAAGAGGCCGCTCAGCGTCTACCGGGCACGATGCgg                               <  1:1561992/150‑1 (MQ=255)
                                                                                                                       tCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTGTCAAACCGCGCCTATGAAGAGGCCGCTCAGCGTCTACCGGGCACGACGCGGc                              <  2:1394224/150‑1 (MQ=255)
                                                                                                                        cAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTGTCAAACCGCGCCTATGAAGAGGCCGCTCAGCGTCTACCGGGCACGATGCGGCa                             <  1:793390/150‑1 (MQ=255)
                                                                                                                        cAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTGTCAAACCGCGCCTATGAAGAGGCCGCTCAGCGTCTACCGGGCACGATGCGGCa                             >  2:1340749/1‑150 (MQ=255)
                                                                                                                        cAGTGGTCGTCACTGGCTGCATCCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTGTCAAACCGCGCCTATGAAGAGGCCGCTCAGCGTCTACCGGGCACGATGCGGCa                             <  1:1668266/150‑1 (MQ=255)
                                                                                                                            ggTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTGTCAAACCGCGCCTATGAAGAGGCCGCTCAGCGTCTACCGGGCACGATGCGGCAGATg                         <  1:189833/150‑1 (MQ=255)
                                                                                                                              tcgtcACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTGTCAAACCGCGCCTATGAAGAGGCCGCTCAGCGTCTACCGGGCACGATGCGGCAGATGGa                       >  2:538211/1‑150 (MQ=255)
                                                                                                                                gtcACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTGTCAAACCGCGCCTATGAAGAGGCCGCTCAGCGTCTACCGGGCACGATGCGGCAGATGGAGt                     >  1:1330593/1‑150 (MQ=255)
                                                                                                                                   aCTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTGTCAAACCGCGCCTATGAAGAGGCCGCTCAGCGTCTACCGGGCACGATGCGGCAGATGGAGTTTa                  >  1:741544/1‑150 (MQ=255)
                                                                                                                                    cTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCTATCAAGCGCATGCTTTGTCAAACCTCGCCTATGAAGAGGCCGCTCAGCGTATACCGTGCACGATGCGGCAGATTGAGTTTAc                 <  2:803799/150‑1 (MQ=255)
                                                                                                                                      ggCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTGTCAAACCGCGCCTATGAAGAGGCCGCTCAGCGTCTACCGGGCACGATGCGGCAGATGGAGTTTACCg               <  2:533548/150‑1 (MQ=255)
                                                                                                                                      ggCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTGTCAAACCGCGCCTATGAAGAGGCCGCTCAGCGTCTACCGGGCACGATGCGGCAGATGGAGTTTACCg               >  1:392479/1‑150 (MQ=255)
                                                                                                                                      ggCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTGTCAAACCGCGCCTATGAAGAGGCCGCTCAGCGTCTACCGGGCACGATGCGGCAGATGGAGTTTACCg               <  2:1067413/150‑1 (MQ=255)
                                                                                                                                          gCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTGTCAAACCGCGCCTATGAAGAGGCCGCTCAGCGTCTACCGGGCACGATGCGGCAGATGGAGTTTACCGTAcc           >  2:53467/1‑150 (MQ=255)
                                                                                                                                          gCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTGTCAAACCGCGCCTATGAAGAGGCCGCTCAGCGTCTACCGGGCACGATGCGGCAGATGGAGTTTACCGTAcc           <  2:591932/150‑1 (MQ=255)
                                                                                                                                          gCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTGTCAAACCGCGCCTATGAAGAGGCCGCTCAGCGTCTACCGGGCACGATGCGGCAGATGGAGTTTACCGTAcc           >  2:59541/1‑150 (MQ=255)
                                                                                                                                          gCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTGTCAAACCGCGCCTATGAAGAGGCCGCTCAGCGTCTACCGGGCACGATGCGGCAGATGGAGTTTACCGTAcc           <  2:1488765/150‑1 (MQ=255)
                                                                                                                                          gCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTGTCAAACCGCGCCTATGAAGAGGCCGCTCAGCGTCTACCGGGCACGATGCGGCAGATGGAGTTTACCGTAcc           >  2:519860/1‑150 (MQ=255)
                                                                                                                                          gCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTGTCAAACCGCGCCTATGAAGAGGCCGCTCAGCGTCTACCGGGCACGATGCGGCAGATGGAGTTTACCGTAcc           <  2:800287/150‑1 (MQ=255)
                                                                                                                                             tGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTGTCAAACCGCGCCTATGAAGAGGCCGCTCAGCGTCTACCGGGCACGATGCGGCa                             >  2:580267/1‑129 (MQ=255)
                                                                                                                                             tGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTGTCAAACCGCGCCTATGAAGAGGCCGCTCAGCGTCTACCGGGCACGATGCGGCa                             <  1:580267/129‑1 (MQ=255)
                                                                                                                                             ggCGGCAATGTTGTACAACGTAGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTGTCAAACCGCGCCTATGAAGAGGCCGCTCAGCGTCTACCGGGCACGATGCGGCAGATGGAGTTTACCGTACCcgg        >  2:128746/2‑150 (MQ=255)
                                                                                                                                              gCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTGTCAAACCGCGCCTATGAAGAGGCCGCTCAGCGTCTACCGGGCACGATGCGGCAGATGGAGTTTACCGTACCcggc       >  2:406233/1‑150 (MQ=255)
                                                                                                                                              gCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTGTCAAACCGCGCCTATGAAGAGGCCGCTCAGCGTCTACCGGGCACGATGCGGCAGATGGAGTTTACCGTACCcggc       <  2:1218818/150‑1 (MQ=255)
                                                                                                                                                   tACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTGTCAAACCGCGCCTATGAAGAGGCCGCTCAGCGTCTACCGGGCACGATGCGGCAGATGGAGTTTACCGTACCCGGCGGTgc  <  1:1154975/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TACCGACGATGATTTATACGATACGGTGATTGGCTACCGTGGCGGCAACTGGATTTATGAGTGGGCCACCCAGGCGATGGTGTGGCAACAAAAAGCCTGTGCGGAAGACGATCCGCAACTCAGTGGTCGTCACTGGCTGCATGCGGCTACGTTGTACAACATTGCCGCCTATCCTCATCTGAAAGGAGATGACCTGGCCGAGCAAGCGCAGGCTTTGTCAAACCGCGCCTATGAAGAGGCCGCTCAGCGTCTACCGGGCACGATGCGGCAGATGGAGTTTACCGTACCCGGCGGTGC  >  NC_000913/256766‑257062

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: