Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,194,609 T→A 7.2% intergenic (+84/+112) yqiK → / ← sibD protein YqiK/small RNA SibD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,194,6090TA 7.2% 273.4 / 8.0 92intergenic (+84/+112)yqiK/sibDprotein YqiK/small RNA SibD
Reads supporting (aligned to +/- strand):  ref base T (46/39);  new base A (4/3);  total (50/42)
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.85e-01

GCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAAGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGCATTGCGTGGGGTCTGGCGCTCATCAACAATCGGGGAGTCAGCAGGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGGCCTTACAGAAGG  >  NC_000913/3194461‑3194757
                                                                                                                                                    |                                                                                                                                                    
gcTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAAGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGgtt                                                                                                                                                     <  1:1570104/150‑1 (MQ=255)
 cTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAAGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGgttg                                                                                                                                                    <  1:594071/150‑1 (MQ=255)
 cTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAAGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGgttg                                                                                                                                                    >  1:1250236/1‑150 (MQ=255)
   cACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAAGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGgttgtt                                                                                                                                                  >  1:1609473/1‑150 (MQ=255)
    aCTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAAGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGttgttg                                                                                                                                                 <  1:1219366/150‑1 (MQ=255)
      tGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAAGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTgca                                                                                                                                               <  1:1250400/150‑1 (MQ=255)
           gcATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAAGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGgt                                                                                                                                                      <  1:1485123/138‑1 (MQ=255)
           gcATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAAGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGgt                                                                                                                                                      >  2:1485123/1‑138 (MQ=255)
            cATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAAGCAGAGTAAACCCTCATCCAAGAATAATCCCTCGACAAAAAAGCCGGCCAATTTTAATATCTGCGTTCCGCAAAAAAGTCCAAATGATCAGCTTGTTGCAGCAatt                                                                                                                                         >  2:1601755/1‑150 (MQ=255)
             aTTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAAGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCAtttt                                                                                                                                        >  1:1347106/1‑150 (MQ=255)
                  cTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAAGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTg                                                                                                                                   >  1:75376/1‑150 (MQ=255)
                     cACCCTTAACCTCAACAACTCCCGTCGAAGAAAAAGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACg                                                                                                                                >  2:297429/1‑150 (MQ=255)
                     cACCCTTAACCTCAACAACTCCCGTCGAAGAAAAAGCAGAGTAAACCCTCTTCCAGAAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGCTCAGGAAGTTGCAGCATTTTGCGTGAcc                                                                                                                                >  2:673033/1‑149 (MQ=255)
                     cACCCTTAACCTCAACAACTCCCGTCGAAGAAAAAGCAGAGTAAACCCTCTTCCAGAAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATCATCAGGAGGTTGCAGCATTTTGCGTGAcc                                                                                                                                >  2:673036/1‑149 (MQ=255)
                       cccTTAACCTCAACAACTCCCGTCGAAGAAAAAGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGc                                                                                                                              >  2:296388/1‑150 (MQ=255)
                        ccTTAACCTCAACAACTCCCGTCGAAGAAAAAGCAGAGTAAACCCTCATCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCt                                                                                                                             >  2:1300369/1‑150 (MQ=255)
                              ccTCAACAACTCCCGTCGAAGAAAAAGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAGGTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGAc                                                                                                                                 <  1:1201940/140‑1 (MQ=18)
                              ccTCAACAACTCCCGTCGAAGAAAAAGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGCTCAGGACGTTGCAGCGTTTTGCGTGAc                                                                                                                                 >  2:1201940/1‑140 (MQ=255)
                                 caacaaCTCCCGTCGAAGAAAAAGCAGAGTAAACGCTCTTCCATCAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAGGTGCAAATGATCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAgg                                                                                                                    <  1:673036/150‑1 (MQ=17)
                                 aaacaacTCCCGTCGAAGAAAAAGCAGATTAAACCCTCTTCCAGAAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAGGTGCAAATGATCAGGACGTTGCAGCGTTTTGCGTGACCGCTCGGGGAAgg                                                                                                                    <  1:673033/149‑1 (MQ=17)
                                  aacaacTCCCGTCGAAGAAAAAGCAGGGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGc                                                                                                                   <  2:122573/150‑1 (MQ=255)
                                  aacaacTCCCGTCGAAGAAAAAGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGc                                                                                                                   <  2:905793/150‑1 (MQ=255)
                                    caacTCCCGTCGAAGAAAAAGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCaa                                                                                                                 <  2:1007510/150‑1 (MQ=255)
                                          ccGTCGAAGAAAAAGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGc                                                                                                           <  1:326050/150‑1 (MQ=255)
                                           cggcGAAGAAAAAGCAGACTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGAGGGCTCGGGGAAGGCAAAACTgga                                                                                                          <  2:1635579/147‑3 (MQ=255)
                                           cGTCGAAGAAAAGGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCa                                                                                                          >  1:655768/1‑150 (MQ=255)
                                              cGAAGAAAAAGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCa                                                                                                       >  1:1182086/1‑150 (MQ=255)
                                              cGAAGAAAAAGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCa                                                                                                       <  2:1260849/150‑1 (MQ=255)
                                               gaagaaAAAGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCag                                                                                                      >  2:270538/1‑150 (MQ=255)
                                                aagaaAAAGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCaga                                                                                                     >  1:187801/1‑150 (MQ=255)
                                                aagaaAAAGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCaga                                                                                                     <  1:732000/150‑1 (MQ=255)
                                                 agaaAAAGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCagag                                                                                                    >  2:237828/1‑150 (MQ=255)
                                                 agaaAAAGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCagag                                                                                                    >  1:1668397/1‑150 (MQ=255)
                                                   aaaaaGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGaa                                                                                                  >  1:114764/1‑150 (MQ=255)
                                                     aaaGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAg                                                                                                <  2:369139/150‑1 (MQ=255)
                                                      aaGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAgg                                                                                                                    >  2:1027497/1‑129 (MQ=255)
                                                      aaGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAgg                                                                                                                    <  1:1027497/129‑1 (MQ=255)
                                                      aaGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGc                                                                                               >  1:33715/1‑150 (MQ=255)
                                                       aGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGCa                                                                                              <  1:869739/150‑1 (MQ=255)
                                                        gCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGCAt                                                                                             >  2:782351/1‑150 (MQ=255)
                                                        gCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGATGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGCAt                                                                                             >  1:490807/1‑150 (MQ=255)
                                                              tAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGCATTGCGTg                                                                                       <  2:1311836/150‑1 (MQ=255)
                                                               aaaCCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGCATTGCGTgg                                                                                      >  1:155867/1‑150 (MQ=255)
                                                                aaCCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGCATTGCGTggg                                                                                     >  1:1277264/1‑150 (MQ=255)
                                                                aaCCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGCATTGCGTggg                                                                                     <  1:1580096/150‑1 (MQ=255)
                                                                 aCCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGCATTGCGTgggg                                                                                    <  1:296388/150‑1 (MQ=255)
                                                                 aCCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGCATTGCGTgggg                                                                                    <  2:1005214/150‑1 (MQ=255)
                                                                  cccTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGCATTGCGTGGGGt                                                                                   >  1:532804/1‑150 (MQ=255)
                                                                  cccTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGCATTGCGTGGGGt                                                                                   >  2:1486153/1‑150 (MQ=255)
                                                                   ccTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGCATTGCGTGGGGTc                                                                                  >  1:55301/1‑150 (MQ=255)
                                                                         ccAGGAATAATCCCTGGCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCTGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGCATTGCGTGGGGTCTGGCGCt                                                                           >  2:612893/1‑150 (MQ=255)
                                                                          cAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGCATTGCGTGGGGTCTGGCGCt                                                                           >  2:1448820/1‑150 (MQ=255)
                                                                           aGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGCATTGCGTGGGGTCTGGCGCtc                                                                          <  2:1358512/150‑1 (MQ=255)
                                                                            ggAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGCATTGCGTGGGGTCTGGCGCtca                                                                         <  2:674276/150‑1 (MQ=255)
                                                                              aataatCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGCATTGCGTGGGGTCTGGCGCtcatc                                                                       <  1:1271575/150‑1 (MQ=255)
                                                                                  atCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGTAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGCATTGCGTGGGGTCTGGCGCTCATcaaca                                                                   >  1:1418114/1‑150 (MQ=255)
                                                                                  atCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGCATTGCGTGGGGTCTGGCGCTCATcaaca                                                                   >  2:685686/1‑150 (MQ=255)
                                                                                        ggCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGGTCGGGGAAGGCAAAACTGCATCAGAGAAAGCATTGCGTGGGGTCTGGCGCTCATCAACAATCggg                                                             <  1:1300369/150‑1 (MQ=255)
                                                                                        ggCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGCATTGCGTGGGGTCTGgc                                                                              <  2:508703/133‑1 (MQ=255)
                                                                                        ggCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGCATTGCGTGGGGTCTGgc                                                                              >  1:508703/1‑133 (MQ=255)
                                                                                          cccAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCTTGACGGCTCGGGGAAGGCAAAAATGCATAAGAGAAAGCATTGCGTGGGGTCTGGCGCTCATCAACAATcctgca                                                           >  1:719798/1‑145 (MQ=255)
                                                                                          cccAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGCATTGCGTGGGGTCTGGCGCTCATCAACAATCGGGGa                                                           >  1:210405/1‑150 (MQ=255)
                                                                                           ccAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGCATTGCGTGGGGTCTGGCGCTCATCAACAATCGGGGAg                                                          >  1:32830/1‑150 (MQ=255)
                                                                                           ccAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGCATTGCGTGGGGTCTGGCGCTCATCAACAATCGGGGAg                                                          >  1:736046/1‑150 (MQ=255)
                                                                                            cAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGCATTGCGTGGGGTCTGGCGCTCATCAACAATCGGGGAGt                                                         >  1:1306258/1‑150 (MQ=255)
                                                                                            cAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGCTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGCATTGAGTGGGGACTGGCGCTCATCAACAATCGGGGAGt                                                         >  2:973725/1‑150 (MQ=255)
                                                                                             aaaaaGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGCATTGCGTGGGGTCTGGCGCTCATCAACAATCGGGGAGTc                                                        >  1:1367608/1‑150 (MQ=255)
                                                                                             aaaaaGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGCATTGCGTGGGGTCTGGCGCTCATCAACAATCGGGGAGTc                                                        <  1:685686/150‑1 (MQ=255)
                                                                                                   ccGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGCATTGCGTGGGGTCTGGCGCTCATCAACAATCGGGGAGTCAGCAgg                                                  >  2:274017/1‑150 (MQ=255)
                                                                                                          ttttttATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGCATTGCGTGGGGTCTGGCGCTCATCAACAATCGGGGAGTCAGCAGGGGGCTGa                                           <  2:409879/150‑1 (MQ=255)
                                                                                                               tataTCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGCATTGCGTGGGGTCTGGCGCTCATCAACAATCGGGGAGTCAGCAGGGGGCTGAAACgg                                      <  1:1570253/150‑1 (MQ=255)
                                                                                                                atatCTGCGTTCCGCTAAAAAGTGCAAATGATCTGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGCATTGCGTGGGGTCTGGCGCTCATCAACAATCGGGGAGTCAGCAGGGGGCTGAAACggg                                     <  1:612893/150‑1 (MQ=255)
                                                                                                                atatCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGCATTGCGTGGGGTCTGGCGCTCATCAACAATCGGGGAGTCAGCAGGGGGCTGAAACggg                                     <  2:1250236/150‑1 (MQ=255)
                                                                                                                  atCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCTTTTTGCGTGACCGCTCGGGGa                                                                                                                       >  1:293589/1‑66 (MQ=255)
                                                                                                                  atCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCTTTTTGCGTGACCGCTCGGGGa                                                                                                                       <  2:299355/66‑1 (MQ=255)
                                                                                                                  atCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCTTTTTGCGTGACCGCTCGGGGa                                                                                                                       <  2:293589/66‑1 (MQ=255)
                                                                                                                  atCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGCATTGCGTGGGGTCTGGCGCTCATCAACAATCGGGGAGTCAGCAGGGGGCTGAAACGGGaa                                   >  2:1145356/1‑150 (MQ=255)
                                                                                                                  atCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACCGCTCGGGGa                                                                                                                       >  1:299355/1‑66 (MQ=255)
                                                                                                                    cTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGCATTGCGTGGGGTCTGGCGCTCATCAACAATCGGGGAGTCAGCAGGGGGCTGAAACGGGAAAg                                 <  2:619295/150‑1 (MQ=255)
                                                                                                                     tGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGCATTGCGTGGGGTCTGGCGCTCATCAACAATCGGGGAGTCAGCAGGGGGCTGAAACGGGAAAGc                                >  1:1115965/1‑150 (MQ=255)
                                                                                                                     tGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGCATTGCGTGGGGTCTGGCGCTCATCAACAATCGGGGAGTCAGCAGGGGGCTGAAACGGGAAAGc                                >  1:1326377/1‑150 (MQ=255)
                                                                                                                      gCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGCATTGCGTGGGGTCTGGCGCTCATCAACAATCGGGGAGTCAGCAGGGGGCTGAAACGGGAAAGcc                               <  2:1326377/150‑1 (MQ=255)
                                                                                                                        gTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGCATTGCGTGGGGTCTGGCGCTCATCAACAATCGGGGAGTCAGCAGGGGGCTGAAACGGGAAAGcccc                             >  1:24540/1‑150 (MQ=255)
                                                                                                                              cTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGCATTGCGTGGGGTCTGGCGCTCATCAACAATCGGGGAGTCAGCAGGGGGCTGAAACGGGAAAGCCCCTCCCga                       <  1:237828/150‑1 (MQ=255)
                                                                                                                                aaaaaGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGCATTGCGTGGGGTCTGGCGCTCATCAACAATCGGGGTGTCAGCAGGGGGCTGAAACGGGAAAGCCCCTCCCgaag                     >  2:103772/1‑150 (MQ=255)
                                                                                                                                aaaaaGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGCATTGCGTGGGGTCTGGCGCTCATCAACAATCGGGGAGTCAGCAGGGGGCTGAAACGGGAAAGCCCCTCCCgaag                     <  2:55301/150‑1 (MQ=255)
                                                                                                                                   aaGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGCATTGCGTGGGGTCTGGCGCTCATCAACAATCGGGGAGTCAGCAGGGGGCTGAAACGGGAAAGCCCCTCCCGaagaag                  <  1:270538/150‑1 (MQ=255)
                                                                                                                                      tGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGCATTGCGTGGGGTCTGGCGCTCATCAACAATCGGGGAGTCAGCAGGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAgggg               <  2:820849/150‑1 (MQ=255)
                                                                                                                                       gCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGGTCGGGGGAGGCAAAACTGCATCAGAGAAAGCATTGCGTGGGGGCTGGCGCTCATCAACAATCGGGGAGTCAGCAGGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAggggg              <  2:1347106/150‑1 (MQ=255)
                                                                                                                                       gCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGCATTGCGTGGGGTCTGGCGCTCATCAACAATCGGGGAGTCAGCAGGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAggggg              <  2:1580354/150‑1 (MQ=255)
                                                                                                                                        cAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGCATTGCGTGGGGTCTGGCGCTCATCAACAATCGGGGAGTCAGCAGGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGGc             >  2:1359893/1‑150 (MQ=255)
                                                                                                                                        cAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGCATTGCGTGGGGTCTGGCGCTCATCAACAATCGGGGAGTCAGCAGGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGGc             >  2:684157/1‑150 (MQ=255)
                                                                                                                                        cAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGCATTGCGTGGGGTCTGGCGCTCATCAACAATCGGGGAGTCAGCAGGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGGc             <  2:1418114/150‑1 (MQ=255)
                                                                                                                                         aaaTGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGCATTGCGTGGGGTCTGGCGCTCATCAACAATCGGGGAGTCAGCAGGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGGcc            >  1:552373/1‑150 (MQ=255)
                                                                                                                                         aaaTGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGCATTGCGTGGGGTCTGGCGCTCATCAACAATCGGGGAGTCAGCAGGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGGcc            >  1:922638/1‑150 (MQ=255)
                                                                                                                                            tGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGCATTGCGTGGGGTCTGGCGCTCATCAACAATCGGGGAGTCAGCAGGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGGCCTTa         <  2:1182086/150‑1 (MQ=255)
                                                                                                                                               tCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGCATTGCGTGGGGTCTGGCGCTCATCAACAATCGGGGAGTCAGCAGGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGGCCTTACAg      <  2:1115965/150‑1 (MQ=255)
                                                                                                                                                  ggTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGCATTGCGTGGGGTCTGGCGCTCATCAACAATCGGGGAGTCAGCAGGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGGCCTTACAGAAg   <  2:1655990/150‑1 (MQ=255)
                                                                                                                                                  ggTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGCATTGCGTGGGGTCTGGCGCTCATCAACAATCGGGGAGTCAGCAGGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGGCCTTACAGAAg   >  2:207084/1‑150 (MQ=255)
                                                                                                                                                   gttgttGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGCATTGCGTGGGGTCTGGCGCTCATCAACAATCGGGGAGTCAGCAGGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGGCCTTACAGAAgg  >  1:633893/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GCTCACTGGCGGCATTGACTTCACCCTTAACCTCAACAACTCCCGTCGAAGAAAAAGCAGAGTAAACCCTCTTCCAGGAATAATCCCTGGCCCAAAAAGCCGGGGATTTTTTATATCTGCGTTCCGCTAAAAAGTGCAAATGATCAGGTTGTTGCAGCATTTTGCGTGACGGCTCGGGGAAGGCAAAACTGCATCAGAGAAAGCATTGCGTGGGGTCTGGCGCTCATCAACAATCGGGGAGTCAGCAGGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGGCCTTACAGAAGG  >  NC_000913/3194461‑3194757

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: