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

GTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTT  >  NC_000913/3260907‑3261205
                                                                                                                                                    |                                                                                                                                                      
gTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCAc                                                                                                                                                       <  1:907123/150‑1 (MQ=255)
  aGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCAcccc                                                                                                                                                    >  1:681225/1‑151 (MQ=255)
          atatTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTa                                                                                                                                             >  2:639489/1‑150 (MQ=255)
          atatTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTa                                                                                                                                             >  1:587612/1‑150 (MQ=255)
           tttccggcgcgtTGTTTTACTCCCCCGCCAGCCAGTCCATGAAGTTATCGGGTCTCCCCATCTTTTTTTCTAAATCCAGCCCTTCGCCCACCAGCGGTGCTTTTTTCGGCCCGACCTGAATCTCCAGCTTCGCGCCCCCCCCGCCGtt                                                                                                                                              <  2:236739/142‑1 (MQ=255)
            atTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAAtt                                                                                                                                          >  1:895090/1‑151 (MQ=255)
               cAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCg                                                                                                                                       >  1:694064/1‑151 (MQ=255)
                   gcgcTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGt                                                                                                                                      >  2:240809/1‑148 (MQ=255)
                            tACGGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCGGTCCATCCCTTTGTCGTAATCCAGCCCGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGCCCTGAATCTTCAGCTTCTCGCCCACCCCGCCGTTAATTGCGTAGAGCAGCGt                                                                                                                            <  1:212354/149‑1 (MQ=39)
                                 ccccGCCAGCAAGCCCATGAAGTGATCGGGGCTGTCCATCACTTTGTCTTAATCCAGCCCGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTAAATCTTCAGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCg                                                                                                                       <  1:97473/147‑1 (MQ=255)
                                   ccGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCCCTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGttt                                                                                                                          <  1:788764/144‑1 (MQ=255)
                                    ccccagccagTCCAGGAAGGGATCGAGGCTGCCCATCCCTTTGTCGTAATCCAGCCCGTCGTCCACCGGCGGTGCTGTTTTCGGCCCGACTTGAATCTTCATCTTCTCTCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTc                                                                                                                        <  2:569870/143‑1 (MQ=255)
                                    cGCCAGCCAGTCCATGAAGTGATCGGGGCTGTCCATCTCTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTTTTTTCGGCCCAACCTGATTCTTCAGCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAgg                                                                                                                  <  1:1697/151‑1 (MQ=38)
                                      ccagccagTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGtt                                                                                                                >  1:206236/1‑151 (MQ=255)
                                       cagccagcCCAGAAAGTAATCGAGGCTGTCCATCTCTTTGTCTAAACCCACCCCTTCGCCCACCGGCGGTGCTGTTTTCGCCCCGCCCTGAATCTTCGCCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTa                                                                                                               <  2:369913/151‑1 (MQ=255)
                                        agccagTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGtt                                                                                                                >  2:684874/1‑149 (MQ=255)
                                             gTCCATGAAGTGATGGAGGCTGTCCATCTCTTTGTCGTAATCCAGCCCGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGt                                                                                                         <  1:129204/151‑1 (MQ=255)
                                               ccATGAAGTGATCGAGGCTGTCCATCCTTTTGTCGAAACCCAGCCCGCCGTCCACCAGGGGGGCTTTTTTGGGCCCGCCCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCGGCGTTTTCGCCAGGTTAg                                                                                                              <  2:129747/144‑1 (MQ=255)
                                                cATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGc                                                                                                       >  2:23018/1‑150 (MQ=255)
                                                     aGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCaa                                                                                                  >  2:700629/1‑150 (MQ=255)
                                                        gaccGAGGCTGTCCATCACTTTGTCGTAATCCAGCCCGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGaa                                                                                              <  1:337383/148‑1 (MQ=255)
                                                        gATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGaa                                                                                              <  1:658669/151‑1 (MQ=255)
                                                        gATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGa                                                                                               >  1:239024/1‑150 (MQ=255)
                                                                   cccATCTCTTTTTCGAACCCCAGCCCGTCTTCCACCAGGGGGGCGTTTTTGGGCCCGCCCTAATTCTTCGGCTCCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCAtt                                                                                       <  2:216117/146‑1 (MQ=255)
                                                                    ccATCATTTTGTCGTAACCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTg                                                                                     <  1:311197/148‑1 (MQ=255)
                                                                    ccATCACTTTGTCGTAATCCAGCCCGCCGTCCATCAGCGGGGCTGTTTTCGGCCCCACCTGATTCTTCACCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTg                                                                                     <  2:1037976/148‑1 (MQ=255)
                                                                     cATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCtt                                                                                  >  1:277548/1‑150 (MQ=255)
                                                                        cccTTTGTCGTAACCCAGCACGTCGTCCACCAGCGGGGCTTTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTAc                                                                                <  2:515814/147‑1 (MQ=255)
                                                                           tttGTCGTAATCCAGCACGTCTTTCATCAGCGGTGCTGGTTTCAGGCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCTCGTGCACCAAAGAACTGCATTTGCTTAcc                                                                               >  2:151056/1‑147 (MQ=255)
                                                                           tttGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATc                                                                           >  2:222864/1‑151 (MQ=255)
                                                                            ttGTCGTAATCCCGCACGTCGTCCATCAGGGGTGCTTTTTTCGCCCCGACCTGAATCTTCAGCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATc                                                                           <  1:57726/150‑1 (MQ=255)
                                                                             tGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCAc                                                                         >  1:191828/1‑151 (MQ=255)
                                                                                 gTAATCCAGCACGTCGTCCATCAGCGGTGCTTTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCAt                                                                      <  1:22834/150‑1 (MQ=255)
                                                                                  tAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTg                                                                    >  1:693546/1‑151 (MQ=255)
                                                                                   aaTCCAGCACGTCGTCCATCGGCGGTGCTGTTTTCGCCCCGCCCTGAATCTTCAGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCAtt                                                                     <  2:239022/149‑1 (MQ=255)
                                                                                   aaTCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTgg                                                                   >  2:980404/1‑151 (MQ=255)
                                                                                      ccAGCCCGTCTTCCATCAGCGGTGCTGTTTGCGCCCCGCCCTGACTTTTCAGTTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGc                                                                 <  1:172161/150‑1 (MQ=255)
                                                                                      ccAGCCCGTCGTCCATCGGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTggg                                                                  <  1:516093/149‑1 (MQ=255)
                                                                                       cAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCt                                                                >  1:726554/1‑150 (MQ=255)
                                                                                        aGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACATGAATCTTCAGCTTCTCGTCCACCCCGCCGGTAATTGCGTAGAGTAGCGTTTTTGCCAGGTTAGCGCGTGCACCAAAGCACTGCATTTGCTTACCAATCACCATTg                                                                    >  2:114894/1‑145 (MQ=255)
                                                                                           acctcccccATCGGGGGTGTTTTTTCCGGCCGCGCCTAAGTCTTCAGCTTCTCTCCCTCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGAc                                                             <  1:910504/142‑1 (MQ=255)
                                                                                            cgtccaccATCAGCGGTATTGTTTCCCGCCGGACCGGAATCTTCAGCTTCGCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACgca                                                          <  2:236188/145‑1 (MQ=255)
                                                                                                tccctccggCGGTGTTGTTTTCGGCCCGACCTAAATCTTCACCTTCTCCTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGAc                                                             <  2:745418/136‑1 (MQ=255)
                                                                                                  ccATCAGCGGTGTTGTTTTCGGCCCGCCCTGAATCTTCAGCTTCTCGCCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACgcagca                                                       <  1:803267/148‑1 (MQ=255)
                                                                                                  ccATCAGCGGTGTTGTTTTCGGCCCGACCTGATTCTTCAGCTCCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACgcag                                                         <  2:895081/146‑1 (MQ=255)
                                                                                                     tCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAAt                                                 >  1:283542/1‑151 (MQ=255)
                                                                                                     tCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAAt                                                 >  1:612255/1‑151 (MQ=255)
                                                                                                      cAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATc                                                >  1:458372/1‑151 (MQ=255)
                                                                                                         cgttgtttttttCGCCCCCCCCTGAATTTTCAGCTCCCCCACCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATc                                                <  1:114489/140‑1 (MQ=255)
                                                                                                          ggTGCTGTTTCCGGCCCGACCTGAATCTCCACCTTCTCGTCCACCCCCCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGATTAGCGCGTGCACCAAAGAACTGCATCTGCTTACCACCCACCATTGGGCTGACGCAGCAGGCACTCGCGt                                            >  2:1017185/1‑151 (MQ=255)
                                                                                                           gTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATcgcg                                             <  2:747775/149‑1 (MQ=255)
                                                                                                               tGTTTCCGCCGCGCCATGAATCTTCAGCTTCTTGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATcgc                                              <  2:231033/144‑1 (MQ=255)
                                                                                                                 ttttCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCCATCGCGTAAtcgtc                                     >  2:760915/1‑151 (MQ=255)
                                                                                                                 ttttCGCCCCGCCCTAAATCTTCAGTTTCCCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAAtcgt                                      <  1:451953/150‑1 (MQ=255)
                                                                                                                  tttCGGCCCGACCTGAATCTTCAGCTTCTCGTGCACCCCGCCGTTAATTGCGCAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGAAATCGCGTAAtcgtcg                                    >  2:372760/1‑151 (MQ=255)
                                                                                                                  tttCGGCCCGACCTAAATCTTCAGCTCCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAAtcgtc                                     <  2:553807/150‑1 (MQ=255)
                                                                                                                   ttCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGc                                   >  2:1042913/1‑151 (MQ=255)
                                                                                                                     cGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTg                                 >  2:235015/1‑151 (MQ=255)
                                                                                                                      ggCCCGCCCTGAATCTTCAGCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCt                                  <  1:745901/149‑1 (MQ=255)
                                                                                                                      gccccGCCCTGAATTTCCAGCTTCTCGACCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTg                                 <  1:664806/148‑1 (MQ=255)
                                                                                                                         ccGACCTGAATCTTCAGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGa                             <  1:971865/151‑1 (MQ=255)
                                                                                                                          cGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGa                             >  2:789501/1‑150 (MQ=255)
                                                                                                                                gAATCTTCAGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCa                        <  2:985249/149‑1 (MQ=255)
                                                                                                                                       cAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATc               >  2:712019/1‑151 (MQ=255)
                                                                                                                                           ttCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATc               >  1:536810/1‑147 (MQ=255)
                                                                                                                                           ttCTCCTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCa              <  2:840210/148‑1 (MQ=255)
                                                                                                                                           tcctcGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACg                   <  2:1055791/141‑1 (MQ=255)
                                                                                                                                             cccGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCGCCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCa              <  1:961738/144‑1 (MQ=255)
                                                                                                                                                gTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTc       >  1:333728/1‑150 (MQ=255)
                                                                                                                                                  ccACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCAtt    >  1:684856/1‑151 (MQ=255)
                                                                                                                                                  ccACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCAtt    >  1:15797/1‑151 (MQ=255)
                                                                                                                                                    aCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCAtttt  >  2:209358/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 < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: