Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,441,661 C→A 6.5% intergenic (‑56/‑103) fabB ← / → mnmC 3‑oxoacyl‑[acyl‑carrier‑protein] synthase I/fused 5‑methylaminomethyl‑2‑thiouridine‑forming enzyme methyltransferase and FAD‑dependent demodification enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,441,6610CA 6.5% 416.9 / 12.6 119intergenic (‑56/‑103)fabB/mnmC3‑oxoacyl‑[acyl‑carrier‑protein] synthase I/fused 5‑methylaminomethyl‑2‑thiouridine‑forming enzyme methyltransferase and FAD‑dependent demodification enzyme
Reads supporting (aligned to +/- strand):  ref base C (59/52);  new base A (3/5);  total (62/57)
Fisher's exact test for biased strand distribution p-value = 4.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 125 total aligned reads displayed.

TCCTTCACGCAGAGATGCCAGGACTTCCTGCTGGTTATTACCGATGCTGGAAACAATGCCCAGGCCAGTAATCACTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCT  >  NC_000913/2441519‑2441808
                                                                                                                                              |                                                                                                                                                   
   ttCACGCAGAGATGCCAGGACTTCCTGCTGGTTATTACCGATGCTGGAAACAATGCCCAGGCCAGTAATCACTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACAAGTGTACGct                                                                                                                                            >  1:568279/1‑148 (MQ=255)
    tCACGCAGAGATGCCAGGACTTCCTGCTGGTTATTACCGATGCTGGAAACAATGCCCAGGCCAGTAATCACTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACAAGTGTACGcta                                                                                                                                           <  2:568279/149‑3 (MQ=255)
    tCACGCAGAGAAGCCAGGACTTCCTGCTGGTTATTACCGATGCTGGAAACAATGCCCAGGCCAGTAATCACTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACAAGTGTACGcta                                                                                                                                           <  3:863418/149‑3 (MQ=255)
      aCGCAGAGATGCCAGGACTTCCTGCTGGTTATTACCGATGCTGGAAACAATGCCCAGGCCAGTAATCACTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACAAGTGTACGctatt                                                                                                                                         <  4:524462/149‑5 (MQ=255)
          agagATGCCAGGACTTCCTGCTGGTTATTACCGATGCTGGAAACAATGCCCAGGCCAGTAATCACTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTacac                                                                                                                                                     >  3:211241/1‑133 (MQ=255)
                gCCAGGACTTCCTGCTGGTTATTACCGATGCTGGAAACAATGCCCAGGCCAGTAATCACTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGAt                                                                                                                               <  2:554237/149‑1 (MQ=255)
                  cAGGACTTCCTGCTGGTTATTACCGATGCTGGAAACAATGCCCAGGCCAGTAATCACTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCa                                                                                                                             >  4:405366/1‑149 (MQ=255)
                   aGGACTTCCGGCTGGTTATTACCGATGCTGGAAACAATGCCCAAGTCAGTAGTCACTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGATTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAg                                                                                                                            >  2:549319/1‑149 (MQ=255)
                       cTTCCTGCTGGTTATTACCGATGCTGGAAACAATGCCCAGGCCAGTAATCACTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTa                                                                                                                                                <  1:135983/125‑1 (MQ=255)
                       cTTCCTGCTGGTTATTACCGATGCTGGAAACAATGCCCAGGCCAGTAATCACTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCAt                                                                                                                        <  4:566700/149‑1 (MQ=255)
                       cTTCCTGCTGGTTATTACCGATGCTGGAAACAATGCCCAGGCCAGTAATCACTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACAAGTGTACGc                                                                                                                                             >  4:712609/1‑128 (MQ=255)
                        ttCCTGCTGGTTATTACCGATGCTGGAAACAATGCCCAGGCCAGTAATCACTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCAtt                                                                                                                       >  4:377310/1‑149 (MQ=255)
                         tCCTGCTGGTTATTACCGATGCTGGAAACAATGCCCAGGCCAGTAATCACTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCAttt                                                                                                                      <  3:595361/149‑1 (MQ=255)
                           ctgctgGTTATTACCGATGCTGGAAACAATGCCCAGGCCAGTAATCACTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTaa                                                                                                                    <  1:191539/149‑1 (MQ=255)
                           ctgctgGTTATTACCGATGCTGGAAACAATGCCCAGGCCAGTAATCACTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGACCAGCCATTTaa                                                                                                                    <  1:489142/149‑1 (MQ=255)
                               tgGTTATTACCGATGCTGGAAACAATGCCCAGGCCAGTAATCACTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATaga                                                                                                                <  3:377310/149‑1 (MQ=255)
                               tgGTTATTACCGATGCTGGAAACAATGCCCAGGCCAGTAATCACTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATaga                                                                                                                >  1:448702/1‑149 (MQ=255)
                                 gTTATTACCGATGCTGGAAACAATGCCCAGGCCAGTAATCACTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAgaga                                                                                                              <  1:347239/149‑1 (MQ=255)
                                  ttattaCCGATGCTGGAAACAATGCCCAGGCCAGTAATCACTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATagag                                                                                                               <  2:311052/147‑1 (MQ=255)
                                  ttattaCCGATGCTGGAAACAATGCCCAGGCCAGTAATCACTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATagag                                                                                                               >  1:311052/1‑147 (MQ=255)
                                   tattaCCGATGCTGGAAACAATGCCCAGGCCAGTAATCACTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGaaa                                                                                                            >  4:425969/1‑149 (MQ=255)
                                   tattaCCGATGCTGGAAACAATGCCCAGGCCAGTAATCACTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGaaa                                                                                                            >  4:243864/1‑149 (MQ=255)
                                        ccGATGCTGGAAACAATGCCCAGGCCAGTAATCACTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGaa                                                                                                             >  3:399318/1‑143 (MQ=255)
                                        ccGATGCTGGAAACAATGCCCAGGCCAGTAATCACTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTgc                                                                                                       <  2:836749/149‑1 (MQ=255)
                                        ccGATGCTGGAAACAATGCCCAGGCCAGTAATCACTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTgc                                                                                                       <  4:101720/149‑1 (MQ=255)
                                             gCTGGAAACAATGCCCAGGCCAGTAATCACTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGc                                                                                                  <  1:536631/149‑1 (MQ=255)
                                                 gAAACAATGCCCAGGCCAGTAATCACTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTa                                                                                              >  1:19745/1‑149 (MQ=255)
                                                          cccAGGCCAGTAATCACTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACtt                                                                                                                                                   <  3:738851/87‑1 (MQ=255)
                                                           ccAGGCCAGTAATCACTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGc                                                                                                                           <  2:637429/110‑1 (MQ=255)
                                                           ccAGGCCAGTAATCACTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGc                                                                                                                           >  1:637429/1‑110 (MQ=255)
                                                              ggCCAGTAATCACTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTaa                                                                                 <  1:780452/149‑1 (MQ=255)
                                                                  aGTAATCACTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATc                                                                             >  2:492165/1‑149 (MQ=255)
                                                                    tAATCACTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGa                                                                           <  2:850210/149‑1 (MQ=255)
                                                                    tAATCACTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGa                                                                           <  3:62312/149‑1 (MQ=255)
                                                                        cACTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCa                                                                       >  3:36028/1‑149 (MQ=255)
                                                                         aCTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCAc                                                                      <  2:37223/149‑1 (MQ=255)
                                                                         aCTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGAGCAGCCTTACACACATCGCTAAGATCGAGCCAc                                                                      >  4:13239/1‑149 (MQ=255)
                                                                           tGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCg                                                                    >  4:774361/1‑149 (MQ=255)
                                                                             cACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTGTCGAACGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCGATTTAATAGAGAAATTTGCGCGGCCTTACACACATCGCTAAGATCGAGCCACCGcc                                                                  >  4:570254/1‑149 (MQ=255)
                                                                                 tttCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTa                                                              >  2:310730/1‑149 (MQ=255)
                                                                                 tttCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTa                                                              >  4:729903/1‑149 (MQ=255)
                                                                                    cattcaATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGa                                                           <  3:243864/149‑1 (MQ=255)
                                                                                     attcaATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGGGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGAc                                                          >  1:542562/1‑149 (MQ=255)
                                                                                     attcaATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGAc                                                          <  2:647415/149‑1 (MQ=255)
                                                                                      ttcaATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACg                                                         >  2:448719/1‑149 (MQ=255)
                                                                                      ttcaATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACg                                                         >  1:634761/1‑149 (MQ=255)
                                                                                       tcaATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGa                                                        >  1:109342/1‑149 (MQ=255)
                                                                                       tcaATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGa                                                        <  3:425969/149‑1 (MQ=255)
                                                                                       tcaATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGa                                                        <  4:99263/149‑1 (MQ=255)
                                                                                        caATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAg                                                       <  2:761226/149‑1 (MQ=255)
                                                                                         aaTACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCt                                                                                                >  2:650336/1‑107 (MQ=255)
                                                                                         aaTACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCt                                                                                                <  1:650336/107‑1 (MQ=255)
                                                                                         aaTACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGAt                                                                              >  2:763335/1‑125 (MQ=255)
                                                                                         aaTACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGAt                                                                              <  1:763335/125‑1 (MQ=255)
                                                                                         aaTACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGt                                                      >  4:655697/1‑149 (MQ=255)
                                                                                          aTACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTa                                                     >  1:227415/1‑149 (MQ=255)
                                                                                          aTACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTa                                                     >  3:98080/1‑149 (MQ=255)
                                                                                           tACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTaa                                                    >  1:928889/1‑149 (MQ=255)
                                                                                           tACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTaa                                                    <  3:199392/149‑1 (MQ=255)
                                                                                            aCCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAAc                                                   >  1:88180/1‑149 (MQ=255)
                                                                                              ctctGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACtt                                                 >  3:206078/1‑149 (MQ=255)
                                                                                              ctctGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACtt                                                 >  1:615303/1‑149 (MQ=255)
                                                                                               tctGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTaca                                                                                            >  4:420340/1‑105 (MQ=255)
                                                                                               tctGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTaca                                                                                            <  3:420340/105‑1 (MQ=255)
                                                                                                 tGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACg                                              >  4:892077/1‑149 (MQ=255)
                                                                                                  gTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGt                                             <  3:40835/149‑1 (MQ=255)
                                                                                                  gTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGt                                             <  3:530569/149‑1 (MQ=255)
                                                                                                       tCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAac                                        >  4:380567/1‑149 (MQ=255)
                                                                                                         gCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAacac                                      >  1:433813/1‑149 (MQ=255)
                                                                                                           acaTAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACacta                                    <  1:492165/149‑1 (MQ=255)
                                                                                                            caTAGAGTAAGTTTCGAATGCATAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACactac                                   <  1:73001/149‑1 (MQ=255)
                                                                                                            caTAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACactac                                   >  3:71198/1‑149 (MQ=255)
                                                                                                             aTAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTaa                                                             <  1:489729/122‑1 (MQ=255)
                                                                                                             aTAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTaa                                                             >  2:489729/1‑122 (MQ=255)
                                                                                                               agagTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTcc                                <  3:774361/149‑1 (MQ=255)
                                                                                                                gagTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCa                               >  2:242453/1‑149 (MQ=255)
                                                                                                                  gTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACaa                                                                                                                                      >  2:589966/1‑44 (MQ=255)
                                                                                                                  gTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACaa                                                                                                                                      <  1:589966/44‑1 (MQ=255)
                                                                                                                  gTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCtt                                                                                               >  3:453503/1‑83 (MQ=255)
                                                                                                                  gTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATa                             <  1:793263/149‑1 (MQ=255)
                                                                                                                   tAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATc                                                                             >  2:896792/1‑100 (MQ=255)
                                                                                                                   tAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATc                                                                             <  1:896792/100‑1 (MQ=255)
                                                                                                                      gTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAAc                         <  1:824912/149‑1 (MQ=255)
                                                                                                                       tttCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAAcc                        >  4:399264/1‑149 (MQ=255)
                                                                                                                            aaTGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCa                   >  4:151347/1‑149 (MQ=255)
                                                                                                                                cacaATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAg                                                       <  2:762009/109‑1 (MQ=255)
                                                                                                                                cacaATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAg                                                       >  1:762009/1‑109 (MQ=255)
                                                                                                                                 acaATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACactac                                   <  4:724829/128‑1 (MQ=255)
                                                                                                                                 acaATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACactac                                   >  3:724829/1‑128 (MQ=255)
                                                                                                                                   aaTAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGATGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGa            >  4:409887/1‑149 (MQ=255)
                                                                                                                                    aTAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGaa           <  2:928889/149‑1 (MQ=255)
                                                                                                                                    aTAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGaa           <  1:310730/149‑1 (MQ=255)
                                                                                                                                      ggCGTACAAGTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTc                                 <  4:48780/124‑1 (MQ=255)
                                                                                                                                        cGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTa       <  2:750164/149‑1 (MQ=255)
                                                                                                                                         gTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTaa      <  4:97339/149‑1 (MQ=255)
                                                                                                                                          tACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAAt     >  2:643251/1‑149 (MQ=255)
                                                                                                                                           acacTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATg    <  3:380567/149‑1 (MQ=255)
                                                                                                                                             acTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTAc                                               <  2:295172/104‑1 (MQ=255)
                                                                                                                                             acTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTAc                                               >  1:295172/1‑104 (MQ=255)
                                                                                                                                             acTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCt  >  1:389726/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                                   
TCCTTCACGCAGAGATGCCAGGACTTCCTGCTGGTTATTACCGATGCTGGAAACAATGCCCAGGCCAGTAATCACTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCT  >  NC_000913/2441519‑2441808

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