Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,441,682 A→C 15.2% intergenic (‑77/‑82) 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,6820AC15.2% 219.7 / 34.0 92intergenic (‑77/‑82)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 A (40/38);  new base C (6/8);  total (46/46)
Fisher's exact test for biased strand distribution p-value = 7.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

AGTAATCACTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCAT  >  NC_000913/2441585‑2441780
                                                                                                 |                                                                                                  
aGTAATCACTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATc                                                                                                  >  2:2394724/1‑100 (MQ=255)
 gTAATCACTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCa                                                                                                 <  1:2566893/100‑1 (MQ=255)
  tAATCACTGCGCGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATGGCTTACACTTGTACGCCGAACAAGTCCGATCAg                                                                                                <  2:1996444/100‑1 (MQ=255)
  tAATCACTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGCTCAg                                                                                                >  1:665510/1‑100 (MQ=255)
  tAATCACTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGACGAACAAGTCCGATCAg                                                                                                >  1:2481685/1‑100 (MQ=255)
   aaTCACTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGc                                                                                               <  2:2189130/100‑1 (MQ=255)
   aaTCACTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGc                                                                                               <  2:2188861/100‑1 (MQ=255)
   aaTCACTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGc                                                                                               <  1:2892336/100‑1 (MQ=255)
   aaTCACTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGc                                                                                               <  1:1930376/100‑1 (MQ=255)
   aaTCACTGCACGTTTCATTCAATACATCTGTAAGTCGCACATAGAGTAAGTTTCGAATTCACAATAACGTACACTTGTACGCCGAACAAGTCCGATCAGc                                                                                               >  2:1504053/1‑100 (MQ=255)
    aTCACTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGCTCAGcc                                                                                              >  1:2830020/1‑100 (MQ=255)
    aTCACTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGcc                                                                                              <  2:2725726/100‑1 (MQ=255)
     tCACTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGCTCAGCCa                                                                                             <  1:2681403/100‑1 (MQ=255)
      cACTGCACGTTTCATTCAATACCTCTGTAAGACGCACATAGAGTAAGTTTCGAATGCTCACTAGCGTGCACTTGTACGCCGAACAAGTCCGATCAGCCAt                                                                                            >  2:2130345/1‑100 (MQ=255)
       aCTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCAtt                                                                                           <  1:1703212/100‑1 (MQ=255)
       aCTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCATACACTTGTACGCCAAACAAGTCCGATCAGCCAtt                                                                                           >  1:725412/1‑100 (MQ=255)
        cTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCAttt                                                                                          <  1:2395358/100‑1 (MQ=255)
          gCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGCGTAAGTGTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTaa                                                                                        <  2:2050321/100‑1 (MQ=255)
          gCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTaa                                                                                        >  1:1860729/1‑100 (MQ=255)
            aCGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGCTCAGCCATTTAATa                                                                                      <  1:839519/100‑1 (MQ=255)
            aCGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATa                                                                                      <  1:2684673/100‑1 (MQ=255)
             cGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATag                                                                                     >  2:999357/1‑100 (MQ=255)
             cGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATag                                                                                     <  1:2757047/100‑1 (MQ=255)
               tttCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGCTCAGCCATTTAATagag                                                                                   >  2:1518522/1‑100 (MQ=255)
                ttcattcaATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAgaga                                                                                  <  2:934231/100‑1 (MQ=255)
                  cattcaATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGaaa                                                                                >  1:565365/1‑100 (MQ=255)
                    ttcaATACCTCTGTAAGTCGCACATAAAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCAGATCAGCCATTTAATAGAGATAtt                                                                              >  2:1912280/1‑100 (MQ=255)
                      caATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTg                                                                            >  1:828631/1‑100 (MQ=255)
                         tACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTgcgc                                                                         <  1:723568/100‑1 (MQ=255)
                          aCCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTGATAGAGAAATTTGCGCa                                                                        >  1:2974605/1‑100 (MQ=255)
                           ccTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGCTCAGCCATTTAATAGAGAAATTTGCGCAg                                                                       <  1:113830/100‑1 (MQ=255)
                           ccTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAg                                                                       <  1:992347/100‑1 (MQ=255)
                           ccTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCAGATCAGCCATTTAATAGAGAAATTTGCGCAg                                                                       >  2:1016020/1‑100 (MQ=255)
                             tctGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGcc                                                                     <  2:565362/100‑1 (MQ=255)
                               tGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCtt                                                                   >  1:1501406/1‑100 (MQ=255)
                               tGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCtt                                                                   >  1:1501505/1‑100 (MQ=255)
                                 tAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTac                                                                 >  2:1781005/1‑100 (MQ=255)
                                 tAAGTCGCACATAGAGTAAGGTTCGAATGCACAATAGCGTATACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTa                                                                  <  2:2312063/99‑1 (MQ=255)
                                   aGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACGTGTACGCCGAACAAGTCCGCTCAGCCATTTAATAGAGAAATTTGCGCAGCCTTacac                                                               <  2:2135691/100‑1 (MQ=255)
                                    gTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTacaca                                                              >  2:1955740/1‑100 (MQ=255)
                                     tCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTacacac                                                             <  2:795684/100‑1 (MQ=255)
                                          caTAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGCAATTTGCGCAGCCTTACACACATCGc                                                        <  2:2748480/100‑1 (MQ=255)
                                           aTAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACGAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGc                                                        <  2:2303126/99‑1 (MQ=255)
                                            tAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTa                                                      >  2:2370408/1‑100 (MQ=255)
                                                   aGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGa                                               <  2:2363635/100‑1 (MQ=255)
                                                   aGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATAAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGa                                               >  2:648908/1‑100 (MQ=255)
                                                    gTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAg                                              <  1:2946686/100‑1 (MQ=255)
                                                     tttCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAg                                              <  2:618434/99‑1 (MQ=255)
                                                       tCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCa                                           >  2:2750789/1‑100 (MQ=255)
                                                         gAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCAcc                                         <  1:5618/100‑1 (MQ=255)
                                                          aaTGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGCTCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCg                                        >  2:2637037/1‑100 (MQ=255)
                                                            tGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGcc                                      <  2:712879/100‑1 (MQ=255)
                                                            tGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGcc                                      >  2:1953270/1‑100 (MQ=255)
                                                             cttgttcggCGTACAAGTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCt                                     >  1:57822/9‑100 (MQ=255)
                                                              cacaATAGCGTACACTTGTACGCCGAACAAGTCCGCTCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTg                                    <  1:2141528/100‑1 (MQ=255)
                                                              cacaATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATATGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCTCCTg                                    >  1:1649383/1‑100 (MQ=255)
                                                                caATAGCGTACACTTGTACGCCGAACAAGTCCGCTCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACACCGCTAAGATCGAGCCACCGCCTGTa                                  <  2:149866/100‑1 (MQ=255)
                                                                caATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTa                                  <  1:358252/100‑1 (MQ=255)
                                                                caATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTa                                  >  1:447731/1‑100 (MQ=255)
                                                                caATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTa                                  >  1:2785958/1‑100 (MQ=255)
                                                                caATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTa                                  >  1:1595308/1‑100 (MQ=255)
                                                                 ttcggCGTACAAGTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTaa                                 >  2:2389103/5‑100 (MQ=255)
                                                                     gatccgacTTGTTCGACGTACAAGTCCGATTAACCATTTAATAGAAAAATTTTCGCATCCTTCGACACATCGCTAAGATCAAGCCGCCCCCTGTAAGAc                              >  2:171229/7‑99 (MQ=255)
                                                                     gCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACg                             <  1:662010/100‑1 (MQ=255)
                                                                     gCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTAAAAACAACACTAAGATCAAGCCACCGCCAGTAAGAc                              >  1:1345249/1‑99 (MQ=255)
                                                                      cGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATAGAGCCACCGCCTGTAAGACGa                            >  1:1643094/1‑100 (MQ=255)
                                                                         acacTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTa                         <  1:237690/100‑1 (MQ=255)
                                                                         acacTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTa                         <  2:1279823/100‑1 (MQ=255)
                                                                         acacTTGTACGACGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTa                         <  2:2481656/100‑1 (MQ=255)
                                                                          cacTTGTACGCCGAACAAGTCCGCTCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTaa                        >  1:1207909/1‑100 (MQ=255)
                                                                          cacTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTaa                        >  1:2102137/1‑100 (MQ=255)
                                                                           acTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACCTCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAAc                       <  2:1616258/100‑1 (MQ=255)
                                                                           acTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAAc                       >  2:2347982/1‑100 (MQ=255)
                                                                               gTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTAc                   >  1:1783734/1‑100 (MQ=255)
                                                                                 aCGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGt                 >  1:534945/1‑100 (MQ=255)
                                                                                  cGCCGAACAAGTCCGCTCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTg                >  1:2016465/1‑100 (MQ=255)
                                                                                  cGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTg                >  1:746062/1‑100 (MQ=255)
                                                                                  cGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTg                >  2:1526123/1‑100 (MQ=255)
                                                                                  cGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACAAACATCGCTAAGATCGAGCCACCACCTGTAAGACGAGTAACTTACGTg                >  1:2412802/1‑100 (MQ=255)
                                                                                   gCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGa               >  2:2332153/1‑100 (MQ=255)
                                                                                      gAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAac            <  1:1588284/100‑1 (MQ=255)
                                                                                      agcgaaattcGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAac            >  2:2341957/10‑100 (MQ=255)
                                                                                       aaCAAGTCCGCTCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAaca           <  1:2558005/100‑1 (MQ=255)
                                                                                         cAAGTCCGCTCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACact         <  2:665506/100‑1 (MQ=255)
                                                                                           aGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACactac       <  2:2102114/100‑1 (MQ=255)
                                                                                            gTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACactact      <  2:1888967/100‑1 (MQ=255)
                                                                                            gTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACactact      >  1:2677401/1‑100 (MQ=255)
                                                                                            gTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCAACCGCCTGTAAGACGAGTAAATTACGTGAAACactac       >  2:2832226/1‑99 (MQ=255)
                                                                                              ccGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTcc    <  1:2278574/100‑1 (MQ=255)
                                                                                              ccGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTcc    <  1:2448994/100‑1 (MQ=255)
                                                                                               cGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGTGACCGCCTGTAAGACGAGTAACTTACGTGAAACATGACTCCa   <  2:287013/100‑1 (MQ=255)
                                                                                                gATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCAt  <  2:746058/100‑1 (MQ=255)
                                                                                                 |                                                                                                  
AGTAATCACTGCACGTTTCATTCAATACCTCTGTAAGTCGCACATAGAGTAAGTTTCGAATGCACAATAGCGTACACTTGTACGCCGAACAAGTCCGATCAGCCATTTAATAGAGAAATTTGCGCAGCCTTACACACATCGCTAAGATCGAGCCACCGCCTGTAAGACGAGTAACTTACGTGAAACACTACTCCAT  >  NC_000913/2441585‑2441780

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