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

GTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGC‑TTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTT  >  NC_000913/3260910‑3261205
                                                                                                                                                  |                                                                                                                                                      
gTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCGCTTTGGCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGC‑TTCTCGTCCAcccc                                                                                                                                                    <  1:862223/150‑1 (MQ=255)
gTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGC‑TTCTCGTCCACCccg                                                                                                                                                   <  2:1046305/151‑1 (MQ=255)
 tGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGC‑TTCTCATCCACCccgc                                                                                                                                                  >  2:1361904/1‑151 (MQ=255)
  gAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGC‑TTCTCGTCCACCccgcc                                                                                                                                                 <  2:1461150/151‑1 (MQ=255)
     tGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGC‑TTCTCGTCCACCCCGCCGt                                                                                                                                               >  2:1005495/1‑150 (MQ=255)
           tCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGC‑TTCTCGTCCACCCCGCCGTTAATTGc                                                                                                                                        >  1:58715/1‑151 (MQ=255)
           tCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGC‑TTCTCGTCCACCCCGCCGTTAATTGc                                                                                                                                        >  2:31866/1‑151 (MQ=255)
           tCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGC‑TTCTCGTCCACCCCGCCGTTAATTGc                                                                                                                                        >  1:1431804/1‑151 (MQ=255)
                gcgcTGTTTTCCGGCCCCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGC‑TTCTCGCCCCCCCCGCCGTTAATTg                                                                                                                                         <  2:142093/145‑1 (MQ=255)
                gcgcTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGC‑TTCTCGTCCACCCCGCCGTTAATTGCGt                                                                                                                                      >  2:1529700/1‑148 (MQ=255)
                gcgcTGATGTACTGCAACGCCAGCCAGTTCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGC‑TTCTCGTCCACCCCGCCGTTAATTGCGt                                                                                                                                      >  1:17887/1‑148 (MQ=255)
                     gATGTCCTGCCCCGCCGGCCATTCAATGAAGTGATCGAGGCTGTCCATCCCTTTGTCTAAATCCACCACGTCCCCCACCGGCGGTGCTGTTTGCGGCCCGCCTGGAATTTTCATT‑TCGTCGCCACCCCCGCCGTTAATTGCGTAGagca                                                                                                                                <  1:1487104/149‑1 (MQ=255)
                      aTGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGC‑TTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCg                                                                                                                             >  1:1161764/1‑151 (MQ=255)
                       tttccTGCCCCCCCACCCCGTCCATGAAGTGAGCGGGGCTGTCCATCTCTTTGTCGTAACCCAGCACTTCTTCCACCGGCGGTGCTGTTTTCGGCCCGCCCGAAATCTTCAGC‑TTCGCGCCCCCCCCGCCGTTAAtt                                                                                                                                          <  1:1160400/133‑1 (MQ=255)
                       tGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGC‑TTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCg                                                                                                                             >  1:170134/1‑150 (MQ=255)
                              cccccccagccagCCCATGAAGTGACCGAGGCTGTCCATCCCTTTGTCGTAATCCACCACGTCGCCCATCAGCGGTGCTGTTTTCGCCCCGCCCTGAATCTTCACC‑TTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGtt                                                                                                                           <  1:616556/140‑1 (MQ=255)
                              ccccccaaccctgTCCATAAGGTGACCGAGGCTGTCCACAATTTTGTCTAAATCCACCACATCGCCCATCAGCGGTGCTGTTTTCGGCCGGCCTTAAATTTTCAGT‑TGCGCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGtttt                                                                                                                         <  2:1219539/135‑1 (MQ=255)
                              ccccGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCTCTTTGTCGTAACCCAGCACGTCGTCCATCAGCGGTGCTTTTTTCGGCCCGACCTGAATCTTCAGC‑TTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTc                                                                                                                        <  2:330172/146‑1 (MQ=255)
                                 cGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGC‑TTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAgg                                                                                                                  >  2:241178/1‑151 (MQ=255)
                                   ccagccagTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGC‑TTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGtt                                                                                                                >  2:1028440/1‑151 (MQ=255)
                                       ccagTGCATAAAGTGCTCGGGGCTGTCCATCACTTTGTCTTAATCCCTCACGCCGTCCATCAGGGGTGTTGTTTTCGCCCCGCCTATATTCTTCAGC‑TTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAgc                                                                                                             <  1:1363514/150‑1 (MQ=255)
                                       ccagTCCATGAATTGATCGGGGCTGTCCATCCCTTTGTCGTACTCCCCCCCGTCGCCCATCGGCGGTCCTGTTTCCGGCCCGCCCGGAATCTTCAGC‑TTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAg                                                                                                              <  1:348898/149‑1 (MQ=255)
                                       ccagTCCATGAAGTGATCGAGGCTGTCCATCTCTTTGTCGTAATCCAGCACGTCGTCCATCGGCGGTGCTTTTTTCGGCCCGACCTAAATCTTCAGC‑TTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAgc                                                                                                             <  2:1192305/150‑1 (MQ=255)
                                       ccagCCCAGGAAGGGATCGAGGCTCTCCATCCTTTTGTCGTAACCCCCCCCGTCGCCAATCAGCGGTCCTTTTTTCGGCCCGCCCTAAATCTTCAGC‑TTCTCTTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTa                                                                                                               <  1:1414725/148‑1 (MQ=255)
                                         agTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGC‑TTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGcgcg                                                                                                          >  1:575435/1‑151 (MQ=255)
                                           tCCATGAAGTGATCGAGGCTGTCCATCCTTTTGTCGTAACCCAGCCCGTCGTCCATCAGCGGTGTTGTTTGCGGCCCGACCTAAATCTTCAGC‑TTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGt                                                                                                         <  1:1352090/150‑1 (MQ=255)
                                              aTGAAGGGATCGGGGCTGTCCATCTTTTTGTCGAAATCCACCACGTCGCCCACCGGCGGTGTTTTTTCCGCCCCGCCTTAAATCTTCAGC‑TTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGt                                                                                                         <  1:1495078/147‑1 (MQ=255)
                                               tGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGC‑TTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCAcc                                                                                                    >  2:1347631/1‑151 (MQ=255)
                                                gAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGC‑TTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCAcc                                                                                                    >  2:1047514/1‑150 (MQ=255)
                                                 aggttaTCGAGGTTGTCCATCTCTTTGTCGTAATCCACCACGTCGCCCATCAGGGTTGCTTTTTGCGGCCCGCCCTGACTCTTCAGT‑TTCGGGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCa                                                                                                   <  1:1291040/145‑1 (MQ=255)
                                                 aaGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGC‑TTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCa                                                                                                   >  2:965371/1‑150 (MQ=255)
                                                 aaGTGATCGAGGCTGCCCATCCCTTTGTCGTAACCCAGCCCTTCTTCCTTCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGC‑TTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCa                                                                                                   <  2:1161828/150‑1 (MQ=255)
                                                 aaGGGACCGGGGCTGTCCATCGTTTTGACGTAACCCAGCCCGTCGTCCATCAGCGGTGCTTTTTCCGGCCCGACCGGAATCTTAGTT‑TTCTCCTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCAc                                                                                                     <  2:1359327/148‑1 (MQ=255)
                                                   gTGATCGGGGCTGTCCATCACTTTGTCGTAACCCAGCACGTCGTCCATCGGCGGTGTTGTTTTCGGCCCGACCTGAATCTTCACC‑TTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAg                                                                                                <  2:1426653/151‑1 (MQ=255)
                                                       tCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGC‑TTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACt                                                                                            >  1:599877/1‑151 (MQ=255)
                                                        cGAGGCTGTCCATCATTTTGTCGTAACCCAGCCCGTCTTCCATCAGCGGTGTTTTTTTCGGCCCGACCTGAATCTTCAGC‑TTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACt                                                                                            <  2:147934/150‑1 (MQ=255)
                                                                tctatctttttGTCGTAATCCCGCACGTCGTCCATCGGCGTTGTTTTTTTCGCCCCGTCCTGAATCTTCAGT‑TTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTg                                                                                     <  2:815902/141‑1 (MQ=255)
                                                                 ccATCATTTTGTCGTAACCCAGCACGTCGTCCATCAGCGGTGTTTTTTTCGGCCCGCCCTGAATCTTCAGC‑TTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCtt                                                                                  <  1:1383313/151‑1 (MQ=255)
                                                                   aTCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGC‑TTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTAc                                                                                >  2:283820/1‑151 (MQ=255)
                                                                      aCTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGC‑TTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCa                                                                              >  2:1336211/1‑150 (MQ=255)
                                                                        tttGTCGTAACCCAGCCCGTCGCCCATCGGCGGTGCTTTTTTCGGCCCGCCCTGAATCTTCAGC‑TTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAAt                                                                            <  2:115150/150‑1 (MQ=255)
                                                                        ttggTCGTAATCCCGCACGTCGTCCCTCAGCGGTGCTGTTTCCGGCCCTACCTGAATCTTCAGC‑CTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAAt                                                                            <  2:572395/147‑1 (MQ=255)
                                                                          tGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGC‑TTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCAc                                                                         >  1:4/1‑151 (MQ=255)
                                                                            tCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGC‑TTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCAcc                                                                        >  2:1378291/1‑150 (MQ=255)
                                                                             cGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGC‑TTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCAt                                                                      >  1:836120/1‑151 (MQ=255)
                                                                             cGTAACCCGGCACGTCGCCCATCAGCGGTGTTTTTTGCGGCCCGCCCTAAATCTTCATC‑TTCCCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCa                                                                       <  1:923525/150‑1 (MQ=255)
                                                                                aaTCCAGCCCGTCGCCCATCAGCGGTGCTGTTTTCGGCCCGCCCTGAATCTTCAGC‑TTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCAt                                                                      <  2:1398419/148‑1 (MQ=255)
                                                                                   ccAGCCCGTCGCCCATCAGGGGTGTTGTTTTCGGCCCGCCCTGAATCTTCAGC‑TTCTCTTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTggg                                                                  <  1:956278/149‑1 (MQ=255)
                                                                                   ccAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGC‑TTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCt                                                                >  2:1489718/1‑151 (MQ=255)
                                                                                    cAGCACGTCGTCCATCCGCTGTTGTGTTTTCGCCCCGCCTTGATTCTCCGGC‑TCCCCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCGCCATTGGGc                                                                 <  1:2886/149‑1 (MQ=255)
                                                                                       cccgtcgtCCATCGGCGGTGTTTTTTCCGCCCCGACCGGAATTTTCACC‑TTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCAta                                                                     <  1:908556/140‑2 (MQ=255)
                                                                                           tcgccCATCAGCGGTGCTGTTTTCGGCCCGCCCTGAATCTTCAGC‑TTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACgca                                                          <  2:323120/149‑1 (MQ=255)
                                                                                                 accAGGGGTGTTGTTTCCGCCCCGA‑CAAAATCCTCAGCTTTCCCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACgca                                                          <  1:835078/141‑1 (MQ=255)
                                                                                                    aGCGGTGCTGTTTTCGGCCCGCCCTAAATCTTCAGC‑TTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATc                                                <  2:819601/150‑1 (MQ=255)
                                                                                                     gggggttttttttCGCCCCGCCATGAATCTTCAGC‑TCCTCGCCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATcg                                               <  1:595320/141‑1 (MQ=255)
                                                                                                     gCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGC‑TTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATc                                                <  1:1410243/149‑1 (MQ=255)
                                                                                                       ggTGTTGTTTCCGGCCCGACCTAAATCTTCACC‑TTCTCCCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGt                                            <  2:17890/151‑1 (MQ=255)
                                                                                                       ggTGCTTTTTTCGGCCCGCCCTAAATCTTCAGC‑TTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAgg                                                     <  1:1191190/142‑1 (MQ=255)
                                                                                                       ggTGCTTTTTTCGGCCCGACCTGAATCTTCAGC‑TTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCa                                                   <  1:1550352/144‑1 (MQ=255)
                                                                                                         tGCTGTTTTCGGCCCGACCTGAATCTTCAGC‑TTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTa                                           >  1:641849/1‑150 (MQ=255)
                                                                                                               ttgcGGCCCGCCCTGAATCTTCAGC‑TTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAAtcgtc                                     <  1:1378868/147‑1 (MQ=255)
                                                                                                                ttCGGCCCGACCTGAATCTTCAGC‑TTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGc                                   >  2:1269567/1‑151 (MQ=255)
                                                                                                                  cGGCCCGACCTGAATCTTCAGC‑TTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTg                                 >  1:1113936/1‑151 (MQ=255)
                                                                                                                     cccGCCCTAAATCTTCAGC‑TTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTg                              <  1:1041037/151‑1 (MQ=255)
                                                                                                                     cccGACCTGAATCTTCAGG‑TTCTCGTCCCCTCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGATGACGCAGCAGGCAATCGCGTAATCGTCGCTGtt                               <  1:520086/150‑1 (MQ=255)
                                                                                                                         cccTGAATCTTCGGT‑TCTTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGt                          <  1:626832/150‑1 (MQ=255)
                                                                                                                          ccTGAATCTTCAGC‑TTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAg                           <  2:579636/149‑1 (MQ=255)
                                                                                                                            ctgatttgtcaatttCTCTTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAGCTGCAGTTGCTTACCAATCCCCTTTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTc                         <  1:887618/137‑1 (MQ=255)
                                                                                                                                tCTTCAGC‑TTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACg                   >  2:470076/1‑151 (MQ=255)
                                                                                                                                 tttccgggttCCCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAg                       <  2:949820/139‑1 (MQ=255)
                                                                                                                                 cTTCAGC‑TTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGc                  >  2:1130219/1‑151 (MQ=255)
                                                                                                                                    cAGC‑TTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTAc                    <  2:905900/146‑1 (MQ=255)
                                                                                                                                     acc‑TTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTAc                    <  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)
                                                                                                                                                  |                                                                                                                                                      
GTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGC‑TTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTT  >  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: ‑