Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 994,637 A→C 6.4% S66A (TCC→GCC)  ssuC ← aliphatic sulfonate ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913994,6370AC 6.4% 217.1 / 4.9 82S66A (TCC→GCC) ssuCaliphatic sulfonate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (46/29);  new base C (2/3);  total (50/32)
Fisher's exact test for biased strand distribution p-value = 3.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.87e-01

CAGCGGGATCAGCGCCAGATGCGGCACGTTGCGCAACATCTGAATTGAGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGTTGCCACACCGCCACGATGCCCACCGGTAAAA  >  NC_000913/994488‑994786
                                                                                                                                                     |                                                                                                                                                     
cAGCGGGATCAGCGCCAGATGCGGCACGTTGCGCAACATCTGAATTGAGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGa                                                                                                                                                       >  2:269266/1‑150 (MQ=255)
    gggATCAGCGCCAGATGCGGCACGTTGCGCAACATCTGAATTGAGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCTCCTGGCGGGCGCTg                                                                                                                                                   <  1:480368/150‑1 (MQ=255)
       aTCAGCGCCAGATGCGGCACGTTGCGCAACATCTGAATTGAGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCTCCTGGCGGGCGCTGAtt                                                                                                                                                <  1:1188566/150‑2 (MQ=255)
       aTCAGCGCCAGATGCGGCACGTTGCGCAACATCTGAATTGAGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATc                                                                                                                                                >  2:979564/1‑150 (MQ=255)
        tCAGCGCCAGATGCGGCACGTTGCGCAACATCTGAATTGAGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCg                                                                                                                                               >  2:648551/1‑150 (MQ=255)
         cAGCGCCAGATGCGGCACGTTGCGCAACATCTGAATTGAGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCTCCTGGCGGGCGCTGAttgg                                                                                                                                              <  1:1042183/150‑4 (MQ=255)
         cAGCGCCAGATGCGGCACGTTGCGCAACATCTGAATTGAGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGc                                                                                                                                              >  2:785249/1‑150 (MQ=255)
           gcgcCAGATGCGGCACGTTGCGCAACATCTGAATTGAGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCa                                                                                                                                            >  1:495400/1‑150 (MQ=255)
            cgcCAGATGCGGCACGTTGCGCAACATCTGAATTGAGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCACCCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCac                                                                                                                                           >  1:453869/1‑149 (MQ=255)
            cgcCAGATGCGGCACGTTGCGCAACATCTGAATTGAGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAg                                                                                                                                           >  2:303823/1‑150 (MQ=255)
             gcCAGATGCGGCACGTTGCGCAACATCTGAATTGAGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCCACATCGACCCCAGCGATCCGCCAATTGAAAAGCCCATCAGCGACCGCCAGGTGCTGATCGCCAGa                                                                                                                                          >  2:191565/1‑150 (MQ=255)
                    gCGGCACGTTGCGCAACATCTGAATTGAGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGcc                                                                                                                                   <  2:1225683/150‑1 (MQ=255)
                       gCACGTTGCGCAACATCTGAATTGAGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCcaca                                                                                                                                >  2:534701/1‑150 (MQ=255)
                          cGTTGCGCAACATCTGAATTGAGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCCACCCCAGCGATCCGCCCATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGtt                                                                                                                             >  2:74797/1‑150 (MQ=255)
                             tGCGCAACATCTGAATTGATGTATACAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCACTCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTcgc                                                                                                                          <  1:138741/150‑1 (MQ=255)
                             tGCGCAACATCTGAATTGAGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTcgc                                                                                                                          <  1:143640/150‑1 (MQ=255)
                               cgcAACATCTGAATTGAGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCACCCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTcgccc                                                                                                                        >  1:258969/1‑149 (MQ=255)
                               cgcAACATCTGAATTGAGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCCGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTcgccg                                                                                                                        >  2:32158/1‑150 (MQ=255)
                                   aCATCTGAATTGAGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTgg                                                                                                                    >  1:496279/1‑150 (MQ=255)
                                   aCATCTGAATTGAGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTgg                                                                                                                    <  1:1237506/150‑1 (MQ=255)
                                        tGAATTGAGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGgag                                                                                                               >  2:516998/1‑150 (MQ=255)
                                             tGAGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGt                                                                                                          >  2:735265/1‑150 (MQ=255)
                                               aGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTcc                                                                                                        <  2:399480/150‑1 (MQ=255)
                                               aGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTcc                                                                                                        <  1:333679/150‑1 (MQ=255)
                                               aGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTcc                                                                                                        >  1:1058682/1‑150 (MQ=255)
                                                ggTATCGAGCAGCAGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCCGCGATCCGCCAATTGAAAAGCCTATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCa                                                                                                       <  2:484813/150‑1 (MQ=255)
                                                ggTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCa                                                                                                       <  1:1240465/150‑1 (MQ=255)
                                                     ccAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACg                                                                                                  >  2:1109676/1‑150 (MQ=255)
                                                       agcagcCGTTCTCCCCAGCGAGATAGCCCGCTACTCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGCAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGcc                                                                                                >  2:623981/1‑150 (MQ=255)
                                                               gtcaccccAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCACCGATCCGGCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGAAGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTcaccac                                                                                        <  1:1002276/146‑1 (MQ=255)
                                                                 ctcCCAAGCGAGATAGCCAGCTAATCAGTCCCAAAATCAACCCCACCGATCCGCCAAATGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTGGCCGCTGGCGGAGAGGGTCCAGCCCGCCGTCACCACTc                                                                                      >  2:150666/1‑150 (MQ=255)
                                                                   ccccAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCCCCACTCCt                                                                                    >  1:647256/1‑150 (MQ=255)
                                                                   ccccAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCt                                                                                    <  1:419529/150‑1 (MQ=255)
                                                                      cAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTcc                                                                                 >  2:680587/1‑150 (MQ=255)
                                                                        gCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCCCAGTTCCCCGCTGGCGGAGAGCGTCCAGAACGCCGTCCCCACTCCTTCCgg                                                                               >  2:1145642/1‑150 (MQ=255)
                                                                        gCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCgg                                                                               <  2:1283830/150‑1 (MQ=255)
                                                                        gCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCgg                                                                               <  1:1113118/150‑1 (MQ=255)
                                                                          gcgaTAGCCCGCTAATCAGTCCCAAAATCAACCCCACCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCCCAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCATCCGGTg                                                                             <  2:1119119/148‑1 (MQ=255)
                                                                          gagaTAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTg                                                                             <  1:1203682/150‑1 (MQ=255)
                                                                           agaTAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAACCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGa                                                                            <  1:1119041/150‑1 (MQ=255)
                                                                            gaTAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGaa                                                                           <  1:724780/150‑1 (MQ=255)
                                                                             aTAGCCCGCTTATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATGAGCGCCCGCCAGTAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACGACTCCTTCCGGTGAAg                                                                          <  2:508868/150‑1 (MQ=255)
                                                                             aTAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAg                                                                          <  2:519425/150‑1 (MQ=255)
                                                                                gCCCGCTAATCAGTCCCAAAATCAACCCCCGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCCCCACTCCTTCCGGGGAAGGcc                                                                       >  1:620822/1‑149 (MQ=255)
                                                                                gCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCa                                                                       >  1:34213/1‑150 (MQ=255)
                                                                                gCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCa                                                                       >  2:480548/1‑150 (MQ=255)
                                                                                gCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCa                                                                       >  1:1155344/1‑150 (MQ=255)
                                                                                 cccGCTAATCAGGCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCGGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCaa                                                                      <  1:1316589/150‑1 (MQ=255)
                                                                                     cTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATa                                                                  <  1:269266/150‑1 (MQ=255)
                                                                                          cAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGt                                                             >  2:477799/1‑150 (MQ=255)
                                                                                              cccAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCTGGTGAAGGCAAAATACGCGTCGAc                                                         >  1:447409/1‑150 (MQ=255)
                                                                                               ccAAAATCAACCCCAGCGCTCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGCGCGTCCAGAACGCCGTCACCACTCCTTCCGGGGAAGGCAAAATACGCGTCGAcc                                                        >  2:304202/1‑149 (MQ=255)
                                                                                               ccAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGAcc                                                        >  2:277297/1‑149 (MQ=255)
                                                                                                cAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGAcag                                                       <  2:647256/150‑1 (MQ=255)
                                                                                                 aaaaTCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGCcagc                                                      >  2:961216/1‑150 (MQ=255)
                                                                                                   aaTCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGAcagcca                                                    <  2:184904/150‑1 (MQ=255)
                                                                                                    aTCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGGGCTGATCGCCAGATGCTGCAACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCGTGTGCAGGCAAAATACGCGTCGAcagccag                                                   >  2:1167636/1‑150 (MQ=255)
                                                                                                    aTCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGAcagccag                                                   >  1:302837/1‑150 (MQ=255)
                                                                                                    aTCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGAcagccac                                                   >  2:479478/1‑149 (MQ=255)
                                                                                                          cccAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACCGCCAGCCAAcc                                             >  2:30258/1‑150 (MQ=255)
                                                                                                          cccAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAAcc                                             >  2:77077/1‑150 (MQ=255)
                                                                                                          cccAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTGCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCCCTCCTTCCGGTGAAGGCATACTACGCGACGCCAGCCAGCCAAcc                                             >  2:215651/1‑150 (MQ=255)
                                                                                                            cAGCGATCCGCCAAGTGAAACGCCAATCAGCGCCCGCCAGGAGCGGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGCGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCga                                           <  1:1175630/150‑1 (MQ=255)
                                                                                                                   ccGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCCACCGAGGAGGcc                                    >  2:730798/1‑150 (MQ=255)
                                                                                                                      ccAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGt                                 >  1:1214279/1‑150 (MQ=255)
                                                                                                                       cAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGtt                                >  2:1126936/1‑150 (MQ=255)
                                                                                                                         aTTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCCGCCAACCGAGGAGGCCAGTTGc                              >  1:972322/1‑150 (MQ=255)
                                                                                                                          ttGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCCCCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGTTGcc                             >  1:497874/1‑150 (MQ=255)
                                                                                                                          ttGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGTTGcc                             <  2:729963/150‑1 (MQ=255)
                                                                                                                           tGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACCGCCAGCCCACCGAGGAGGCCAGTTGccc                            >  1:1435/1‑149 (MQ=255)
                                                                                                                             aaaaGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGTTGCcaca                          >  1:387874/1‑150 (MQ=255)
                                                                                                                              aaaGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGTTGCCacac                         <  1:648551/150‑1 (MQ=255)
                                                                                                                               aaGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGTTGCCACAcc                        <  2:495400/150‑1 (MQ=255)
                                                                                                                                aGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGTTGCCACACCg                       >  1:1047991/1‑150 (MQ=255)
                                                                                                                                 ggcgaTCAGCTCCTGGCGGGCGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCACCCGAGGAGGCCAGTTGCCACACCGc                      >  2:647886/5‑150 (MQ=255)
                                                                                                                                 gCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGTTGCCACCCCGc                      >  2:36688/1‑150 (MQ=255)
                                                                                                                                   cgaTCAGCTCCTGGCGGGCGCTGATCGCCAGATGCTGCCACAGTTCACCGCTGGCGGAGCGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGTTGCCACACCGccc                    >  2:1135920/3‑149 (MQ=255)
                                                                                                                                        aGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGCCAAAAGACGCGTCGACAGCCCGCCAACCGAGGAGGCCCGGTGCCACCCCGCCACGATg               >  2:784463/1‑150 (MQ=255)
                                                                                                                                              cgccAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGTTGCCACACCGCCACGATGCCCAcc         <  2:557747/150‑1 (MQ=255)
                                                                                                                                                ccAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGTTGCCACACCGCCACGATGCCCACCgg       >  2:1273107/1‑150 (MQ=255)
                                                                                                                                                ccAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGTTGCCACACCGCCACGATGCCCACCgg       <  1:516998/150‑1 (MQ=255)
                                                                                                                                                ccAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGTTGCCACACCGCCACGATGCCCACCgg       >  1:969670/1‑150 (MQ=255)
                                                                                                                                                  aGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGTTGCCACACCGCCACGATGCCCACCGGTa     <  2:118877/150‑1 (MQ=255)
                                                                                                                                                     aGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGTTGCCACACCGCCACGATGCCCACCGGTaaaa  >  1:837209/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CAGCGGGATCAGCGCCAGATGCGGCACGTTGCGCAACATCTGAATTGAGGTATCCAGCAGCCGTTCTCCCCAGCGAGATAGCCCGCTAATCAGTCCCAAAATCAACCCCAGCGATCCGCCAATTGAAAAGCCAATCAGCGCCCGCCAGGAGCTGATCGCCAGATGCTGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGTTGCCACACCGCCACGATGCCCACCGGTAAAA  >  NC_000913/994488‑994786

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: