Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,324,838 T→A 7.6% Y31* (TAT→TAA) ‡ yciQ → enhancer of membrane protein expression; putative inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,324,8380TA 7.6% 292.9 / 10.1 87Y31* (TAT→TAA) ‡yciQenhancer of membrane protein expression; putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (44/35);  new base A (3/4);  total (47/40)
Fisher's exact test for biased strand distribution p-value = 6.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.72e-01

GGCAGAAACCGACAACGGTTATTGATCTTACCGATGATACGCCGGTCGTGGTGCGTGAAGGCGTAGGTGATGTGAAGCCTTTCTTATAACAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATGGATGCAGCCTGATGGGGACGCTGCACTGCTGAACTATCATATTGTTGGCGTCACCCGTGATGGTATTCCTGAACCCTGGCATCTTGACTGGCATATCGGGTTAATGAGTATTGTCGT  >  NC_000913/1324630‑1325092
                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                              
ggCAGAAACCGACAACGGTTATTGATCTTACCGATGATACGCCGGTCGTGGTGCGTGAAGGCGTAGGTGATGTGAAGCCTTTCTTATAACAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCAACTCTGCAACCAGATGGGGTCCTCGACATCaaa                                                                                                                                                                                              >  2:389493/1‑275 (MQ=255)
     aaaCCGACAACGCTTATTGATCTTACCGATGATACGCCGGTCGTGATGCGTGAAGGCGTAGGTGATGTGAAGCCTTTCTTATAACAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTAACATCTCTAAGTTATGCAGAAAACACGGAGATCCATTCTTATGAAGAAGGGAACTCGCTCTTTGATTTTAAAGCCACTCTGCAACCAGATGGGGTTCTCgcactctgtctct                                                                                                                                                                                          >  4:50113/1‑262 (MQ=255)
         cGACAACGGTTATTGATCTTACCGATGATACGCCGGTCGTGGTGCGTGAAGGCGTAGGTGATGTGAAGCCTTTCTTATAACAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATggg                                                                                                                                                                                                             <  2:50113/251‑1 (MQ=255)
                            tACCGATGATACGCCGGTCGTGGTGCGTGAAGGCGTAGGTGATGTGAAGCCTTTCTTATAACAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACAt                                                                                                                                                                                                  >  4:367211/1‑243 (MQ=255)
                                gatgatACGCCGGTCGTGGTGCGTGAAGGCGTAGGTGATGTGAAGCCTTTCTTATAACAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCaaa                                                                                                                                                                                              <  2:367216/243‑1 (MQ=255)
                                  tgatACGCCGGTCGTGGTGCGTGAAGGCGTAGGTGATGTGAAGCCTTTCTTATAACAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGAGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCCCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGa                                                                                                                                                            >  2:145975/1‑275 (MQ=255)
                                          cGGTCGTGGTGCGTGAAGGCGTAGGTGATGTGAAGCCTTTCTTATAACAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCAc                                                                                                                                                    <  2:384172/275‑1 (MQ=255)
                                            gTCGTGGTGCGTGAAGGCGTAGGTGATGTGAAGCCTTTCTTATAACAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATNTAGACAACACGGAGACTCCTTTTTATGAAGAAGGGCTCTCGCTCTTTGATGTTCAAGCCACCCTGCAACAAGATNGGGTTCTCGACCTCAAAGAAAACATTCATTNNCAGTCGCGAAACCAGCATATTaactatga                                                                                                                                                  >  3:62553/1‑269 (MQ=255)
                                                tggtgCGTGAAGGCGTAGGTGATGTGAAGCCTTTCTTATAACAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGt                                                                                                                                                                                                                                      <  1:49440/187‑1 (MQ=255)
                                                   tgCGTGAAGGCGTAGGTGATGGGAAGCCTTTCTTATAACAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTa                                                                                                                                            <  1:62553/274‑1 (MQ=255)
                                                              gTAGGTGATGTGAAGCCTTTCTTATAACAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATga                                                                                                                                                                                                                                                              <  8:192613/149‑1 (MQ=255)
                                                              gTAGGGGATGTGAAGCCTTTCATATAACAAATTGCGTTCCACGGATGGAAGACTATGGCAAGAAAGTTGCGCTGCATTTTGCTGTTGATAGAAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTGTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTAc                                                                                                                                 <  2:357016/274‑1 (MQ=255)
                                                               tAGGTGATGTGAAGCCTTTCTTATAACAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATgaa                                                                                                                                                                                                                                                             >  6:189383/1‑149 (MQ=255)
                                                                aGGTGATGTGAAGCCTTTCTTATAACAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCa                                                                                                                               <  2:441077/274‑1 (MQ=255)
                                                                 ggTGATGTGAAGCCTTTCTTATAACAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATTCTCTCTAAGTTATGCAGAAAACACTGAGATTCCTTCTTATGAAGTAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGACAATATTTATTTTCAGGCGTGAAAACAGCAGATTAAACAATGATTTTAACGGGATTTACCaa                                                                                                                              >  4:41879/1‑273 (MQ=255)
                                                                  gTGATGTGAAGCCTTTCTTATAACAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATgaagaa                                                                                                                                                                                                                                                          <  5:15045/149‑1 (MQ=255)
                                                                      tgtgAAGCCTTTCTTATAACAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCCGATTAAGCCAGGCTTTTCACGTTATTTCACACGACTag                                                                                                                        >  4:355245/1‑274 (MQ=255)
                                                                        tgAAGCCTTTCTTATAACAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGAtc                                                                                                                                                                                                                                                    <  5:231614/149‑1 (MQ=255)
                                                                         gAAGCCTTTCTTATAACAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCATGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTCATAAGAAggg                                                                                                                                                                                                                                                       >  5:23293/1‑145 (MQ=255)
                                                                         gAAGCCTTTCTTATAACAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTCATAAGAAggg                                                                                                                                                                                                                                                       <  6:23293/145‑1 (MQ=255)
                                                                          aaGCCTTTCTTATAACAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCACTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGAtctc                                                                                                                                                                                                                                                  <  6:196461/149‑1 (MQ=255)
                                                                            gcgtttCTTATAAGAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATGGATg                                                                                                                   <  2:53532/271‑1 (MQ=255)
                                                                                   ttATAACAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTCATAAGAAGGGAtctccgtgt                                                                                                                                                                                                                                             <  5:455318/145‑6 (MQ=255)
                                                                                   ttATAACAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTCATAAGAAGGGAtctccgtgt                                                                                                                                                                                                                                             >  6:455318/1‑140 (MQ=255)
                                                                                      taACAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGt                                                                                                                                                                                                                                      >  7:415453/1‑149 (MQ=255)
                                                                                           aaTTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAg                                                                                                                                                                                                                                 >  8:59084/1‑149 (MQ=255)
                                                                                           aaTTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATcagc                                                                                                                                                               <  2:325532/215‑1 (MQ=255)
                                                                                            aTTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATgggg                                                                                                                                                                                                            >  4:543189/1‑169 (MQ=255)
                                                                                             ttGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGTGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGcc                                                                                                                                                                                                                               >  7:400816/1‑149 (MQ=255)
                                                                                               gCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCCGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCAc                                                                                                                                                                                                                             <  8:400816/149‑1 (MQ=255)
                                                                                                cGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCt                                                                                                                                                                                                        <  2:543194/169‑1 (MQ=255)
                                                                                                  ttCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTTTGGATGCAGCCTGATGGGGACGCTGCAct                                                                                            >  2:268909/1‑275 (MQ=255)
                                                                                                     cACGGATGGAAGATTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTACCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTCATaagaag                                                                                                                                                                                                                                                         <  5:232268/115‑1 (MQ=255)
                                                                                                     cACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATGGATGCAGCCTGATGGGGACGCTGCActgct                                                                                         >  2:341559/1‑275 (MQ=255)
                                                                                                     cACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATGGATGCAGCCTGATGAGGACGCTGCActgct                                                                                         >  2:486148/1‑275 (MQ=255)
                                                                                                     cACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTACCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTCATaagaag                                                                                                                                                                                                                                                         >  6:232268/1‑115 (MQ=255)
                                                                                                      aCGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCaa                                                                                                                                                                                                                      >  8:62502/1‑149 (MQ=255)
                                                                                                               aGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATggg                                                                                                                                                                                                             >  5:8831/1‑149 (MQ=255)
                                                                                                                    aTGGCAGGGAAGTTTCGCTGCCTTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAAATTACACACGGCTTGTATCGTGATTTACCACGACTATGGATGCAGCCGTATGGGGCAGCCGCACTGCTGCATTATAATAttg                                                                           >  4:66007/1‑274 (MQ=255)
                                                                                                                    aTGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCa                                                                                                                                                                                                                   >  4:362535/1‑138 (MQ=255)
                                                                                                                       gCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGa                                                                                                                                                                                                     >  6:408797/1‑149 (MQ=255)
                                                                                                                        cAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATg                                                                                                                                                                                                               <  2:362540/138‑1 (MQ=255)
                                                                                                                         aGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACa                                                                                                                                                                                                   <  6:396276/149‑1 (MQ=255)
                                                                                                                           ggAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGAttt                                                                                                                                   >  4:207986/1‑211 (MQ=255)
                                                                                                                            gAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATGGATGCAGCCTGATGGGGACGCTGCACTGCTGAACTATCATATTGTTGGCGTc                                                                   <  2:355250/274‑1 (MQ=255)
                                                                                                                               gTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCa                                                                                                                               <  2:207989/211‑1 (MQ=255)
                                                                                                                                  tCGCTGCATTTTTCTGTTGATAGCAGGGCTGTTTGTATCATCGCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATGGATGCAGCCTGATGGGGACGCTGCACTGCTGAACTATCATATTGTTGGCGTCACCCGTg                                                            <  2:66008/275‑1 (MQ=255)
                                                                                                                                                tGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAg                                                                                                                                                                            >  5:69453/1‑149 (MQ=255)
                                                                                                                                                  ttGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGgc                                                                                                                                                                          >  7:6161/1‑149 (MQ=255)
                                                                                                                                                    gATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGgcgc                                                                                                                                                                        >  6:157230/1‑149 (MQ=255)
                                                                                                                                                       aGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGaa                                                                                                                                                                     >  6:126057/1‑149 (MQ=255)
                                                                                                                                                                 tttGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGAt                                                                                                                                                           <  5:408797/149‑1 (MQ=255)
                                                                                                                                                                       tcatcTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCa                                                                                                                                                     <  8:205778/149‑1 (MQ=255)
                                                                                                                                                                           ctctAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATGGa                                                                                                                     >  3:136690/1‑177 (MQ=255)
                                                                                                                                                                           ccctAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTGTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATGGATGCAGCCTGATGCGGACGCTGCACTTCTGAACTATCATATTGTTCGCGTCACCCGTGATGGCATTCCTTAACCCTGGCATCTTCACTGGCATATCggc                   >  4:308752/3‑274 (MQ=255)
                                                                                                                                                                            tctAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCt                                                                                                                                                >  5:129728/1‑149 (MQ=255)
                                                                                                                                                                             ctAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCtt                                                                                                                                               <  8:415453/149‑1 (MQ=255)
                                                                                                                                                                              tAAGTTATGCAGTAAACACGGATATTCCTTATTAGTAAGAAGGGATCTCGCTCTTTGATGTTTAAGGCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATGGATGCAGCCTGATGGGGACGCTGCACTGCTGAACTATCATATTGTTGGCGTCACCCGTGATGGTATTCCTGAACCCTGGCATCTTGACTGGCATATCGGGTTa                <  4:341555/275‑1 (MQ=255)
                                                                                                                                                                              gtagtttagCAGAAAAACCGGGAATCCCTTCTTAAGAAGAAGGGTTCTCGCGATTTGATGATGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATAAAAGAAAATATTCATTTTCAGGCGCGAAATGAGCGGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATGGATGCAGCCTGATGGGGACGCTGCACTGCTGAACTATCATATTGTTGGCGTCACCCGTGATGGTATTCCTGAACCCTGGCATCTTGACTGGCATATCGGGTTa                <  4:389488/267‑1 (MQ=255)
                                                                                                                                                                               aaGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATGGATGCa                                                                                                                 <  1:136690/177‑1 (MQ=255)
                                                                                                                                                                               aaGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATGGATGCAGCCTGATGGGGACGCTGCACTGCTGAACTATCATATTGTTGGCGTCACCCGTGATgg                                                        <  2:67979/234‑1 (MQ=255)
                                                                                                                                                                                aGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATGGATGCAGCCTGATGGGGACGCTGCACTGCTGAACTATCATATTGTTGGCGTCACCCGTGATGGTATTCCTGAACCCTGGCATCTTGACTGGCATATCGGGTTAAt              <  2:308756/275‑1 (MQ=255)
                                                                                                                                                                                  ttATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATc                                                                                                                                          <  6:129728/149‑1 (MQ=255)
                                                                                                                                                                                     tGCAGTAAACACGGAGCTTCCCTCGTATGAAGAAGGGATCTCGCCCTTTGATGTTGAAGCCACATTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATACATTTTCAGGCGCGTAATCAGCAGactt                                                                                                                                                         >  4:522963/1‑128 (MQ=255)
                                                                                                                                                                                     tGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATGGATGCAGCCTGATGGGGACGc                                                                                                  >  4:403651/1‑186 (MQ=255)
                                                                                                                                                                                     tGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATGGATGCAGCCTGATGGGGACGCTGCACTGCTGAACTATCATATTGTTGGCGTCACCCGGGAAGGTAGTCCTGAACCCTGCCATCTTGACTGGCATAGCGGGTTACTCAGtt         >  2:13307/1‑274 (MQ=255)
                                                                                                                                                                                        aGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGAtt                                                                                                                                    >  6:151161/1‑149 (MQ=255)
                                                                                                                                                                                         gAAAACACGGATATCCCTGCTTATGAAGAAGTGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATTGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCa                                                                                                                                                     <  2:522968/131‑1 (MQ=255)
                                                                                                                                                                                         gAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGAttt                                                                                                                                   >  5:92602/1‑149 (MQ=255)
                                                                                                                                                                                         gAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATGGATGCAGCCTGATGGGGACGCTGCa                                                                                              <  2:403656/186‑1 (MQ=255)
                                                                                                                                                                                            aaCACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATGGATGCAGCCTGATGGGGACGCTGCACTGCTGAACTATCATATTGTTGGCGTCACCCGTGATGGTATTCCTGAACCCTGGCATCTTGACTGGCATATCGGGTTAATGAGTATTGTCGt  >  2:157599/1‑275 (MQ=255)
                                                                                                                                                                                                  gagaTCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATgggg                                                                                                                                                                                                            >  4:190792/1‑67 (MQ=255)
                                                                                                                                                                                                    gaTCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATc                                                                                                                                          >  3:128466/1‑131 (MQ=255)
                                                                                                                                                                                                    gaTCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTAt                                                                                                                        <  6:8831/149‑1 (MQ=255)
                                                                                                                                                                                                     aTCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATg                                                                                                                       >  7:255475/1‑149 (MQ=255)
                                                                                                                                                                                                     aTCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATg                                                                                                                       <  5:189383/149‑1 (MQ=255)
                                                                                                                                                                                                      tCCCTTCTTCATAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATgg                                                                                                                      >  5:150344/1‑149 (MQ=255)
                                                                                                                                                                                                      tCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCt                                                                                                                                                                                                        <  2:190794/67‑1 (MQ=255)
                                                                                                                                                                                                       cccTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGa                                                                                                                                      <  7:329041/132‑1 (MQ=255)
                                                                                                                                                                                                       cccTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGa                                                                                                                                      >  8:329041/1‑132 (MQ=255)
                                                                                                                                                                                                       cccTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATGGa                                                                                                                     <  6:69453/149‑1 (MQ=255)
                                                                                                                                                                                                       cccTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATGGa                                                                                                                     >  5:148344/1‑149 (MQ=255)
                                                                                                                                                                                                        ccTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAATCACGGCTTTTATCGTGATTTACCACGACTATGGAt                                                                                                                    >  8:310550/1‑149 (MQ=255)
                                                                                                                                                                                                        ccTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGa                                                                                                                                      <  1:128466/131‑1 (MQ=255)
                                                                                                                                                                                                            cttATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATGGATGCAg                                                                                                                >  7:50627/1‑149 (MQ=255)
                                                                                                                                                                                                             ttATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATGGATGCAGc                                                                                                               >  5:73863/1‑149 (MQ=255)
                                                                                                                                                                                                                aTGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATTTTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATGGATGCAGCCTGATGGGGACGCTGCACTGCTGAAc                                                                                     >  4:125084/1‑173 (MQ=255)
                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                              
GGCAGAAACCGACAACGGTTATTGATCTTACCGATGATACGCCGGTCGTGGTGCGTGAAGGCGTAGGTGATGTGAAGCCTTTCTTATAACAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATGGATGCAGCCTGATGGGGACGCTGCACTGCTGAACTATCATATTGTTGGCGTCACCCGTGATGGTATTCCTGAACCCTGGCATCTTGACTGGCATATCGGGTTAATGAGTATTGTCGT  >  NC_000913/1324630‑1325092

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: