Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 176,979 C→T 8.7% intergenic (+25/+22) erpA → / ← yadS iron‑sulfur cluster insertion protein ErpA/conserved inner membrane protein YadS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913176,9790CT 8.7% 202.1 / 7.7 69intergenic (+25/+22)erpA/yadSiron‑sulfur cluster insertion protein ErpA/conserved inner membrane protein YadS
Reads supporting (aligned to +/- strand):  ref base C (30/33);  new base T (3/3);  total (33/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.33e-01

TGAGCCTGCAATATCTGGTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGAAAATGTGT  >  NC_000913/176830‑177128
                                                                                                                                                     |                                                                                                                                                     
tGAGCCTGCAATATCTGGTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACTTGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGTATAAAc                                                                                                                                                       <  2:136333/150‑1 (MQ=255)
tGAGCCTGCAATATCTGGTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAAc                                                                                                                                                       <  1:615027/150‑1 (MQ=255)
     cTGCAATATCTGGTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATcc                                                                                                                                                  >  2:935875/1‑150 (MQ=255)
     cTGCAATATCTGGTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATcc                                                                                                                                                  <  1:1531494/150‑1 (MQ=255)
      tGCAATATCTGGTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCg                                                                                                                                                 >  2:974800/1‑150 (MQ=255)
      tGCAATATCTGGTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCg                                                                                                                                                 <  2:1172707/150‑1 (MQ=255)
      tGCAATATCTGGTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCg                                                                                                                                                 <  1:219381/150‑1 (MQ=255)
        cAATATCTGGTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGTTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGc                                                                                                                                               <  1:1125022/150‑1 (MQ=255)
        cAATATCTGGTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAGTTTATCCGGc                                                                                                                                               <  2:925508/150‑1 (MQ=255)
         aaTATCTGGTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTAGCTAATCTGTTGTTCATGTGCCGGATAGTTTATCCGGca                                                                                                                                              <  1:1191669/150‑1 (MQ=255)
          atatCTGGTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAGTTTATCCGGcac                                                                                                                                             <  1:20361/150‑1 (MQ=255)
          atatCTGGTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGcac                                                                                                                                             <  1:670601/150‑1 (MQ=255)
           tatCTGGTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGcaca                                                                                                                                            <  1:273090/150‑1 (MQ=255)
             tCTGGTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACAtc                                                                                                                                          <  1:904737/150‑1 (MQ=255)
                ggTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACAtcttc                                                                                                                                       >  1:1500629/1‑150 (MQ=255)
                ggTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACAtcttc                                                                                                                                       <  2:247202/150‑1 (MQ=255)
                  tCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTtctc                                                                                                                                     <  1:935875/150‑1 (MQ=255)
                    ggcggTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTtctctc                                                                                                                                   <  1:598552/150‑1 (MQ=255)
                     gcggTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTctctct                                                                                                                                  <  2:1104711/150‑1 (MQ=255)
                           ccGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACg                                                                                                                            >  2:1449342/1‑150 (MQ=255)
                              ttGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGccc                                                                                                                         >  2:1086849/1‑150 (MQ=255)
                              ttGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGGGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGccc                                                                                                                         <  2:875976/150‑1 (MQ=255)
                                  ttATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATtc                                                                                                                     <  1:1526703/150‑1 (MQ=255)
                                   tataCCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATtct                                                                                                                    <  1:882920/150‑1 (MQ=255)
                                    ataCCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATtctc                                                                                                                   >  1:1603834/1‑150 (MQ=255)
                                      aCCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCAt                                                                                                                 >  1:328339/1‑150 (MQ=255)
                                       ccGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATc                                                                                                                >  2:585957/1‑150 (MQ=255)
                                           aGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAgc                                                                                                            >  1:1249097/1‑150 (MQ=255)
                                              tCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCa                                                                                                         <  2:757718/150‑1 (MQ=255)
                                               cTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGTATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCGGCGCaa                                                                                                        <  2:428254/150‑1 (MQ=255)
                                                tGGAAGGTTCTCGTTTCATCGTGACCAATCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCTGCGCaaa                                                                                                       <  2:299016/150‑1 (MQ=255)
                                                 ggAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAAc                                                                                                      >  1:1570663/1‑150 (MQ=255)
                                                   aaGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGt                                                                                                    <  2:951726/150‑1 (MQ=255)
                                                       ttCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGt                                                                                                <  1:1333834/150‑1 (MQ=255)
                                                           cGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCt                                                                                            >  1:1454112/1‑150 (MQ=255)
                                                                 aacGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGCCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGAt                                                                                      >  2:1516358/3‑150 (MQ=255)
                                                                    gTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGcc                                                                                   >  1:847425/1‑150 (MQ=255)
                                                                        ccAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACg                                                                               <  1:46429/150‑1 (MQ=255)
                                                                           aCCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAAt                                                                            >  2:760050/1‑150 (MQ=255)
                                                                             ccGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAg                                                                          >  2:309699/1‑150 (MQ=255)
                                                                                    cgAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAg                                                                   >  2:1269252/1‑150 (MQ=255)
                                                                                          gCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAAt                                                             <  1:939704/150‑1 (MQ=255)
                                                                                          gCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAAt                                                             >  1:1613284/1‑150 (MQ=255)
                                                                                               tgcggttgcggtTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATa                                                        <  2:1570663/150‑1 (MQ=255)
                                                                                                gcggttgcggttCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAg                                                       <  1:1219252/150‑1 (MQ=255)
                                                                                                 cggttgcggttCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGc                                                      <  1:1449342/150‑1 (MQ=255)
                                                                                                   gttgcggttCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGt                                                    >  1:814440/1‑150 (MQ=255)
                                                                                                   gttgcggttCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGt                                                    >  2:215799/1‑150 (MQ=255)
                                                                                                         gttCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGAc                                              <  2:328339/150‑1 (MQ=255)
                                                                                                           tcttcCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCa                                            <  1:309699/150‑1 (MQ=255)
                                                                                                              tcCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGc                                         <  2:1609163/150‑1 (MQ=255)
                                                                                                                   tAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCcatc                                    >  1:23611/1‑150 (MQ=255)
                                                                                                                    aGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCcatca                                   >  2:1027919/1‑150 (MQ=255)
                                                                                                                       aTCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATAc                                >  2:124748/1‑150 (MQ=255)
                                                                                                                         cTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTg                              >  1:1229497/1‑150 (MQ=255)
                                                                                                                         cTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTg                              >  1:1285081/1‑150 (MQ=255)
                                                                                                                           aaTCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGc                            <  1:1269252/150‑1 (MQ=255)
                                                                                                                            aTCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCt                           >  1:431489/1‑150 (MQ=255)
                                                                                                                             tctcTTGATCATCTGCCGGATTCACTATCCGTCACATCTTCTCGGTGCAACGGCCCTTTTCATCCTGCGCAAACGTCGGTAGCTTAATATGCCAACGAATAGCCGCCATCCGACTCAATAGCGACACGACCATGCCCATCATACTGGCTg                          <  2:688682/150‑1 (MQ=255)
                                                                                                                             tCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCGATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTg                          >  1:741591/1‑150 (MQ=255)
                                                                                                                             tCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTg                          <  2:1613284/150‑1 (MQ=255)
                                                                                                                               tgttgtTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGtt                        >  2:486762/1‑150 (MQ=255)
                                                                                                                                 gagaagaTGTGCCGGATAGTTTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTc                      >  1:342750/7‑150 (MQ=255)
                                                                                                                                 gagaagaTGTGCCGGATAGTTTATCCGGCACATCTTCTCTATTCAACGCCCATTCTCCTCCAGCGCAAACGTCGGTAGCTTAAGATGGCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTAGCTGttta                      >  1:1467260/7‑149 (MQ=255)
                                                                                                                                  tgttCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTcc                     <  2:1225894/150‑1 (MQ=255)
                                                                                                                                    ttCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAg                   <  2:1454112/150‑1 (MQ=255)
                                                                                                                                     tCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGt                  <  1:1086849/150‑1 (MQ=255)
                                                                                                                                        tgtgCCGGATAGTTTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGt               >  2:146771/1‑150 (MQ=255)
                                                                                                                                         gtgCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCCTCCAGCGCAAACGTCGGTAGCTTAAGATGCCACCGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCTTGCCCATCATACTGGCTGTTTCCAGTGGTa              >  2:896377/1‑150 (MQ=255)
                                                                                                                                           gCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCTGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACg            >  1:623325/1‑150 (MQ=255)
                                                                                                                                                 tAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGAAAAt      >  1:999398/1‑150 (MQ=255)
                                                                                                                                                     cTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCGATCATACTGTCTGTTTCCAGTGGTATGGAAAATgtgt  <  2:1249097/150‑1 (MQ=255)
                                                                                                                                                     cTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGAAAATgtgt  >  2:582238/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGAGCCTGCAATATCTGGTCGGCGGTTCCGTTGATTATACCGAAGGTCTGGAAGGTTCTCGTTTCATCGTGACCAACCCGAACGCGAAAAGCACCTGCGGTTGCGGTTCTTCCTTTAGTATCTAATCTGTTGTTCATGTGCCGGATAAACTATCCGGCACATCTTCTCTCTTCAACGCCCATTCTCATCCAGCGCAAACGTCGGTAGCTTAAGATGCCAACGAATAGCCGCCAGCCGAATCAATAGCGTCACGACCATGCCCATCATACTGGCTGTTTCCAGTGGTACGGAAAATGTGT  >  NC_000913/176830‑177128

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: