Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,261,0550AC44.2% 126.4 / 24.7 77V455G (GTG→GGG) tdcEpyruvate formate‑lyase 4/2‑ketobutyrate formate‑lyase
Reads supporting (aligned to +/- strand):  ref base A (38/5);  new base C (0/34);  total (38/39)
Fisher's exact test for biased strand distribution p-value = 7.98e-17
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTT  >  NC_000913/3260905‑3261204
                                                                                                                                                       |                                                                                                                                                     
aTGTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCa                                                                                                                                                       <  2:1014741/151‑1 (MQ=255)
    aGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCAcccc                                                                                                                                                   >  2:256836/1‑151 (MQ=255)
     gTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCccg                                                                                                                                                  >  2:43849/1‑151 (MQ=255)
             tatTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAAt                                                                                                                                          >  1:75677/1‑151 (MQ=255)
             tatTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAAt                                                                                                                                          >  2:301958/1‑151 (MQ=255)
              atTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAAtt                                                                                                                                         >  2:828814/1‑151 (MQ=255)
                 cAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCg                                                                                                                                      >  1:712095/1‑151 (MQ=255)
                    cgcgcTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGt                                                                                                                                     >  1:1169467/1‑149 (MQ=255)
                     gcgcTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGt                                                                                                                                     >  2:375584/1‑148 (MQ=255)
                     gcgcTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGt                                                                                                                                     >  2:757261/1‑148 (MQ=255)
                       gcttttttcctcccccGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGCCGTCCATCAGCGGGGCTGTTTTCGGCCCGACCAGAATCTTCAGCTTCTCGTCCCCCCCGCCGTTAATTGCGTaga                                                                                                                                  <  1:241289/135‑1 (MQ=255)
                         tGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGagc                                                                                                                                >  2:793387/1‑149 (MQ=255)
                            tGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGt                                                                                                                           >  1:338529/1‑151 (MQ=255)
                             gTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGtt                                                                                                                          >  2:318386/1‑151 (MQ=255)
                                cgcccccGCCACCCAGTCCATGAAGTGAGCGAGGCTGTCCATCACTTTGTCGTAACCCAGCCCGCCGCCCACCAGCGGTGCTGTTTTCGGCCCGACCTGATTCTTCAGCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGtttt                                                                                                                        <  1:243639/145‑1 (MQ=255)
                                ccccccccccagccagTCCAGGAAGTGATCGGGGCTGCCCATCTCTTTGTCGTAATCCAGCACGTCGTCCACCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGTTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGttt                                                                                                                         <  1:114529/141‑1 (MQ=255)
                                cTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTc                                                                                                                       >  2:472339/1‑151 (MQ=255)
                                     ccgcgcgccagCCCGTGAAGTAGTCGAGGCTGCCCATCTTTTTGTCGAAACCCACCCGGTCTTCCATCAGTGGTGTTGTTTTCGGCCCGCCCTGAATTTTCAGTTTCCCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGc                                                                                                                     <  2:455832/142‑1 (MQ=255)
                                     ccGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGGGGTGCTGTTTTCGGCCCGCCCTGAATCTTCAGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTc                                                                                                                       <  2:1237141/146‑1 (MQ=255)
                                     ccGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGcc                                                                                                                    >  1:918025/1‑149 (MQ=255)
                                        ccagccagTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGtt                                                                                                               >  1:400892/1‑151 (MQ=255)
                                          acccagTCCATGAAGATATGAAGGCTCTCCATCCTTTTGTATTAATCCACCGGGTCTCCAATCGGAGGTGGTGTTTTGGGCCCAACCTAATTCTTCGCCTCCGCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAg                                                                                                             <  2:616134/149‑1 (MQ=255)
                                            ccagTCCATGAAGTGATCGAGGCTGTCCATCATTTTGTCGTAACCCAGCCCGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGCCCTGAATCTTCAGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAg                                                                                                             <  1:1160977/149‑1 (MQ=255)
                                             cagccaaTGAAGTGATCGGGGCTGTCCATCACTTTGTCGTACCCCAGCCCGTCTTCCATCAGCGGTGCTGTTTTCGCCCCGCCTTGAATCTTCACCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAgcgc                                                                                                          <  1:783971/145‑1 (MQ=255)
                                               gTCCATGAAGGGATCGAGGCTGTCCATCACTTTGTCGTAACCCAGCCCGTCGTCCATCGGCGGGGCTGTTTTCGGCCCGACCTGAATCTTCAGTTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGcgcg                                                                                                         <  1:528575/150‑1 (MQ=255)
                                                    tGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCAcc                                                                                                   >  2:786698/1‑151 (MQ=255)
                                                    tGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCAcc                                                                                                   >  1:1148753/1‑151 (MQ=255)
                                                      aacggagcGAGGCTGTCCACCACTTTGTCGTAAGCCCGCTCTTCGCCCTGCGGCGGTTGTTTTTTCCCCCCACCTAAACTCTTCAACTCTCCCTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGt                                                                                                        <  2:34543/137‑1 (MQ=255)
                                                       agggttgaaGGCTGCCTATCGCTTTGTCGTAATCAACCAGGTCGTCCTCCGGCGGTGCGTTTTTCGCCCTTCCCGAACTATACTGTTGTTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCaaa                                                                                                <  2:493127/143‑1 (MQ=255)
                                                          gATCGAGGCTTTCCATCACTTTGTCTTAACCCAGCACGTCGCCCATCGGCGGTGCTGTTTTCGGCCCGACTTGAATTTTCAGCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGa                                                                                              <  2:1012120/150‑1 (MQ=255)
                                                             cGAGGCTGTCCATCATTTTGTCGTAATCCAGCCCGCCGTCCATCAGCGGTGCTGTTTTCGGCCCGCCCTGAATCTTCAGCTTCTCGCCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACt                                                                                           <  2:561098/150‑1 (MQ=255)
                                                              gAGGCTGTCCTTCATTTTGTCGTAATCCAGCTCGCCGTCCATCAGGGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCGCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCaa                                                                                                 <  1:700298/143‑1 (MQ=255)
                                                                ggttttCCATCACTTTGTCGTAACCCAGCACGTCGCCCATCGGCGGTGCTTTTTCCGGCCCGCCCTGAATCTTCAGCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCa                                                                                        <  1:769046/145‑1 (MQ=255)
                                                                ggCTTTCAATCTTTTTGTCGTAACCCAGCACGTCGCCCATCAGCGGTGCTGTTTTCGGCCCGCCCTGAATCTCCAGTTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCa                                                                                        <  2:1232693/150‑1 (MQ=255)
                                                                  cTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCAttt                                                                                     >  1:235017/1‑151 (MQ=255)
                                                                      ccATCTTTTTGTCGTAATCCAGCACGTCGCCCACCAGCGGTGTTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCTCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGc                                                                                   <  1:831602/149‑1 (MQ=255)
                                                                      ccATCACTTTGTCGTAATCCAGCACGTCGTCCACCAGCGGTGCTTTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCt                                                                                  <  2:430370/150‑1 (MQ=255)
                                                                      ccATCACTTTGTCGTAATCCAGCACGTCGCCCATCAGCGGGGTTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCAttt                                                                                     <  1:163332/147‑1 (MQ=255)
                                                                          cgctttttagaaaTCCAGCACGCCGCCCAGCAGCGGTGCTGTTTTCGGCCGTACCTGACCTCACGGCTTCGCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTAc                                                                               <  2:945034/138‑1 (MQ=255)
                                                                          cattttGTCAAAACCCAGCCCGGCGTCCACCGGGGGGGCTGTGTTCGGCCCGCCCTAAATCTTCGGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTa                                                                                <  1:301740/145‑1 (MQ=255)
                                                                             tttGTCGAAACCCAGCACGCCCTCCACCAGCGGTGCTTTTTGCGGCCCGCCCTGAATCTTCACCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCa                                                                             <  1:829058/148‑1 (MQ=255)
                                                                                  cGTAACCCAGCACGTCGCCCACCAGCGGTGCTTTTTTCGGCCCGACCGAAATTTTCAGTTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCAcc                                                                       <  2:438035/149‑1 (MQ=255)
                                                                                     accccaccgcgtcgtCCTGCGGGGGTGTTTTTTTCGCCCCGCCCTGAATCTTCAGTTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCa                                                                         <  1:476796/135‑1 (MQ=255)
                                                                                     aacccAGCCCGTCGTCCATCAGGGGTGCTGTTTCCGCCCCGCCCTAAATCTTCAGCTTCCCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCAc                                                                        <  2:712097/142‑1 (MQ=255)
                                                                                      aTCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGc                                                                >  2:745351/1‑151 (MQ=255)
                                                                                        ccAGCACGTCGCCCATCGGCGGTGTTGTTTTCGGCCCGACCTGAATCTCCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGc                                                                <  1:244941/150‑1 (MQ=255)
                                                                                         cAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGCCCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTg                                                              <  2:717521/151‑1 (MQ=255)
                                                                                                  gtcCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACgcagc                                                       >  2:169939/1‑149 (MQ=255)
                                                                                                    ccATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGc                                                   >  1:1130731/1‑151 (MQ=255)
                                                                                                    ccATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGc                                                   >  2:1252562/1‑151 (MQ=255)
                                                                                                      aTCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCa                                                  >  2:972740/1‑150 (MQ=255)
                                                                                                        cggCGGTGTTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAAt                                                <  2:1061480/148‑1 (MQ=255)
                                                                                                        cAGGGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATc                                               <  2:824989/151‑1 (MQ=255)
                                                                                                         aGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATcg                                              >  2:114276/1‑151 (MQ=255)
                                                                                                          gCGGTGTTGTTTTCGGCCCGCCCTGAATCTTCAGCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATcg                                              <  2:828299/150‑1 (MQ=255)
                                                                                                          gCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATcgc                                             >  2:464002/1‑151 (MQ=255)
                                                                                                                  gTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAAtcgt                                     >  2:717246/1‑151 (MQ=255)
                                                                                                                   ttttCGGCCCGCCCTGACTCTTCAGCTTCCCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAAtcg                                      <  2:891397/149‑1 (MQ=255)
                                                                                                                    tttCGCCCGCACCTAAATCTCCAGCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCCCCATTGGGCTGACGCAGCAGGCAATCGCGTAAtcgtc                                    <  2:8188/150‑1 (MQ=255)
                                                                                                                        ggCCCGACCTGAATCTTCAGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGc                                  <  2:872438/148‑1 (MQ=255)
                                                                                                                        ggCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGt                               >  2:1002007/1‑151 (MQ=255)
                                                                                                                          cccGCCCTGAATCTTCAGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGt                               <  2:387262/149‑1 (MQ=255)
                                                                                                                           ccGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTg                             >  2:542472/1‑150 (MQ=255)
                                                                                                                              aCCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGt                         >  1:2281/1‑151 (MQ=255)
                                                                                                                                cTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTc                        >  2:1249510/1‑150 (MQ=255)
                                                                                                                                    aTCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTAc                   >  2:580387/1‑151 (MQ=255)
                                                                                                                                      tCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACg                  >  1:1178410/1‑151 (MQ=255)
                                                                                                                                       cTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGc                 >  1:1056453/1‑151 (MQ=255)
                                                                                                                                        ttCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGc                 >  2:71964/1‑150 (MQ=255)
                                                                                                                                         tCAGCTTCCCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACg                  <  2:1108781/148‑1 (MQ=255)
                                                                                                                                          caccTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACg                  <  2:739165/144‑1 (MQ=255)
                                                                                                                                          cAGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATc              <  2:1169477/151‑1 (MQ=255)
                                                                                                                                          cAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATc              >  2:470273/1‑151 (MQ=255)
                                                                                                                                           aGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCAt               <  2:1577/149‑1 (MQ=255)
                                                                                                                                           aGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATc              >  2:810517/1‑150 (MQ=255)
                                                                                                                                            ctttCGCGCCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGt                     <  2:1105122/140‑1 (MQ=255)
                                                                                                                                               tctcGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATc              >  1:121943/1‑146 (MQ=255)
                                                                                                                                                      cACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCAttt  >  1:639730/1‑151 (MQ=255)
                                                                                                                                                       |                                                                                                                                                     
ATGTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTT  >  NC_000913/3260905‑3261204

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