Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,956,973:1 +A 6.1% coding (848/894 nt) ilvY ← DNA‑binding transcriptional dual regulator IlvY

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,956,9731.A 6.1% 593.2 / 33.7 147coding (848/894 nt)ilvYDNA‑binding transcriptional dual regulator IlvY
Reads supporting (aligned to +/- strand):  ref base . (65/73);  new base A (3/6);  total (68/79)
Fisher's exact test for biased strand distribution p-value = 5.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.80e-01

Only 100 of 147 total aligned reads displayed.

CCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGCTGTTTTCCAGCACCACTTCCGGCAACAATGCCACGCCACAGCCGAGTGC  >  NC_000913/3956827‑3957121
                                                                                                                                                   |                                                                                                                                                    
 cACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATg                                                                                                                                                  >  2:430538/1‑150 (MQ=255)
 cACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATg                                                                                                                                                  <  2:373288/150‑1 (MQ=255)
     acgaAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGc                                                                                                                                              <  2:60401/150‑1 (MQ=255)
        aaaCCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGTGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGc                                                                                                                                           <  2:551438/150‑1 (MQ=255)
        aaaCCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGc                                                                                                                                           >  2:1848495/1‑150 (MQ=255)
         aaacAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCt                                                                                                                                          <  2:1237536/147‑1 (MQ=255)
         aaCCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCt                                                                                                                                          >  2:1234222/1‑150 (MQ=255)
         aaCCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCt                                                                                                                                          >  2:1459241/1‑150 (MQ=255)
            cAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTACATGAGCCGCtaat                                                                                                                                       <  1:351217/150‑4 (MQ=255)
            cAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTACATGAGCCGCtaat                                                                                                                                       <  1:419499/150‑4 (MQ=255)
             aaTAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCttttt                                                                                                                                      >  1:1253346/1‑150 (MQ=255)
               tAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTACATGAGCCGCtaattga                                                                                                                                    <  1:1362811/150‑7 (MQ=255)
                aaCCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTACATGAGCCGCtaattgag                                                                                                                                   <  1:413706/150‑8 (MQ=255)
                aaCCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTg                                                                                                                                   >  1:1818461/1‑150 (MQ=255)
                 aCCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTgc                                                                                                                                  <  2:2406151/150‑1 (MQ=255)
                 aCCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTATgc                                                                                                                                  <  1:1603033/150‑1 (MQ=255)
                  cccGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTgcg                                                                                                                                 <  1:141063/150‑1 (MQ=255)
                  cccGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTgcg                                                                                                                                 <  2:1253346/150‑1 (MQ=255)
                   ccgtagTTAAGGTTCTTTTTGGCGGTTTAGAGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTgcgc                                                                                                                                >  2:1057551/6‑150 (MQ=255)
                    cGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGca                                                                                                                               >  2:2153638/1‑150 (MQ=255)
                     ggCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGcac                                                                                                                              <  1:274764/150‑1 (MQ=255)
                     ggCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGcac                                                                                                                              <  1:1149871/150‑1 (MQ=255)
                        gTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGcacaca                                                                                                                           <  1:1672166/150‑1 (MQ=255)
                        gTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGcacaca                                                                                                                           <  1:1338015/150‑1 (MQ=255)
                          tCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACg                                                                                                                         <  1:584898/150‑1 (MQ=255)
                           cAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGc                                                                                                                        >  2:1301433/1‑150 (MQ=255)
                            aGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGcc                                                                                                                       >  2:2088074/1‑150 (MQ=255)
                            aGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGcc                                                                                                                       <  1:467118/150‑1 (MQ=255)
                             ggTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGccc                                                                                                                      >  1:1307723/1‑150 (MQ=255)
                                tCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGGCTCAATTAGCGGCTCATGTAGCCACTTCTTCTGCGCACACACGCCCAGc                                                                                                                   <  1:180044/150‑1 (MQ=255)
                                   ttttGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCa                                                                                                                <  2:962364/150‑1 (MQ=255)
                                       ggCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACg                                                                                                            >  2:982100/1‑150 (MQ=255)
                                        gCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACgg                                                                                                           >  1:890739/1‑150 (MQ=255)
                                           ggTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGt                                                                                                        >  1:55161/1‑150 (MQ=255)
                                            gTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGCTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGtt                                                                                                       >  1:238900/1‑150 (MQ=255)
                                               aGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTtc                                                                                                    >  1:1045753/1‑150 (MQ=255)
                                                gggAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGcacaca                                                                                                                           >  2:1174981/1‑126 (MQ=255)
                                                gggAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGcacaca                                                                                                                           <  1:1174981/126‑1 (MQ=255)
                                                gggAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTtct                                                                                                   >  2:2077165/1‑150 (MQ=255)
                                                gggAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTtct                                                                                                   <  2:2005429/150‑1 (MQ=255)
                                                 ggAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTtctc                                                                                                  >  1:2363657/1‑150 (MQ=255)
                                                  gAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGcac                                                                                                                              >  2:1705502/1‑121 (MQ=255)
                                                  gAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCa                                                                                                 >  2:1816358/1‑150 (MQ=255)
                                                  gAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCa                                                                                                 >  1:875406/1‑150 (MQ=255)
                                                            tGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGct                                                                                       >  1:104346/1‑150 (MQ=255)
                                                            tGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGct                                                                                       <  1:2294363/150‑1 (MQ=255)
                                                             gATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGctc                                                                                      <  2:728634/150‑1 (MQ=255)
                                                                    cttcgcttATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAatc                                                                               >  2:2027849/1‑150 (MQ=255)
                                                                       cgcttATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAatcatc                                                                            <  2:1586947/150‑1 (MQ=255)
                                                                        gcttATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACCCGCCCAGCTCAAACGGCGTTTTTTCATCGCTGCGCTCTAAAatcatcc                                                                           >  2:967487/1‑149 (MQ=255)
                                                                        gcttATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAAtcatca                                                                           <  1:1848495/150‑1 (MQ=255)
                                                                        gcttATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAAtcatca                                                                           <  1:2249728/150‑1 (MQ=255)
                                                                        gcttATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCACTCTAAAAtcatca                                                                           >  1:2010980/1‑150 (MQ=255)
                                                                          ttATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCAcg                                                                         <  2:32745/150‑1 (MQ=255)
                                                                            aTCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCGGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCAcgcg                                                                       >  1:203975/1‑150 (MQ=255)
                                                                             tCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCgg                                                                      <  2:2294940/150‑1 (MQ=255)
                                                                             tCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCgg                                                                      <  2:1988421/150‑1 (MQ=255)
                                                                              cAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGt                                                                     <  2:1763720/150‑1 (MQ=255)
                                                                              cAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGt                                                                     >  2:1074299/1‑150 (MQ=255)
                                                                              cAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGt                                                                     <  2:2222811/150‑1 (MQ=255)
                                                                              cAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGt                                                                     >  1:216286/1‑150 (MQ=255)
                                                                                ggCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTa                                                                   >  1:1990430/1‑150 (MQ=255)
                                                                                 gCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTAc                                                                  >  1:68395/1‑150 (MQ=255)
                                                                                   cTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGc                                                                <  1:823509/150‑1 (MQ=255)
                                                                                    tACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCa                                                               <  1:201126/150‑1 (MQ=255)
                                                                                     aCCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCAc                                                              <  1:816919/150‑1 (MQ=255)
                                                                                     aCCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCAc                                                              <  2:1321627/150‑1 (MQ=255)
                                                                                     aCCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCAc                                                              >  2:646106/1‑150 (MQ=255)
                                                                                      cccGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCAcc                                                             <  1:721272/150‑1 (MQ=255)
                                                                                             gACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCgg                                                      >  1:280039/1‑150 (MQ=255)
                                                                                             gACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTATAAAATCATCACGCGTTTAAGCACCTGTTCgg                                                      <  2:100478/150‑1 (MQ=255)
                                                                                              aCAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCggg                                                     <  2:738736/150‑1 (MQ=255)
                                                                                                aaCGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGc                                                   <  2:890739/150‑1 (MQ=255)
                                                                                                 aCGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGCt                                                  <  2:2363657/150‑1 (MQ=255)
                                                                                                   gTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGCTGt                                                <  1:2169964/150‑1 (MQ=255)
                                                                                                     cATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTACATGAGCCGCtaa                                                                                                                                        <  1:72171/60‑3 (MQ=37)
                                                                                                         tgtgGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGCTGTTTTCCa                                          <  1:182435/150‑1 (MQ=255)
                                                                                                            ggTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGCTGTTTTCCAGca                                       >  1:999403/1‑150 (MQ=255)
                                                                                                            ggTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGCTGTTTTCCAGca                                       <  1:393635/150‑1 (MQ=255)
                                                                                                             gTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGCTGTTTTCCAGcac                                      <  2:1885933/150‑1 (MQ=255)
                                                                                                             gTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGCTGTTTTCCAGcac                                      <  2:560350/150‑1 (MQ=255)
                                                                                                                 ggCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGCTGTTTTCCAGCACCACt                                  >  2:422846/1‑150 (MQ=255)
                                                                                                                 ggCAAAATCTTCCAGAATCCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGCTGTTTTCCAGCACCACt                                  >  2:1184850/1‑150 (MQ=255)
                                                                                                                  gCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGCTGTTTTCCAGCACCACtt                                 >  1:677839/1‑150 (MQ=255)
                                                                                                                   cAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGCTGTTTTCCAGCACCACTTc                                >  1:1116706/1‑150 (MQ=255)
                                                                                                                       aTCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGCTGTTTTCCAGCACCACTTCCGGc                            <  1:1459241/150‑1 (MQ=255)
                                                                                                                        tCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGCTGTTTTCCAGCACCACTTCCGGca                           <  2:1538533/150‑1 (MQ=255)
                                                                                                                          ttCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGCTGTTTTCCAGCACCACTTCCGGcaac                         <  1:1216737/150‑1 (MQ=255)
                                                                                                                              aGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGCTGTTTTCCAGCACCACTTCCGGCAACAATg                     <  1:1301244/150‑1 (MQ=255)
                                                                                                                                aaTGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGCTGTTTTCCAGCACCACTTCCGGCAACAATgcc                   <  1:1057551/150‑1 (MQ=255)
                                                                                                                                aaTGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGCTGTTTTCCAGCACCACTTCCGGCAACAATgcc                   <  1:2077165/150‑1 (MQ=255)
                                                                                                                                    ccTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGCTGTTTTCCAGCACCACTTCCGGCAACAATgccacgc               >  1:1671785/1‑150 (MQ=255)
                                                                                                                                     cTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGCTGTTTTCCAGCACCACTTCCGGCAACAATgccacgcc              <  1:2442810/150‑1 (MQ=255)
                                                                                                                                     cTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGCTGTTTTCCAGCACCACTTCCGGCAACAATgccacgcc              >  1:467759/1‑150 (MQ=255)
                                                                                                                                       gaaaaaGCGGCTACATGAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGCTGTTTTCCAGCACCACTTCCGGCAACAATgccacgccac            >  2:1524174/6‑150 (MQ=255)
                                                                                                                                       cAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCTGCTCAAACGGCGTTTTCTGATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGCTGTTTTCCAGCACCACTTCCGGGAACAATgccaggccac            <  2:1506542/150‑1 (MQ=255)
                                                                                                                                         aaaaGCGGCTACATGAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGCTGTTTTCCAGCACCACTTCCGGCAACAATGCCACGCCACAg          >  1:804155/4‑150 (MQ=255)
                                                                                                                                         aTTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGCTGTTTTCCAGCACCACTTCCGGCAACAATGCCACGCCACAg          >  1:1953201/1‑150 (MQ=255)
                                                                                                                                                gCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGCTGTTTTCCAGCACCACTTCCGGCAACAATGCCACGCCACAGCCGAGTg   <  1:1234222/150‑1 (MQ=255)
                                                                                                                                                 cTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTATCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGCTGTTTTCCAGCACCACTTCCGGCAACAATGCCACGCCACAGCCGAGTGc  >  2:1296425/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CCACGACGAAACCAATAACCCGGCGTTCAGGTTCTTTTTGGCGGGTTAGGGAAAAATGCCTGATAGCGCTTCGCTTATCAGGCCTACCCGCGCGACAACGTCATTTGTGGTTCGGCAAAATCTTCCAGAATGCCTCAATTAGCGGCTCATGTAGCCGCTTTTTCTGCGCACACACGCCCAGCTCAAACGGCGTTTTCTCATCGCTGCGCTCTAAAATCATCACGCGGTTACGCACCGGTTCGGGGCTGTTTTCCAGCACCACTTCCGGCAACAATGCCACGCCACAGCCGAGTGC  >  NC_000913/3956827‑3957121

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: