Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,122,332 A→T 5.5% intergenic (‑45/+48) hslV ← / ← ftsN peptidase component of the HslUV protease/essential cell division protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,122,3320AT 5.5% 340.3 / 5.8 104intergenic (‑45/+48)hslV/ftsNpeptidase component of the HslUV protease/essential cell division protein
Reads supporting (aligned to +/- strand):  ref base A (55/43);  new base T (2/4);  total (57/47)
Fisher's exact test for biased strand distribution p-value = 4.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 112 total aligned reads displayed.

CTTTTTCACGTTGCCTTTCATTACGGTATTGCCCAACGTGGCCTGACCATCACCAGCGATGACCACATGGCCGTTACGGCGTACGCTTACTATAGTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGCTGTCTGCGTTCTCTTTGCCTTTCACCGGGCCAATGA  >  NC_000913/4122186‑4122480
                                                                                                                                                  |                                                                                                                                                    
cTTTTTCCCGCTGCGTTTCATTCCGGTTTTGGCCTACGTGGCCTGTCCATCACCAGCGATGACCACTTGGCCGTTTCGGCGTACGCTTACTATAGTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACaaaa                                                                                                                                                    <  4:340242/149‑1 (MQ=255)
cTTTTTCACGTTGCCTTTCATTACGGTATTGCCCAACGTGGCCTGACCATCACCAGCGATGACCACATGGCCGTTACGGCGTACGCTTACTATAGTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACaaaa                                                                                                                                                    <  5:90410/149‑1 (MQ=255)
cTTTTTCACGTTGCCTTTCATTACGGTATTGCCCAACGTGGCCTGACCATCACCAGCGATGACCACATGGCCGTTACGGCGTACGCTTACTATAGTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACaaaa                                                                                                                                                    <  8:417730/149‑1 (MQ=255)
cTTTTTCACGTTGCCTTTCATTACGGTATTGCCCAACGTGGCCCGACCATCACCAGCGATGACCACATGGCCGTTACGGCGTACGCTTACTATAGTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACaaaa                                                                                                                                                    <  7:395507/149‑1 (MQ=255)
  ttttCACGTTGCCTTTCATTACGGTATTGCCCAACGTGGCCTGACCATCACCAGCGATGACCACATGGCCGTTACGGCGTACGCTTACTATAGTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAg                                                                                                                                                   >  5:404590/1‑148 (MQ=255)
        cggtgaaaggcaaagagaacgcagacagcaccctcaatcggttgaagatggcgggtcatacaaactgcattcggctcgccgccgggggttgaaaccctcaaaatcccccccatctataattgcattatgCCCCGTACAAAAGTACgggg                                                                                                                                            >  3:15652‑M1/130‑149 (MQ=255)
        cGTTGCCTTTCATTACGGTATTGCCCAACGTGGCCTGACCATCACCAGCGATGACCACATGGCCGTTACGGCGTACGCTTACTATAGTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACgggg                                                                                                                                            <  8:6309/149‑1 (MQ=255)
        cGTTGCCTTTCATTACGGTATTGCCCAACGTGGCCTGACCATCACCAGCGATGACCACATGGCCGTTACGGCGTACGCTTACTATAGTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACgggg                                                                                                                                            <  6:404590/149‑1 (MQ=255)
         gTTGCCTTTCATTACGGTATTGCCCAACGTGGCCTGACCATCACCAGCGATGACCACATGGCCGTTACGGCGTACGCTTACTATAGTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACgggg                                                                                                                                            >  2:11354/1‑148 (MQ=255)
         gTTGCCTTTCATTACGGTATTGCCCAACGTGGCCTGACCATCACCAGCGATGACCACATGGCCGTTACGGCGTACGCTTACTATAGTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACgggg                                                                                                                                            >  8:438595/1‑148 (MQ=255)
         gTTGCCTTTCATTACGGTATTGCCCAACGTGGCCTGACCATCACCAGCGATGACCACATGGCCGTTACGGCGTACGCTTACTATAGTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACgggg                                                                                                                                            >  8:335406/1‑148 (MQ=255)
         gTTGCCTTTCATTACGGTATTGCCCAACGTGGCCTGACCATCACCAGCGATGACCACATGGCCGTTACGGCGTACGCTTACTATAGTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGc                                                                                                                                           <  1:48368/149‑1 (MQ=255)
             ccTTTCATTACGGTATTGCCCAACGTGGCCTGACCATCACCAGCGATGACCACATGGCCGTTACGGCGTACGCTTACTATAGTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAGCCCCGTACTTTTGTACggggtttgt                                                                                                                                       <  2:63951/149‑6 (MQ=255)
              cTTTCATTACGGTATTGCCCAACGTGGCCTGACCATCACCAGCGATGACCACATGGCCGTTACGGCGTACGCTTACTATAGTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAAt                                                                                                                                      >  3:390708/1‑149 (MQ=255)
                   aTTACGGTATTGCCCCACGGGGCCTGACCATCACCCGCGATGACCACATGGCCGTTACGGCGGACGCTTACTATAGTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCaa                                                                                                                                  <  4:228076/148‑1 (MQ=255)
                    ttACGGTATTGCCCAACGTGGCCTGACCATCACCAGCGATGACCACATGGCCGTTACGGCGTACGCTTACTATAGTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCaa                                                                                                                                  >  7:442714/1‑147 (MQ=255)
                    ttACGGTATTGCCCAACGTGGCCTGACCATCACCAGCGATGACCACATGGCCGTTACGGCGTACGCTTACTATAGTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCaa                                                                                                                                  <  8:442714/147‑1 (MQ=255)
                    ttACGGTATTGCCCAACGTGGCCTGACCATCACCAGCGATGACCACATGGCCGTTACGGCGTACGCTTACTATAGTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAAtt                                                                                                                                >  7:277907/1‑149 (MQ=255)
                    ttACGGTATTGCCCAACGTGGCCTGACCATCACCAGCGATGACCACATGGCCGTTACGGCGTACGCTTACTATAGTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAAt                                                                                                                                 >  6:156374/1‑148 (MQ=255)
                       cGGTATTGCCCAACGTGGCCTGACCATCACCAGCGATGACCACATGGCCGTTACGGCGTACGCTTACTATAGTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATtat                                                                                                                              <  1:11354/148‑1 (MQ=255)
                                        gCCTGCCCATAACCAGCGATGACAACATGGCCGTTACGGCGTACTCTTAATATAGTTGTCAAGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGa                                                                                                            <  5:134000/149‑1 (MQ=255)
                                        gCCTGACCATCACCAGCGATTACCACATGGCCGTTACGGCGTACGCTTACTATAGTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGa                                                                                                            <  5:82413/149‑1 (MQ=255)
                                        gCCTGACCATCACCAGCGATGACCACATGGCCGTTACGGCGTACGCTTACTATAGTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGa                                                                                                            <  8:277907/149‑1 (MQ=255)
                                         ccTGACCATCACCAGCGCTGACCACATGGCCGTTACGGCGTACGCTTACTATAGTTGTCACGAGCTGACACCTTGGTTACGAATACAGAGTACAACCCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAg                                                                                                           >  8:435575/1‑149 (MQ=255)
                                          cTGACCATCACCAGCGATGACCACATGGCCGTTACGGCGTACGCTTACTATAGTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACaaaa                                                                                                                                                    >  7:382549/1‑107 (MQ=255)
                                          cTGACCATCACCAGCGATGACCACATGGCCGTTACGGCGTACGCTTACTATAGTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACaaaa                                                                                                                                                    <  8:382549/107‑1 (MQ=255)
                                          cTGACCATCACCAGCGATGACCACATGGCCGTTACGGCGTACGCTTACTATAGTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAgg                                                                                                          <  7:125395/149‑1 (MQ=255)
                                          cTGACCATCACCAGCGATGACCACATGGCCGTTACGGCGTACGCTTACTATAGTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAgg                                                                                                          <  7:335406/149‑1 (MQ=255)
                                          cTGACCATCACCAGCGATGACCACATGGCCGTTACGGCGTACGCTTACTATAGTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAgg                                                                                                          >  1:300840/1‑149 (MQ=255)
                                               cATCACCAGCGATGACCACATGGCCGTTACGGCGTACGCTTACTATAGTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGttt                                                                                                      >  1:324156/1‑148 (MQ=255)
                                               cATCACCAGCGATGACCACATGGCCGTTACGGCGTACGCTTACTATAGTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTc                                                                                                     <  4:416529/149‑1 (MQ=255)
                                                aTCACCAGCGATGACCACATGGCCGTAACGGCGTACGCTTACTATAGTTGACACCAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATCCAATAATAGATGGGGGGTATTTTGAGGGATTc                                                                                                     >  8:232192/1‑148 (MQ=255)
                                                 ttgaagatggcgggtcatacaaactgcattcggctcgccgccgggggttgaaaccctcaaaatcccccccatctataattgcattatgCCCCGTACTTTTGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGtt                                                                                                       >  8:270273‑M1/89‑145 (MQ=255)
                                                 ttgaagatggcgggtcatacaaactgcattcggctcgccgccgggggttgaaaccctcaaaatcccccccatctataattgcattatgCCCCGTACTTTTGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGtt                                                                                                       <  7:270273‑M1/57‑1 (MQ=255)
                                                 tCACCAGCGATGACCACATGGCCGTTACGGCGTACGCTTACTATAGTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCaa                                                                                                   >  1:453897/1‑149 (MQ=255)
                                                    ccAGCGATGACCACATGGCCGTTACGGCGTACGCTTACTATAGTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAAccc                                                                                                >  3:245205/1‑149 (MQ=255)
                                                    ccAGCGATGACCACATGGCCGTTACGGCGTAAGCTTACTATAGTTGTCACGAGCTGACCCCTTTGTTACGAATACAGAGTACAAACACCGTACAAAAGTACCGGGCATAATGCAATTATAGATTGGGGGGATTTTGAGGGTTTCAAccc                                                                                                >  3:231318/1‑149 (MQ=255)
                                                     cAGCGATGACCACATGGCCGTTACGGCGTACGCTTACTATAGTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCCAcccc                                                                                               >  6:434721/1‑149 (MQ=255)
                                                         gATGACCACATGGCCGTTACGGCGTACGCTTACTATAGTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCcggc                                                                                           <  2:324156/149‑1 (MQ=255)
                                                              ccACATGGCCGTTACGGCGTACGCTTACTATAGTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAg                                                                                                           <  7:78278/128‑1 (MQ=255)
                                                               cgcaTGGCCGTTACGGCGTACGCTTACTATAGTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAg                                                                                     >  2:454867/3‑149 (MQ=255)
                                                               cacaTGGCCGTTACGGCGTACGCTTACTATAGTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATATATGGGGGGGATTTTGAGGGTTTCAACCCCCGgcggcggg                                                                                     >  2:92624/1‑147 (MQ=255)
                                                               cacaTGGCCGTTACGGCGTACGCTTACTATAGTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAg                                                                                     >  5:172158/1‑149 (MQ=255)
                                                               cacaTGGCCGTTACGGCGTACGCTTACTATAGTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAg                                                                                     >  3:305587/1‑149 (MQ=255)
                                                                        gTTACGGCGTACGCTTACTATAGTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATg                                                                              >  1:454484/1‑147 (MQ=255)
                                                                        gTTACGGCGTACGCTTACTATAGTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCa                                                                            <  4:305587/149‑1 (MQ=255)
                                                                            cGGCGTACGCTTACTATAGTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGtt                                                                         <  4:390708/148‑1 (MQ=255)
                                                                             ggCGTACGCTTACTATAGTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTg                                                                       >  3:392030/1‑149 (MQ=255)
                                                                              gCGTACGCTTACTATAGTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGt                                                                      <  4:392030/149‑1 (MQ=255)
                                                                                gTACGCTTACTATAGTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTAt                                                                    <  1:454867/149‑1 (MQ=255)
                                                                                         cTATAGTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCa                                                           >  2:294757/1‑149 (MQ=255)
                                                                                          tataGTTGTCACGAGCTGACCCCTTGTTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGa                                                                                                            <  1:458834/99‑1 (MQ=255)
                                                                                          tataGTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGa                                                                                                            >  2:458834/1‑99 (MQ=255)
                                                                                           ataGTTGTCACGAGCTGCCCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATc                                                         >  4:12126/1‑149 (MQ=255)
                                                                                            taGTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTa                                                                     >  6:354605/1‑136 (MQ=255)
                                                                                            taGTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTa                                                                     <  5:354605/136‑1 (MQ=255)
                                                                                               ttGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATgggg                                                                                                                      <  1:463085/84‑1 (MQ=255)
                                                                                               ttGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATgggg                                                                                                                      >  2:463085/1‑84 (MQ=255)
                                                                                               ttGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTc                                                      >  1:396936/1‑148 (MQ=255)
                                                                                                      gAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCt                                                        <  1:3436/139‑1 (MQ=255)
                                                                                                      gAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCt                                                        >  2:3436/1‑139 (MQ=255)
                                                                                                      gAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGAtt                                              <  5:156374/149‑1 (MQ=255)
                                                                                                      gAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGAtt                                              <  2:396936/149‑1 (MQ=255)
                                                                                                        gCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGa                                            <  2:300840/149‑1 (MQ=255)
                                                                                                         cTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAg                                           <  4:245205/149‑1 (MQ=255)
                                                                                                         cTGACCACTTGGTTACGAATACAGATTACAAACCCCGTACAAAATTACGTGGCATAATGCAATTATAGATGGGGGTGAGTTTGATGGTTTCAACCCCCGGCGGCGAGACGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAg                                           <  6:373517/149‑1 (MQ=255)
                                                                                                          tGACCCCTTGGTTGCGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAgg                                          >  4:239347/1‑149 (MQ=255)
                                                                                                           gACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCCCCCGAGTGAggg                                         >  4:73038/1‑149 (MQ=255)
                                                                                                           gACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAggg                                         >  5:111763/1‑149 (MQ=255)
                                                                                                            tcggcTTGGTTACTAATACAGAGTACAAACCCCGTAAAAAAGTACGGGGCATAATGAAATTATAGATGGGTGGGATTTTGATGGTTTCAACCCCcggc                                                                                           <  6:402882/94‑1 (MQ=255)
                                                                                                            aCCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAgg                                          <  2:454484/147‑1 (MQ=255)
                                                                                                              cccTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGc                                      >  2:226031/1‑149 (MQ=255)
                                                                                                               ccTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGCt                                     >  5:74898/1‑149 (MQ=255)
                                                                                                                cTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGCTg                                    >  6:198602/1‑149 (MQ=255)
                                                                                                                     ttACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGCTGTCTGc                               >  7:220173/1‑149 (MQ=255)
                                                                                                                     ttACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGCTGTCTGc                               >  3:189605/1‑149 (MQ=255)
                                                                                                                       aCGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGCTGTCTGCg                              <  1:402367/148‑1 (MQ=255)
                                                                                                                        cGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGCTGTCTGCGtt                            <  4:213860/149‑1 (MQ=255)
                                                                                                                        cGAATACAGAGTACAAAACCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGCTGTCTGCGtt                            <  6:36546/149‑1 (MQ=255)
                                                                                                                         gAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGCTGTCTGCGTtc                           >  5:317354/1‑149 (MQ=255)
                                                                                                                         gAATACAGAGTACAAACCCAGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGCTGTCTGCGTtc                           <  6:172158/149‑1 (MQ=255)
                                                                                                                            taattgcattatgCCCCGTACTTTTGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGCTGTCTGCGTTctct                        >  3:400088‑M1/14‑149 (MQ=255)
                                                                                                                               agagTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGCTGTCTGCGTTctct                        >  1:295525/1‑146 (MQ=255)
                                                                                                                               agagTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGCTGTCTGCGTTctct                        <  2:295525/146‑1 (MQ=255)
                                                                                                                               agagTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGCTGTCTGCGTTCTCTTTg                     <  8:220173/149‑1 (MQ=255)
                                                                                                                                   attatgCCCCGTACTTTTGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGCTGTCTGCGTTCTCTTTGCCtt                 <  7:246715‑M1/143‑1 (MQ=255)
                                                                                                                                     cAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGCTGTCTGCGTtct                          <  7:421582/138‑1 (MQ=255)
                                                                                                                                     cAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGCTGTCTGCGTtct                          >  8:421582/1‑138 (MQ=255)
                                                                                                                                     cAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGCTGTCTGCGTTCTCTTTGCCTTTc               <  7:435575/149‑1 (MQ=255)
                                                                                                                                      aaaCCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGCGCCGAATgaa                                                                            >  2:363372/1‑85 (MQ=255)
                                                                                                                                      aaaCCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGCTGTCTGCGTTCTCTTTg                     >  1:136481/1‑142 (MQ=255)
                                                                                                                                      aaaCCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGCTGTCTGCGTTCTCTTTGCCTTTCa              >  7:97869/1‑149 (MQ=255)
                                                                                                                                           CCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGCGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGCTGTCTGCGTTCTCTTTGCCTTTCACCggg         <  1:400204‑M1/149‑1 (MQ=255)
                                                                                                                                           CCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGCTGTCTGCGTTCTCttt                      >  8:247075‑M1/1‑136 (MQ=255)
                                                                                                                                           CCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGCTGTCTGCGTTCTCttt                      <  7:247075‑M1/136‑1 (MQ=255)
                                                                                                                                           CCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGCTGTCTGCGTTCTCTTTGCCTTTCACCggg         <  4:114306‑M1/149‑1 (MQ=255)
                                                                                                                                             GTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGCTGTCTGCGTTCTCTTTGCCTTTCACCGGGcc       >  2:444009‑M1/1‑149 (MQ=255)
                                                                                                                                                CAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTGGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGCTGTCTGCGTTCTCTTTGCCTTTCACCGGGCCaa     >  4:7233‑M1/1‑148 (MQ=255)
                                                                                                                                                 aaaaGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGCTGTCTGCGTTCTCTTTGCCTTTCACCGGGCCAAt    >  2:119856/1‑148 (MQ=255)
                                                                                                                                                  aaaGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGCTGTCTGCGTTCTCTTTGCCTTTCACCGGGCCAATGa  <  5:434721/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CTTTTTCACGTTGCCTTTCATTACGGTATTGCCCAACGTGGCCTGACCATCACCAGCGATGACCACATGGCCGTTACGGCGTACGCTTACTATAGTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGCTGTCTGCGTTCTCTTTGCCTTTCACCGGGCCAATGA  >  NC_000913/4122186‑4122480

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