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)
           gatcgTATTTATTAAATACAAATTACCGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATTAAAAAACGcccc                                                                                         <  1:1270656/100‑1 (MQ=255)
           gatcgTATTTATTAAATACAAATTACCGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATTAAAAAACGccc                                                                                          >  2:1011853/1‑99 (MQ=255)
            atcgTATTTATTAAATACAAATTACCGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATTAAAAAACGCCCCg                                                                                        >  2:208134/1‑100 (MQ=255)
            atcgTATTTATTAAATACAAATTACCGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATTAAAAAACGCCCCg                                                                                        >  1:3290006/1‑100 (MQ=255)
            atcgTATTTATTAAATACAAATTACCGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATTAAAAAACGCCCCg                                                                                        <  1:109945/100‑1 (MQ=255)
               gTATTTATTAAATACAAATTACCGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATTAAAAAACGCCCCGTaa                                                                                     <  1:3075337/100‑1 (MQ=255)
                 atttattAAATACAAATTACCGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATTAAAAAAACGCCCCGTAat                                                                                    <  1:2634659/100‑1 (MQ=255)
                   ttattaAATACAAATTACCGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATTAAAAAACGCCCCGTAAtata                                                                                 >  2:1009739/1‑100 (MQ=255)
                     attaAATACAAATTACCGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATTAAAAAACGCCCCGTAATATAAc                                                                               >  2:612657/1‑100 (MQ=255)
                     attaAATACAAATTACCGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATTAAAAAACGCCCCGTAATATAAc                                                                               >  2:1087557/1‑100 (MQ=255)
                        aaaTACAAATTACCGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATTAAAAAACGCCCCGTAATATAACTCa                                                                            >  2:2727388/1‑100 (MQ=255)
                         aaTACAAATTACCGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATTAAAAAAACGCCCCGTAATATAACTCa                                                                            >  1:1029060/1‑100 (MQ=255)
                              aaaTTACCGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATTAAAAAACGCCCCGTAATATAACTCAGACTaa                                                                      <  1:1087559/100‑1 (MQ=255)
                               aaTTACCGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATTAAAAAACGCCCCGTAATATAACTCAGACTAAt                                                                     <  2:1633932/100‑1 (MQ=255)
                                aTTACCGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATTAAAAAACGCCCCGTAATATAACTCAGACTAATc                                                                    <  1:2478620/100‑1 (MQ=255)
                                aTTACCGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATTAAAAAACGCCCCGTAATATAACTCAGACTAATc                                                                    <  2:1719106/100‑1 (MQ=255)
                                aTTACCGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATTAAAAAACGCCCCGTAATATAACTCAGACTAATc                                                                    <  1:2868170/100‑1 (MQ=255)
                                   aCCGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATTAAAAAACGCCCCGTAATATAACTCAGACTAATCAtt                                                                 <  1:2075618/100‑1 (MQ=255)
                                   aCCGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAAAAAAAAAAAAAAAAACGCCCCGTAATATAACTCAGACTAATCAt                                                                  <  2:2570328/99‑1 (MQ=255)
                                    ccGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATTAAAAAACGCCCCGTAATATAACTCAGACTAATCATTa                                                                >  2:948111/1‑100 (MQ=255)
                                    ccGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATTAAAAAACGCCCCGTAATATAACTCAGACTAATCATTa                                                                <  1:1079134/100‑1 (MQ=255)
                                    ccGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATTAAAAAACGCCCCGTAATATAACTCAGACTAATCATTa                                                                >  1:2621453/1‑100 (MQ=255)
                                    ccGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATTAAAAAACGCCCCGTAATATAACTCAGACTAAACATTa                                                                >  1:1414114/1‑100 (MQ=255)
                                     cGATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATTAAAAAACGCCCCGTAATATAACTCAGACTAATCATTaa                                                               <  1:1369718/100‑1 (MQ=255)
                                      gATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATTAAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAg                                                              >  1:3266374/1‑100 (MQ=255)
                                      gATATTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGAAAAAAAAAAAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAg                                                              <  2:343281/100‑1 (MQ=255)
                                        tatTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATTAAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGcc                                                            >  1:2756972/1‑100 (MQ=255)
                                        tatTTAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATTAAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGcc                                                            >  1:5000/1‑100 (MQ=255)
                                           ttAACCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATTAAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTAc                                                         >  2:2875383/1‑100 (MQ=255)
                                            cgatcTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATTAAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACa                                                        <  2:433319/96‑1 (MQ=255)
                                             aaCCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATTAAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACAt                                                       <  1:2916737/100‑1 (MQ=255)
                                              aCCTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATTAAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGTCTACAtt                                                      <  1:1799607/100‑1 (MQ=255)
                                               ccTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATTAAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTg                                                     <  2:3000596/100‑1 (MQ=255)
                                               ccTTATAATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATTAAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTg                                                     >  2:2133174/1‑100 (MQ=255)
                                                cTTATAATTACAAATATTATATAAAAGCAAATATATGTAAAGCGTTTCATTAAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACAtt                                                      >  2:555185/1‑98 (MQ=255)
                                                  tataATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATTAAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACaagatcg                                                  >  1:433320/1‑94 (MQ=255)
                                                  tataATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATTAAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGt                                                  <  2:745920/100‑1 (MQ=255)
                                                    taATTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATTAAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAg                                                <  2:5000/100‑1 (MQ=255)
                                                     aaTTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATTAAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAgg                                               >  1:980109/1‑100 (MQ=255)
                                                     aaTTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATTAAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAgg                                               >  2:135701/1‑100 (MQ=255)
                                                     aaTTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATTAAAAAACGCCCCGTAATATAAATCAGACTAATCATTAAGCCTACATTGCGTAg                                                >  2:321257/1‑99 (MQ=255)
                                                     aaTTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATTAAAAAACGCCCCGTAACATAACTCAGACTAATCATTAAGCCTACATTGCGTAgg                                               <  2:1829713/100‑1 (MQ=255)
                                                      aTTACAATTATTTTATTAATGCAAATCTATGTAAAGCGGGGCATTCAATAACGCCCCGAAATATAACTCAGACTAATCATTAAGCCTACATTACGTAGGc                                              >  2:2656640/1‑100 (MQ=255)
                                                      aTTACAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATTAAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTaaga                                              >  1:2160655/1‑97 (MQ=255)
                                                         aCAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATTAAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTAt                                           >  1:1449899/1‑100 (MQ=255)
                                                         aCAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATTAAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTAt                                           >  2:90062/1‑100 (MQ=255)
                                                          cAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATTAAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTAtt                                          >  1:2320686/1‑100 (MQ=255)
                                                          cAATTATTTTATTAATGCAAATATATGTAAAGCGGGGCATTAAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTAtt                                          >  1:1489799/1‑100 (MQ=255)
                                                           aaTTATTTTATTAATGCAAATATATGTAAAGCGGGGCATTAAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTAttt                                         >  2:1139967/1‑100 (MQ=255)
                                                                  ttattaATGCAAATATATGTAAAGCGGGGCATTAAAAAACGCCCCGTAATATAACTCGGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGatttta                                  >  2:2389441/1‑100 (MQ=255)
                                                                  ttattaATGCAAATATATGTAAAGCGGGGCATTAAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGatttta                                  <  1:1195092/100‑1 (MQ=255)
                                                                  ttattaATGCAAATATATGTAAAGCGGGGAATTAAAAAACGCCCCGTAATATAAATCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGatttt                                   >  1:2974605/1‑99 (MQ=255)
                                                                   tattaATGCAAATATATGTAAAGCGGGGCATTAAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGattttat                                 >  1:251292/1‑100 (MQ=255)
                                                                   tattaATGCAAATATATGTAAAGCGGGGCATTAAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGattttat                                 >  1:3317747/1‑100 (MQ=255)
                                                                   tattaATGCAAATATATGTAAAGCGGGGCATTAAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGattttat                                 >  1:1912643/1‑100 (MQ=255)
                                                                   tattaATGCAAATATATGTAAAGCGGGGCATTAAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGattttat                                 >  2:1025349/1‑100 (MQ=255)
                                                                    attaATGCAAATATATGTAAAGCGGGGCATTAAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGattttatt                                >  2:3767/1‑100 (MQ=255)
                                                                    attaATGCAAATATATGTAAAGCGGGGCATTAAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGattttatt                                <  1:1139970/100‑1 (MQ=255)
                                                                      taATGCAAATATATGTAAAGCGGGGCATTAAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGattttatttt                              <  1:1492998/100‑1 (MQ=255)
                                                                      taATGCAAATATATGTAAAGCGGGGCATTAAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGattttatttt                              <  1:588153/100‑1 (MQ=255)
                                                                       aaTGCAAATATATGTAAAGCGGGGCATTAAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGattttatagat                             >  2:1385342/1‑96 (MQ=255)
                                                                       aaTGCAAATATATGTAAAGCGGGGCATTAAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGATTTTATTTTg                             >  2:2348054/1‑100 (MQ=255)
                                                                        aTGCAAATATATGTAAAGCGGGGCATTAATAAACGCCCCGTAATATAACTCAGACTAATCATTAGGTCTACATTGCGTAGGCTATTTGATTTTATTTTGc                            <  2:2050677/100‑1 (MQ=255)
                                                                        aTGCAAATATATGTAAAGCGGGGCATTAAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGATTTTATTTTGc                            >  1:162607/1‑100 (MQ=255)
                                                                         tGCAAATATATGTAAAGCGGGGCATTAAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGATTTTATTTTGcc                           >  2:179826/1‑100 (MQ=255)
                                                                         tGCAAATATATGTAAAGCGGGGCATTAAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGATTTTATTTTGcc                           >  2:3166168/1‑100 (MQ=255)
                                                                         tGCAAATATATGTAAAGCGGGGCATTAAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGATTTTATTTTGcc                           <  1:2133179/100‑1 (MQ=255)
                                                                          gCAAATATATGTAAAGCGGGGCATTAAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGATTTTATTTTGccc                          >  2:3120432/1‑100 (MQ=255)
                                                                          gCAAATATATGTAAAGCGGGGCATTAAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGATTTTATTTTGccc                          >  1:1193884/1‑100 (MQ=255)
                                                                          gCAAATATATGTAAAGCGGGGCATTAAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGATTTTATTTTGccc                          <  2:2320679/100‑1 (MQ=255)
                                                                            aaaTATATGTAAAGCGGGGCATTAAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGATTTTATTTTGCCCAg                        <  1:458428/100‑1 (MQ=255)
                                                                            aaaTATATGTAAAGCGCGGCATTAAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGATTTTATTTTGCCCAg                        <  2:1989543/100‑1 (MQ=255)
                                                                             aaTATATGTAAAGCGGGGCATTAAAAAACGCCCCGTAATATAACTCAGACTAACCATTAAGCCTACATTGCGTAGGCTATTTGATTTTATTTTGCCCAGa                       <  2:1885888/100‑1 (MQ=255)
                                                                               tatatGTAAAGCGGGGCATTAAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGATTTTATTTTGCCCAGAAt                     <  1:2180362/100‑1 (MQ=255)
                                                                               tatatGTAAAGCGGGGCATTAAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGATTTTATTTTGCCCAGAAt                     <  1:628301/100‑1 (MQ=255)
                                                                               tatatGTAAAGCGGGGCATTAAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGATTTTATTTTGCCCAGAAt                     <  1:901151/100‑1 (MQ=255)
                                                                               atcttGTAAAGCGGGGCATTAAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGATTTTATTTTGCCCAGAAt                     <  1:712825/96‑1 (MQ=255)
                                                                                 tatGTAAAGCGGGGCATTAAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGATTTTATTTTGCCCAGAAtgc                   >  1:468622/1‑100 (MQ=255)
                                                                                  atGTAAAGCGGGGCATTAAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGATTTTATTTTGCCCAGAAtgct                  >  1:2898577/1‑100 (MQ=255)
                                                                                   tGTAAAGCGGGGCATTAAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGATTTTATTTTGCCCAGAAtagat                 >  2:712825/1‑96 (MQ=255)
                                                                                    gTAAAGCGGGGCATTAAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGATTTTATTTTGCCCAGAAtgctgc                >  1:3072777/1‑100 (MQ=255)
                                                                                    gTAAAGCGGGGCATTAAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGATTTTATTTTGCCCAGAAtgctgc                >  1:115190/1‑100 (MQ=255)
                                                                                      ctaGCGGGGCATTAAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGATTTTATTTTGCCCAGAATGCTGCtt              <  2:863817/98‑1 (MQ=255)
                                                                                         gCGGGGCATTAAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGATTTTATTTTGCCCAGAATGCTGCTTTgg           <  2:1489796/100‑1 (MQ=255)
                                                                                          cGGGGCATTAAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGATTTTATTTTGCCCAGAATGCTGCTTTGgc          <  2:2324115/100‑1 (MQ=255)
                                                                                           ggggCATTAAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGATTTTATTTTGCCCAGAATGCTGCTTTGgcg         >  2:3301360/1‑100 (MQ=255)
                                                                                             ggCATTAAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGTTTTTATTTTTCTTTTTTTGCTGCTTTGgcgcg       >  2:177059/1‑100 (MQ=255)
                                                                                              gCATTAAAAAACGCCCCGTAATATAACTCAGACTAATCATTAAGCCTACATTGCGTAGGCTATTTGATTTTATTTTGCCCAGAATGCTGCTTTGgcgcgc      >  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: