Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,121,962 C→A 8.5% G109V (GGT→GTT)  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,9620CA 8.5% 307.0 / 13.7 106G109V (GGT→GTT) hslVpeptidase component of the HslUV protease
Reads supporting (aligned to +/- strand):  ref base C (47/50);  new base A (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 = 8.68e-01

Only 100 of 109 total aligned reads displayed.

GTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCTTTGACCAGATGGCCCTGATGCATTTCCAGTTTACGTTCAAACAGT  >  NC_000913/4121814‑4122111
                                                                                                                                                    |                                                                                                                                                     
gTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTAcc                                                                                                                                                      <  2:186998/150‑1 (MQ=255)
gTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGGGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTAcc                                                                                                                                                      <  2:868912/150‑1 (MQ=255)
   gCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGt                                                                                                                                                   >  1:807429/1‑150 (MQ=255)
   gCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGt                                                                                                                                                   <  2:657545/150‑1 (MQ=255)
   gCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGGTAACGGt                                                                                                                                                   <  1:1126773/150‑1 (MQ=255)
    ccggCAATATACAACGCCATTTCAGCAATATCACGGGCGCTAAGGTCAGTGTTTTCTCACCGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGGCACCGTTACCGGtg                                                                                                                                                  <  2:1033572/147‑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)
         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)
                       tttCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGttt                                                                                                                               <  2:1390534/150‑1 (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)
                              aaTTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGAGATGCAGttt                                                                                                                               >  1:207024/1‑143 (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                                                                                                               <  2:1364326/150‑1 (MQ=255)
                                       ggCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCcag                                                                                                               >  1:608769/1‑150 (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)
                                                  cAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCcg                                                                                                                         <  1:580057/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)
                                                                  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)
                                                                          ccgcAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGTTGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTAc                                                                            >  2:1605526/1‑150 (MQ=255)
                                                                          ccgcAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTAc                                                                            >  2:255484/1‑150 (MQ=255)
                                                                          ccgcAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTAc                                                                            >  1:1450579/1‑150 (MQ=255)
                                                                          ccgcAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTAc                                                                            >  1:1319404/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                                                                                            <  2:291518/127‑1 (MQ=255)
                                                                                 cTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTgcgc                                                                                            >  1:291518/1‑127 (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)
                                                                                           gggCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGc                                                           <  1:997619/150‑1 (MQ=255)
                                                                                             gccgccGGAGCTGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTc                                                         >  2:243852/1‑150 (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)
                                                                                                        cGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCtt                                              >  1:1349190/1‑150 (MQ=255)
                                                                                                          aTAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCTTTg                                            <  1:1586705/150‑1 (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)
                                                                                                                                        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)
                                                                                                                                                   aCCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCTTTGACCAGATGGCCCTGATGCATTTCCAGTTTACGTTCAAACAg   <  2:126300/150‑1 (MQ=255)
                                                                                                                                                    ccGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCAAGTCTTTTGCCAGCTCAATGGCGGCTTTGCCCAGATGGCCATGATGCATTTCCAGTTTACGTTCAAACAGt  >  2:883205/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GTCGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCTTTGACCAGATGGCCCTGATGCATTTCCAGTTTACGTTCAAACAGT  >  NC_000913/4121814‑4122111

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: