Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,601,632 A→G 5.2% F152S (TTC→TCC)  ftsX ← putative ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,601,6320AG 5.2% 299.3 / 4.5 95F152S (TTC→TCC) ftsXputative ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (46/44);  new base G (3/2);  total (49/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.68e-01

AATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACTTTCTCCACGCCTTGCTCGGCCTGCAACTGTGCCACCACGCCCGCAGCAGCGTCATCGTCCAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGT  >  NC_000913/3601484‑3601780
                                                                                                                                                    |                                                                                                                                                    
aaTGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGaa                                                                                                                                                     >  5:149142/1‑149 (MQ=255)
      gTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCAc                                                                                                                                                >  6:89730/1‑149 (MQ=255)
           tCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGt                                                                                                                                           >  6:309948/1‑149 (MQ=255)
             tGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGc                                                                                                                                         >  5:260670/1‑149 (MQ=255)
               cGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGc                                                                                                                                         <  8:21003/147‑1 (MQ=255)
               cGTGATACGATCACGCAGCGTATTCAGTGAATCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGc                                                                                                                                         <  2:79043/147‑1 (MQ=255)
                 tGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCt                                                                                                                                     <  8:268193/149‑1 (MQ=255)
                       gATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCAc                                                                                                                                 >  3:994/1‑147 (MQ=255)
                       gATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACga                                                                                                                               <  7:120264/149‑1 (MQ=255)
                       gATCACGCAGCGAATTCAGTGATTCCGTCCCCTTGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACga                                                                                                                               <  5:238581/149‑1 (MQ=255)
                         tCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGGACTCAcc                                                                                                                                               <  6:28657/131‑1 (MQ=255)
                         tCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGGACTCAcc                                                                                                                                               >  5:23286/1‑131 (MQ=255)
                         tCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGGACTCAcc                                                                                                                                               >  5:28657/1‑131 (MQ=255)
                         tCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGGACTCAcc                                                                                                                                               <  6:23286/131‑1 (MQ=255)
                             gCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGc                                                                                                                                         <  6:307643/133‑1 (MQ=255)
                             gCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGc                                                                                                                                         >  5:307643/1‑133 (MQ=255)
                                gCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAg                                                                                                                      <  6:185148/149‑1 (MQ=255)
                                    aTTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCa                                                                                                                  >  6:103111/1‑149 (MQ=255)
                                     ttCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCCGACCAGTTACGCAACTCACCCCGTGCGTCTTCACGAGAACGATAGTTCAc                                                                                                                 <  8:174912/149‑1 (MQ=255)
                                     ttCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCAc                                                                                                                 <  5:309948/149‑1 (MQ=255)
                                      tCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACt                                                                                                                >  4:298687/1‑149 (MQ=255)
                                       cAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACtt                                                                                                               >  8:166810/1‑149 (MQ=255)
                                        aGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACttt                                                                                                              <  1:15829/149‑1 (MQ=255)
                                        aGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGAGCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGACAGATAGATCACttt                                                                                                              >  4:193665/1‑149 (MQ=255)
                                         gTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACTTtc                                                                                                             <  6:308119/149‑1 (MQ=255)
                                           gtttCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCTCTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACTTtctc                                                                                                           <  2:201313/147‑1 (MQ=255)
                                                  tCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACTTTCTCCACGCCt                                                                                                    >  4:353857/1‑149 (MQ=255)
                                                     ccTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACTTTCTCCACGCCTTGc                                                                                                 >  8:29839/1‑149 (MQ=255)
                                                      cTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACTTTCTCCACGCCTTGCt                                                                                                >  5:304926/1‑149 (MQ=255)
                                                      cTGGAAATCGAGTTTCGGCATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACTTTCTCCACGCCTTGc                                                                                                 >  3:56201/1‑148 (MQ=255)
                                                       tGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACTTTCTCCACGCCTTGCTc                                                                                               >  4:273111/1‑149 (MQ=255)
                                                         gAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACTTTCTCCACGCCTTGCTc                                                                                               <  6:260670/147‑1 (MQ=255)
                                                          aaaTCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACTTTCTCCACGCCTTGCTCGGc                                                                                            <  5:377413/149‑1 (MQ=255)
                                                            aTCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACTTTCTCCACGCCTTGCTCGGCCt                                                                                          <  7:166810/149‑1 (MQ=255)
                                                             tCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACTTTCTCCACGCCTTGCTCGGCCTg                                                                                         <  8:323325/149‑1 (MQ=255)
                                                                 gTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACTTTCTCCACGCCTTGCTCGGCCTGCAAc                                                                                     >  5:198814/1‑149 (MQ=255)
                                                                 gTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACTTTCTCCACGCCTTGCTCGGCCTGCAAc                                                                                     >  1:367012/1‑149 (MQ=255)
                                                                   ttCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACTTTCTCCACGCCTTGCTCGGCCTGCAACtg                                                                                   >  7:7748/1‑149 (MQ=255)
                                                                    tCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACTTTCTCCCCGCCTTGCTCGGCCTGCAACtgt                                                                                  >  5:293574/1‑149 (MQ=255)
                                                                    tCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACTTTCTCCACGCCTTGCTCGGCCTGCAACtgt                                                                                  >  5:4333/1‑149 (MQ=255)
                                                                     cGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACTTTCTCCACGCCTTGCTCGGCCTGCAACtgt                                                                                  <  5:89730/148‑1 (MQ=255)
                                                                     cGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACTTTCTCCACGCCTTGCTCGGCCTGCAACtgt                                                                                  <  8:292748/148‑1 (MQ=255)
                                                                     cGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACCTTCTCCACGCCTTGCTCGGCCTGCAACtgt                                                                                  <  5:103111/148‑1 (MQ=255)
                                                                        gATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACTTTCTCCACGCCTTGCTCGGCCTGCAACTGTGcca                                                                              <  4:994/149‑1 (MQ=255)
                                                                           caccacCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACTTTCTCCACGCCTTGCTCGGCCTGCAACTGTGccacca                                                                           >  6:214309/1‑149 (MQ=255)
                                                                                ccGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACTTTCTCCACGCCTTGCTCGGCCTGCAACTGTGCCACCACGccc                                                                      >  2:95012/1‑149 (MQ=255)
                                                                                 cGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACTTTCTCCACGCCTTGCTCGGCCTGCAACTGTGCCACCACGCCCg                                                                     <  6:4333/149‑1 (MQ=255)
                                                                                   ccACTGCCGGAAGCGGGTTTTCTTCCCGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACTTTCTCCACGCCTTGCTCGGCCTGCAACTGTGCCACCACGCCCg                                                                     >  6:178694/1‑147 (MQ=255)
                                                                                    cACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTAACGGAACTCACCCAGTGCGTCTACACGAGAAAGATAGTTCACTTACTCCACGCCTTGCTCGGCCTGCAACTGTGCCACCACGCCCgcag                                                                  >  3:195691/1‑149 (MQ=255)
                                                                                     aCTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGaa                                                                                                                            <  7:373498/90‑1 (MQ=255)
                                                                                     aCTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGaa                                                                                                                            >  8:373498/1‑90 (MQ=255)
                                                                                     aCTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACTTTCTCCACGCCTTGCTCGGCCTGCAACTGTGCCACCACGCCCgcag                                                                  >  3:31266/1‑148 (MQ=255)
                                                                                      cTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACTTTCTCCACGCCTTGCTCGGCCTGCAACTGTGCCACCACGCCCgcagca                                                                <  3:353857/149‑1 (MQ=255)
                                                                                       tGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACTTTCTCCACGCCTTGCTCGGCCTGCAACTGTGCCACCACGCCCGcagcag                                                               >  4:184938/1‑149 (MQ=255)
                                                                                          cGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACTTTCTCCACGCCTTGCTCGGCCTGCAACTGTGCCACCACGCCCGCAGCAGCGt                                                            >  1:216278/1‑149 (MQ=255)
                                                                                           ggAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACTTTCTCCACGCCTTGCTCGGCCTGCAACTGTGCCACCACGCCCGCAGCAGCggc                                                           <  6:318058/149‑3 (MQ=255)
                                                                                            gAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGAt                                                                                                                        >  4:289192/1‑87 (MQ=255)
                                                                                            gAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGAt                                                                                                                        <  3:289192/87‑1 (MQ=255)
                                                                                            gAAGCGGGTTTTCTTCCAACATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACTTTCTCCACGCCTTGCTCGGCCTGCAACTGTGCCACCACGCCCGCAGCAGCGTCa                                                          >  4:234775/1‑149 (MQ=255)
                                                                                             aaGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACTTTCTCCACGCCTTGCTCGGCCt                                                                                          >  4:316460/1‑116 (MQ=255)
                                                                                             aaGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACTTTCTCCACGCCTTGCTCGGCCt                                                                                          <  3:316460/116‑1 (MQ=255)
                                                                                                 cgggtttCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACTTTCTCCACGCCTTGCTCGGCCTGCAACTGTGCCACCACGCCCGCAGCAGCGTCATCGTc                                                     <  6:149142/145‑1 (MQ=255)
                                                                                                            cAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACTTTCTCCACGCCTTGCTCGGCCTGCAACTGTGCCACCACGCCCGCAGCAGCGTCATCGTCCAGCGttttt                                           <  7:29839/148‑1 (MQ=255)
                                                                                                             aGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACTTTCTCCACGCCTTGCTCGGCCTGCAACTGTGCCACCACGCCCGCAGCAGCGTCATc                                                        <  2:372799/134‑1 (MQ=255)
                                                                                                             aGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACTTTCTCCACGCCTTGCTCGGCCTGCAACTGTGCCACCACGCCCGCAGCAGCGTCATc                                                        >  1:372799/1‑134 (MQ=255)
                                                                                                              gCATATCCAGCGCACCACCAAAACCAGAACAGTTACGGGACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACTTTCTCCACGCCTTGCTCGGCCTGCAACTGTGCCCCCACGCCCGCAGCAGCGTCATCGTCCAGCGGTTTTTg                                         >  2:136665/1‑148 (MQ=255)
                                                                                                               cATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACTTTCTCCACGCCTTGCTCGGCCTGCAACTGTGCCACCACGCCCGCAGCAGCGTCATCGTCCAGCGTTTTTTGCa                                       <  3:334843/149‑1 (MQ=255)
                                                                                                               cATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACTTTCTCCACGCCTTGCTCGGCCTGCAACTGTGCCACCACGCCCGCAGCAGCGTCATCGTCCAGCGTTTTTTGCa                                       <  2:367012/149‑1 (MQ=255)
                                                                                                               cATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACTTTCTCCACGCCTTGCTCGGCCTGCAACTGTGCCACCACGCCCGCAGCAGCGTCATCGTCCAGCGTTTTTTGCa                                       <  6:304926/149‑1 (MQ=255)
                                                                                                                   tCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACTTTCTCCACGCCTTGCTCGGCCTGCAACTGTGCCACCACGCCCGCAGCAGCGTCATCGTCCAGc                                                 <  3:303290/135‑1 (MQ=255)
                                                                                                                   tCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACTTTCTCCACGCCTTGCTCGGCCTGCAACTGTGCCACCACGCCCGCAGCAGCGTCATCGTCCAGc                                                 >  4:303290/1‑135 (MQ=255)
                                                                                                                      aGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACTTTCTCCACGCCTTGCTCGGCCTGCAACTGTGCCACCACGCCCGCAGCAGCGTCATCGTCCAGCGTTTTTTGCAGATaaa                                 >  4:353637/1‑148 (MQ=255)
                                                                                                                       gcgcACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACTTTCTCCACGCCTTGCTCGGCCTGCAACTGTGCCACCACGCCCGCAGCAGCGTCATCGTCCAGCGTTTTTTGCAGATaaa                                 >  2:112370/1‑147 (MQ=255)
                                                                                                                       gcgcACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACTTTCTCCACGCCTTGCTCGGCCTGCAACTGTGCCACCACGCCCGCAGCAGCGTCATCGTCCAGCGTTTTTTGCAGATaaa                                 >  2:96229/1‑147 (MQ=255)
                                                                                                                       gcgcACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACTTTCTCCACGCCTTGCTCGGCCTGCAACTGTGCCACCACGCCCGCAGCAGCGTCATCGTCCAGCGTTTTTTGCAGATAAACa                               <  4:56201/149‑1 (MQ=255)
                                                                                                                       gcgcACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACTTTCTCCACGCCTTGCTCGGCCTGCAACTGTGCCACCACGCCCGCAGCAGCGTCATCGTCCAGCGTTTTTTGCAGATAAACa                               <  1:96229/149‑1 (MQ=255)
                                                                                                                          caccacCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACTTTCTCCACGCCTTGCTCGGCCTGCAACTGTGCCACCACGCCCGCAGCAGCGTCATCGTCCAGCGTTTTTTGCAGATAAACAGt                             >  3:220416/1‑148 (MQ=255)
                                                                                                                           cccaccaAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACTTTCTCCACGCCTTGCTCGGCCTGCAACTGTGCCACCACGCCCGCAGCAGCGGCATCGTCCAGCGTTTTTTGCAGATATACAGt                             <  5:327725/146‑1 (MQ=255)
                                                                                                                           accaccAAAACCAGACCAGTTACGGACCTCACCCAGTGCGTATTCACGAGAAAGCTAGTTCACTTTCTCCACGCCTTGCTCGGCCTGCAACTGTCCCACCACGCCCGCAGCAGCGTCATCGTCCCGCGTTTTTTGCAGCTAACCAgg                             >  6:327725/1‑146 (MQ=255)
                                                                                                                             caccaAAACCAGCCCAGGTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACTTTCTCCACGCCTTGCTCGGCCTGCAACTGTGCCACCACGCCCGCAGCAGCGTCATCGTCCAGCGTTTTTTGCAGATAAACAGTGAtt                         >  1:231876/1‑149 (MQ=255)
                                                                                                                             caccaAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACTTTCTCCACGCCTTGCTCGGCCTGCAACTGTGCCACCACGCCCGCAGCAGCGTCATCGTCCAGCGTTTTTTGCAGATAAACa                               <  2:265577/143‑1 (MQ=255)
                                                                                                                             caccaAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACTTTCTCCACGCCTTGCTCGGCCTGCAACTGTGCCACCACGCCCGCAGCAGCGTCATCGTCCAGCGTTTTTTGCAGATAAACa                               >  1:265577/1‑143 (MQ=255)
                                                                                                                             caccaAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACTTTCTCCACGCCTTGCTCGGCCTGCAACTGTGCCACCACGCCCGCAGCAGCGTCATCGTCCAGCGTTTTTTGCAGATAAACAGTGAtt                         >  1:68916/1‑149 (MQ=255)
                                                                                                                              accaAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACTTTCTCCACGCCTTGCTCGGCCTGCAACTGTGCCACCACGCCCGCAGCAGCGTCATCGTCCAGCGTTTTTTGCAGATAAACAGTGAttt                        <  8:7748/149‑1 (MQ=255)
                                                                                                                               ccaAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACTTTCTCCACGCCTTGCTCGGCCTGCAACTGTGCCACCACGCCCGCAGCAGCGTCATCGTCCAGCGTTTTTTGCAGATAAACAGTGATTTg                       >  5:245091/1‑149 (MQ=255)
                                                                                                                                 aaaaCCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACTTTCTCCACGCCTTGCTCGGCCTGCAACTGTGCCACCACGCCCGCAGCAGCGTCATCGTCCAGCGTTTTTTGCAGATAAACAGTGATTTGCg                     <  5:35716/149‑1 (MQ=255)
                                                                                                                                      cagaccagTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACTTTCTCCACGCCTTGCTCGGCCTGCAACTGTGCCACCACGCCCGCAGCAGCGTCATCGTCCAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGAc                >  6:347750/1‑149 (MQ=255)
                                                                                                                                        gaccagTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACTTTCTCCACGCCTTGCTCGGCCTGCAACTGTGCCACCACGCCCGCAGCAGCGTCATCGTCCAGCGTTTTTTg                                         <  7:358345/122‑1 (MQ=255)
                                                                                                                                        gaccagTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACTTTCTCCACGCCTTGCTCGGCCTGCAACTGTGCCACCACGCCCGCAGCAGCGTCATCGTCCAGCGTTTTTTg                                         >  8:358345/1‑122 (MQ=255)
                                                                                                                                         accagTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACTTTCTCCACGCCTTGCTCGGCCTGCAACTGTGCCACCACGCCCGCAGCAGCGTCATCGTCCAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGa             <  4:31266/149‑1 (MQ=255)
                                                                                                                                         accagTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACTTTCTCCACGCCTTGCTCGGCCTGCAACTGTGCCACCACGCCCGCAGCAGCGTCATCGTCCAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGa             <  3:273111/149‑1 (MQ=255)
                                                                                                                                            agTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACTTTCTCCACGCCTTGCTCGGCCTGCAACTGTGCCACCACGCCCGCAGCAGCGTCATCGTCCAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGataa          <  3:184938/149‑1 (MQ=255)
                                                                                                                                             gTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACTTTCTCCACGCCTTGCTCGGCCTGCAACTGTGCCACCACGCCCGCAGCAGCGTCATCGTCCAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGataat         <  3:140529/149‑1 (MQ=255)
                                                                                                                                              ttACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACTTTCTCCACGCCTTGCTCGGCCTGCAACTGTGCCACCACGCCCGCAGCAGCGTCATCGTCCAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGataata        >  4:77111/1‑149 (MQ=255)
                                                                                                                                              ttACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACTTTCTCCACGCCTTGCTCGGCCTGCAACTGTGCCACCACGCCCGCAGCAGCGTCATCGTCCAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGataat         >  5:147532/1‑148 (MQ=255)
                                                                                                                                              ttACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACTTTCTCCACGCCTTGCTCGGCCTGCAACTGTGCCACCACGCCCGCAGCAGCGTCATCGCCAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATAc       >  2:261619/1‑149 (MQ=255)
                                                                                                                                                    aaCTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACTTTCTCCACGCCTTGCTCGGCCTGCaa                                                                                      >  2:254803/1‑65 (MQ=255)
                                                                                                                                                    aaCTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACTTTCTCCACGCCTTGCTCGGCCTGCaa                                                                                      <  1:254803/65‑1 (MQ=255)
                                                                                                                                                    aaCTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACTTTCTCCACGCCTTGCTCGGCCTGCAACTGTGCCACCACGCCCGCAGCAGCGTCATCGTCCAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGt  >  3:77506/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AATGCCGTTAATCTGCGTGATACGATCACGCAGCGTATTCAGTGATTCCGTCCCCTGGAAATCGAGTTTCGGGATCACCACCGCCACTGCCGGAAGCGGGTTTTCTTCCAGCATATCCAGCGCACCACCAAAACCAGACCAGTTACGGAACTCACCCAGTGCGTCTTCACGAGAAAGATAGTTCACTTTCTCCACGCCTTGCTCGGCCTGCAACTGTGCCACCACGCCCGCAGCAGCGTCATCGTCCAGCGTTTTTTGCAGATAAACAGTGATTTGCGGTGACGGATAATACTGCGT  >  NC_000913/3601484‑3601780

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