Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,213,187:1 +A 5.8% intergenic (+28/‑47) rrfE → / → yjaA 5S ribosomal RNA of rrnE operon/stress‑induced protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,213,1871.A 5.8% 970.6 / 51.3 226intergenic (+28/‑47)rrfE/yjaA5S ribosomal RNA of rrnE operon/stress‑induced protein
Reads supporting (aligned to +/- strand):  ref base . (100/113);  new base A (6/7);  total (106/120)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 228 total aligned reads displayed.

CTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTGCCTTCAGTAAC  >  NC_000913/4213043‑4213329
                                                                                                                                                 |                                                                                                                                              
cTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCgg                                                                                                                                            >  4:1353146/1‑149 (MQ=255)
cTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCgg                                                                                                                                            <  2:675149/149‑1 (MQ=255)
 tGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCTATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCggg                                                                                                                                           <  2:144197/149‑1 (MQ=255)
 tGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATACGG‑CCggg                                                                                                                                           <  4:881443/149‑1 (MQ=255)
 tGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCggg                                                                                                                                           <  2:405988/149‑1 (MQ=255)
 tGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCggg                                                                                                                                           >  3:478048/1‑149 (MQ=255)
 tGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCggg                                                                                                                                           <  3:45596/149‑1 (MQ=255)
 tGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCggg                                                                                                                                           <  4:74587/149‑1 (MQ=255)
  ggcggcCGTAGCGCAGTGGATCCACCTGACCACCTGCCGAACTCAGATGTGTAACGCCGTCGCGCAGAAGGTAGTGTGGGGTCTCCCCATGCGAGACTAGGGAACTGCCAGTAATCAAATTAGAAAAACCCCGATCCATAAGG‑CCgggg                                                                                                                                          >  4:751160/1‑149 (MQ=14)
     ggcAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGG‑CCTTATGGACCGGGGttt                                                                                                                                       <  4:1107693/149‑1 (MQ=14)
        aGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGtttttt                                                                                                                                    >  3:701723/1‑149 (MQ=255)
        aGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAAGTAGAAAAACCCCGG‑CCTTATGGACCGGGGtttttc                                                                                                                                    <  3:807031/149‑2 (MQ=14)
         gTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGG‑CCTTATGGACCGGGGtttttct                                                                                                                                   <  1:1625342/149‑3 (MQ=14)
           aGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCa                                                                                                                                 <  4:1044225/149‑1 (MQ=255)
            gcgcGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCat                                                                                                                                >  2:1290189/1‑149 (MQ=255)
                ggtggtCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCa                                                                                                                            >  2:499046/1‑149 (MQ=255)
                   ggtCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAatt                                                                                                                         <  2:1076286/149‑1 (MQ=255)
                       ccACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTAttt                                                                                                                     >  4:1501235/1‑149 (MQ=255)
                         aCCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGc                                                                                                                   >  1:185099/1‑149 (MQ=255)
                           cTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCat                                                                                                                 <  1:1648182/149‑1 (MQ=255)
                             gACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCatga                                                                                                               >  3:1634444/1‑149 (MQ=255)
                                cccATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCAtgatga                                                                                                            >  3:1463596/1‑149 (MQ=255)
                                 ccATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAat                                                                                                                          <  2:1216013/134‑1 (MQ=255)
                                 ccATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTAttt                                                                                                                     >  4:705414/1‑139 (MQ=255)
                                 ccATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGaa                                                                                                           <  3:529504/149‑1 (MQ=255)
                                 ccATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGaa                                                                                                           <  1:1013943/149‑1 (MQ=255)
                                  cATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTCTTGCATATCAATTATTTGCATGATGAAg                                                                                                          <  3:1326706/149‑1 (MQ=255)
                                   aTGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAgg                                                                                                         <  2:1498693/149‑1 (MQ=255)
                                    tGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCCTGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAggg                                                                                                        <  4:1689855/149‑1 (MQ=255)
                                       cGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAAt                                                                                                     <  2:177855/149‑1 (MQ=255)
                                            tCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCAt                                                                                                >  2:855846/1‑149 (MQ=255)
                                            tCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCAt                                                                                                >  3:123054/1‑149 (MQ=255)
                                              aGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCa                                                                                                                                 <  1:1579126/114‑1 (MQ=255)
                                              aGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGt                                                                                              >  2:1629691/1‑149 (MQ=255)
                                              aGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGt                                                                                              >  1:1545526/1‑149 (MQ=255)
                                                aaGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCa                                                                                            >  1:668407/1‑149 (MQ=255)
                                                  gTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCatgat                                                                                                              >  1:1638980/1‑129 (MQ=255)
                                                     aaaCGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTc                                                                                        <  1:755070/148‑1 (MQ=255)
                                                     aaaCGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTc                                                                                        >  2:1106445/1‑148 (MQ=255)
                                                     aaaCGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTc                                                                                        <  1:1106445/148‑1 (MQ=255)
                                                     aaaCGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTc                                                                                        >  2:755070/1‑148 (MQ=255)
                                                      aaCGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGttt                                                                                                                                       >  4:1614698/1‑100 (MQ=255)
                                                      aaCGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTg                                                                                      >  4:279582/1‑149 (MQ=255)
                                                                      tGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATtcgtc                                                                      <  2:1004649/149‑1 (MQ=255)
                                                                      tGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATtcgtc                                                                      <  2:1682677/149‑1 (MQ=255)
                                                                      tGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATtcgtc                                                                      <  4:1579690/149‑1 (MQ=255)
                                                                                gggTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAAt                                                            <  3:1551851/149‑1 (MQ=255)
                                                                                gggTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAAt                                                            >  3:1363810/1‑149 (MQ=255)
                                                                                    ctcCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACt                                                        <  4:478048/149‑1 (MQ=255)
                                                                                       cccATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGtt                                                     >  4:34528/1‑149 (MQ=255)
                                                                                       cccATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGtt                                                     >  4:910540/1‑149 (MQ=255)
                                                                                       cccATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGtt                                                     >  1:264250/1‑149 (MQ=255)
                                                                                        ccATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAATACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTc                                                    <  4:681098/149‑1 (MQ=255)
                                                                                        ccATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTc                                                    <  1:1684810/149‑1 (MQ=255)
                                                                                        ccATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATCACTGTTc                                                    <  4:1018622/149‑1 (MQ=255)
                                                                                            gCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGTGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCACGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGa                                                <  4:1266060/149‑1 (MQ=255)
                                                                                               agagTAGGGAACTGCCCGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGTGGTTTTTTGCATATCCATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAAGTACTGTACGCGATCt                                             <  2:370321/149‑1 (MQ=255)
                                                                                               agagTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCt                                             <  3:1501235/149‑1 (MQ=255)
                                                                                                  gTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGa                                          <  2:264250/149‑1 (MQ=255)
                                                                                                   tAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGaa                                         >  1:729374/1‑149 (MQ=255)
                                                                                                    aGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGaaa                                        <  3:1086972/149‑1 (MQ=255)
                                                                                                    aGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGaaa                                        <  1:859623/149‑1 (MQ=255)
                                                                                                    aGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGaaa                                        >  2:868738/1‑149 (MQ=255)
                                                                                                    aGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGaaa                                        <  1:1290189/149‑1 (MQ=255)
                                                                                                     gggAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAg                                       >  2:1639209/1‑149 (MQ=255)
                                                                                                      ggAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGc                                      >  4:1147136/1‑149 (MQ=255)
                                                                                                      ggAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGc                                      <  3:579014/149‑1 (MQ=255)
                                                                                                      ggAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGc                                      >  1:1230249/1‑149 (MQ=255)
                                                                                                         aCTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCaaa                                   <  3:664321/149‑1 (MQ=255)
                                                                                                          cTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTACTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAAc                                  <  4:826461/149‑1 (MQ=255)
                                                                                                             ccAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAAc                                  >  3:283335/1‑146 (MQ=255)
                                                                                                               aGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGaa                             <  3:242472/149‑1 (MQ=255)
                                                                                                               aGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGaa                             >  3:623707/1‑149 (MQ=255)
                                                                                                               aGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGaa                             >  1:538029/1‑149 (MQ=255)
                                                                                                               aGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGaa                             <  3:796818/149‑1 (MQ=255)
                                                                                                               aGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGaa                             <  1:139720/149‑1 (MQ=255)
                                                                                                               aGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGaa                             >  2:563425/1‑149 (MQ=255)
                                                                                                                 gCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAgt                           <  1:1220018/149‑1 (MQ=255)
                                                                                                                  cATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAgtg                          >  1:816007/1‑149 (MQ=255)
                                                                                                                  cATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAgtg                          >  1:608187/1‑149 (MQ=255)
                                                                                                                   aTCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAgtgt                         >  1:1407084/1‑149 (MQ=255)
                                                                                                                    tCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTc                        <  3:1181231/149‑1 (MQ=255)
                                                                                                                    tCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTc                        >  4:1715444/1‑149 (MQ=255)
                                                                                                                        aTTAGAAAAACCCCGGACCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGga                    >  4:1177915/1‑149 (MQ=255)
                                                                                                                          tAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGgaga                  <  3:885539/149‑1 (MQ=255)
                                                                                                                           aGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGAc                 <  4:1529116/149‑1 (MQ=255)
                                                                                                                           aGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAACTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGAc                 >  3:1609982/1‑149 (MQ=255)
                                                                                                                            gAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAg                      <  1:964470/143‑1 (MQ=255)
                                                                                                                            gAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAg                      >  2:964470/1‑143 (MQ=255)
                                                                                                                             aaaacCCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATCTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGCCGc               <  2:1703793/149‑1 (MQ=255)
                                                                                                                             aaaaaCCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGc               >  2:629009/1‑149 (MQ=255)
                                                                                                                             aaaaaCCCCGG‑CCTTATGGACCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGAc                 >  3:116540/1‑147 (MQ=255)
                                                                                                                               aaaCCTCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTg             >  4:1426389/1‑149 (MQ=255)
                                                                                                                                 aCCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAg                                                                                                          <  3:1539951/54‑1 (MQ=255)
                                                                                                                                   cccGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAAt                                                                 >  1:713977/1‑93 (MQ=255)
                                                                                                                                   cccGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGATGCTGCCtt         >  2:1035569/1‑149 (MQ=255)
                                                                                                                                   cccGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTGCCtt         >  1:465721/1‑149 (MQ=255)
                                                                                                                                      ggTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTGCCTTCAg      <  4:1182519/149‑1 (MQ=255)
                                                                                                                                             aaGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGTGAATCTCATGTCAGTTCTGTATATCCAAATtcgtc                                                                      <  2:1317337/78‑1 (MQ=255)
                                                                                                                                             aaGG‑CCGGGGTTTTGTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATtcgtc                                                                      >  1:1317337/1‑78 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
CTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGGTCTCCCCATGCGAGAGTAGGGAACTGCCAGGCATCAAATTAGAAAAACCCCGGTCCATAAGG‑CCGGGGTTTTTTGCATATCAATTATTTGCATGATGAAGGGAATCTCATGTCAGTTCTGTATATCCAAATTCGTCGTAATCAAATTACTGTTCGCGATCTTGAAAGCAAACGTGAAGTGTCAGGAGACGCTGCCTTCAGTAAC  >  NC_000913/4213043‑4213329

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑