Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,206,822 A→G 5.7% K376K (AAA→AAG yehM → uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,206,8220AG 5.7% 222.0 / 2.5 70K376K (AAA→AAGyehMuncharacterized protein
Reads supporting (aligned to +/- strand):  ref base A (32/34);  new base G (2/2);  total (34/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAAGCCCTGAACGCGCCGTTGCCGTGGAGCTATCGCGGCGTTATTCATCCCGATACCGATCCGATTCTGCTAACGTTGATAGACACATTAGCGGGTGACGGATTCGGTAAACTTGCCCCTTCTACGCCACAACCGCCTCTGCCAAAAGATGTCACCTGCGAACTGGAACGTACCGCAATCTCCCTTCCGGCGGAGCTTACCTTAAATCGCTTTACCCCCGATGGGCTGGCGCAAAGTCAGGTGTTACATCGGCTGGCAATACTGGAGATCCCTGGGATTGTACGCCAGCAGGG  >  NC_000913/2206676‑2206968
                                                                                                                                                  |                                                                                                                                                  
gaagCCCTGAACGCGCCGTTGCCGTGGAGCTATCGCGGCGTTATTCATCCCGATACCGATCCGATTCTGCTAACGTTGATAGACACATTAGCGGGTGACGGATTCGGTAAACTTGCCCCTTCTACGCCACAACCGCCTCTGCCAAAAGa                                                                                                                                                  <  6:164477/149‑1 (MQ=255)
   gCCCTGAACGCGCCGTTGCCGTGGAGCTATCGCGGCGTTATTCATCCCGATACCGATCCGATTCTGCTAACGTTGATAGACACATTAGCGGGTGACGGATTCGGTAAACTTGCCCCTTCTACGCCACAACCGCCTCTGCCAAAAGATGt                                                                                                                                               <  7:53474/149‑1 (MQ=255)
    cccTGAACGCGCCGTTGCCGTGGAGCTATCGCGGCGTTATTCATCCCGATACCGATCCGATTCTGCTAACGTTGATAGACACATTAGCGGGTGACGGATTCGGTAAACTTGCCCCTTCTACGCCACAACCGCCTCTGCCAAAAGATGTc                                                                                                                                              <  7:215538/149‑1 (MQ=255)
     ccTGAACGCGCCGTTGCCGTGGAGCTATCGCGTCGTTATTCATCCCGATACCGATCCGATTCTGCTAACGTTGATAGACACATTAGCGGGTGACGGATTCGGTAAACTTGCCCCTTCTACGCCACAACCGCCTCTGCCAAAAGATGTCa                                                                                                                                             <  7:46759/149‑1 (MQ=255)
            gcgcCGTTGCCGTGGAGCTATCGCGGCGTTATTCATCCCGATACCGATCCGATTCTGCTAACGTTGATAGACACATTAGCGGGTGACGGATTCGGTAAACTTGCCCCTTCTACGCCACAACCGCCTCTGCCAAAAGATGTCACCTGCg                                                                                                                                       >  8:151788/1‑148 (MQ=255)
             cgcCGTTGCCGTGGAGCTATCGCGGCGTTATTCATCCCGATACCGATCCGATTCTGATAACGTTGATAGACACATTAGCGGGTGACGGATTCGGTAAACTTGCCCCTTCTACGCCACAACCGCCTCTGCCAAAAGATGTCACCTGCGaa                                                                                                                                     <  5:350141/149‑1 (MQ=255)
                 gTTGCCGTGGAGCTATCGCGGCGTTATTCATCCCGATACCGATCCGATTCTGCTAACGTTGATAGACACATTAGCGGGTGACGGATTCGGTAAACTTGCCCCTTCTACGCCACAACCGCCTCTGCCAAAAGATGTCACCTGCGAACTgg                                                                                                                                 <  2:380162/149‑1 (MQ=255)
                       gTGGAGCTATAGCGGCGTTATTCATCCCGATACCGATCCGATTCTGCTAACGTTGATAGACACATTAGCGGGTGACGGATTCGGTAAACTTGCCCCTTCTACGCCACAACCGCCTCTGCCAAAAGATGTCACCTGCGAACTGGAACGTa                                                                                                                           >  5:192161/1‑149 (MQ=255)
                        tGGAGCTATCGCGGCGTTATTCATCCCGATACCGATCCGATTCTGCTAACGTTGATAGACACATTAGCGGGTGTCGGATTCGGTAAACTTGCCCCTTCTACGCCACAACCGCCTCTGCCAAAAGATGTCACCTGCGAACTGGAACGTAc                                                                                                                          <  4:311136/149‑1 (MQ=255)
                          gAGCTATCGCGGCGTTATTCATCCCGATACCGATCCGATTCTGCTAACGTTGATAGACACATTAGCGGGTGACGGATTCGGTAAACTTGCCCCTTCTACGCCACAACCGCCTCTGCCAAAAGATGTCACCTGCGAACTGGAACGTACCg                                                                                                                        <  6:22445/149‑1 (MQ=255)
                             cTATCGCGGCGTTATTCATCCCGATACCGATCCGATTCTGCTAACGTTGATAGACACATTAGCGGGTGACGGATTCGGTAAACTTGCCCCTTCTACGCCTCAACCGCCTCTGCCAAAAGATGTCACCTGCGAACTGGAACGTACCGCaa                                                                                                                     >  2:96939/1‑149 (MQ=255)
                                tCGCGGCGTTATTCATCCCGATACCGATCCGATTCTGCTAACGTTGATAGACACATTAGCGGGTGACGGATTCGGTAAACTTGCCCCTTCTACGCCACACCCGCCTCTGCCAAAAGATGTCACCTGCGCACTGGCACGTACCGCAAtct                                                                                                                  >  4:43895/1‑149 (MQ=255)
                                  gcggcgTTATTCATCCCGATACCGATCCGATTCTGCTAACGTTGATAGACACATTAGCGGGTGACGGATTCGGTAAACTTGCCCCTTCTACGCCACAACCGCCTCTGCCAAAAGATGTCACCTGCGAACTGGAACGTACCGCAATCTcc                                                                                                                >  1:2151/1‑149 (MQ=255)
                                  gcggcgTTATTCATCCCGATACCGATCCGATTCTGCTAACGTTGATAGACACATTAGCGGGTGACGGATTCGGTAAACTTGCCCCTTCTACGCCACAACCGCCTCTGCCAAAAGATGTCACCTGCGAACTGGAACGTACCGCAATCTcc                                                                                                                >  3:93718/1‑149 (MQ=255)
                                  gcggcgTTATTCATCCCGATACCGATCCGATTCTGCTAACGTTGATAGACACATTAGCGGGTGACGGATTCGGTAAACTTGCCCCTTCTACGCCACAACCGCCTCTGCCAAAAGATGTCACCTGCGAACTGGAACGTACCGCAATCTcc                                                                                                                <  7:151788/149‑1 (MQ=255)
                                    ggcgTTATTCATCCCGATACCGATCCGATTCTGCTAACGTTGATAGACACATTAGCGGGTGACGGATTCGGTAAACTTGCCCCTTCTCCGCCACAACCGCCTCTGCCAAAAGATGTCACCTGCGAACTGGAACGTACCGAAATCTCCCt                                                                                                              >  4:211374/1‑149 (MQ=255)
                                        ttATTCATCCCGATACCGATCCGATTCTGCTAACGTTGATAGACACATTAGCGGGTGACGGATTCGGTAAACTTGCCCCTTCTACGCCACAACCGCCTCTGCCAAAAGATGTCACCTGCGAACTGGAACGTACCGCAATCTCCCTTCcg                                                                                                          <  6:192161/149‑1 (MQ=255)
                                         tATTCATCCCGATACCGATCCGATTCTGCTAACGTTGATAGACACATTAGCGGGTGACGGATTCGGTAAACTTGCCCCTTCTACGCCACAACCGCCTCTGCCAAAAGATGTCACCTGCGAACTGGAACGTACCGCAATCTCCCTTCcgg                                                                                                         >  5:98229/1‑149 (MQ=255)
                                              aTCCCGATACCGATCCGATTCTGCTAACGTTGATAGACACATTAGCGGGTGACGGATTCGGTAAACTTGCCCCTTCTACGCCACAACCGCCTCTGCCAAAAGATGTCACCTGCGAACTGGAACGTACCGCAATCTCCCTTCCGGCGGAg                                                                                                    >  6:28546/1‑149 (MQ=255)
                                               tCCCGATACCGATCCGATTCTGCTAACGTTGATAGACACATTAGCGGGTGACGGATTCGGTAAACTTGCCCCTTCTACGCCACAACCGCCTCTGCCAAAAGATGTCACCTGCGAACTGGAACGTACCGCAATCTCCCTTCCGGCGGAGc                                                                                                   <  5:28546/149‑1 (MQ=255)
                                                   gATACCGATCCGATTCTGCTAACGTTGATAGACACATTAGCGGGTGACGGATTCGGTAAACTTGCCCCTTCTACGCCACAACCGCCTCTGCCAAAAGATGTCACCTGCGAACTGGAACGTACCGCAATCTCCCTTCCGGCGGAGCTTAc                                                                                               <  8:341318/149‑1 (MQ=255)
                                                        cgatccgatTCTGCTAACGTTGATAGACACATTAGCGGGTGACGGATTCGGTACACTTGCCCCTTCTACGCCACAACCGCCTCTGCCAAAAGATGTCACCTGCGAACTGGAACGTACCGCAATCTCCCTTCCGGCGGAGCTTACCTTa                                                                                           >  5:169231/1‑148 (MQ=255)
                                                         gatccgatTCTGCTAACGTTGATAGACACATTAGCGGGTGACGGATTCGGTAAACTTGCCCCTTCTACGCCACAACCGCCTCTGCCAAAAGATGTCACCTGCGAACTGGAACGTACCGCAATCTCCCTTCCGGCGGAGCTTACCTTaaa                                                                                         >  5:57025/1‑149 (MQ=255)
                                                         gatccgatTCTGCTAACGTTGATAGACACATTAGCGGGTGACGGATTCGGTAAACTTGCCCCTTCTACGCCACAACCGCCTCTGCCAAAAGATGTCACCTGCGAACTGGAACGTACCGCAATCTCCCTTCCGGCGGAGCTTACCTTaaa                                                                                         >  3:297667/1‑149 (MQ=255)
                                                              gatTCTGCTAACGTTGATAGACACATTAGCGGGTGACGGATTCGGTAAACTTGCCCCTTCTACGCCACAACCGCCTCTGCCAAAAGATGTCACCTGCGAACTGGAACGTACCGCAATCTCCCTTCCGGCGGAGCTTACCTTAAATCGCt                                                                                    <  1:303881/149‑1 (MQ=255)
                                                               atTCTGCTAACGTTGATAGACACATTCGCGGGGGACGGATTCGGTAAACTTGCCCCTTCTACGCCACAACCGCCTCTGCCAAAAGATGTCACCTGCGAACTGGAACGTACCGCAATCTCCCTTCCGGCGGAGCTTACCTTAAATCGCtt                                                                                   <  6:134319/149‑1 (MQ=255)
                                                               atTCTGCTAACGTTGATAGACACATTAGCGGGTGACGGATTCGGTAAACTTGCCCCTTCTACGTCACAACCGCCTCTGCCAAAAGATGTCACCTGCGAACTGGAACGTACCGCAATCTCCCTTCCGGCGGAGCTTACCTTAAATCGCtt                                                                                   >  4:137700/1‑149 (MQ=255)
                                                                   tGCTAACATTGATAGACACATTAGCGGGTGACGGATTCGGTAAACTTGCCCCTTCTACGCCACAACCGCCTCTGCCAAAAGATGTCTCCTGCGAACTGGAACGTACCGCAATCTCCCTTCCGGCGGAGCTTACCTTAAATCGCTTTAcc                                                                               >  6:309630/1‑149 (MQ=255)
                                                                     cTAACGTTGATAGACACATTAGCGGGTGACGGATTCGGTAAACTTGCCCCTTCTACGCCACAACCGCCTCTGCCAAAAGATGTCACCTGCGAACTGGAACGTACCGCAATCTCCCTTCCGGCGGAGCTTACCTTAAATCGCTTTAc                                                                                <  1:60273/146‑1 (MQ=255)
                                                                      tAACGTTGATAGACACATTAGCGGGTGACGGATTCGGTAAACTTGCCCCTTCTACGCCACAACCGCCTCTGCCAAAAGATGTCACCTGCGAACTGGAACGTACCGCAATCTCCCTTCCGGCGGAGCTTACCTTAAATCGCtt                                                                                   <  8:243726/142‑1 (MQ=255)
                                                                      tAACGTTGATAGACACATTAGCGGGTGACGGATTCGGTAAACTTGCCCCTTCTACGCCACAACCGCCTCTGCCAAAAGATGTCACCTGCGAACTGGAACGTACCGCAATCTCCCTTCCGGCGGAGCTTACCTTAAATCGCtt                                                                                   >  7:243726/1‑142 (MQ=255)
                                                                             gATAGACACATTAGCGGGTGACGGATTCGGTAAACTTGCCCCTTCTACGCCACAACCGCCTCTGCCAAAAGATGTCACCTGCGAACTGGAACGTACCGCAATCTCCCTTCCGGCGGAGCTTACCTTAAATCGCTTTACCCCCGATggg                                                                      <  4:297667/148‑1 (MQ=255)
                                                                              aTAGACACATTAGCGGGTGACGGATTCGGTAAACTTGCCCCTTCTACGCCACAACCGCCTCTGCCAAAAGATGTCACCTGCGAACTGGAACGTACCGCAATCTCCCTTCCGGCGGAGCTTACCTTAAATCGCTTTACCCCCGATGGGCt                                                                    >  2:384935/1‑149 (MQ=255)
                                                                               tAGACACATTAGCGGGTGACGGATTCGGTAAACTTGCCCCTTCTACGCCACAACCGCCTCTGCCAAAAGATGTCACCTGCGAACTGGAACGTACCGCAATCTCCCTTCCGGCGGAGCTTACCTTAAATCGCTTTACCCCCGATGGGCTg                                                                   >  4:172042/1‑149 (MQ=255)
                                                                                 gACACATTAGCGGGTGACGGATTCGGTAAACTTGCCCCTTCTACGCCACAACCGCCTCTGCCAAAAGATGTCACCTGCGAACTGGAACGTACCGCAATCTCCCTTCCGGCGGAGCTTACCTTAAATCGCTTTACCCCCGATGGGCTGgc                                                                 <  4:100758/149‑1 (MQ=255)
                                                                                    acaTTAGCGGGTGACGGATTCGGTAAACTTGCCCCTTCTACGCCACAACCGCCTCTGCCAAAAGATGTCACCTGCGAACTGGAACGTACCGCAATCTCCCTTCCGGCGGAGCTTACCTTAAATCGCTTTACCCCCGATGGGCTGgcgc                                                               >  4:258307/1‑148 (MQ=255)
                                                                                      aTTAGCGGGTGACGGATTCGGTAAACTTGCCCCTTCTACGCCACATCCGCCTCTGCCAAAAGATGTCACCTGCGAACTGGAACGTACCGCAATCTCCCTTCCGGCGGAGCTTACCTTAAATCGCTTTACCCCCGATGGGCTGGCGCaaa                                                            <  8:356886/149‑1 (MQ=255)
                                                                                      aTTAGCGGGTGACGGATTCGGTAAACTTGCCCCTTCTACGCCACAACCGCCTCTGCCAAAGGATGTCCCCTGCGAACTGGAACGTACCGCAATCTCCCTTCCGGCGGAGCTTACCTTAAATCGCTTTACCCCCGATGGGCTGGCGCa                                                              >  8:278422/1‑147 (MQ=255)
                                                                                      aTTAGCGGGTGACGGATTCGGTAAACTTGCCCCTTCTACGCCACAACCGCCTCTGCCAAAGGATGTCACCTGCGAACTGGAACGTACCGCAATCTCCCTTCCGGCGGAGCTTACCTTAAATCGCTTTACCCCCGATGGGCTGGCGCa                                                              >  8:278424/1‑147 (MQ=255)
                                                                                      aTTAGCGGGTGACGGATTCGGTAAACTTGCCCCTTCTACGCCACAACCGCCTCTGCCAAAAGATGTCACCTGCGAACTGGAACGTACCGCAATCTCCCTTCCGGCGGAGCTTACCTTAAATCGCTTTACCCCCGATGGGCTGGCGCaaa                                                            <  4:320086/149‑1 (MQ=255)
                                                                                      aTTAGCGGGTGACGGATTCGGTAAACTTGCCCCTTCTACGCCACAACCGCCTCTGCCAAAAGATGTCACCTGCGAACTGGAACGTACCGCAATCTCCCTTCCGGCGGAGCTTACCTTAAATCGCTTTACCCCCGATGGGCTGGCGCaaa                                                            <  4:314768/149‑1 (MQ=255)
                                                                                       ttAGCGGGTGACGGATTCGGTAAACTTGCCCCTTCTACGCCACAACCGCCTCTGCCAAAAGATGTCACCTGCGAACTGGAACGTACCGCAATCTCCCTTCCGGCGGAGCTTACCTTAAATCGCTTTACCCCCGATGGGCTGGCGCAAAg                                                           >  3:325174/1‑149 (MQ=255)
                                                                                            gggTGACGGATTCGGTAAACTTGCCCCTTCTACGCCACAACCGCCTCTGCCAAAAGATGTCACCTGCGAACTGGAACGTACCGCAATCTCCCTTCCGGCGGAGCTTa                                                                                                >  1:155872/1‑107 (MQ=255)
                                                                                            gggTGACGGATTCGGTAAACTTGCCCCTTCTACGCCACAACCGCCTCTGCCAAAAGATGTCACCTGCGAACTGGAACGTACCGCAATCTCCCTTCCGGCGGAGCTTa                                                                                                <  2:155872/107‑1 (MQ=255)
                                                                                            gggTGACGGATTCGGTAAACTTGCCCCTTCTACGCCACAACCGCCTCTGCCAAAAGATGTCACCTGCGAACTGGAACGTACCGCAATCTCCCTTCCGGCGGAGCTTACCTTAAATCGCTTTACCCCCGATGGGCTGGCGCAAAGTCAgg                                                      <  2:152947/149‑1 (MQ=255)
                                                                                                 aCGGATTCGGTAAACTTGCCCCTTCTACGCCACAACCGCCTCTGCCAAAAGATGTCACCTGCGAACTGGAACGTACCGCAATCTCCCTTCCGGCGGAGCTTACCTTAAATCGCTTTACCCCCGATGGGCTGGCGCAAAGTCAGGTGTTa                                                 <  4:325174/149‑1 (MQ=255)
                                                                                                    gATTCGGTAAACTTGCCCCTTCTACGCCACAACCGCCTCTGCCAAAAGATGTCACCTGCGAACTGGAACGTACCGCAATCTCCCTTCCGGCGGAGCTTACCTTAAATCGCTTTACCCCCGATGGGCTGGCGCAAAGTCAGGTGTTACAt                                              <  1:96939/149‑1 (MQ=255)
                                                                                                        cGGTAAACTTGCCCCTTCTACGCAACAACCGCCTCTGCCAAAAGATGTCACCTGCTAACTGGAAAGTACAGCAATATCCCTTCCGGCGGCGCTTACCTTAAATCGATTTACACCAGATCGGCTGGCGAAAAGTCAGGTGTTACATCGGc                                          >  5:389473/1‑149 (MQ=255)
                                                                                                         ggTAAACTTGCCCCTTCTACGCCCCAACCGCCTCTGCCAAAAGATGTCACCTGCGCACTGGAACGTACCGCAATCTCCCTTCCGGCGGAGCTTACCTTAAATCGCTTTACCCCCGATGGGCTGGCGCAAAGTCAGGTGTTACATCGGCt                                         >  2:375366/1‑149 (MQ=255)
                                                                                                          gTAAACTTGCCCCTTCTACGCCACAACCGCCTCTGCCAAAGGATGTCACCTGCGAACTGGAACGTACCGCAATCTCCCTTCCGGCGGAGCTTACCTTAAATCGCTTTACCCCCGATGGGCTGGCGCAAAGTCAGGTGTTACATCGGCTg                                        <  7:278424/149‑1 (MQ=255)
                                                                                                          gTAAACTTGCCCCTTCTACGCCACAACCGCCTCTGCCAAAGGATGTCACCTGCGAACTGGAACGTACCGCAATCTCCCTTCCGGCGGAGCTTACCTTAAATCGCTTTACCCCCGATGGGCTGGCGCAAAGTCAGGTGTTACATCGGCTg                                        <  7:278422/149‑1 (MQ=255)
                                                                                                             aaCTTGCCCCTTCTACGCCACAACCGCCTCTGCCAAAAGATGTCACCTGCGAACTGGAACGTACCGCAATCTCCCTTCCGGCGGAGCTTACCTTAAATCGCTTTACCCCCGATGGGCTGGCGCAAAGTCAGGTGTTACATCGGCTGGc                                      >  4:314581/1‑148 (MQ=255)
                                                                                                                    cccTTCTACGCCACAACCGCCTCTGCCAAAAGATGTCACCTGCGAACTGGAACGTACCGCAATCTCCCTTCCGGCGGAGCTTACCTTAAATCGCTTTACCCCCGATGGGCTGGCGCAAAGTCAGGTGTTACATCGGCTGGCAATACTg                               >  2:264974/1‑148 (MQ=255)
                                                                                                                    cccTTCTACGCCACAACCGCCTCTGCCAAAAGATGTCACCTGCGAACTGGAACGTACCGCAATCTCCCTTCCGGCGGAGCTTACCTTAAATCGATTTACCCCCGATGGGCTGGCGCAAAGTCAGGTGTTACATCGGCAGGCAATACTg                               >  6:189161/1‑148 (MQ=255)
                                                                                                                     ccTTCTACGCCACAACCGCCTCTGCCAAAAGATGTCACCTGCGAACTGGAACGTACCGCAATCTCCCTTCCGGCGGAGCTTACCTTAAATCGCTTTACCCCCGATGGGCTGGCGCAAAGTCAGGTGTTACATCGGCTGGCAATACTGga                             <  3:258307/149‑1 (MQ=255)
                                                                                                                     ccTTCTACGCCACAACCGCCTCTGCCAAAAGATGTCACCTGCGAACTGGAACGTACCGCAATCTCCCTTCCGGCGGAGCTTACCTTAAATCGCTTTACCCCCGATGGGCTGGCGCAAAGTCAGGTGTTACATCGGCTGGCAATACTGga                             <  1:384935/149‑1 (MQ=255)
                                                                                                                     ccTTCTACGCCACAACCGCCTCTGCCAAAAGATGTCACCTGCGAACTGGAACGTACCGCAATCTCCCTTCCGGCGGAGCTTACCTTAAATCGCTTTACCCCCGATGGGCTGGCGCAAAGTCAGGTGTTACATCGGCTGGCAATACTGga                             <  6:57025/149‑1 (MQ=255)
                                                                                                                         cTACGCCACAACCGCCTCTGCCAAAAGATGTCACCTGCGAACTGGAACGTACCGCAATCTCCCTTCCGGCGGAGCTTACCTTAAATCGCTTTACCCCCGATGGGCTGGCGCAAAGTCAGGTGTTACATCGGCTGGCAATACTGGAGATc                         >  2:343897/1‑149 (MQ=255)
                                                                                                                             gCCACAACCGCCTCTGCCAAAAGATGTCACCTGCGAACTGGAACGTACCGCAATCTCCCTTCCGGCGGAGCTTACCTTAAATCGCTTTACCCCCGATGGGCTGGCGCAAAGTCAGGTGTTACATCGGCTGGCAATACTGGAGATCCCTg                     >  3:253669/1‑149 (MQ=255)
                                                                                                                                acaACCGCCTCTGCCAAAAGATGTCACCTGCGAACTGGAACGTACCGCAATCTCCCTTCCGGCGGAGCTTACCTTAAATCGCTTTACCCCCGATGGGCTGGCGCAAAGTCAGGTGTTACATCGGCTGGCAATACTGGAGATCCCTGGGa                  <  1:343897/149‑1 (MQ=255)
                                                                                                                                 caACCGCCTCTGCCAAAAGATGTCACCTGCGAACTGGAACGTACCGCAATCTCCCTTCCGGCGGAGCTTACCTTAAATCGCTTTACCCCCGATGGGCTGGCGCAAAGTCAGGTGTTACATCGGCTGGCAATACTGGAGATCCCTGGGAt                 <  6:98229/149‑1 (MQ=255)
                                                                                                                                    ccGCCTCTGCCAAAAGATGTCACCTGCGAACTGGAACGTACCGCAATCTCCCTTCCGGCGGAGCTTACCTTAAATCGCTTTACCCCCGATGGGCTGGCGCAAAGTCAGGTGTTACATCGGCTGGCAATACTGGAGATCCCTGGGATTg               >  2:271271/1‑148 (MQ=255)
                                                                                                                                     cGCCTCTGCCAAAAGATGTCACCTGCGAACTGGAACGTACCGCAATCTCCCTTCCGGCGGAGCTTACCTTAAATCGCTTTACCCCCGATGGGCTGGCGCAAAGTCAGGTGTTACATCGGCTGGCAATACTGGAGATCCCTGGGATTGTa             <  3:172042/149‑1 (MQ=255)
                                                                                                                                       ccTCTGCCAAAAGATGTCACCTGCGAACTGGAACGTACCGCAATCTCCCTTCCGGCGGAGCTTACCTTAAATCGCTTTACCCCCGATGGGCTGGCGCAAAGTCAGGTGTTACATCGGCTGGCAATACTGGAGATCCCTGGGATTGTa             >  5:207504/1‑147 (MQ=255)
                                                                                                                                       ccTCTGCCAAAAGATGTCACCTGCGAACTGGAACGTACCGCAATCTCCCTTCCGGCGGAGCTTACCTTAAATCGCTTTACCCCCGATGGGCTGGCGCAAAGTCAGGTGTTACATCGGCTGGCAATACTGGAGATCCCTGGGATTGTa             <  6:207504/147‑1 (MQ=255)
                                                                                                                                        ctctGCCAAAAGATGTCACCTGCGAACTGGAACGTACCGCAATCTCCCTTCCGGCGGAGCTTACCTTAAATCGCTTTACCCCCGATGGGCTGGCGCAAAGTCAGGTGTTACATCGGCTGGCAATACTGGAGATCCCTGGGATTGTACGc          >  4:318535/1‑149 (MQ=255)
                                                                                                                                        ctctGCCAAAAGATGTCACCTGCGAACTGGAACGTACCGCAATCTCCCTTCCGGCGGAGCTTACCTTAAATCGCTTTACCCCCGATGGGCTGGCGCAAAGTCAGGTGTTACATCGGCTGGCAATACTGGAGATCCCTGGGATTGTACGc          >  6:254999/1‑149 (MQ=255)
                                                                                                                                           tGCCAAAAGATGTCACCTGCGAACTGGAACGTACCGCAATCTCCCTTCCGGCGGAGCTTACCTTAAATCGCTTTACCCCCGATGGGCTGGCGCAAAGTCAGGTGtt                                                  >  2:90021/1‑106 (MQ=255)
                                                                                                                                           tGCCAAAAGATGTCACCTGCGAACTGGAACGTACCGCAATCTCCCTTCCGGCGGAGCTTACCTTAAATCGCTTTACCCCCGATGGGCTGGCGCAAAGTCAGGTGtt                                                  <  1:90021/106‑1 (MQ=255)
                                                                                                                                             ccAAAAGATGTCACCTGCGAACTGGAACGTACCGCAATCTCCCTTCCGGCGGAGCTTACCTTAAATCGCTTTACCCCCGATGGGCTGGCGCAAAGTCAGGTGTTACATCGGCTGGCAATACTGGAGATCCCTGGGATTGTACGCcagca     <  2:2151/149‑1 (MQ=255)
                                                                                                                                               aaaaGATGTCACCTGCGAACTGGAACGTACCGCAATCTCCCTTCCGGCGGAGCTTACCTTAAATCGCTTTACCCCCGATGGGCTGGCGCAAAGTCAGGTGTTACATCGGCTGGCAATACTGga                             <  1:256941/123‑1 (MQ=255)
                                                                                                                                               aaaaGATGTCACCTGCGAACTGGAACGTACCGCAATCTCCCTTCCGGCGGAGCTTACCTTAAATCGCTTTACCCCCGATGGGCTGGCGCAAAGTCAGGTGTTACATCGGCTGGCAATACTGga                             >  2:256941/1‑123 (MQ=255)
                                                                                                                                                aaaGATGTCACCTGCGAACTGGAACGTACCGCAATCTCCCTTCCGGCGGAGCTTACCTTAAATCGCTTTACCCCCGATGGGCTGGCGCAAAGGCAGGTGTTACATCGGCTGGCAATACTGGAGATCCCTGGGATTGTACGCCAGCAggg  <  3:43895/149‑1 (MQ=255)
                                                                                                                                                aaaGATGTCACCTGCGAACTGGAACGTACCGCAATCTCCCTTCCGGCGGAGCTTACCTTAAATCGCGTTACCCCCGATGGGCTGGCGCAAAGTCAGGTGTTACATCGGCTGGCAATACTGGAGATCCCTGGGATTGTACGCCAGCAggg  <  6:169231/149‑1 (MQ=255)
                                                                                                                                                 aaGATGTCACCTGCGAACTGGAACGTACCGCAATCTCCCTTCCGGCGGAGCTTACCTTAAATCGCTTTACCCCCGATGGGCTGGCGCAAAGTCAGGTGTTACATCGGCTGGCAATACTGGAGATCCCTGGGATTGTACGCCAGCAggg  >  8:327733/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GAAGCCCTGAACGCGCCGTTGCCGTGGAGCTATCGCGGCGTTATTCATCCCGATACCGATCCGATTCTGCTAACGTTGATAGACACATTAGCGGGTGACGGATTCGGTAAACTTGCCCCTTCTACGCCACAACCGCCTCTGCCAAAAGATGTCACCTGCGAACTGGAACGTACCGCAATCTCCCTTCCGGCGGAGCTTACCTTAAATCGCTTTACCCCCGATGGGCTGGCGCAAAGTCAGGTGTTACATCGGCTGGCAATACTGGAGATCCCTGGGATTGTACGCCAGCAGGG  >  NC_000913/2206676‑2206968

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