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

TTTCATCATCTGCACTTTCCGCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCG‑C‑TTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCGGGCGCGGGGAAAGGGACTCAGG  >  NC_000913/497035‑497223
                                                                                                     |                                                                                         
tttcATCATCTGCACTTTCCGCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCG‑c                                                                                            >  1:2920375/1‑100 (MQ=255)
 ttcATCATCTGCACTTTCCGCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCG‑C‑t                                                                                          >  1:786480/1‑100 (MQ=255)
   catcatCTGCACTTTCCGCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCG‑C‑ttt                                                                                        >  2:314243/1‑100 (MQ=255)
     tcatcTGCACTTTCCGTAAATTATCTCGCCATTAACTGTTTCAGCCCTAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCG‑C‑ttttt                                                                                      <  1:937676/100‑1 (MQ=255)
     tcatcTGCACTTTCCGCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCG‑C‑ttttt                                                                                      <  1:1660373/100‑1 (MQ=255)
     tcagcTGCACTTTCCGCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCG‑C‑ttttt                                                                                      <  2:394797/100‑1 (MQ=255)
       atcTGCACTTTCCGCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCG‑C‑ttttttc                                                                                    >  1:2853759/1‑98 (MQ=255)
       atcTGCACTTTCCGCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCG‑C‑ttttttc                                                                                    >  1:2335142/1‑98 (MQ=255)
        tcTGCACTTTCCGCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCG‑C‑ttttttcc                                                                                   <  2:2888656/100‑4 (MQ=255)
          tGCACTTTCCGCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCG‑C‑ttttttcaaa                                                                                 >  2:2290586/1‑95 (MQ=255)
           gCACTTTCCGCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTTGTATCGTTTATCG‑C‑ttttttcaaa                                                                                 <  2:130901/99‑6 (MQ=255)
           gCACTTTCCGCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCG‑CTTTTTTCaaaa                                                                                 >  2:381203/1‑100 (MQ=255)
           gCACTTTCCGCAAATTATCTCGCAATTAACCGTTTCAGCCCCAGGTTCCTTTCTTGAGTCAATCACATGTTGTTGGTATCGCTTATCG‑C‑ttttttcacaa                                                                                >  2:99117/1‑94 (MQ=255)
              cTTTCCGCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCG‑CTTTTTTCAAAAAAt                                                                              >  2:2344396/1‑100 (MQ=255)
              cTTTCCGCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCG‑C‑TTTATTaaaaaaat                                                                             <  2:2930066/100‑3 (MQ=255)
               tttCCGCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCG‑CTTTTTTCAAAAAAtt                                                                             <  2:505558/100‑1 (MQ=255)
                ttCCGCAAATTATTTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCG‑CTTTTTTCAAAAAATTc                                                                            <  2:1464086/100‑1 (MQ=255)
                 tCCGCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCG‑CTTTTTTCAAAAAATTCg                                                                           <  2:1123194/100‑1 (MQ=255)
                 tCCGCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCG‑CTTTTTTCAAAAAATTCg                                                                           <  1:1396475/100‑1 (MQ=255)
                 tCCGCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCG‑CTTTTTTCAAAAAATTCg                                                                           <  1:2917301/100‑1 (MQ=255)
                  ccGCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTAAAACAATCGCCAGATGGAAGTATCGTTTAACG‑CTTTTTTCAAAAAATTCGa                                                                          >  1:348093/1‑100 (MQ=255)
                   cGCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCG‑CTTTTTTCAAAAAATTCGac                                                                         >  2:2506716/1‑100 (MQ=255)
                    gCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCG‑CTTTTTTCAAAAAATTCGaca                                                                        >  1:166960/1‑100 (MQ=255)
                    gCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCG‑CTTTTTTCAAAAAATTCGac                                                                         <  2:2737095/99‑1 (MQ=255)
                    gCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGATTATCG‑C‑TTTTTCAAAAAATTCGacac                                                                       >  1:2024893/1‑100 (MQ=255)
                       aaTTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTAGGTGGTATCGTTTATCG‑CTTTTTTCAAAAAATTCGACACAt                                                                     >  1:919653/1‑100 (MQ=255)
                         ttATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCG‑CTTTTTTCAAAAAATTCGACACAtt                                                                    >  1:2405267/1‑99 (MQ=255)
                           aTCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCG‑CTTTTTTCAAAAAATTCGACACATTTTa                                                                 <  1:1926007/100‑1 (MQ=255)
                            tctcGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGAT‑TTTTTAAAAAAATTCGACACATTTTa                                                                 <  2:1051108/99‑1 (MQ=255)
                            tctcGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCG‑CTTTTTTCAAAAAATTCGACACATTTTaa                                                                >  2:496473/1‑100 (MQ=255)
                             ctcGCCATTAAGCGTTTCGGCCCCAGGTGCCTTTCTTGAGGCACTCGCCTGGTGGGGGTATCTTTTATCG‑CTTTTTTCAAAAAATTCGACACATTTTAAg                                                               >  2:1341714/1‑100 (MQ=255)
                              tcGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCG‑CTTTTTTCAAAAAATTCGACACATTTTAAgg                                                              <  2:824773/100‑1 (MQ=255)
                               cGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCG‑CTTTTTTCAAAAAATTCGACACATTTTAAggg                                                             >  2:33488/1‑100 (MQ=255)
                                  cATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCG‑CTTTTTTCAAAAAATTCGACACATTTTAAGGGGAt                                                          >  1:571163/1‑100 (MQ=255)
                                  cATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCG‑CTTTTTTCAAAAAATTCGACACATTTTAAGGGGAt                                                          >  2:928116/1‑100 (MQ=255)
                                       aCCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCG‑CTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCg                                                     <  2:2011209/100‑1 (MQ=255)
                                       aCCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCG‑CTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCg                                                     >  2:1011975/1‑100 (MQ=255)
                                       aCCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCG‑CTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCg                                                     >  1:2701471/1‑100 (MQ=255)
                                        ccGTTTCAGCCTCAGTTGCCTTTCTTGAGGCAATGTCATGTTGGTGGTATCGTTTATCG‑CTTTTTTCATAAAATTCGACATTTTTTAAGGGGATTTTCGc                                                    <  2:2664431/100‑1 (MQ=255)
                                        ccGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCG‑CTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGc                                                    <  2:2558605/100‑1 (MQ=255)
                                        ccGTTTCAGCCACAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATTG‑CTTTTTTCAAAAAATTGGTCACATTTTAAGGGGAtttt                                                       >  2:2262309/1‑97 (MQ=255)
                                         cttttCATCCCCATGTTCCTTTCTTAAGTCAATTGCCTGTTTGTGGTATCGTTTATCG‑CTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGc                                                    <  1:1358631/97‑1 (MQ=255)
                                          gTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCTCCTGTTGGTGGTATCGTTTATCG‑CTTTTTTCAAAAATTTCGACACATTTTAAGGGGATTTTCGCaa                                                  >  2:1974799/1‑100 (MQ=255)
                                          gTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCG‑CTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCaa                                                  <  1:2236924/100‑1 (MQ=255)
                                           tttCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCG‑CTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCaa                                                  <  2:1784996/99‑1 (MQ=255)
                                             tCAGCCCCAGGTGCCTTTCTTGAGGCAAGCGCCTGTTGGTGGTATCGTTTATCG‑CTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGc                                               <  2:759833/100‑1 (MQ=255)
                                              cAGCTCCAGGTGCCTTTCTTGAGGCAATCGCCTTTTGGTGTTATCGTTTATCG‑CTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCg                                              <  1:2685994/100‑1 (MQ=255)
                                              cAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTTGTGGTATCGTTTATCG‑CTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCg                                              <  1:1613440/100‑1 (MQ=255)
                                              cAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCG‑CTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCg                                              <  2:1842778/100‑1 (MQ=255)
                                                gCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCG‑CTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTa                                            <  2:365536/100‑1 (MQ=255)
                                                 ccccAGGTGCCTTTCTTGAGGCAATCGCCTGTTGTTGGTATCGTTTATCG‑CTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTAt                                           <  2:1592515/100‑1 (MQ=255)
                                                      ggTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCG‑CTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCAtcc                                      <  1:2042927/100‑3 (MQ=255)
                                                      ggTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCG‑CTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTc                                      <  2:3021746/100‑1 (MQ=255)
                                                         gCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCG‑CTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGc                                   >  2:1591718/1‑100 (MQ=255)
                                                          cttttCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCG‑CTTTTTTCAAAAAATTCGACACAGTTGAAGGGGATTTTCGCAATGCGTATCATTCTGCt                                  <  1:1286534/98‑1 (MQ=255)
                                                          ccTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCG‑CTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCt                                  >  1:1681371/1‑100 (MQ=255)
                                                           cTTTCTTGAGGCAATCGCCTGTTGGTGTTATCGTTTATCT‑CTTTTTTCAAAAAATTCGACACATTTTAAGGGAATTTTCGCAATGCGTATCATTCTGCtt                                 >  1:970508/1‑100 (MQ=255)
                                                              tCTTTAGGCAATCGCCTGTTGGTGGTATCGTTTATCG‑CTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGgc                              <  2:377101/100‑1 (MQ=255)
                                                              tCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCG‑CTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGgc                              <  1:2460354/100‑1 (MQ=255)
                                                              tCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCG‑CTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGgc                              <  2:2223345/100‑1 (MQ=255)
                                                              tCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCG‑CTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATGCTGCTTGgc                              <  1:2359899/100‑1 (MQ=255)
                                                                 tGAGGCAATCGCCTGTTGGTGGTATCGTTTATCG‑CTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATAATTCTGCTTGGCGCt                           >  2:859585/1‑100 (MQ=255)
                                                                    ggCAATCGCCTGTTGGTGGTATCGTTTATCG‑CTTTTTTTAAAATTTTCGACACATTTTACGGGGCTTCTCGCAATGCGTATCATTCTGCCTGCCGCTcc                         >  2:2416554/1‑99 (MQ=255)
                                                                    ggCAATCGCCTGTTGGTGGTATCGTTTATCG‑CTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCg                        >  2:2117184/1‑100 (MQ=255)
                                                                    ggCAATCGCCTGTTGGTGGTATCGTTTATCG‑CTTTTTTCAAAAAATTCGACACATTATAAGGGGATTTTCGCAAAGCGTATCATTCTGCTTGGAGCTCCg                        >  2:2267752/1‑100 (MQ=255)
                                                                      cAATCGCCTGTTGGTGGTATCGTTTATCG‑CTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATTCGTATCATTCTGCTTGgc                              >  2:636158/1‑92 (MQ=255)
                                                                      cAATCGCCTGTTGGTGGTATCGTTTATCG‑CTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCggg                      <  1:2626300/100‑1 (MQ=255)
                                                                      cAATCGCCTGTTGGTGGTATCGTTTATCG‑CTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCggg                      <  1:234766/100‑1 (MQ=255)
                                                                      cAATCGCCTGTTGGTGGTATCGTTTATCG‑CTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCggg                      >  2:765467/1‑100 (MQ=255)
                                                                        aTCGCCTGTTGGTGGTATCGTTTATCG‑CTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCGGgcg                    >  1:2573617/1‑100 (MQ=255)
                                                                         tCGCCTGTTGGTGGTATCGTTTATCG‑CTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCGGgcgc                   <  2:2920343/100‑1 (MQ=255)
                                                                           gCCTGTTGGTGGTATCGTTTATCG‑CTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCGGGCGCgg                 >  2:2152388/1‑100 (MQ=255)
                                                                            tcTGTTGGTGGTATCGTTTATCG‑CTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCGGGCGCggg                <  1:2210708/99‑1 (MQ=255)
                                                                            ccTGTTGGTGGTATCGTTTATCG‑CTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCGGGCGCggg                <  1:1407309/100‑1 (MQ=255)
                                                                             cTGTTGGTGGTATTGTTTATCG‑CTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCGGGCGCgg                 <  1:319049/98‑1 (MQ=255)
                                                                              tGTTGGTGGTATCGTTTATCG‑CTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCGGGCGCGGGGa              <  1:2214443/100‑1 (MQ=255)
                                                                               gTTGGTGGTATCGTTTATCG‑CTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCGGGCGCgggaga             >  2:2210685/1‑97 (MQ=255)
                                                                               gTTGGTGGTATCGTTTATCG‑CTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCGGGCGCGGGGaa             <  1:314246/100‑1 (MQ=255)
                                                                                  ggtggtATCGTTTATCG‑CTTTTTTTATTTTTTTTTACACATTTTAAGGGGATTTTCGCACTGCGTATCATTCTGCTTGGCGCTCTGGGCGCGCGGAAAgg          >  2:1461085/1‑100 (MQ=255)
                                                                                  ggtggtATCGTTTATCG‑CTTTTTTTATTTTTTTTTAAACATTTTAAGGGGATATTAACAATGCGTATCATTCTGCTTGGCGCtttggtc                     >  2:1217382/1‑83 (MQ=255)
                                                                                  ggtggtATCGTTTATCG‑CTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCGGGCGCGGGGAAAgg          <  2:2693446/100‑1 (MQ=255)
                                                                                  ggtggtATCGTTTATCG‑CTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCGGGCGCGGGGAAAgg          >  2:1449859/1‑100 (MQ=255)
                                                                                  ggtggtATCGTTTATCG‑CTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCGCCGGGCGCGGGGAAAgg          <  1:381206/100‑1 (MQ=255)
                                                                                   gtggtATCGTTTATCG‑CTTTTTTTATATAATTCAAAACATATTAAGGGGATTTTCACAATGCGTATCATTCTGCTTGGCGCTCCGGGCGCGGAGAAAggg         >  2:999733/1‑100 (MQ=255)
                                                                                   gtggtATCGTTTATCG‑CTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCGGGCGCGGGGAAAggg         <  2:1681351/100‑1 (MQ=255)
                                                                                       tATCGTTTATCG‑CTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCGGGCGCGGGGAAAGGGACTc     >  1:360645/1‑100 (MQ=255)
                                                                                       tATCGTTTATCG‑CTTTTTTCAAAAAATGCGACACATTATAAGGAGATTTTCGCAATGCGTATCATTCTGCTTGCCGCTCCGGGCGCGGGGAAAGGGACTc     <  1:1341726/100‑1 (MQ=255)
                                                                                         tCGTTTATCG‑CTTTTTTCAAAAAATTCGACACATTTTAAGGGGTTTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCGGGCGCGGGGAAAGGGACTCAg   <  1:1787551/100‑1 (MQ=255)
                                                                                         tCGTTTATCG‑CTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCGGGCGCGGGGAAAGGGACTCAg   <  1:1949122/100‑1 (MQ=255)
                                                                                          cGTTTATCG‑CTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCGGGCGCGGGGAAAGGGACTCAgg  <  1:765471/100‑1 (MQ=255)
                                                                                                     |                                                                                         
TTTCATCATCTGCACTTTCCGCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCG‑C‑TTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCGGGCGCGGGGAAAGGGACTCAGG  >  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: ‑