Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 497,139 (T)5→6 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.T89.9% 256.4 / 22.2 79intergenic (+140/‑41)htpG/adkprotein refolding molecular co‑chaperone Hsp90, Hsp70‑dependent; heat‑shock protein; ATPase/adenylate kinase
Reads supporting (aligned to +/- strand):  ref base . (3/5);  new base T (34/37);  total (37/42)
Fisher's exact test for biased strand distribution p-value = 7.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTCATCATCTGCACTTTCCGCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCGGGCGCGGGGAAAGGGACTCAGGCTCAG  >  NC_000913/497035‑497228
                                                                                                    |                                                                                                 
tttcATCATCTGCACTTTCCGCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGc                                                                                                    >  1:1775181/1‑100 (MQ=255)
tttcATCATCTGCACTTTCCGCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGc                                                                                                    <  2:1327863/100‑1 (MQ=255)
  tcatcaTCTGCACTTTCCGCAAATTATCTCGCCATTATCCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCtt                                                                                                <  1:222851/100‑1 (MQ=255)
   catcatCTGCACTTTCCGCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCttt                                                                                              >  2:2280793/1‑100 (MQ=255)
    atcatcTGCACTTTCCGCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCtttt                                                                                             <  2:396887/100‑1 (MQ=255)
      catcTGCACTTTCCGCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCtttttt                                                                                           >  1:2124362/1‑99 (MQ=255)
      catcTGCACTTTCCGCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCtttttt                                                                                           >  2:1804708/1‑99 (MQ=255)
      catcTGCACTTTCCGCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCtttttt                                                                                           <  1:1264571/100‑2 (MQ=255)
      catcTGCACTTTCCGCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCtttttt                                                                                           >  2:855233/1‑99 (MQ=255)
       atcTGCACTTTCCGCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCttttttc                                                                                          <  2:15054/100‑3 (MQ=255)
        tcTGCACTTTCCGCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCttttttca                                                                                         >  2:2274073/1‑97 (MQ=255)
          tGCACTTTCCGCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCttttttcaaa                                                                                       >  2:187978/1‑95 (MQ=255)
          tGCACTTTCCGCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCttttttcaaa                                                                                       >  1:1082546/1‑95 (MQ=255)
          tGCACTTACCGCAAATTATCTCGCCATTAACCGTTTCAGCCTCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTATCTTCaaa                                                                                       >  2:1278143/1‑100 (MQ=255)
           gCACTTTCCGCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCaaaa                                                                                      >  2:1214220/1‑100 (MQ=255)
             aCTTTCCGCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTCAAAAAAt                                                                                   >  1:2020819/1‑100 (MQ=255)
              cTTTCCGCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAAt                                                                                   >  2:717431/1‑100 (MQ=255)
              cTTTCCGCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAAt                                                                                   <  1:484979/100‑1 (MQ=255)
                ttCCGCAAATTATCTTGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTc                                                                                 <  1:2302102/100‑1 (MQ=255)
                 tCCGCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCg                                                                                <  1:717432/100‑1 (MQ=255)
                  ccGCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTTAAAAAATTCGa                                                                               >  2:995580/1‑100 (MQ=255)
                   cGCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGac                                                                              >  2:126782/1‑100 (MQ=255)
                   cGCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGac                                                                              <  1:1436245/100‑1 (MQ=255)
                   cGCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGAGGTATCGTTTATCGCTTTTTTCAAAAAATTCGac                                                                              >  2:1236162/1‑100 (MQ=255)
                    gCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGaca                                                                             <  1:126784/100‑1 (MQ=255)
                    gCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGaca                                                                             >  1:1910560/1‑100 (MQ=255)
                    gCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGATTTTTAAAAAAAATTCGac                                                                              <  2:2124355/99‑1 (MQ=255)
                        aTTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACAtt                                                                         <  1:1718466/100‑1 (MQ=255)
                        aTTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACAtt                                                                         <  1:34764/100‑1 (MQ=255)
                         ttATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACAttt                                                                        <  2:2212596/100‑1 (MQ=255)
                         ttATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTATTTCAAAAAATTCGACACAttt                                                                        <  2:695076/100‑1 (MQ=255)
                          tATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACAtttt                                                                       <  2:598049/100‑1 (MQ=255)
                           aTCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTa                                                                      <  2:779512/100‑1 (MQ=255)
                           aTCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTa                                                                      <  2:533704/100‑1 (MQ=255)
                           aTCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTa                                                                      >  2:1088818/1‑100 (MQ=255)
                            tctcGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAACTTCGACACAttt                                                                        >  2:2338897/1‑97 (MQ=255)
                             ctcGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAg                                                                    >  1:508164/1‑100 (MQ=255)
                               cGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAggg                                                                  <  1:1214223/100‑1 (MQ=255)
                                gCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAgggn                                                                 <  2:1910553/100‑2 (MQ=255)
                                 ccATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCCGTTGGTGGTAACGTTTATCGCTTTTTTCAAAACATTCGACACATTTTAAGGGGa                                                                >  2:67818/1‑100 (MQ=255)
                                  cATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGAt                                                               <  2:1775176/100‑1 (MQ=255)
                                    ttAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGAttt                                                             <  2:508162/100‑1 (MQ=255)
                                     tAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGAtttt                                                            >  1:348123/1‑100 (MQ=255)
                                     tAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGAtttt                                                            >  1:1523982/1‑100 (MQ=255)
                                      aaCCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACTTTTTAAGGGGATTTTc                                                           >  2:2233586/1‑100 (MQ=255)
                                      aaCCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTc                                                           <  1:855235/100‑1 (MQ=255)
                                      aaCCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTc                                                           <  1:2280801/100‑1 (MQ=255)
                                        ccGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGc                                                         >  2:820499/1‑100 (MQ=255)
                                        ccGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGc                                                         <  1:242027/100‑1 (MQ=255)
                                         cGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCa                                                        <  1:2006138/100‑1 (MQ=255)
                                          gTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCaa                                                       >  1:1943813/1‑100 (MQ=255)
                                            ttCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATg                                                     <  1:1180300/100‑1 (MQ=255)
                                             tCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGc                                                    <  2:1674899/100‑1 (MQ=255)
                                              cAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCg                                                   <  1:1088820/100‑1 (MQ=255)
                                                gCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTa                                                 >  2:1451645/1‑100 (MQ=255)
                                                 ccccAGTTCCCTTTTTTGAGGAAATCGCCTGTTGTTGTTATCTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCg                                                   <  2:258727/96‑1 (MQ=255)
                                                  cccAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATc                                               <  1:1804713/100‑1 (MQ=255)
                                                  cccAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATc                                               >  2:2143893/1‑100 (MQ=255)
                                                       gTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCt                                          <  1:2182181/100‑1 (MQ=255)
                                                            tttCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTg                                     <  1:216510/100‑1 (MQ=255)
                                                             ttCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTgg                                    <  1:2167288/100‑1 (MQ=255)
                                                              cgatctggCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGgcg                                  <  1:1848076/94‑1 (MQ=255)
                                                                 tGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCt                                >  1:818819/1‑100 (MQ=255)
                                                                    ggCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCg                             >  1:15411/1‑100 (MQ=255)
                                                                    ggCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCg                             >  1:1722635/1‑100 (MQ=255)
                                                                    ggCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTCGGCGCTCCg                             <  1:301597/100‑1 (MQ=255)
                                                                    ggCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGgcgagatcg                            >  2:1848071/1‑94 (MQ=255)
                                                                       aaTCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCGGgc                          <  1:224396/100‑1 (MQ=255)
                                                                        aTCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTTTTTTTTTTTTACACATTTTAAGGGGATTTTCTCAATGCGTATCATTCTGCTTAGCGCTCCGGgc                          >  2:396451/1‑99 (MQ=255)
                                                                        aTCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCTTTTTTTTTTACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCGGgcg                         >  2:1804629/1‑100 (MQ=255)
                                                                        aTCGCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCGGgc                          >  2:582903/1‑99 (MQ=255)
                                                                           gCCTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCGGGCGCgg                      <  1:1545201/100‑1 (MQ=255)
                                                                           gCCTGTTGGTGGTATCGTTTATCGCTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCGGGCGCggg                     <  2:2020812/100‑1 (MQ=255)
                                                                            ccTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCGGGCGCggg                     <  2:15410/100‑1 (MQ=255)
                                                                            ccTGTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCGGGCGCggg                     <  1:1723946/100‑1 (MQ=255)
                                                                               gTTGGTGGTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCGGGCGCGGGGaa                  >  2:1421027/1‑100 (MQ=255)
                                                                                 tggtggTATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCGGGCGCAGGGAAAg                <  1:2079326/100‑1 (MQ=255)
                                                                                  tctggtATCGTTTATCGCTTTTTTCAAAAAATTCGACACATGTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCGGGCGCGGGGAAAgg               <  2:853073/98‑1 (MQ=255)
                                                                                    tggtATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCGGGCGCGGGGAAAGGGa             >  2:1838635/1‑100 (MQ=255)
                                                                                     ggtATCGTTTATCGCTTTTTTTTTTTTTTTTGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCGGGCGCGGGGAAAGGGa             >  2:1316435/1‑99 (MQ=255)
                                                                                     ggtATCGTTTATCGCTTTTTTTTTTTTTTTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGTGCTCCGGGCGCGGGGAAAGGGAc            >  2:1202233/1‑100 (MQ=255)
                                                                                     ggtATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCGGGCGCGGGGAAAggag             >  1:853075/1‑97 (MQ=255)
                                                                                     ggtATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCGGGCGCGGGGAAAGGGAc            <  1:1236165/100‑1 (MQ=255)
                                                                                     ggtATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCGGGCGCGGGGAAAGGGAc            <  2:1082544/100‑1 (MQ=255)
                                                                                     ggtATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCGGGCGCGGGGAAAGGGAc            >  2:1956626/1‑100 (MQ=255)
                                                                                      gtATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCGGGCGCGGGGAAAGGGACt           >  2:1332788/1‑100 (MQ=255)
                                                                                      gtATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCGGGCGCGGGGAAAGGGACt           >  2:1334695/1‑100 (MQ=255)
                                                                                      gtATCGTTTATCGCTTTTTTCAAAAAATTCGACACATCTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCGGGCGCGGGGAAAGGGACt           >  2:791380/1‑100 (MQ=255)
                                                                                       tATCGTTTATCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCGGGCGCGGGGAAAGGGACTc          >  2:1411417/1‑100 (MQ=255)
                                                                                        aTCGTTTATCGCTTTTTTCAAAAAATTCGCCACATTTTAAGGGGATTTTCGCAACGCGTATCATTCTGCTTGGCGCTCCGGGCGCGGGGAACGGGACTc          >  2:458274/1‑99 (MQ=255)
                                                                                            tttATCGCTTATTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCGGGCGCGGGGAAAGGGACTCAGGc      <  1:1804634/99‑1 (MQ=255)
                                                                                               aTCGCTTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCGGGCGCGGGGAAAGGGACTCAGGCTCAg  <  1:187980/100‑1 (MQ=255)
                                                                                                    |                                                                                                 
TTTCATCATCTGCACTTTCCGCAAATTATCTCGCCATTAACCGTTTCAGCCCCAGGTGCCTTTCTTGAGGCAATCGCCTGTTGGTGGTATCGTTTATCGCTTTTTCAAAAAATTCGACACATTTTAAGGGGATTTTCGCAATGCGTATCATTCTGCTTGGCGCTCCGGGCGCGGGGAAAGGGACTCAGGCTCAG  >  NC_000913/497035‑497228

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