Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,261,0550AC40.0% 102.0 / 10.0 62V455G (GTG→GGG) tdcEpyruvate formate‑lyase 4/2‑ketobutyrate formate‑lyase
Reads supporting (aligned to +/- strand):  ref base A (29/7);  new base C (0/24);  total (29/33)
Fisher's exact test for biased strand distribution p-value = 7.62e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTT  >  NC_000913/3260907‑3261205
                                                                                                                                                    |                                                                                                                                                      
gTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCAc                                                                                                                                                       <  2:163361/150‑1 (MQ=255)
   gTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGACAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTCCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCAACccg                                                                                                                                                   <  1:121753/151‑1 (MQ=255)
   gTGAATGATATTCAGCGCGCTGATGTACTACACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCccg                                                                                                                                                   <  2:54919/151‑1 (MQ=255)
    tGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCccgc                                                                                                                                                  >  2:242725/1‑151 (MQ=255)
      aaTGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCccgccg                                                                                                                                                >  2:804462/1‑151 (MQ=255)
       aTGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGt                                                                                                                                               >  2:950673/1‑151 (MQ=255)
          atatTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTa                                                                                                                                             >  2:822335/1‑150 (MQ=255)
          atatTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTa                                                                                                                                             >  2:816947/1‑150 (MQ=255)
           tatTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAAt                                                                                                                                           >  2:975330/1‑151 (MQ=255)
            atTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAAtt                                                                                                                                          >  1:299148/1‑151 (MQ=255)
               cAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCg                                                                                                                                       >  1:839802/1‑151 (MQ=255)
                         ttttcctccACCCCCAGCCAGCCCATGAAGTGACGGAGCCTGTCCATCACTTTGTCGTAACCCACCCCGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGCCCTGAATCTTCAGTTTCTCGACCCCCCCGCCGTTAATTGc                                                                                                                                        <  1:914112/132‑1 (MQ=255)
                           gTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGtt                                                                                                                           >  2:144841/1‑151 (MQ=255)
                             aCTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGtttt                                                                                                                         >  1:614061/1‑151 (MQ=255)
                              cgccccccccacccagTCCATGAGGTGAGCGGGGTTGTCCATCTTTTGTTCAAAATCCAGCTCGTCGTCCACCGGGGGTGCTGTTTTCGGCCCGACCTGATTCTTCAGCTTCTCCCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGttt                                                                                                                          <  1:117310/141‑1 (MQ=255)
                                      cgagcttccccATGAGGTGATCGTGGCTGGCCTTCTTTTTGCCGTATTCCACCCCCTCGCCCATAAGGGGTGCTGTTTTCGCCCCGCCCTGATTCTTCACCTTCGCTTCCACGCCGCCGTTAATTGCGTAGagcagc                                                                                                                              <  1:1283972/128‑1 (MQ=255)
                                      ccagccagTCCAGAAGGTTATCGGGGCTCTCCACCTCTTTTACTTAACCCTCGCCGTGGTCCACAAGGGGGGCTGTTTCGGGCCCTACTTGAGTCTTCATTTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGtt                                                                                                                           <  1:753129/140‑1 (MQ=255)
                                       cagccagTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGtt                                                                                                                >  1:434108/1‑150 (MQ=255)
                                          ccagTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCCGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAgcg                                                                                                            >  1:917589/1‑151 (MQ=255)
                                          ccagTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAgcg                                                                                                            >  2:827767/1‑151 (MQ=255)
                                          ccagTCCATGAAGTGACGGGGGCTGCCCTTCCTTTTGTCGAAACCCAGCACGTCCTCCACCAGGGGTGCTTTTTCCGGCCCGCCCGGATTCTTCAGCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAgcg                                                                                                            <  1:1251274/151‑1 (MQ=255)
                                               caaggAAGTGATCGAGGCTGTCCATCACTTTGTCGTAACCCAGCACGTCGACCATCAGCGGTGCTGTTTTCGCCCCGACGTAAATCTTCAGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGc                                                                                                       <  1:639785/147‑1 (MQ=255)
                                                cATGAAGTGTACGAGCCTCTCCATCGCTTTTCCGAAATCCACCGCGTCGGCCGTGGGGGGGGGTGTTTTGGGCCCGCCCAGACCCTTCAGCTTCTCGACCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCa                                                                                                      <  1:602916/151‑1 (MQ=255)
                                                       tGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCaa                                                                                                  >  2:364036/1‑148 (MQ=255)
                                                          tCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACt                                                                                            >  2:678403/1‑151 (MQ=255)
                                                           caaGGCTGTCCATCATTTTCCCTTAATCCTCCACATCGCCTATCGCGGGCGTTGTTTTGGGCCCTTCCTGAATGTTGCGCTTCTAGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGa                                                                                               <  2:266624/145‑1 (MQ=255)
                                                            gAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGc                                                                                          >  2:668536/1‑151 (MQ=255)
                                                              ggCTGTCCATCACTTTGTCGTAATCCAGCCCGTCGTCCATCAGCGGTGTTTTTTTCGGCCCGCCCTGAATCTTCAGTTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCa                                                                                         <  2:1081538/150‑1 (MQ=255)
                                                                        ctttttGTCGTAACCCAGCACGTCTTCCATCAGCGGTGTTTTTTTCGGCCCGCCCTGAATCTTCAGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTa                                                                                 <  1:332391/145‑1 (MQ=255)
                                                                         aCTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCa                                                                              >  2:91230/1‑150 (MQ=255)
                                                                            ttctCGTAACCCAGCCCTTCGTCCATCAGCGGTGCTTTTTCCGCCCCGCCCTGATTCTTCATCTCCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCa                                                                              <  1:135328/144‑1 (MQ=255)
                                                                             tGTCGTAACCCAGCACGTCGTCCATCGGCGGTGCTTTTTTCGGCCCGCCTTAAATCTTCATCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCa                                                                          <  2:755816/150‑1 (MQ=255)
                                                                                cgaaacccAGCCCGTCGTCCATCAGCGGTGCTGTTTTCGCCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCAt                                                                      <  1:731951/145‑1 (MQ=255)
                                                                                 gtacccccgCATGTCGCCCATGAGCGGTGTTTTTTTCGGCCCGCCGTTAATCTTCTACTCCTCGCCCGCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCAcc                                                                        <  1:929944/140‑1 (MQ=255)
                                                                                   aaTCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTgg                                                                   >  2:1016718/1‑151 (MQ=255)
                                                                                     tCCAGCAAGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGAACTGAATCTTCAGCTTCTCGTCCACCCCGACGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGATAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGc                                                                 >  2:879370/1‑151 (MQ=255)
                                                                                         gcccgtcctCCACGAGCGGTGTTGTTTTCGACCGCCCCTAAATCTTCATCTCCTCGCCCGCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGc                                                                 <  1:1240546/144‑1 (MQ=255)
                                                                                          cACGTCGTCCATCAGCGGTGCTGTTTGCGGCCCGACCTGAATCTTCAGCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTg                                                               <  1:975273/148‑1 (MQ=255)
                                                                                                  ccATCAGCGTTGCTGTTTTCGGCCCGCCCTGATTCTTCAGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACgcag                                                         <  2:515984/146‑1 (MQ=255)
                                                                                                   caccgggggTGCTGTTTCCGCCCCGACCTAAATCTTCGGTTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCa                                                   <  2:1054806/144‑1 (MQ=255)
                                                                                                          ggTGTTGTTTTCGGCCCGCCCTGAATCTTCAGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATc                                                <  1:337274/147‑1 (MQ=255)
                                                                                                           gTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGt                                            >  2:72349/1‑150 (MQ=255)
                                                                                                             gtttttttCGGCCCGACCTGAATCTTCAGCTCCTCCCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTaa                                          <  1:22496/146‑1 (MQ=255)
                                                                                                                 ttttCGGCCCGCCCTCATTCTTCACCTCCCCTCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAAt                                         <  1:1003059/147‑1 (MQ=255)
                                                                                                                 ttttCGGCCCGACCTGATTCTTCAGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAAtc                                        <  2:428097/148‑1 (MQ=255)
                                                                                                                 tttgCGGCCCGACTGGAATCTTCAGCTTCTCCCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAAtcg                                       <  1:668512/149‑1 (MQ=255)
                                                                                                                  tttCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAAtcgtc                                     >  1:767405/1‑150 (MQ=255)
                                                                                                                  ttccgcccggcccacaaTCTTCCTCTTCCCCCCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAAtcgtc                                     <  1:234916/135‑1 (MQ=255)
                                                                                                                     cGGCCGGGCTTGATTCTCAAGCTTCCGCTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAAtcgtcg                                    <  2:454051/148‑1 (MQ=255)
                                                                                                                     cGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTg                                 >  2:693092/1‑151 (MQ=255)
                                                                                                                      ggCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGt                                >  1:1270035/1‑151 (MQ=255)
                                                                                                                        cccgccgcgAATCTCAACCTTCCCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCt                                  <  1:168743/139‑1 (MQ=255)
                                                                                                                        cccGTCCTGAATCTTCATCTTCCCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGtt                               <  1:754448/150‑1 (MQ=255)
                                                                                                                             ccTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTc                         >  2:365366/1‑151 (MQ=255)
                                                                                                                             ccTGAATCTTCAGCTCCTCCTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAg                           <  1:550976/149‑1 (MQ=255)
                                                                                                                               taaaTCTTCAGCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTc                         <  1:452389/147‑1 (MQ=255)
                                                                                                                                  aTCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTAc                    >  2:1233885/1‑151 (MQ=255)
                                                                                                                                  aTCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTAc                    >  1:348305/1‑151 (MQ=255)
                                                                                                                                          cTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATc               >  2:464774/1‑148 (MQ=255)
                                                                                                                                               cGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGt        >  1:919944/1‑150 (MQ=255)
                                                                                                                                                gTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTc       >  2:984064/1‑150 (MQ=255)
                                                                                                                                                   cACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCAttt   <  1:78787/151‑1 (MQ=255)
                                                                                                                                                    aCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCAtttt  >  1:516442/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
GTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTT  >  NC_000913/3260907‑3261205

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