Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,261,0550AC53.9% 94.5 / 121.6 76V455G (GTG→GGG) tdcEpyruvate formate‑lyase 4/2‑ketobutyrate formate‑lyase
Reads supporting (aligned to +/- strand):  ref base A (30/5);  new base C (0/41);  total (30/46)
Fisher's exact test for biased strand distribution p-value = 2.51e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.76e-14

GTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTT  >  NC_000913/3260910‑3261205
                                                                                                                                                  |                                                                                                                                                      
gTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCGCTTTGGCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCAcccc                                                                                                                                                    <  1:862223/150‑1 (MQ=255)
gTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCccg                                                                                                                                                   <  2:1046305/151‑1 (MQ=255)
 tGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCATCCACCccgc                                                                                                                                                  >  2:1361904/1‑151 (MQ=255)
  gAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCccgcc                                                                                                                                                 <  2:1461150/151‑1 (MQ=255)
     tGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGt                                                                                                                                               >  2:1005495/1‑150 (MQ=255)
           tCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGc                                                                                                                                        >  1:58715/1‑151 (MQ=255)
           tCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGc                                                                                                                                        >  2:31866/1‑151 (MQ=255)
           tCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGc                                                                                                                                        >  1:1431804/1‑151 (MQ=255)
                gcgcTGTTTTCCGGCCCCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGCCCCCCCCGCCGTTAATTg                                                                                                                                         <  2:142093/145‑1 (MQ=255)
                gcgcTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGt                                                                                                                                      >  2:1529700/1‑148 (MQ=255)
                gcgcTGATGTACTGCAACGCCAGCCAGTTCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGt                                                                                                                                      >  1:17887/1‑148 (MQ=255)
                     gATGTCCTGCCCCGCCGGCCATTCAATGAAGTGATCGAGGCTGTCCATCCCTTTGTCTAAATCCACCACGTCCCCCACCGGCGGTGCTGTTTGCGGCCCGCCTGGAATTTTCATTTCGTCGCCACCCCCGCCGTTAATTGCGTAGagca                                                                                                                                <  1:1487104/149‑1 (MQ=255)
                      aTGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCg                                                                                                                             >  1:1161764/1‑151 (MQ=255)
                       tttccTGCCCCCCCACCCCGTCCATGAAGTGAGCGGGGCTGTCCATCTCTTTGTCGTAACCCAGCACTTCTTCCACCGGCGGTGCTGTTTTCGGCCCGCCCGAAATCTTCAGCTTCGCGCCCCCCCCGCCGTTAAtt                                                                                                                                          <  1:1160400/133‑1 (MQ=255)
                       tGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCg                                                                                                                             >  1:170134/1‑150 (MQ=255)
                              cccccccagccagCCCATGAAGTGACCGAGGCTGTCCATCCCTTTGTCGTAATCCACCACGTCGCCCATCAGCGGTGCTGTTTTCGCCCCGCCCTGAATCTTCACCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGtt                                                                                                                           <  1:616556/140‑1 (MQ=255)
                              ccccccaaccctgTCCATAAGGTGACCGAGGCTGTCCACAATTTTGTCTAAATCCACCACATCGCCCATCAGCGGTGCTGTTTTCGGCCGGCCTTAAATTTTCAGTTGCGCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGtttt                                                                                                                         <  2:1219539/135‑1 (MQ=255)
                              ccccGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCTCTTTGTCGTAACCCAGCACGTCGTCCATCAGCGGTGCTTTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTc                                                                                                                        <  2:330172/146‑1 (MQ=255)
                                 cGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAgg                                                                                                                  >  2:241178/1‑151 (MQ=255)
                                   ccagccagTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGtt                                                                                                                >  2:1028440/1‑151 (MQ=255)
                                       ccagTGCATAAAGTGCTCGGGGCTGTCCATCACTTTGTCTTAATCCCTCACGCCGTCCATCAGGGGTGTTGTTTTCGCCCCGCCTATATTCTTCAGCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAgc                                                                                                             <  1:1363514/150‑1 (MQ=255)
                                       ccagTCCATGAATTGATCGGGGCTGTCCATCCCTTTGTCGTACTCCCCCCCGTCGCCCATCGGCGGTCCTGTTTCCGGCCCGCCCGGAATCTTCAGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAg                                                                                                              <  1:348898/149‑1 (MQ=255)
                                       ccagTCCATGAAGTGATCGAGGCTGTCCATCTCTTTGTCGTAATCCAGCACGTCGTCCATCGGCGGTGCTTTTTTCGGCCCGACCTAAATCTTCAGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAgc                                                                                                             <  2:1192305/150‑1 (MQ=255)
                                       ccagCCCAGGAAGGGATCGAGGCTCTCCATCCTTTTGTCGTAACCCCCCCCGTCGCCAATCAGCGGTCCTTTTTTCGGCCCGCCCTAAATCTTCAGCTTCTCTTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTa                                                                                                               <  1:1414725/148‑1 (MQ=255)
                                         agTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGcgcg                                                                                                          >  1:575435/1‑151 (MQ=255)
                                           tCCATGAAGTGATCGAGGCTGTCCATCCTTTTGTCGTAACCCAGCCCGTCGTCCATCAGCGGTGTTGTTTGCGGCCCGACCTAAATCTTCAGCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGt                                                                                                         <  1:1352090/150‑1 (MQ=255)
                                              aTGAAGGGATCGGGGCTGTCCATCTTTTTGTCGAAATCCACCACGTCGCCCACCGGCGGTGTTTTTTCCGCCCCGCCTTAAATCTTCAGCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGt                                                                                                         <  1:1495078/147‑1 (MQ=255)
                                               tGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCAcc                                                                                                    >  2:1347631/1‑151 (MQ=255)
                                                gAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCAcc                                                                                                    >  2:1047514/1‑150 (MQ=255)
                                                 aggttaTCGAGGTTGTCCATCTCTTTGTCGTAATCCACCACGTCGCCCATCAGGGTTGCTTTTTGCGGCCCGCCCTGACTCTTCAGTTTCGGGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCa                                                                                                   <  1:1291040/145‑1 (MQ=255)
                                                 aaGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCa                                                                                                   >  2:965371/1‑150 (MQ=255)
                                                 aaGTGATCGAGGCTGCCCATCCCTTTGTCGTAACCCAGCCCTTCTTCCTTCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCa                                                                                                   <  2:1161828/150‑1 (MQ=255)
                                                 aaGGGACCGGGGCTGTCCATCGTTTTGACGTAACCCAGCCCGTCGTCCATCAGCGGTGCTTTTTCCGGCCCGACCGGAATCTTAGTTTTCTCCTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCAc                                                                                                     <  2:1359327/148‑1 (MQ=255)
                                                   gTGATCGGGGCTGTCCATCACTTTGTCGTAACCCAGCACGTCGTCCATCGGCGGTGTTGTTTTCGGCCCGACCTGAATCTTCACCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAg                                                                                                <  2:1426653/151‑1 (MQ=255)
                                                       tCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACt                                                                                            >  1:599877/1‑151 (MQ=255)
                                                        cGAGGCTGTCCATCATTTTGTCGTAACCCAGCCCGTCTTCCATCAGCGGTGTTTTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACt                                                                                            <  2:147934/150‑1 (MQ=255)
                                                                tctatctttttGTCGTAATCCCGCACGTCGTCCATCGGCGTTGTTTTTTTCGCCCCGTCCTGAATCTTCAGTTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTg                                                                                     <  2:815902/141‑1 (MQ=255)
                                                                 ccATCATTTTGTCGTAACCCAGCACGTCGTCCATCAGCGGTGTTTTTTTCGGCCCGCCCTGAATCTTCAGCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCtt                                                                                  <  1:1383313/151‑1 (MQ=255)
                                                                   aTCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTAc                                                                                >  2:283820/1‑151 (MQ=255)
                                                                      aCTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCa                                                                              >  2:1336211/1‑150 (MQ=255)
                                                                        tttGTCGTAACCCAGCCCGTCGCCCATCGGCGGTGCTTTTTTCGGCCCGCCCTGAATCTTCAGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAAt                                                                            <  2:115150/150‑1 (MQ=255)
                                                                        ttggTCGTAATCCCGCACGTCGTCCCTCAGCGGTGCTGTTTCCGGCCCTACCTGAATCTTCAGCCTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAAt                                                                            <  2:572395/147‑1 (MQ=255)
                                                                          tGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCAc                                                                         >  1:4/1‑151 (MQ=255)
                                                                            tCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCAcc                                                                        >  2:1378291/1‑150 (MQ=255)
                                                                             cGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCAt                                                                      >  1:836120/1‑151 (MQ=255)
                                                                             cGTAACCCGGCACGTCGCCCATCAGCGGTGTTTTTTGCGGCCCGCCCTAAATCTTCATCTTCCCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCa                                                                       <  1:923525/150‑1 (MQ=255)
                                                                                aaTCCAGCCCGTCGCCCATCAGCGGTGCTGTTTTCGGCCCGCCCTGAATCTTCAGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCAt                                                                      <  2:1398419/148‑1 (MQ=255)
                                                                                   ccAGCCCGTCGCCCATCAGGGGTGTTGTTTTCGGCCCGCCCTGAATCTTCAGCTTCTCTTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTggg                                                                  <  1:956278/149‑1 (MQ=255)
                                                                                   ccAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCt                                                                >  2:1489718/1‑151 (MQ=255)
                                                                                    cAGCACGTCGTCCATCCGCTGTTGTGTTTTCGCCCCGCCTTGATTCTCCGGCTCCCCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCGCCATTGGGc                                                                 <  1:2886/149‑1 (MQ=255)
                                                                                       cccgtcgtCCATCGGCGGTGTTTTTTCCGCCCCGACCGGAATTTTCACCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCAta                                                                     <  1:908556/140‑2 (MQ=255)
                                                                                           tcgccCATCAGCGGTGCTGTTTTCGGCCCGCCCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACgca                                                          <  2:323120/149‑1 (MQ=255)
                                                                                                 accAGGGGTGTTGTTTCCGCCCCGACAAAATCCTCAGCTTTCCCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACgca                                                          <  1:835078/141‑1 (MQ=255)
                                                                                                    aGCGGTGCTGTTTTCGGCCCGCCCTAAATCTTCAGCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATc                                                <  2:819601/150‑1 (MQ=255)
                                                                                                     gggggttttttttCGCCCCGCCATGAATCTTCAGCTCCTCGCCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATcg                                               <  1:595320/141‑1 (MQ=255)
                                                                                                     gCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATc                                                <  1:1410243/149‑1 (MQ=255)
                                                                                                       ggTGTTGTTTCCGGCCCGACCTAAATCTTCACCTTCTCCCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGt                                            <  2:17890/151‑1 (MQ=255)
                                                                                                       ggTGCTTTTTTCGGCCCGCCCTAAATCTTCAGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAgg                                                     <  1:1191190/142‑1 (MQ=255)
                                                                                                       ggTGCTTTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCa                                                   <  1:1550352/144‑1 (MQ=255)
                                                                                                         tGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTa                                           >  1:641849/1‑150 (MQ=255)
                                                                                                               ttgcGGCCCGCCCTGAATCTTCAGCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAAtcgtc                                     <  1:1378868/147‑1 (MQ=255)
                                                                                                                ttCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGc                                   >  2:1269567/1‑151 (MQ=255)
                                                                                                                  cGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTg                                 >  1:1113936/1‑151 (MQ=255)
                                                                                                                     cccGCCCTAAATCTTCAGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTg                              <  1:1041037/151‑1 (MQ=255)
                                                                                                                     cccGACCTGAATCTTCAGGTTCTCGTCCCCTCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGATGACGCAGCAGGCAATCGCGTAATCGTCGCTGtt                               <  1:520086/150‑1 (MQ=255)
                                                                                                                         cccTGAATCTTCGGTTCTTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGt                          <  1:626832/150‑1 (MQ=255)
                                                                                                                          ccTGAATCTTCAGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAg                           <  2:579636/149‑1 (MQ=255)
                                                                                                                            ctgatttgtcaatttCTCTTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAGCTGCAGTTGCTTACCAATCCCCTTTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTc                         <  1:887618/137‑1 (MQ=255)
                                                                                                                                tCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACg                   >  2:470076/1‑151 (MQ=255)
                                                                                                                                 tttccgggttCCCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAg                       <  2:949820/139‑1 (MQ=255)
                                                                                                                                 cTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGc                  >  2:1130219/1‑151 (MQ=255)
                                                                                                                                    cAGCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTAc                    <  2:905900/146‑1 (MQ=255)
                                                                                                                                     accTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTAc                    <  1:136313/143‑1 (MQ=255)
                                                                                                                                         ttCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATc               >  2:610212/1‑147 (MQ=255)
                                                                                                                                              gTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTc       >  1:1411425/1‑150 (MQ=255)
                                                                                                                                               tCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCAt     >  1:1339118/1‑151 (MQ=255)
                                                                                                                                                  aCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCAtttt  >  1:1481060/1‑151 (MQ=255)
                                                                                                                                                  |                                                                                                                                                      
GTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTT  >  NC_000913/3260910‑3261205

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