Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,122,333 A→T 5.8% intergenic (‑46/+47) 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,3330AT 5.8% 331.4 / 8.3 104intergenic (‑46/+47)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 = 3.97e-01

Only 100 of 114 total aligned reads displayed.

CTTTTTCACGTTGCCTTTCATTACGGTATTGCCCAACGTGGCCTGACCATCACCAGCGATGACCACATGGCCGTTACGGCGTACGCTTACTATAGTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGCTGTCTGCGTTCTCTTTGCCTTTCACCGGGCCAATGAC  >  NC_000913/4122186‑4122481
                                                                                                                                                   |                                                                                                                                                    
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)
         gTTGCCTTTCATTACGGTATTGCCCAACGTGGCCTGACCATCACCAGCGATGACCACATGGCCGTTACGGCGTACGCTTACTATAGTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACgggg                                                                                                                                             >  2:11354/1‑148 (MQ=255)
         gTTGCCTTTCATTACGGTATTGCCCAACGTGGCCTGACCATCACCAGCGATGACCACATGGCCGTTACGGCGTACGCTTACTATAGTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACgggg                                                                                                                                             >  8:438595/1‑148 (MQ=255)
             ccTTTCATTACGGTATTGCCCAACGTGGCCTGACCATCACCAGCGATGACCACATGGCCGTTACGGCGTACGCTTACTATAGTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAGCCCCGTACTTTTGTACggggtttgt                                                                                                                                        <  2:63951/149‑6 (MQ=255)
              cTTTCATTACGGTATTGCCCAACGTGGCCTGACCATCACCAGCGATGACCACATGGCCGTTACGGCGTACGCTTACTATAGTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAAt                                                                                                                                       >  3:390708/1‑149 (MQ=255)
                   aTTACGGTATTGCCCCACGGGGCCTGACCATCACCCGCGATGACCACATGGCCGTTACGGCGGACGCTTACTATAGTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCaa                                                                                                                                   <  4:228076/148‑1 (MQ=255)
                    ttACGGTATTGCCCAACGTGGCCTGACCATCACCAGCGATGACCACATGGCCGTTACGGCGTACGCTTACTATAGTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCaa                                                                                                                                   <  8:442714/147‑1 (MQ=255)
                    ttACGGTATTGCCCAACGTGGCCTGACCATCACCAGCGATGACCACATGGCCGTTACGGCGTACGCTTACTATAGTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCaa                                                                                                                                   >  7:442714/1‑147 (MQ=255)
                       cGGTATTGCCCAACGTGGCCTGACCATCACCAGCGATGACCACATGGCCGTTACGGCGTACGCTTACTATAGTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATtat                                                                                                                               <  1:11354/148‑1 (MQ=255)
                                ccAACGTGGCCTGACCATCACCAGCGATGACCACATGGCCGTTACGGCGGACGCTTACTATAGTTGGCACGAGCTGCCCCCTTGGTTACGAATACAGATTACAAACCCCGTACAAAAGTACGTGGCATAATGCAATGATAGATgggtgg                                                                                                                     >  5:373517/1‑149 (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                                                                                                           >  1:300840/1‑149 (MQ=255)
                                          cTGACCATCACCAGCGATGACCACATGGCCGTTACGGCGTACGCTTACTATAGTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAgg                                                                                                           <  7:125395/149‑1 (MQ=255)
                                          cTGACCATCACCAGCGATGACCACATGGCCGTTACGGCGTACGCTTACTATAGTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAgg                                                                                                           <  7:335406/149‑1 (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)
                                                   aCCAGCGATGACCACATGGCCGTTACGGCGTACGCTTACTATAGTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAAcc                                                                                                  <  7:438595/149‑1 (MQ=255)
                                                    ccAGCGATGACCACATGGCCGTTACGGCGTACGCTTACTATAGTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAAccc                                                                                                 >  3:245205/1‑149 (MQ=255)
                                                     cAGCGATGACCACATGGCCGTTACGGCGTACGCTTACTATAGTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCCAcccc                                                                                                >  6:434721/1‑149 (MQ=255)
                                                     cAGCGATGACCACATGGCCGTTACGGCGTACGCTTACTATAGTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAAcccc                                                                                                >  2:402367/1‑149 (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)
                                                                  aTGGCCGTTACGGCGTACGCTTACTATAGTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGcc                                                                                    >  8:297611/1‑148 (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)
                                                                                gTACGCTTACTATAGTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTAt                                                                     <  1:454867/149‑1 (MQ=255)
                                                                                       tACTATAGTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGc                                                              >  6:220335/1‑149 (MQ=255)
                                                                                         cTATAGTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCa                                                            >  2:294757/1‑149 (MQ=255)
                                                                                          tataGTTGTCACGAGCTGACCCCTTGTTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGa                                                                                                             <  1:458834/99‑1 (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)
                                                                                                    aCGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGa                                                 >  3:213860/1‑149 (MQ=255)
                                                                                                      gAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCt                                                         >  2:3436/1‑139 (MQ=255)
                                                                                                      gAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCt                                                         <  1:3436/139‑1 (MQ=255)
                                                                                                      gAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGAtt                                               <  5:156374/149‑1 (MQ=255)
                                                                                                      gAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGAtt                                               <  2:396936/149‑1 (MQ=255)
                                                                                                        gCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGa                                             <  2:300840/149‑1 (MQ=255)
                                                                                                        ctcaaaatcccccccatctataattgcattatgCCCCGTACTTTTGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGAtt                                               <  4:400088‑M1/114‑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                                >  3:189605/1‑149 (MQ=255)
                                                                                                                     ttACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGCTGTCTGc                                >  7:220173/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)
                                                                                                                         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                           >  8:421582/1‑138 (MQ=255)
                                                                                                                                     cAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGCTGTCTGCGTtct                           <  7:421582/138‑1 (MQ=255)
                                                                                                                                     cAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGCTGTCTGCGTTCTCTTTGCCTTTc                <  7:435575/149‑1 (MQ=255)
                                                                                                                                      aaaCCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGCGCCGAATgaa                                                                             >  2:363372/1‑85 (MQ=255)
                                                                                                                                      aaaCCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGCTGTCTGCGTTCTCTTTg                      >  1:136481/1‑142 (MQ=255)
                                                                                                                                      aaaCCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGCTGTCTGCGTTCTCTTTg                      <  2:136481/142‑1 (MQ=255)
                                                                                                                                      aaaCCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGCTGTCTGCGTTCTCTTTGCCTTTCa               >  7:97869/1‑149 (MQ=255)
                                                                                                                                           CCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGCTGTCTGCGTTCTCttt                       >  8:247075‑M1/1‑136 (MQ=255)
                                                                                                                                           CCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGCTGTCTGCGTTCTCttt                       <  7:247075‑M1/136‑1 (MQ=255)
                                                                                                                                           CCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGCTGTCTGCGTTCTCTTTGCCTTTCACCggg          <  1:444009‑M1/149‑1 (MQ=255)
                                                                                                                                           CCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGCTGTCTGCGTTCTCTTTGCCTTTCACCggg          <  4:114306‑M1/149‑1 (MQ=255)
                                                                                                                                             GTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGCTGTCTGCGTTCTCTTTGCCTTTCACCGGGcc        >  2:444009‑M1/1‑149 (MQ=255)
                                                                                                                                                CAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGCTGTCTGCGTTCTCTTTGCCTTTCACCGGGCCaa      >  7:232192‑M1/1‑148 (MQ=255)
                                                                                                                                                CAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTGGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGCTGTCTGCGTTCTCTTTGCCTTTCACCGGGCCaa      >  4:7233‑M1/1‑148 (MQ=255)
                                                                                                                                                 aaaaGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGCTGTCTGCGTTCTCTTTGCCTTTCACCGGGCCAAt     >  2:119856/1‑148 (MQ=255)
                                                                                                                                                  aaaGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGCTGTCTGCGTTCTCTTTGCCTTTCACCGGGCCAATGa   <  5:434721/149‑1 (MQ=255)
                                                                                                                                                   aaGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGCTGTCTGCGTTCTCTTTGCCTTTCACCGGGCCAATGAc  <  3:12126/149‑1 (MQ=255)
                                                                                                                                                   aaGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGCTGTCTGCGTTCTCTTTGCCTTTCACCGGGCCAATGAc  >  1:398964/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CTTTTTCACGTTGCCTTTCATTACGGTATTGCCCAACGTGGCCTGACCATCACCAGCGATGACCACATGGCCGTTACGGCGTACGCTTACTATAGTTGTCACGAGCTGACCCCTTGGTTACGAATACAGAGTACAAACCCCGTACAAAAGTACGGGGCATAATGCAATTATAGATGGGGGGGATTTTGAGGGTTTCAACCCCCGGCGGCGAGCCGAATGCAGTTTGTATGACCCGCCATCTTCAACCGATTGAGGGTGCTGTCTGCGTTCTCTTTGCCTTTCACCGGGCCAATGAC  >  NC_000913/4122186‑4122481

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