Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,441,662 T→G 5.3% intergenic (‑57/‑102) fabB ← / → mnmC beta‑ketoacyl‑[acyl carrier protein] synthase I/fused 5‑methylaminomethyl‑2‑thiouridine‑forming 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,6620TG 5.3% 317.2 / 7.2 95intergenic (‑57/‑102)fabB/mnmCbeta‑ketoacyl‑[acyl carrier protein] synthase I/fused 5‑methylaminomethyl‑2‑thiouridine‑forming methyltransferase and FAD‑dependent demodification enzyme
Reads supporting (aligned to +/- strand):  ref base T (52/38);  new base G (2/3);  total (54/41)
Fisher's exact test for biased strand distribution p-value = 6.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.90e-01

TGAACGTCCTTCACGCAGAGATGCCAGGACTTCCTGCTGGTTATTACCGATGCTGGAAACAATGCCCAGGCCAGTAATCACTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTG  >  NC_000913/2441513‑2441809
                                                                                                                                                     |                                                                                                                                                   
tGAACGTCCTTCACGCAGAGATGCCAGGACTTCCTGCTGGTTATTACCGATGCTGGAAACAATGCCCAGGCCAGTAATCACTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACt                                                                                                                                                     <  1:627487/150‑1 (MQ=255)
tGAACGTCCTTCACGCAGAGATGCCAGGACTTCCTGCTGGTTATTACCGATGCTGGAAACAATGCCCAGGCCAGTAATCACTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACt                                                                                                                                                     <  1:186483/150‑1 (MQ=255)
  aaCGTCCTTCACGCAGAGATGCCAGGACTTCCTGCTGGTTATTACCGATGCTGGAAACAATGCCCAGGCCAGTAATCACTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTg                                                                                                                                                   <  1:1215230/150‑1 (MQ=255)
        cTTCACGCAGAGATGCCAGGACTTCCTGCTGGTTATTACCGATGCTGGAAACAATGCCCAGGCCAGTAATCACTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGcc                                                                                                                                             <  2:435667/150‑1 (MQ=255)
            aCGCAGAGATGCCAGGACTTCCTGCTGGTTATTACCGATGCTGGAAACAATGCCCAGGCCAGTAATCACTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAAc                                                                                                                                         <  1:586645/150‑1 (MQ=255)
              gCAGAGATGCCAGGACTTCCTGCTGGTTATTACCGATGCTGGAAACAATGCCCAGGCCAGTAATCACTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACAGTTGTACGCCGAACaa                                                                                                                                       >  2:277194/1‑150 (MQ=255)
              gCAGAGATGCCAGGACTTCCTGCTGGTTATTACCGATGCTGGAAACAATGCCCAGGCCAGTAATCACTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACAAGTGTACGctattgtg                                                                                                                                       <  2:850339/150‑8 (MQ=255)
              gCAGAGATGCCAGGACTTCCTGCTGGTTATTACCGATGCTGGAAACAATGCCCAGGCCAGTAATCACTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACAAGTGTACGctattgtg                                                                                                                                       <  1:1316628/150‑8 (MQ=255)
                  agaTGCCAGGACTTCCTGCTGGTTATTACCGATGCTGGAAACAATGCCCAGGCCAGTAATCACTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTcc                                                                                                                                   <  1:1096/150‑1 (MQ=255)
                  agaTGCCAGGACTTCCTGCTGGTTATTACCGATGCTGGAAACAATGCCCAGGCCAGTAATCACTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTcc                                                                                                                                   <  2:1104060/150‑1 (MQ=255)
                     tGCCAGGACTTCCTGCTGGTTATTACCGATGCTGGAAACAATGCCCAGGCCAGTAATCACTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGAt                                                                                                                                <  1:585619/150‑1 (MQ=255)
                     tGCCAGGACTTCCTGCTGGTTATTACCGATGCTGGAAACAATGCCCAGGCCAGTAATCACTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGAt                                                                                                                                >  1:253046/1‑150 (MQ=255)
                     tGCCAGGACTTCCTGCTGGTTATTACCGATGCTGGAAACAATGCCCAGGCCAGTAATCACTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGAt                                                                                                                                >  1:1236327/1‑150 (MQ=255)
                      gCCAGGACTTCCTGCTGGTTATTACCGATGCTGGAAACAATGCCCAGGCCAGTAATCACTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATc                                                                                                                               <  1:585598/150‑1 (MQ=255)
                      gCCAGGACTTCCTGCTGGTTATTACCGATGCTGGAAACAATGCCCAGGCCAGCAATCACTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATc                                                                                                                               >  1:1371976/1‑150 (MQ=255)
                         aGGACTTCCTGCTGGTTATTACCGATGCTGGAAACAATGCCCAGGCCAGTAATCACTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGc                                                                                                                            <  1:823411/150‑1 (MQ=255)
                              ttCCTGCTGGGTATTACCGATGCTGCAAACAATGCCCAGGCCAGTAATCACTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCAttt                                                                                                                       <  2:139850/150‑1 (MQ=255)
                                ccTGCTGGTTATTACCGATGCTGGAAACAATGCCCAGGCCAGTAATCACTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCTAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTaa                                                                                                                     >  2:831906/1‑150 (MQ=255)
                                   gctgGTTATTACCGATGCTGTAAACAATGCCCAGGCCAGTAATCACTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATag                                                                                                                  <  2:928623/150‑1 (MQ=255)
                                   gctgGTTATTACCGATGCTGGAAACAATGCCCAGGCCAGTAATCACTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATag                                                                                                                  >  2:315185/1‑150 (MQ=255)
                                   gctgGTTATTACCGATGCTGGAAACAATGCCCAGGCCAGTAATCACTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATag                                                                                                                  >  2:1127256/1‑150 (MQ=255)
                                      ggTTATTACCGATGCTGGAAACAATGCCCAGGCCAGTAATCACTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAgaga                                                                                                               >  2:146411/1‑150 (MQ=255)
                                            taCCGATGCTGGAAACAATGCCCAGGCCAGTAATCACTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTg                                                                                                         >  1:1045280/1‑150 (MQ=255)
                                                gATGCTGGAAACAATGCCCAGGCCAGTAATCACTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCa                                                                                                     <  2:1045280/150‑1 (MQ=255)
                                                   gCTGGAAACAATGCCCAGGCCAGTAATCACTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGcc                                                                                                  >  2:1033255/1‑150 (MQ=255)
                                                    cTGGAAACAATGCCCAGGCCAGTAATCACTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCt                                                                                                 >  2:799626/1‑150 (MQ=255)
                                                     tGGAAACAATGCCCAGGCCAGTAATCACTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTgcgc                                                                                                      >  2:976329/1‑144 (MQ=255)
                                                     tGGAAACAATGCCCAGGCCAGTAATCACTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTgcgc                                                                                                      <  1:976329/144‑1 (MQ=255)
                                                         aaCAATGCCCAGGCCAGTAATCACTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAGTAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTacac                                                                                            <  1:313261/150‑1 (MQ=255)
                                                           cAATGCCCAGGCCAGTAATCACTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTacacac                                                                                          >  1:163880/1‑150 (MQ=255)
                                                              tGCCCAGGCCAGTAATCACTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATc                                                                                       >  2:175650/1‑150 (MQ=255)
                                                               gCCCAGGCCAGTAATCACTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCg                                                                                      >  2:692534/1‑150 (MQ=255)
                                                               gCCCAGGCCAGTAATCACTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTCACACATCGc                                                                                     >  1:325953/1‑150 (MQ=255)
                                                                cccAGGCCAGTAATCACTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGc                                                                                     >  1:339870/1‑150 (MQ=255)
                                                                cccAGGCCAGTAATCACTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGc                                                                                     >  1:1174150/1‑150 (MQ=255)
                                                                  cAGGCCAGTAATCACTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTa                                                                                   <  1:1237117/150‑1 (MQ=255)
                                                                       cAGTAATCACTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATc                                                                              <  2:1468480/150‑1 (MQ=255)
                                                                          tAATCACTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAg                                                                           >  1:1123382/1‑150 (MQ=255)
                                                                          tAATCACTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAg                                                                           >  1:562043/1‑150 (MQ=255)
                                                                           aaTCACTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGc                                                                          >  1:223372/1‑150 (MQ=255)
                                                                              cACTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCAc                                                                       >  1:1023795/1‑150 (MQ=255)
                                                                               aCTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACCCACATCGCTAAGATCGAGCCAcc                                                                      >  1:633494/1‑150 (MQ=255)
                                                                                cTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCAGAACAAGTCAGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCg                                                                     >  2:1502159/1‑150 (MQ=255)
                                                                                  gCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGcc                                                                   <  2:339870/150‑1 (MQ=255)
                                                                                     cGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGt                                                                <  1:484007/150‑1 (MQ=255)
                                                                                     cGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGt                                                                >  1:926246/1‑150 (MQ=255)
                                                                                      gTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTa                                                               >  2:498004/1‑150 (MQ=255)
                                                                                      gTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTa                                                               >  2:630387/1‑150 (MQ=255)
                                                                                       tttCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTaa                                                              <  2:253046/150‑1 (MQ=255)
                                                                                            ttcaATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGa                                                         <  1:1033255/150‑1 (MQ=255)
                                                                                            ttcaATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGa                                                         >  1:340091/1‑150 (MQ=255)
                                                                                             tcaATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAg                                                        <  2:223372/150‑1 (MQ=255)
                                                                                             tcaATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAg                                                        >  1:1530882/1‑150 (MQ=255)
                                                                                             tcaATACCTCTGTAACACGCACATAGAGTAAGATTCGAATGCACAATAGCGTACAGTTGTACGCGGAACAAGTCCGTTCAGCCATTTAATAGCGAAATTTGCGCAGCCTTACACACATCGCTAATATCGAGCCACCGCCTGTAAGACGAg                                                        <  2:86460/150‑1 (MQ=255)
                                                                                              caATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGt                                                       >  1:1394014/1‑150 (MQ=255)
                                                                                               aaTACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTa                                                      >  2:695847/1‑150 (MQ=255)
                                                                                                  aCCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTa                                                                                   <  1:1307090/118‑1 (MQ=255)
                                                                                                  aCCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTa                                                                                   >  2:1307090/1‑118 (MQ=255)
                                                                                                   ccTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACtt                                                  >  2:1287424/1‑150 (MQ=255)
                                                                                                     tctGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTAc                                                >  2:1430904/1‑150 (MQ=255)
                                                                                                       tGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTa                                                      >  2:1086806/1‑142 (MQ=255)
                                                                                                       tGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTa                                                      <  1:1086806/142‑1 (MQ=255)
                                                                                                        gTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTg                                             >  2:561867/1‑150 (MQ=255)
                                                                                                          aaGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGaa                                           <  1:498004/150‑1 (MQ=255)
                                                                                                          aaGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGaa                                           <  1:615705/150‑1 (MQ=255)
                                                                                                           aGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGaaa                                          >  1:445436/1‑150 (MQ=255)
                                                                                                            gTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAac                                         >  2:1152304/1‑150 (MQ=255)
                                                                                                            gTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAac                                         >  2:476762/1‑150 (MQ=255)
                                                                                                            gTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAac                                         >  1:714383/1‑150 (MQ=255)
                                                                                                              cGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAacac                                       >  1:23554/1‑150 (MQ=255)
                                                                                                              cGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGACAAATTTGCGCAGCCTTACACACATCGCTAAGCTCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAaccc                                       >  2:678600/1‑148 (MQ=255)
                                                                                                               gCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACact                                      >  1:136790/1‑150 (MQ=255)
                                                                                                                cacaTAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACacta                                     >  2:196179/1‑150 (MQ=255)
                                                                                                                   aTAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTc                                  >  1:1105335/1‑150 (MQ=255)
                                                                                                                      gagTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCAt                               <  1:315185/150‑1 (MQ=255)
                                                                                                                       agTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATa                              <  1:1152304/150‑1 (MQ=255)
                                                                                                                       agTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATa                              <  1:348158/150‑1 (MQ=255)
                                                                                                                       agTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATa                              <  1:820050/150‑1 (MQ=255)
                                                                                                                       agTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATa                              <  1:825305/150‑1 (MQ=255)
                                                                                                                         tAAGTTTCTAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACa                            <  1:831906/150‑1 (MQ=255)
                                                                                                                             tttCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCt                        <  1:1127256/150‑1 (MQ=255)
                                                                                                                              ttCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTg                       <  2:543950/150‑1 (MQ=255)
                                                                                                                                cGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCCTACAACCTGcc                     >  1:97853/1‑150 (MQ=255)
                                                                                                                                  aaTGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCaa                   <  1:955454/150‑1 (MQ=255)
                                                                                                                                      cacaATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTc               >  2:1335207/1‑150 (MQ=255)
                                                                                                                                       acaATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCg              >  1:927243/1‑150 (MQ=255)
                                                                                                                                        caATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGa             >  2:195093/1‑150 (MQ=255)
                                                                                                                                         ttcggCGTACAAGTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGcc                                                                   <  2:1397709/91‑1 (MQ=255)
                                                                                                                                         ttcggCGTACAAGTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGcc                                                                   >  1:1397709/5‑95 (MQ=255)
                                                                                                                                         ttcggCGTACAAGTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGaa            >  1:1322785/5‑150 (MQ=255)
                                                                                                                                         aaTAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGaa            <  1:175650/150‑1 (MQ=255)
                                                                                                                                          aTAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAAt           <  1:799626/150‑1 (MQ=255)
                                                                                                                                           tAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAAtt          >  2:1072001/1‑150 (MQ=255)
                                                                                                                                             gCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTa        >  2:992231/1‑150 (MQ=255)
                                                                                                                                                 acacTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGc    <  1:692534/150‑1 (MQ=255)
                                                                                                                                                  cacTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCt   <  1:855754/150‑1 (MQ=255)
                                                                                                                                                   acTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTg  >  1:651156/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
TGAACGTCCTTCACGCAGAGATGCCAGGACTTCCTGCTGGTTATTACCGATGCTGGAAACAATGCCCAGGCCAGTAATCACTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCATACAACCTGCCAACCTCGAATTTAATGCTG  >  NC_000913/2441513‑2441809

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