Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 954,813 (A)7→6 intergenic (‑347/+59) focA ← / ← ycaO formate channel/ribosomal protein S12 methylthiotransferase accessory factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913954,8070A.95.7% 493.6 / 14.0 115intergenic (‑341/+65)focA/ycaOformate channel/ribosomal protein S12 methylthiotransferase accessory factor
Reads supporting (aligned to +/- strand):  ref base A (2/3);  new base . (64/46);  total (66/49)
Fisher's exact test for biased strand distribution p-value = 6.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 130 total aligned reads displayed.

GTAACAATATCGATCGTATTTATTAAATACAAATTACCGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATTAAAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGATTTTATTTTGCCCAGAATGCTGCTTTGGCGCGCTGCA  >  NC_000913/954708‑954906
                                                                                                   |                                                                                                   
gTAACAATATCGATCGTATTTATTAAATACAAATTACCGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATTa                                                                                                     <  1:3095950/100‑1 (MQ=255)
gTAACAATATCGATCGTATTTATTAAATACAAATTACCGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATTa                                                                                                     <  2:1951408/100‑1 (MQ=255)
    cAATATCGATCGTATTTATTAAATACAAATTACCGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATTaaaaa                                                                                                 <  1:1659539/100‑1 (MQ=255)
    cAATATCGATCGTATTTATTAAATACAAATTACCGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATTaaaaa                                                                                                 <  2:521641/100‑1 (MQ=255)
    cAATATCGATCGTATTTATTAAATACAAATTACCGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATTaaaaa                                                                                                 >  2:2000767/1‑100 (MQ=255)
     aaTATCGATCGTATTTATTAAATACAAATTACCGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATTaaaaaa                                                                                                >  2:1194776/1‑100 (MQ=255)
     aaTATCGATCGTATTTATTAAATACAAATTACCGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATTaaaaaa                                                                                                >  2:1195089/1‑100 (MQ=255)
      atatCGATCGTATTTATTAAATACAAATTACCGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATTaaaaaac                                                                                               <  2:2932932/100‑2 (MQ=255)
       tatCGATCGTATTTATTAAATACAAATTACCGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATTaaaaaacg                                                                                              >  1:1633934/1‑98 (MQ=255)
          cgctcgTATTTATTAAATACAAATTACCGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATTaaaaaacgcc                                                                                            <  1:2465773/96‑5 (MQ=255)
           gatcgTATTTATTAAATACAAATTACCGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATTaaaaaacgc                                                                                             >  2:2064089/1‑94 (MQ=255)
           gatcgTATTTATTAAATACAAATTACCGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATT‑AAAAAACGcccc                                                                                         <  1:1270656/100‑1 (MQ=255)
           gatcgTATTTATTAAATACAAATTACCGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATT‑AAAAAACGccc                                                                                          >  2:1011853/1‑99 (MQ=255)
            atcgTATTTATTAAATACAAATTACCGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATT‑AAAAAACGCCCCg                                                                                        >  2:208134/1‑100 (MQ=255)
            atcgTATTTATTAAATACAAATTACCGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATT‑AAAAAACGCCCCg                                                                                        >  1:3290006/1‑100 (MQ=255)
            atcgTATTTATTAAATACAAATTACCGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATT‑AAAAAACGCCCCg                                                                                        <  1:109945/100‑1 (MQ=255)
               gTATTTATTAAATACAAATTACCGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATT‑AAAAAACGCCCCGTaa                                                                                     <  1:3075337/100‑1 (MQ=255)
                 atttattAAATACAAATTACCGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATTAAAAAAACGCCCCGTAat                                                                                    <  1:2634659/100‑1 (MQ=255)
                   ttattaAATACAAATTACCGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAAtata                                                                                 >  2:1009739/1‑100 (MQ=255)
                     attaAATACAAATTACCGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAAc                                                                               >  2:612657/1‑100 (MQ=255)
                     attaAATACAAATTACCGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAAc                                                                               >  2:1087557/1‑100 (MQ=255)
                        aaaTACAAATTACCGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCa                                                                            >  2:2727388/1‑100 (MQ=255)
                         aaTACAAATTACCGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATTAAAAAAACGCCCCGTAATATAACTCa                                                                            >  1:1029060/1‑100 (MQ=255)
                              aaaTTACCGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCAGACTaa                                                                      <  1:1087559/100‑1 (MQ=255)
                               aaTTACCGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAAt                                                                     <  2:1633932/100‑1 (MQ=255)
                                aTTACCGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAATc                                                                    <  1:2478620/100‑1 (MQ=255)
                                aTTACCGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAATc                                                                    <  2:1719106/100‑1 (MQ=255)
                                aTTACCGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAATc                                                                    <  1:2868170/100‑1 (MQ=255)
                                   aCCGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAATCAtt                                                                 <  1:2075618/100‑1 (MQ=255)
                                   aCCGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAAAAAAAAA‑AAAAAAAACGCCCCGTAATATAACTCAGACTAATCAt                                                                  <  2:2570328/99‑1 (MQ=255)
                                    ccGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAATCATTa                                                                >  2:948111/1‑100 (MQ=255)
                                    ccGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAATCATTa                                                                <  1:1079134/100‑1 (MQ=255)
                                    ccGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAATCATTa                                                                >  1:2621453/1‑100 (MQ=255)
                                    ccGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAAACATTa                                                                >  1:1414114/1‑100 (MQ=255)
                                     cGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAATCATTaa                                                               <  1:1369718/100‑1 (MQ=255)
                                      gATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAg                                                              >  1:3266374/1‑100 (MQ=255)
                                      gATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAG‑AAAAAAAAAAAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAg                                                              <  2:343281/100‑1 (MQ=255)
                                        tatTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGcc                                                            >  1:2756972/1‑100 (MQ=255)
                                        tatTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGcc                                                            >  1:5000/1‑100 (MQ=255)
                                           ttAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTAc                                                         >  2:2875383/1‑100 (MQ=255)
                                            cgatcTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACa                                                        <  2:433319/96‑1 (MQ=255)
                                             aaCCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACAt                                                       <  1:2916737/100‑1 (MQ=255)
                                              aCCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGTCTACAtt                                                      <  1:1799607/100‑1 (MQ=255)
                                               ccTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTg                                                     <  2:3000596/100‑1 (MQ=255)
                                               ccTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTg                                                     >  2:2133174/1‑100 (MQ=255)
                                                cTTATAATTACAAATATTATATAAAAGCAAATATATGTAAAGCGTTTCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACAtt                                                      >  2:555185/1‑98 (MQ=255)
                                                  tataATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACaagatcg                                                  >  1:433320/1‑94 (MQ=255)
                                                  tataATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGt                                                  <  2:745920/100‑1 (MQ=255)
                                                    taATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAg                                                <  2:5000/100‑1 (MQ=255)
                                                     aaTTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAgg                                               >  1:980109/1‑100 (MQ=255)
                                                     aaTTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAgg                                               >  2:135701/1‑100 (MQ=255)
                                                     aaTTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAAATCAGACTAATCATTAAGCCTACATTGCGTAg                                                >  2:321257/1‑99 (MQ=255)
                                                     aaTTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAACATAACTCAGACTAATCATTAAGCCTACATTGCGTAgg                                               <  2:1829713/100‑1 (MQ=255)
                                                      aTTACAATTATTTTATTAATGCAAATCTATGTAAAGCGGGGCATT‑CAATAACGCCCCGAAATATAACTCAGACTAATCATTAAGCCTACATTACGTAGGc                                              >  2:2656640/1‑100 (MQ=255)
                                                      aTTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTaaga                                              >  1:2160655/1‑97 (MQ=255)
                                                         aCAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTAt                                           >  1:1449899/1‑100 (MQ=255)
                                                         aCAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTAt                                           >  2:90062/1‑100 (MQ=255)
                                                          cAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTAtt                                          >  1:2320686/1‑100 (MQ=255)
                                                          cAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTAtt                                          >  1:1489799/1‑100 (MQ=255)
                                                           aaTTATTTTATTAATGCAAATATATGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTAttt                                         >  2:1139967/1‑100 (MQ=255)
                                                                  ttattaATGCAAATATATGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCGGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGatttta                                  >  2:2389441/1‑100 (MQ=255)
                                                                  ttattaATGCAAATATATGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGatttta                                  <  1:1195092/100‑1 (MQ=255)
                                                                  ttattaATGCAAATATATGTAAAGCGGGGAATT‑AAAAAACGCCCCGTAATATAAATCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGatttt                                   >  1:2974605/1‑99 (MQ=255)
                                                                   tattaATGCAAATATATGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGattttat                                 >  1:251292/1‑100 (MQ=255)
                                                                   tattaATGCAAATATATGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGattttat                                 >  1:3317747/1‑100 (MQ=255)
                                                                   tattaATGCAAATATATGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGattttat                                 >  1:1912643/1‑100 (MQ=255)
                                                                   tattaATGCAAATATATGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGattttat                                 >  2:1025349/1‑100 (MQ=255)
                                                                    attaATGCAAATATATGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGattttatt                                >  2:3767/1‑100 (MQ=255)
                                                                    attaATGCAAATATATGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGattttatt                                <  1:1139970/100‑1 (MQ=255)
                                                                      taATGCAAATATATGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGattttatttt                              <  1:1492998/100‑1 (MQ=255)
                                                                      taATGCAAATATATGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGattttatttt                              <  1:588153/100‑1 (MQ=255)
                                                                       aaTGCAAATATATGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGattttatagat                             >  2:1385342/1‑96 (MQ=255)
                                                                       aaTGCAAATATATGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGATTTTATTTTg                             >  2:2348054/1‑100 (MQ=255)
                                                                        aTGCAAATATATGTAAAGCGGGGCATT‑AATAAACGCCCCGTAATATAACTCAGACTAATCATTAGGTCTACATTGCGTAGGCTATTTGATTTTATTTTGc                            <  2:2050677/100‑1 (MQ=255)
                                                                        aTGCAAATATATGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGATTTTATTTTGc                            >  1:162607/1‑100 (MQ=255)
                                                                         tGCAAATATATGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGATTTTATTTTGcc                           >  2:179826/1‑100 (MQ=255)
                                                                         tGCAAATATATGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGATTTTATTTTGcc                           >  2:3166168/1‑100 (MQ=255)
                                                                         tGCAAATATATGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGATTTTATTTTGcc                           <  1:2133179/100‑1 (MQ=255)
                                                                          gCAAATATATGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGATTTTATTTTGccc                          >  2:3120432/1‑100 (MQ=255)
                                                                          gCAAATATATGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGATTTTATTTTGccc                          >  1:1193884/1‑100 (MQ=255)
                                                                          gCAAATATATGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGATTTTATTTTGccc                          <  2:2320679/100‑1 (MQ=255)
                                                                            aaaTATATGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGATTTTATTTTGCCCAg                        <  1:458428/100‑1 (MQ=255)
                                                                            aaaTATATGTAAAGCGCGGCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGATTTTATTTTGCCCAg                        <  2:1989543/100‑1 (MQ=255)
                                                                             aaTATATGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAACCATTAAGCCTACATTGCGTAGGCTATTTGATTTTATTTTGCCCAGa                       <  2:1885888/100‑1 (MQ=255)
                                                                               tatatGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGATTTTATTTTGCCCAGAAt                     <  1:2180362/100‑1 (MQ=255)
                                                                               tatatGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGATTTTATTTTGCCCAGAAt                     <  1:628301/100‑1 (MQ=255)
                                                                               tatatGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGATTTTATTTTGCCCAGAAt                     <  1:901151/100‑1 (MQ=255)
                                                                               atcttGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGATTTTATTTTGCCCAGAAt                     <  1:712825/96‑1 (MQ=255)
                                                                                 tatGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGATTTTATTTTGCCCAGAAtgc                   >  1:468622/1‑100 (MQ=255)
                                                                                  atGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGATTTTATTTTGCCCAGAAtgct                  >  1:2898577/1‑100 (MQ=255)
                                                                                   tGTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGATTTTATTTTGCCCAGAAtagat                 >  2:712825/1‑96 (MQ=255)
                                                                                    gTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGATTTTATTTTGCCCAGAAtgctgc                >  1:3072777/1‑100 (MQ=255)
                                                                                    gTAAAGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGATTTTATTTTGCCCAGAAtgctgc                >  1:115190/1‑100 (MQ=255)
                                                                                      ctaGCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGATTTTATTTTGCCCAGAATGCTGCtt              <  2:863817/98‑1 (MQ=255)
                                                                                         gCGGGGCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGATTTTATTTTGCCCAGAATGCTGCTTTgg           <  2:1489796/100‑1 (MQ=255)
                                                                                          cGGGGCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGATTTTATTTTGCCCAGAATGCTGCTTTGgc          <  2:2324115/100‑1 (MQ=255)
                                                                                           ggggCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGATTTTATTTTGCCCAGAATGCTGCTTTGgcg         >  2:3301360/1‑100 (MQ=255)
                                                                                             ggCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGTTTTTATTTTTCTTTTTTTGCTGCTTTGgcgcg       >  2:177059/1‑100 (MQ=255)
                                                                                              gCATT‑AAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGATTTTATTTTGCCCAGAATGCTGCTTTGgcgcgc      >  2:1709216/1‑100 (MQ=255)
                                                                                                   aaaaaaaCGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGATTTTATTTTGCCCAGAATGCTGCTTTGGCGCGCTGCa  <  2:1029059/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
GTAACAATATCGATCGTATTTATTAAATACAAATTACCGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATTAAAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGATTTTATTTTGCCCAGAATGCTGCTTTGGCGCGCTGCA  >  NC_000913/954708‑954906

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