Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913811,7210TG40.9% 126.3 / 18.9 90V67G (GTG→GGG) bioCmalonyl‑ACP O‑methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base T (8/44);  new base G (35/1);  total (45/45)
Fisher's exact test for biased strand distribution p-value = 1.20e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGG  >  NC_000913/811582‑811871
                                                                                                                                            |                                                                                                                                                      
gAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGGACGTCACGCGCAGGGGACgggccta                                                                                                                                              >  2:1073543/1‑145 (MQ=255)
gAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAg                                                                                                                                             <  2:702507/151‑1 (MQ=255)
 aGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGGACGGCACTCGCCGGTGaacgccgt                                                                                                                                               >  1:61170/1‑141 (MQ=255)
  gCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGAGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGAt                                                                                                                                           <  1:421807/151‑1 (MQ=255)
    aaCATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGAtct                                                                                                                                         <  2:63019/151‑1 (MQ=255)
     aCATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGGCACGCGCAGGGGACGGCCTTAGAtct                                                                                                                                         >  1:458061/1‑150 (MQ=255)
     aCATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAAAGTCACGCGCCGGGGACGGccct                                                                                                                                               >  1:987502/1‑142 (MQ=255)
       aTGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGGTCAGCGCAGGGGACGGCCCTAGAtctct                                                                                                                                       >  1:984257/1‑150 (MQ=255)
       aTGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGGCCCGCGCAGGGGGAGGCCCTAg                                                                                                                                             >  2:62614/1‑144 (MQ=255)
              tCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGGGACGGCCTTAGATCTCTCgccgcc                                                                                                                                >  1:356591/1‑150 (MQ=255)
                 aCAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTcgccgc                                                                                                                                 <  2:353754/146‑1 (MQ=255)
                 aCAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGGCCCGCGCCGGGGAGGGCCTTAGATCTCTCGCCGCCaa                                                                                                                              >  2:14849/1‑149 (MQ=255)
                   aGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGGACGTCACGCGCAGGGGACGGCCTTAGATTTCTCCCCGCCCATg                                                                                                                            >  1:267311/1‑149 (MQ=255)
                     cgcCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGGACGGCCCCCGCCGGGGACGGCCCTCGATCTCTCcccccccacgc                                                                                                                           >  2:993259/1‑143 (MQ=255)
                       ccAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGtt                                                                                                                      <  2:549928/151‑1 (MQ=255)
                       ccAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGt                                                                                                                       <  1:968605/150‑1 (MQ=255)
                           agTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGGACGGGCCGCGCAGGGGGCGGCCCCTGAGATCTCTCCGCCCAAACTTGTTCacg                                                                                                                  >  2:532792/1‑149 (MQ=255)
                            gTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGGACGTCACGCGCCGGGGGCGGCCTTTGATATCTCGCCGCCAAAGCTTGGTGAgg                                                                                                                  >  2:723653/1‑150 (MQ=255)
                             tGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGGACGGCACGCGCCGGGGACGGCCTTTGATCTCCCGCCGCCAATGCTTGTTCAgg                                                                                                                  >  2:335544/1‑149 (MQ=255)
                             tGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGAAGGCCCTCGATCTCTCGCCGCCAATGCTTGTTCAGGc                                                                                                                 >  1:703791/1‑150 (MQ=255)
                              gCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCAc                                                                                                               <  2:233389/151‑1 (MQ=255)
                               cTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGGACGCCACGCGCAGGGGACGGCCCTAGATCTCTCGCCGCCCAAGCTTGTTCAgg                                                                                                                  >  2:1069850/1‑147 (MQ=255)
                                tGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACCCCCACGTACTGGAGGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCAc                                                                                                               <  2:1260964/149‑1 (MQ=255)
                                       ttACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACCCGCACGGGGAGGCCCTATATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGaa                                                                                                        >  2:973466/1‑149 (MQ=255)
                                        tACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAg                                                                                                       <  1:404573/149‑1 (MQ=255)
                                             gCAATGCTTCCACAGCGTAAATACACCCCCGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTTACGCGCAGGAGACGGCCTTAGAGTTCTCGCCGCCGATGGTTGTGCCGGCACGCCAGGAGGGTg                                                                                                   >  2:39256/1‑148 (MQ=255)
                                             gCAATGCTTCCACAGCGTAAATACACCCCCGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCg                                                                                                <  2:712814/151‑1 (MQ=255)
                                             gCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAGAGTCCCCGGCCGGGGACGGCCCTACGACTCTTGCCCCCAAAGCTTGTTCTAGCACGcccg                                                                                                          >  2:188219/1‑139 (MQ=255)
                                               aaTGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCa                                                                                              <  1:47479/151‑1 (MQ=255)
                                               aaTGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGGCACGCGCCGGGGACGGCCTTAGATCTCTCGCCGCCCATGCTTGTTTAGGCACGCCAAAAGGATGcc                                                                                                 >  2:509377/1‑148 (MQ=255)
                                                aTGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGGGACGGCCTTAAATCTCTCTCCCCCCAAGCCTGTTCAGGGACGCCAGAAAGATggc                                                                                                 >  2:752836/1‑145 (MQ=255)
                                                 tGCTTCCACAGCGTAAATACACCCACGTCCTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAg                                                                                             <  2:913639/150‑1 (MQ=255)
                                                  gCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGGACGCCACGCGCAGGGGACGGCCTTAGATCTCTCGCCGCCCCTGCTTGTTCAGGCACGCCACAAGGGTGCCGCAGa                                                                                            >  2:725914/1‑150 (MQ=255)
                                                  gCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGGGAAGGCCCTAGAACTCTCGCCGCCAATGCTTTGTCAGGCACGCCCGAAGGATggc                                                                                                 >  2:330681/1‑143 (MQ=255)
                                                  gCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAAAGTCCCGCGCAGGGGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTTAGGCACGCCAGAAGGATGCCGCAg                                                                                             >  2:61382/1‑149 (MQ=255)
                                                    ttCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGGACGTCCCGCGCCGGGGACGGCCCTAGATCTCTCCCCGCCAAAGCTTTTTAAGGCCCGCCCGAAAGGAGCCCCAGAc                                                                                           >  1:1271375/1‑149 (MQ=255)
                                                         caGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTAt                                                                                     <  2:751045/150‑1 (MQ=255)
                                                          aGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCt                                                                                   <  2:604552/151‑1 (MQ=255)
                                                           gCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTg                                                                                  <  1:945921/151‑1 (MQ=255)
                                                           gCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTg                                                                                  <  2:103421/151‑1 (MQ=255)
                                                           gCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGGGAAGGCCTTAGATCTCTCCCCGCCCATGCTTTTTCAGGCACGCCAGAAGGATGCCGCAGACCATTAttt                                                                                   >  1:60046/1‑148 (MQ=255)
                                                               aaaTACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCgg                                                                              <  2:307162/151‑1 (MQ=255)
                                                                  tACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGGGACGGCCCTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCggggg                                                                           >  1:1178863/1‑149 (MQ=255)
                                                                   acacCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGgag                                                                           <  1:1031742/150‑1 (MQ=255)
                                                                    cgcccACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGgag                                                                           <  1:40515/147‑1 (MQ=255)
                                                                       ccACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATc                                                                      <  1:883562/151‑1 (MQ=255)
                                                                           gTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAAt                                                                  <  1:962671/151‑1 (MQ=255)
                                                                           gTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGGGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTTAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGGGATATCGaa                                                                   >  1:1075081/1‑150 (MQ=255)
                                                                           gTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGGGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCCGACCATTATCTGGCGGGGGGTATCGaa                                                                   >  1:341962/1‑150 (MQ=255)
                                                                               tGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAGCCGCCCCCTCCCGAGACAGCCCTACATCTCTCTCCGACAAAGCTTGGTCAGGGACGGCCGAAGGCTGGCGGATAACATTATATGGGGGGGGAACTCGCATccc                                                               >  1:306306/1‑150 (MQ=255)
                                                                                ggACGCGGGTTGTGGACCTGGCTGGATGGGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCt                                                              <  2:721833/150‑1 (MQ=255)
                                                                                ggACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCCGGTGACGGCCCTAGATCTCTCGCCGCCCATGCTTGTTCAGGGACGCCAGAAGGATGGCGCAGACCATTATCTTGCGGGGGGTAACGAATccc                                                               >  2:835679/1‑149 (MQ=255)
                                                                                ggACGCGCGTTGTGGACCTGGCTGGAGGGGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTTTCGGCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCt                                                              <  1:214419/150‑1 (MQ=255)
                                                                                  aCGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGGACGTCACGCGCAGGGGACGGCCCTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGAACATTATCTGGCGGGGGAGATCGAATCCCCGcc                                                           >  1:305873/1‑151 (MQ=255)
                                                                                   cgcgGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGcc                                                           <  2:503551/150‑1 (MQ=255)
                                                                                    gcgGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGGACGTCCCGCGCCGGGGACGGCCTTAGATCCCCCGCCCCCCAAGCCTGTTCAGGCACCCCCGAAGGGCGCCGCCGAACCCTTTCTGGGGGGGGGTATTGaaagcccgcc                                                           >  2:62160/1‑141 (MQ=255)
                                                                                    gcgGGTTGTGGACCTGGCTGGAGGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCg                                                          <  1:1276517/150‑1 (MQ=255)
                                                                                      gggTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGGGACGGCCTTAGATCTCTCGCCGCCCATGCTTGTTCAGGCACGCCCGAAGGATGCCGCAGACCATTATCTGGCGGGGGATATCGAATCCCTGCCGtt                                                        >  1:988456/1‑150 (MQ=255)
                                                                                      gggTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAAAGGCACGCGCCGGGGACGGCCCTAGATCTCTCGCCGCCCAAGCTTGGTTAGGCACGCCCGAAGGGTGCCGCAGACCATTATCTGGCGGGGGACATCGAATCCCCGcc                                                           >  2:451457/1‑147 (MQ=255)
                                                                                        gTTGTGGCCCTGGCTGGATGAGCCGCCACTGGCGGGAGCGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGc                                                     <  1:562041/151‑1 (MQ=255)
                                                                                        gTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAAAGTCCCCCGGAGGGGAAGGCCCTTGATCCCCCCCCCCCACTCACTGTTCAGGCACCGCCCAACGATGTCGCCGACCCCTAACTGGGGGGAGGTatcccatcccccgc                                                           >  2:1218785/1‑141 (MQ=255)
                                                                                        gTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAAAGGCACGCGCAGGGGACGGCCCTAGATCTCTCGCCGCCCAAGCTTGTTTAGGGGCGCCAGAAGGAAGCCGCCGACCATTATCTGGCGGGGGATATTGAATTCCTGCCCTTAGc                                                     >  1:1159699/1‑151 (MQ=255)
                                                                                          tgtgGACCTGTTTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCg                                                    <  2:652530/150‑1 (MQ=255)
                                                                                            gggcccgggctggATGGGCCGCCCCTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACt                                                 <  2:1216447/144‑1 (MQ=255)
                                                                                             ggcccTGGCTGGATGGGCCGCCACTGGCGGGAACGTCACGCGCAGGGGACGGCCTTAGACCTCTCGCCGCAAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGAc                                                  <  2:426908/146‑1 (MQ=255)
                                                                                                 ctggctggATGAGCCGCCACTGGCGGGAACGGCACGCGCAGGGGACGGCCCTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCCGAAGGATGCCGCAGACCCTTATCTGGCGGGGGATATCGAAACCCCGCCGGTAGCGAccg                                                >  2:104394/1‑145 (MQ=255)
                                                                                                  tggctggATGAGCCGCCACTGGCGGGAACCTCACGCGCAGGGGAAGGCCCTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGGAGCCGCAGACCATTATCTGGCGGGGGATATCGAATCCCTGCCCGTAGCGACTGCGAc                                            >  2:1156854/1‑150 (MQ=255)
                                                                                                   ggctgaATGGGCCGCCACTGGCGGGAACGTCCCGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCTGCTGCGACg                                           <  2:1159647/150‑1 (MQ=255)
                                                                                                      tggAGGTGCCCCCACGGGCGGAAACGTCCCGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATCCCGCAGACCATTATCTGGCGGGAAATTTCGAATCCCTGCCGTTAGCGACTGCGACGtt                                         <  2:378439/149‑1 (MQ=255)
                                                                                                       ggATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGtt                                         <  2:1262609/148‑1 (MQ=255)
                                                                                                       ggATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCg                                       <  2:621834/150‑1 (MQ=255)
                                                                                                            ggcccccACTGGCGGAAACGTCACGCGAAGGTGCCGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAGGGATGCCGCAGACCATTAGCTGCCGGGAGATTTCGAATCCCTGCCGTTAGCGACTGCGACGTTCg                                       <  2:686698/144‑1 (MQ=255)
                                                                                                             gccgccACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTg                                 <  1:463073/150‑1 (MQ=255)
                                                                                                              ccgccACTGGCGGGGACGTCACGCGCCGGGGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACCCCAGAAGGATGCCGCAGACCATTATCTGGCGGGGGGTATCGAATCCCTGCCGTTAGCGACTGCGACGGTCGGTCtt                                  >  2:802427/1‑148 (MQ=255)
                                                                                                               cgccACTGGCGGGAACGTCACGCGCAGGTGACGGCCCTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCCGAAGGATGCCGCAGACCATTATCTGGCGGGGGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCAt                              >  1:279786/1‑151 (MQ=255)
                                                                                                                gccACTGGCGGGAACGCCACGCGCAGGGGACGGCCTTAGATTTCTCGCCGCCAATGCTTGTTTAGGCACCCCCGAAGGATGCCGCAGACCATTATCTCGCGGGGGATATCGGATCCCTGCCGTTAGCGACTGCGACGGTCGATCTTGcct                              >  2:808494/1‑148 (MQ=255)
                                                                                                                  cACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGa                           <  1:840823/151‑1 (MQ=255)
                                                                                                                  cACGGGCGGGAACGTCACGCGAAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGCCCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGa                           <  2:120994/151‑1 (MQ=255)
                                                                                                                        cGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAAtc                     >  1:562962/1‑151 (MQ=255)
                                                                                                                         gggAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAAtct                    <  2:10787/151‑1 (MQ=255)
                                                                                                                          ggAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAAtctc                   <  1:835710/151‑1 (MQ=255)
                                                                                                                          gaaaCGTCACGCGCAGGTGACGGCCTTAGATCTCCCCCCGCCAATGCTTGTTCAGGCACCCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAAtctc                   <  2:916076/149‑1 (MQ=255)
                                                                                                                            aaCGTCCCGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAAtctc                   <  2:770229/149‑1 (MQ=255)
                                                                                                                            aaCGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCgc                 >  2:297137/1‑151 (MQ=255)
                                                                                                                                gTCCCGCGCAGTTGACGGCTTTAGACCTCTCGCCGCCAAGGCTTGTGCAGGCAGGCCGGAAGGATGCCGCGGACCATTATCTGGCGGGGGATATCGAATCCCTGCCGTTAGCGGCTGCGACGTTCGATCTTGCATGGAGCAATCTCgca                <  2:532977/149‑1 (MQ=255)
                                                                                                                                  cACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCgcagtg             <  1:1218043/150‑1 (MQ=255)
                                                                                                                                   aCGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCgcagtgca           <  1:5452/151‑1 (MQ=255)
                                                                                                                                     gcgcAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCgcagtgcag          >  2:308008/1‑150 (MQ=255)
                                                                                                                                     gcgcAGGTGACGGCCTTAATTCTCTCGCCGCCAATGCTGGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCgcagtgcagt         <  2:388053/151‑1 (MQ=255)
                                                                                                                                           gTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCg   >  1:897690/1‑151 (MQ=255)
                                                                                                                                            tGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCgg  >  1:202470/1‑151 (MQ=255)
                                                                                                                                            |                                                                                                                                                      
GAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCGGGAACGTCACGCGCAGGTGACGGCCTTAGATCTCTCGCCGCCAATGCTTGTTCAGGCACGCCAGAAGGATGCCGCAGACCATTATCTGGCGGGAGATATCGAATCCCTGCCGTTAGCGACTGCGACGTTCGATCTTGCATGGAGCAATCTCGCAGTGCAGTGGTGCGG  >  NC_000913/811582‑811871

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