Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,121,959 T→G 8.5% N110T (AAC→ACC)  hslV ← peptidase component of the HslUV protease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,121,9590TG 8.5% 277.5 / 16.8 106N110T (AAC→ACC) hslVpeptidase component of the HslUV protease
Reads supporting (aligned to +/- strand):  ref base T (47/50);  new base G (3/6);  total (50/56)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.58e-01

Only 100 of 108 total aligned reads displayed.

AAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCTTTGACCAGATGGCCCTGATGCATTTCCAGTTTACGTTCAAAC  >  NC_000913/4121810‑4122108
                                                                                                                                                     |                                                                                                                                                     
aaaTGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGGGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGt                                                                                                                                                       <  1:510707/150‑1 (MQ=255)
    gTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTAcc                                                                                                                                                   <  2:186998/150‑1 (MQ=255)
    gTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGGGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTAcc                                                                                                                                                   <  2:868912/150‑1 (MQ=255)
       gCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGt                                                                                                                                                >  1:807429/1‑150 (MQ=255)
       gCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGGTAACGGt                                                                                                                                                <  1:1126773/150‑1 (MQ=255)
        ccggCAATATACAACGCCATTTCAGCAATATCACGGGCGCTAAGGTCAGTGTTTTCTCACCGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGGCACCGTTACCGGtg                                                                                                                                               <  2:1033572/147‑1 (MQ=255)
        ccTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGtg                                                                                                                                               >  1:1688318/1‑150 (MQ=255)
          tGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGtgat                                                                                                                                             <  2:684558/150‑1 (MQ=255)
          tGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGtgat                                                                                                                                             >  2:1214356/1‑150 (MQ=255)
           gCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGGTAACGGtgacg                                                                                                                                            <  2:463839/150‑3 (MQ=255)
            cAATATCCACCGCCTTGTCAGCAATTGCACGGCCGCGAAGTTAAGTGTTGTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTCAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGGTAACGGtgacgt                                                                                                                                           <  2:836801/150‑4 (MQ=255)
            cAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACCGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGACCCGATAGCAATAAGATCGTTATCTGGCTGCCCCACGTCACCGTTCCCGAtgatga                                                                                                                                           >  2:1155349/1‑150 (MQ=255)
             aaTATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTgatgat                                                                                                                                          <  2:375667/150‑1 (MQ=255)
             aaTATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGGTAACGGtgacgtg                                                                                                                                          <  2:1280614/150‑5 (MQ=255)
              atatCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGGTAACGGtgacgtgg                                                                                                                                         <  2:757730/150‑6 (MQ=255)
              aaatCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACATCACCGGTAACGGtgacgtgg                                                                                                                                         <  1:616014/148‑6 (MQ=255)
                 tCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGc                                                                                                                                      <  2:1141510/150‑1 (MQ=255)
                 tCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCACCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGc                                                                                                                                      >  1:156218/1‑150 (MQ=255)
                  ccAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCACAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGCTCGTTTTCTGGCGGCACCACCTCACCGTAACCGGGGCTGATAcgcg                                                                                                                                     >  1:257033/1‑150 (MQ=255)
                       gCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCa                                                                                                                                >  1:616059/1‑150 (MQ=255)
                        ccTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGTGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAg                                                                                                                               >  1:1248024/1‑150 (MQ=255)
                        ccTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAg                                                                                                                               >  2:1586705/1‑150 (MQ=255)
                               aGCAATTTCACGGGCGCTAAGTTCACTGTTGTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGGCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATc                                                                                                                        <  1:1047190/150‑1 (MQ=255)
                                 cAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCcg                                                                                                                      >  2:1331979/1‑150 (MQ=255)
                                 cAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTATGGCTGCACAACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCAtacc                                                                                                                      >  2:781747/1‑147 (MQ=255)
                                  aaTTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGTGATGCAGttt                                                                                                                            <  2:207024/143‑1 (MQ=255)
                                   cgttCCCGGGCGCTACGTGCAGTGGTTTCTAACCGCGCGCGCGCCGCAGCCGGGGCGTAAGGGGCGCCGGAGCCGCTCGCCAGAAGATCGTGGTCTGGCGGCACCCCGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCcgcg                                                                                                                    <  1:112011/148‑1 (MQ=255)
                                   aTTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTAACCGGTGATGATAAGCGATGCAGTATCATCcgcg                                                                                                                    >  1:348113/1‑150 (MQ=255)
                                    tttCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGa                                                                                                                   >  1:1418344/1‑150 (MQ=255)
                                      tcccGGGCGCTAAGTTGAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACt                                                                                                                 <  1:768213/147‑1 (MQ=255)
                                         cGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGcc                                                                                                              >  2:1020520/1‑150 (MQ=255)
                                           ggCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCcag                                                                                                            >  1:608769/1‑150 (MQ=255)
                                           ggCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCcag                                                                                                            <  2:1364326/150‑1 (MQ=255)
                                               cTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGt                                                                                                        >  2:426153/1‑150 (MQ=255)
                                                 aaGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGc                                                                                                      <  2:1688318/150‑1 (MQ=255)
                                                     tCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCa                                                                                                  <  2:645723/150‑1 (MQ=255)
                                                     tCAGAGTTTTCTAACAGCGCGCGCGCCGCAGCCAGGGCGAAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGACACCGTTACCGGTGATGATAAGCGAAGCAGATTCATCAGCGACTGCCAGCAGAGCTTCa                                                                                                  >  1:335347/1‑150 (MQ=255)
                                                      cAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCcg                                                                                                                      >  2:580057/1‑129 (MQ=255)
                                                      cAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCcg                                                                                                                      >  2:586677/1‑129 (MQ=255)
                                                      cAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCcg                                                                                                                      <  1:586677/129‑1 (MQ=255)
                                                      cAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCaa                                                                                                 >  2:1373863/1‑150 (MQ=255)
                                                         ggtTTTCTAACAGAGCGCGGGACGGAGGCTGGGCGTAAGGGCCGCAGGAGCCTATTGCAATAAGATCGTTTTCTGGCTGCACCACTTCACCGTTACCGGTTATGATAAGCGATGCAGTTTCATCCGCGACTGGCAGCAGTGATTCAAGtt                                                                                              <  1:962722/149‑1 (MQ=255)
                                                          gtTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGttt                                                                                             <  1:1373863/150‑1 (MQ=255)
                                                           ttttCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTg                                                                                            <  1:260694/150‑1 (MQ=255)
                                                            tttCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTgc                                                                                           <  2:75633/150‑1 (MQ=255)
                                                             ttCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTgcg                                                                                          >  2:525949/1‑150 (MQ=255)
                                                                 aaCAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCgcagc                                                                                      <  1:415033/150‑1 (MQ=255)
                                                                      cgcgcgcgcCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGGTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGGGCTTCAAGTTTGCGCAGCATGCg                                                                                 >  1:204122/1‑150 (MQ=255)
                                                                        cgcgcgcCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGAt                                                                               <  1:1214356/150‑1 (MQ=255)
                                                                        cgcgcgcCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGAt                                                                               <  1:1687977/150‑1 (MQ=255)
                                                                         gcgcggCGCAGCCTGTGCGTAAGGGCCGCCTTAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCTGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATc                                                                              <  2:489552/150‑1 (MQ=255)
                                                                         gcgcgcCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATc                                                                              >  1:939092/1‑150 (MQ=255)
                                                                             gccgcAGCCGGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTa                                                                          <  1:1337979/150‑1 (MQ=255)
                                                                              ccgcAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGTTGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTAc                                                                         >  2:1605526/1‑150 (MQ=255)
                                                                              ccgcAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTAc                                                                         >  2:255484/1‑150 (MQ=255)
                                                                              ccgcAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTAc                                                                         >  1:1319404/1‑150 (MQ=255)
                                                                              ccgcAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTAc                                                                         >  1:1450579/1‑150 (MQ=255)
                                                                                 cAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGcc                                                                      <  1:590694/150‑1 (MQ=255)
                                                                                  aGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCa                                                                     >  1:126300/1‑150 (MQ=255)
                                                                                    ccTGGGCGTCCGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCCCGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGt                                                                   <  2:257033/150‑1 (MQ=255)
                                                                                     cTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTgcgc                                                                                         >  1:291518/1‑127 (MQ=255)
                                                                                     cTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTgcgc                                                                                         <  2:291518/127‑1 (MQ=255)
                                                                                        ggCGTACGGGCCGCCGGAGCCGATAGCAATAAGCTCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATCAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCttt                                                               <  1:237832/150‑1 (MQ=255)
                                                                                        ggCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCttt                                                               <  1:232802/150‑1 (MQ=255)
                                                                                           gTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGc                                                            >  2:1623568/1‑150 (MQ=255)
                                                                                              aGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAg                                                         <  2:772757/150‑1 (MQ=255)
                                                                                              aGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAg                                                         >  1:1200653/1‑150 (MQ=255)
                                                                                               gggCCGCCTGAGCCTATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGAGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGc                                                        <  2:1185353/150‑1 (MQ=255)
                                                                                               gggCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGc                                                        <  2:1418344/150‑1 (MQ=255)
                                                                                                     ccgGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAAcg                                                  >  1:228438/1‑150 (MQ=255)
                                                                                                       ggAGCCGATAGCAATAAGTTCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCACcggc                                                >  2:839274/1‑150 (MQ=255)
                                                                                                       ggAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATg                                                                                   <  2:1408117/115‑1 (MQ=255)
                                                                                                       ggAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATg                                                                                   >  1:1408117/1‑115 (MQ=255)
                                                                                                       ggAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAAcggc                                                <  1:1623568/150‑1 (MQ=255)
                                                                                                           ccGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACAGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCt                                            >  1:137542/1‑150 (MQ=255)
                                                                                                            cGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCtt                                           >  1:909937/1‑150 (MQ=255)
                                                                                                              aTAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCTTTg                                         <  1:1586705/150‑1 (MQ=255)
                                                                                                                      aaGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCTTTGACCAGATg                                 >  2:122774/1‑150 (MQ=255)
                                                                                                                       aGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCTTTGACCAGATgg                                <  2:608769/150‑1 (MQ=255)
                                                                                                                        gATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCTTTGACCAGATGGc                               <  1:426153/150‑1 (MQ=255)
                                                                                                                            cttttCTGGCTGCACCACGTCACCGTTACCGGGGATGATACGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTGTTGCCAGCTCAACGGCGGCTTTGACCAGATGGCCCTg                           <  1:1202998/149‑1 (MQ=255)
                                                                                                                              tttCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCATGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCTTTGACCAGATGGCCCTGAt                         >  2:675302/1‑150 (MQ=255)
                                                                                                                               ttCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCTTTGACCAGATGGCCCTGATg                        >  2:964264/1‑150 (MQ=255)
                                                                                                                               ttCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCTTTGACCAGATGGCCCTGATg                        <  1:525949/150‑1 (MQ=255)
                                                                                                                                tCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTATCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCTTTGACCAGATGGCCCTGATGc                       <  2:230639/150‑1 (MQ=255)
                                                                                                                                    gCTGCACCACGGCACCGTTACCGGTGATGATAAGAGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCTTTGACCAGATGGCCCTGATGCAttt                   <  2:951100/150‑1 (MQ=255)
                                                                                                                                      tGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCACCGGCGGCTTTGACCAGATGGCCCTGATGCATTTcc                 >  2:1463394/1‑150 (MQ=255)
                                                                                                                                      tGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCTTTGACCAGATGGCCCTGATGCATTTcc                 <  2:204122/150‑1 (MQ=255)
                                                                                                                                       gccccacGTCACCGTTACCGGTGATGATAAGCCCGGCAGATTCCTCCGCGACTGCCAGCCGTGAATCAAGAATGCGCACCAGTCGATCGGTACGCCAGTCTATTGCCAGCTCAACGGCGGCTTTGACCAGCTGGCCCTGATGCATTACCa                >  2:1089695/4‑150 (MQ=255)
                                                                                                                                        cttatcatCACCGGTAACGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCTTTGACCAGATGGCCCTGATGCATTTCCAg               >  1:903425/8‑150 (MQ=255)
                                                                                                                                            tcatCACCGGTAACGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCTTTGACCAGATGGCCCTGATGCATTTCCAGTTTa           >  1:1256666/4‑150 (MQ=255)
                                                                                                                                             cGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCTTTGACCAGATGGCCCTGATGCATTTCCAGTTTAc          <  1:122774/150‑1 (MQ=255)
                                                                                                                                              atCACCGGTAACGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCATTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGAGGCTATGACCAGATGGCCCTGATGCTTTTCCAGTTTACg         >  2:1501217/2‑150 (MQ=255)
                                                                                                                                               tCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCTTTGACCAGATGGCCCTGATGGATTTCCAGTTTACGt        <  1:255484/150‑1 (MQ=255)
                                                                                                                                               gctccGTTCCCGGTGATGATAAGCGATGCAGTTTCATCGGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGAATGCGATCGGTACGCCAGTCTTTTGCCTGCTCATCGGCGGCTTTGACCAGATGGCCCTGATGCATTTCCAGTTTACGt        <  1:641494/147‑1 (MQ=255)
                                                                                                                                                cACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCTTTGACCAGATGGCCCTGATGCATTTCCAGTTTACGtt       >  2:499488/1‑150 (MQ=255)
                                                                                                                                                cACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCTTTGACCAGATGGCCCAGATGCATTTCCAGTTTACGtt       >  1:163203/1‑150 (MQ=255)
                                                                                                                                                  ccGTTTCCGGTGGTGATAAGCGCTGCAGTTTCATCCGCGACTGCCAGCAGTGCGTCAAGTTTGCGCAGGATGCGATCGGTTCGGCCGTGTTTTGCCAGCTCACCGGCGGCTTTGCCCAGATGGCCCTGCTGCATTTCCAGTTTACGGTCa     >  2:1361277/1‑150 (MQ=255)
                                                                                                                                                    gTTACCGCTTATGATAAGCGATTCAGTATCATCCTCTACTTCCAGCAGGGCTTCAAGTTTGAGCTGCATGCGTTCGGTACGCCAGTCTTTTGTCAGCTCAACGGCGTCTTTGACCAGATGGCCCTGATTCATTTCCAGTGTACGGTCaaa   <  2:846088/150‑1 (MQ=255)
                                                                                                                                                     ttACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCTTTGACCAGATGGCCCTGATGCATTTCCAGTTTACGTTCAAAc  <  1:1020520/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AAATGTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCTTTGACCAGATGGCCCTGATGCATTTCCAGTTTACGTTCAAAC  >  NC_000913/4121810‑4122108

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: