Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 497,139:1 (T)5→6 100% intergenic (+145/‑36) htpG → / → adk protein refolding molecular co‑chaperone Hsp90, Hsp70‑dependent; heat‑shock protein; ATPase/adenylate kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913497,1341.T96.2% 249.7 / 2.8 78intergenic (+140/‑41)htpG/adkprotein refolding molecular co‑chaperone Hsp90, Hsp70‑dependent; heat‑shock protein; ATPase/adenylate kinase
Reads supporting (aligned to +/- strand):  ref base . (1/2);  new base T (33/42);  total (34/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTCATCATCTGCACTTTCCGCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCGGGCGCGGGGAAAGGGACTCAGG  >  NC_000913/497035‑497223
                                                                                                     |                                                                                         
tttcATCATCTGCACTTTCCGCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGc                                                                                            >  1:2920375/1‑100 (MQ=255)
 ttcATCATCTGCACTTTCCGCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCt                                                                                          >  1:786480/1‑100 (MQ=255)
   catcatCTGCACTTTCCGCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCttt                                                                                        >  2:314243/1‑100 (MQ=255)
     tcatcTGCACTTTCCGTAAATTATCTCGCCATTAACTGTTTCAGCCCTAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCttttt                                                                                      <  1:937676/100‑1 (MQ=255)
     tcatcTGCACTTTCCGCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCttttt                                                                                      <  1:1660373/100‑1 (MQ=255)
     tcagcTGCACTTTCCGCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCttttt                                                                                      <  2:394797/100‑1 (MQ=255)
       atcTGCACTTTCCGCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCttttttc                                                                                    >  1:2853759/1‑98 (MQ=255)
       atcTGCACTTTCCGCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCttttttc                                                                                    >  1:2335142/1‑98 (MQ=255)
        tcTGCACTTTCCGCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCttttttcc                                                                                   <  2:2888656/100‑4 (MQ=255)
          tGCACTTTCCGCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCttttttcaaa                                                                                 >  2:2290586/1‑95 (MQ=255)
           gCACTTTCCGCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTTGTATCGTTTATCGCttttttcaaa                                                                                 <  2:130901/99‑6 (MQ=255)
           gCACTTTCCGCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCaaaa                                                                                 >  2:381203/1‑100 (MQ=255)
           gCACTTTCCGCAAATTATCTCGCAATTAACCGTTTCAGCCCCAGGTTCCTTTCTTGAGTCAATCACATGTTGTTGGTATCGCTTATCGCttttttcacaa                                                                                >  2:99117/1‑94 (MQ=255)
              cTTTCCGCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAAt                                                                              >  2:2344396/1‑100 (MQ=255)
              cTTTCCGCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTATTaaaaaaat                                                                             <  2:2930066/100‑3 (MQ=255)
               tttCCGCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAAtt                                                                             <  2:505558/100‑1 (MQ=255)
                ttCCGCAAATTATTTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTc                                                                            <  2:1464086/100‑1 (MQ=255)
                 tCCGCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCg                                                                           <  2:1123194/100‑1 (MQ=255)
                 tCCGCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCg                                                                           <  1:1396475/100‑1 (MQ=255)
                 tCCGCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCg                                                                           <  1:2917301/100‑1 (MQ=255)
                  ccGCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTAAAACAATCGCCAGATGGAAGTATCGTTTAACGCTTTTTTCAAAAAATTCGa                                                                          >  1:348093/1‑100 (MQ=255)
                   cGCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGac                                                                         >  2:2506716/1‑100 (MQ=255)
                    gCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGaca                                                                        >  1:166960/1‑100 (MQ=255)
                    gCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGac                                                                         <  2:2737095/99‑1 (MQ=255)
                    gCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGATTATCGCTTTTTCAAAAAATTCGacac                                                                       >  1:2024893/1‑100 (MQ=255)
                       aaTTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTAGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACAt                                                                     >  1:919653/1‑100 (MQ=255)
                         ttATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACAtt                                                                    >  1:2405267/1‑99 (MQ=255)
                           aTCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTa                                                                 <  1:1926007/100‑1 (MQ=255)
                            tctcGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGATTTTTTAAAAAAATTCGACACATTTTa                                                                 <  2:1051108/99‑1 (MQ=255)
                            tctcGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTaa                                                                >  2:496473/1‑100 (MQ=255)
                             ctcGCCATTAAGCGTTTCGGCCCCAGGTGCCTTTCTTGAGGCACTCGCCTGGTGGGGGTATCTTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAg                                                               >  2:1341714/1‑100 (MQ=255)
                              tcGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAgg                                                              <  2:824773/100‑1 (MQ=255)
                               cGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAggg                                                             >  2:33488/1‑100 (MQ=255)
                                  cATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGAt                                                          >  1:571163/1‑100 (MQ=255)
                                  cATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGAt                                                          >  2:928116/1‑100 (MQ=255)
                                       aCCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCg                                                     <  2:2011209/100‑1 (MQ=255)
                                       aCCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCg                                                     >  2:1011975/1‑100 (MQ=255)
                                       aCCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCg                                                     >  1:2701471/1‑100 (MQ=255)
                                        ccGTTTCAGCCTCAGTTGCCTTTCTTGAGGCAATGTCATGTTGGTGGTATCGTTTATCGCTTTTTTCATAAAATTCGACATTTTTTAAGGGGATTTTCGc                                                    <  2:2664431/100‑1 (MQ=255)
                                        ccGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGc                                                    <  2:2558605/100‑1 (MQ=255)
                                        ccGTTTCAGCCACAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATTGCTTTTTTCAAAAAATTGGTCACATTTTAAGGGGAtttt                                                       >  2:2262309/1‑97 (MQ=255)
                                         cttttCATCCCCATGTTCCTTTCTTAAGTCAATTGCCTGTTTGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGc                                                    <  1:1358631/97‑1 (MQ=255)
                                          gTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCTCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAATTTCGACACATTTTAAGGGGATTTTCGCaa                                                  >  2:1974799/1‑100 (MQ=255)
                                          gTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCaa                                                  <  1:2236924/100‑1 (MQ=255)
                                           tttCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCaa                                                  <  2:1784996/99‑1 (MQ=255)
                                             tCAGCCCCAGGTGCCTTTCTTGAGGCAAGCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGc                                               <  2:759833/100‑1 (MQ=255)
                                              cAGCTCCAGGTGCCTTTCTTGAGGCAATCGCCTTTTGGTGTTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCg                                              <  1:2685994/100‑1 (MQ=255)
                                              cAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTTGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCg                                              <  1:1613440/100‑1 (MQ=255)
                                              cAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCg                                              <  2:1842778/100‑1 (MQ=255)
                                                gCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTa                                            <  2:365536/100‑1 (MQ=255)
                                                 ccccAGGTGCCTTTCTTGAGGCAATCGCCTGTTGTTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTAt                                           <  2:1592515/100‑1 (MQ=255)
                                                      ggTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCAtcc                                      <  1:2042927/100‑3 (MQ=255)
                                                      ggTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTc                                      <  2:3021746/100‑1 (MQ=255)
                                                         gCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGc                                   >  2:1591718/1‑100 (MQ=255)
                                                          cttttCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACAGTTGAAGGGGATTTTCGCAATGCGTATCATTCTGCt                                  <  1:1286534/98‑1 (MQ=255)
                                                          ccTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCt                                  >  1:1681371/1‑100 (MQ=255)
                                                           cTTTCTTGAGGCAATCGCCTGTTGGTGTTATCGTTTATCTCTTTTTTCAAAAAATTCGACACATTTTAAGGGAATTTTCGCAATGCGTATCATTCTGCtt                                 >  1:970508/1‑100 (MQ=255)
                                                              tCTTTAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGgc                              <  2:377101/100‑1 (MQ=255)
                                                              tCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGgc                              <  1:2460354/100‑1 (MQ=255)
                                                              tCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGgc                              <  2:2223345/100‑1 (MQ=255)
                                                              tCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATGCTGCTTGgc                              <  1:2359899/100‑1 (MQ=255)
                                                                 tGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATAATTCTGCTTGGCGCt                           >  2:859585/1‑100 (MQ=255)
                                                                    ggCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTTAAAATTTTCGACACATTTTACGGGGCTTCTCGCAATGCGTATCATTCTGCCTGCCGCTcc                         >  2:2416554/1‑99 (MQ=255)
                                                                    ggCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCg                        >  2:2117184/1‑100 (MQ=255)
                                                                    ggCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTATAAGGGGATTTTCGCAAAGCGTATCATTCTGCTTGGAGCTCCg                        >  2:2267752/1‑100 (MQ=255)
                                                                      cAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATTCGTATCATTCTGCTTGgc                              >  2:636158/1‑92 (MQ=255)
                                                                      cAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCggg                      <  1:2626300/100‑1 (MQ=255)
                                                                      cAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCggg                      <  1:234766/100‑1 (MQ=255)
                                                                      cAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCggg                      >  2:765467/1‑100 (MQ=255)
                                                                        aTCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCGGgcg                    >  1:2573617/1‑100 (MQ=255)
                                                                         tCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCGGgcgc                   <  2:2920343/100‑1 (MQ=255)
                                                                           gCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCGGGCGCgg                 >  2:2152388/1‑100 (MQ=255)
                                                                            tcTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCGGGCGCggg                <  1:2210708/99‑1 (MQ=255)
                                                                            ccTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCGGGCGCggg                <  1:1407309/100‑1 (MQ=255)
                                                                             cTGTTGGTGGTATTGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCGGGCGCgg                 <  1:319049/98‑1 (MQ=255)
                                                                              tGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCGGGCGCGGGGa              <  1:2214443/100‑1 (MQ=255)
                                                                               gTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCGGGCGCgggaga             >  2:2210685/1‑97 (MQ=255)
                                                                               gTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCGGGCGCGGGGaa             <  1:314246/100‑1 (MQ=255)
                                                                                  ggtggtATCGTTTATCGCTTTTTTTATTTTTTTTTACACATTTTAAGGGGATTTTCGCACTGCGTATCATTCTGCTTGGCGCTCTGGGCGCGCGGAAAgg          >  2:1461085/1‑100 (MQ=255)
                                                                                  ggtggtATCGTTTATCGCTTTTTTTATTTTTTTTTAAACATTTTAAGGGGATATTAACAATGCGTATCATTCTGCTTGGCGCtttggtc                     >  2:1217382/1‑83 (MQ=255)
                                                                                  ggtggtATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCGGGCGCGGGGAAAgg          <  2:2693446/100‑1 (MQ=255)
                                                                                  ggtggtATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCGGGCGCGGGGAAAgg          >  2:1449859/1‑100 (MQ=255)
                                                                                  ggtggtATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCGCCGGGCGCGGGGAAAgg          <  1:381206/100‑1 (MQ=255)
                                                                                   gtggtATCGTTTATCGCTTTTTTTATATAATTCAAAACATATTAAGGGGATTTTCACAATGCGTATCATTCTGCTTGGCGCTCCGGGCGCGGAGAAAggg         >  2:999733/1‑100 (MQ=255)
                                                                                   gtggtATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCGGGCGCGGGGAAAggg         <  2:1681351/100‑1 (MQ=255)
                                                                                       tATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCGGGCGCGGGGAAAGGGACTc     >  1:360645/1‑100 (MQ=255)
                                                                                       tATCGTTTATCGCTTTTTTCAAAAAATGCGACACATTATAAGGAGATTTTCGCAATGCGTATCATTCTGCTTGCCGCTCCGGGCGCGGGGAAAGGGACTc     <  1:1341726/100‑1 (MQ=255)
                                                                                         tCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGTTTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCGGGCGCGGGGAAAGGGACTCAg   <  1:1787551/100‑1 (MQ=255)
                                                                                         tCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCGGGCGCGGGGAAAGGGACTCAg   <  1:1949122/100‑1 (MQ=255)
                                                                                          cGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCGGGCGCGGGGAAAGGGACTCAgg  <  1:765471/100‑1 (MQ=255)
                                                                                                     |                                                                                         
TTTCATCATCTGCACTTTCCGCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCGGGCGCGGGGAAAGGGACTCAGG  >  NC_000913/497035‑497223

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