Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,439,845 T→A 16.1% intergenic (+52/+27) ytfK → / ← ytfL DUF1107 family protein/UPF0053 family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,439,8450TA16.1% 281.3 / 43.0 111intergenic (+52/+27)ytfK/ytfLDUF1107 family protein/UPF0053 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (40/53);  new base A (8/10);  total (48/63)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.11e-01

Only 100 of 113 total aligned reads displayed.

TCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGAT  >  NC_000913/4439746‑4439943
                                                                                                   |                                                                                                  
tCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGAt                                                                                                    >  1:694893/1‑100 (MQ=255)
 ccGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGAtt                                                                                                   >  1:950543/1‑100 (MQ=255)
    aaGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAg                                                                                                <  1:1719739/100‑1 (MQ=255)
    aaGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAg                                                                                                <  2:2108503/100‑1 (MQ=255)
     aGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAgg                                                                                               <  2:2198200/100‑1 (MQ=255)
      gTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGa                                                                                              >  1:1542764/1‑100 (MQ=255)
      gTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGAATCAGGa                                                                                              <  1:728824/100‑1 (MQ=255)
      gTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAAACGGCTCCTGAATCAGGa                                                                                              >  2:1351192/1‑100 (MQ=255)
       ttAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAg                                                                                             <  2:109231/100‑1 (MQ=255)
       ttAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAg                                                                                             <  1:2277531/100‑1 (MQ=255)
       ttAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAg                                                                                             <  2:154909/100‑1 (MQ=255)
       ttAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAg                                                                                             <  2:2218716/100‑1 (MQ=255)
        tAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGc                                                                                            >  2:1022207/1‑100 (MQ=255)
         aaCCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGcc                                                                                           >  1:1563584/1‑100 (MQ=255)
         aaCCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGAATCAGGAGcc                                                                                           <  2:2191330/100‑1 (MQ=255)
         aaCCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGAATCAGGAGcc                                                                                           <  1:2587527/100‑1 (MQ=255)
          aCCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCg                                                                                          >  2:916142/1‑100 (MQ=255)
          aCCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAAACGGCTCCTGAATCAGGAGCCg                                                                                          <  1:2513370/100‑1 (MQ=255)
            cGTCAGGTTATGCGTCTGCAAACATAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGtt                                                                                        <  2:463670/100‑1 (MQ=255)
            cGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGTAGTAATAAGACGGCTCCTGAATCAGGAGCCGtc                                                                                        <  2:1705517/100‑2 (MQ=255)
             gTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAATGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGAATCAGGAGCCGtcn                                                                                       <  1:1065569/100‑3 (MQ=255)
             gTCAGGCTATGCGTCTGCAAACAGAGATGGCTTTACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTg                                                                                       >  2:1719261/1‑100 (MQ=255)
               cAGGTTATGCGTCTGCAAAAAGAGATCGCTTAAACAAAGTGATATGCAGTAAAAAGTGCTATGCAATAATAAGTCGGCTCCTGAATCAGGAGCCGtt                                                                                        <  2:2434760/97‑1 (MQ=255)
                 ggTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGAATCAGGAGCCGtcttatt                                                                                   <  1:1053891/100‑7 (MQ=255)
                   ttATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGAATCAGGAGCCGtcttattac                                                                                 <  1:1210473/100‑9 (MQ=38)
                      tGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTg                                                                              >  2:2062016/1‑100 (MQ=255)
                      tGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTAATTCAGGAGCCGTTGATGTTTCt                                                                               >  1:2272140/1‑99 (MQ=255)
                        cGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTggg                                                                            >  1:1670964/1‑100 (MQ=255)
                          tCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTggggg                                                                          <  1:837657/100‑1 (MQ=255)
                           cTGCAAACAGGGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATGAGACGGCTCCTGATTCAGGTGCCGTTGATGTTTCTggggg                                                                          <  2:687051/99‑1 (MQ=255)
                             gCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAACCGTTGATATATCTGAGGGTTa                                                                       >  2:2473022/1‑100 (MQ=255)
                              cAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTAc                                                                      >  1:1440785/1‑100 (MQ=255)
                               aaaCAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTAcg                                                                     <  1:350523/100‑1 (MQ=255)
                                aaCAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGCGCCGTTGATGTTTCTGGGGGTTAcgc                                                                    >  2:2129800/1‑100 (MQ=255)
                                aaCAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTAcgc                                                                    >  1:226416/1‑100 (MQ=255)
                                aaCAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTAcgc                                                                    >  2:578832/1‑100 (MQ=255)
                                aaCAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTAcgc                                                                    >  1:2182004/1‑100 (MQ=255)
                                    gagaTGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACg                                                                >  2:1274950/1‑100 (MQ=255)
                                         ggCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTttc                                                           >  1:2432932/1‑100 (MQ=255)
                                          gCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTttct                                                          <  1:1130116/100‑1 (MQ=255)
                                          gCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTttct                                                          <  2:1300781/100‑1 (MQ=255)
                                           cTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTttctt                                                         >  1:1876559/1‑100 (MQ=255)
                                           cTTAACCAAAGTGCTATGCAGTAAAAAGCGCTATGCAGTAATAAGACGGCTCCTGCTTCAGGAGCCGTTGCTGTTTCTGGGGGTTACGCGACGCTttct                                                          >  1:2498993/1‑99 (MQ=255)
                                            ttAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTttcttc                                                        <  1:1799927/100‑1 (MQ=255)
                                            ttAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTttcttc                                                        <  1:147547/100‑1 (MQ=255)
                                              aaCCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTcttctt                                                      <  2:1805450/100‑1 (MQ=255)
                                               aCCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTtcttt                                                     >  1:162048/1‑100 (MQ=255)
                                               aCCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTtcttt                                                     <  2:1291888/100‑1 (MQ=255)
                                                  aaagtgctatgcagtaaaaagtgctatgcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTtctttatc                                                  <  2:300402/100‑1 (MQ=255)
                                                    agtgctatgcagtaaaaagtgctatgcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTtctttatctn                                                <  2:1670955/100‑2 (MQ=255)
                                                       gctatgcagtaaatagtgctatgcagtaaTAAGTCGGCTCCTTATTCAGGAGCCGTTGATGTTTCTGTGGGTTACGCGACGCTTTCTTCTTTATCTTTAg                                             <  2:834460/100‑1 (MQ=255)
                                                       gctatgcagtaaaaagtgctatgcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAg                                             <  1:1577461/100‑1 (MQ=255)
                                                         tatgcagtaaaaagtgctatgcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCg                                           >  2:929559/1‑100 (MQ=255)
                                                         tatgcagtaaaaagtgctatgcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCg                                           <  1:855240/100‑1 (MQ=255)
                                                            gcagtaaaaagtgctatgcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCa                                        <  2:1772559/100‑1 (MQ=255)
                                                            gcagtaaaaagtgctatgcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCa                                        <  2:1038802/100‑1 (MQ=255)
                                                             cagtaaaaagtgctatgcagtaaGAAGACGGCTCCTTATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAg                                       >  2:521461/1‑100 (MQ=255)
                                                               gtaaaaagtgctatgcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACTCGACGCTTTCTTCTTTATCTTTAGCGTCAGGc                                     <  1:2467384/100‑1 (MQ=255)
                                                               gtaaaaagtgctatgcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGc                                     >  1:663243/1‑100 (MQ=255)
                                                                   aaagtgctatgcactaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGtt                                 >  1:1769653/1‑100 (MQ=255)
                                                                    aagtgctatgcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGttt                                >  1:220211/1‑100 (MQ=255)
                                                                     agtgctatgcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGtttt                               >  2:1680248/1‑100 (MQ=255)
                                                                         ctatgcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCg                           <  2:462780/100‑1 (MQ=255)
                                                                           atgcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGaa                         >  1:1908497/1‑100 (MQ=255)
                                                                           atgcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGaa                         >  1:2286910/1‑100 (MQ=255)
                                                                             gcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAg                       <  1:1128556/100‑1 (MQ=255)
                                                                             gcagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAg                       >  1:121790/1‑100 (MQ=255)
                                                                              cagtaaTAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAgg                      <  1:1101780/100‑1 (MQ=255)
                                                                                 taataaGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCg                   <  2:938078/100‑1 (MQ=255)
                                                                                   catcaaCGGCTCCTGAATCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGGCAGGCAGTTTTGGCGAAAGGGCGGt                 >  2:87036/6‑100 (MQ=39)
                                                                                   ataaGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGt                 >  2:1718025/1‑100 (MQ=255)
                                                                                    taaGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTg                <  1:150303/100‑1 (MQ=255)
                                                                                    atcaaCGGCTCCTGAATCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTg                >  1:3495/5‑100 (MQ=255)
                                                                                     aaGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTgg               <  1:762412/100‑1 (MQ=255)
                                                                                      aGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGc              >  2:2053710/1‑100 (MQ=255)
                                                                                      aGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGc              >  2:2053325/1‑100 (MQ=255)
                                                                                      aGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGc              >  1:96060/1‑100 (MQ=255)
                                                                                      aGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGACAGGCAGTTTTGGCGaaa                        >  2:704813/1‑90 (MQ=255)
                                                                                       gACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGcc             >  1:273279/1‑100 (MQ=255)
                                                                                       gACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGcc             <  1:29022/100‑1 (MQ=255)
                                                                                       aaCGGCTCCTGAATCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGcc             >  2:1102582/2‑100 (MQ=255)
                                                                                        aCGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCt            <  2:1177545/100‑1 (MQ=255)
                                                                                        aCGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTCTGGCGAAAGGGCGGTGGCCt            <  1:628286/100‑1 (MQ=255)
                                                                                        aCGGCTCCTGAATCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCt            >  2:2240246/1‑100 (MQ=255)
                                                                                           gCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGc         <  1:605362/100‑1 (MQ=255)
                                                                                           gCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGc         <  2:684232/100‑1 (MQ=255)
                                                                                           gCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGc         <  2:805941/100‑1 (MQ=255)
                                                                                           gCTCCTGAATCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCGTCTTTATCTTTAGCGTCCGCCAGTTTTGGCGCAAGGGCGGTGGCCTTGc         >  2:496221/1‑100 (MQ=255)
                                                                                            cTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCt        >  1:1164071/1‑100 (MQ=255)
                                                                                             tCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGTAAGGGCGGTGGCCTTGCTg       <  2:118209/100‑1 (MQ=255)
                                                                                             tCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTg       <  2:1596103/100‑1 (MQ=255)
                                                                                             tCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTg       <  1:870392/100‑1 (MQ=255)
                                                                                               cTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTc     <  2:2289354/100‑1 (MQ=255)
                                                                                               cTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTc     <  1:1776557/100‑1 (MQ=255)
                                                                                               cTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTc     >  1:875511/1‑100 (MQ=255)
                                                                                                tGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCg    <  1:1790663/100‑1 (MQ=255)
                                                                                                tGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCg    >  2:2289078/1‑100 (MQ=255)
                                                                                                 gATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGa   >  1:125887/1‑100 (MQ=255)
                                                                                                  aTTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGAt  <  2:1619500/100‑1 (MQ=255)
                                                                                                  aTTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGAt  >  1:873300/1‑100 (MQ=255)
                                                                                                   |                                                                                                  
TCCGAAGTTAACCGTCAGGTTATGCGTCTGCAAACAGAGATGGCTTAACCAAAGTGCTATGCAGTAAAAAGTGCTATGCAGTAATAAGACGGCTCCTGATTCAGGAGCCGTTGATGTTTCTGGGGGTTACGCGACGCTTTCTTCTTTATCTTTAGCGTCAGGCAGTTTTGGCGAAAGGGCGGTGGCCTTGCTGTCGAT  >  NC_000913/4439746‑4439943

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