Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,339,624 A→G K348R (AAG→AGG)  yciM → putative heat shock protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,339,6240AG99.6% 762.1 / ‑5.4 237K348R (AAG→AGG) yciMputative heat shock protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (126/110);  total (126/111)
Fisher's exact test for biased strand distribution p-value = 4.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

Only 100 of 255 total aligned reads displayed.

TCCATAAGTTAATGGATTACCACTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAAGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCCGTCTTGTCGGGCCTGGT  >  NC_000913/1339534‑1339713
                                                                                          |                                                                                         
tCCATAAGTTAATGGATTACCACTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAggg                                                                                         >  2:1268603/1‑93 (MQ=255)
tCCATAAGTTAATGGATTACCACTTAAATGAAGCGAAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAggg                                                                                         >  1:4676583/1‑93 (MQ=255)
 ccATAAGTAAAAGGATTTCTAAATAAATGAAGCGGAAGAAGGGCCTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGt                                                                                        <  2:7495474/93‑1 (MQ=255)
   aTAAGTTAATGGATTACCACTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTag                                                                                      >  2:2897106/1‑92 (MQ=255)
   aTAAGTTAATGGATTACCACTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTAc                                                                                      >  2:5928262/1‑93 (MQ=255)
    tAAGTTAATGGATTACCACTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACg                                                                                     <  2:4643214/93‑1 (MQ=255)
    tAAGTTAATGGATTACCACTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACg                                                                                     >  1:9522307/1‑93 (MQ=255)
     aaGTTAATGGATTACCACTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACgt                                                                                    <  2:7179762/93‑1 (MQ=255)
      aGTTAATGGATTAGCACTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACgta                                                                                   <  2:9499177/93‑1 (MQ=255)
      aGTTAATGGATTACCACTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACgta                                                                                   >  1:8007477/1‑93 (MQ=255)
      aGTTAATGGATTACCACTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACgta                                                                                   >  1:5672639/1‑93 (MQ=255)
      aGTTAATGGATTACCAATTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACTTGGTTGGCGAGAGAGTACgta                                                                                   <  2:4725084/93‑1 (MQ=255)
       gTTAATGGATTACCACTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCTTGACATGGTTGGCGAGAGGGTACgtag                                                                                  >  1:8663832/1‑93 (MQ=255)
       gTTAATGGATGACCACTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACgtag                                                                                  <  2:8897858/93‑1 (MQ=255)
        ttAATGGATTACCACTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACgtagt                                                                                 >  2:2078059/1‑93 (MQ=255)
         tAATGGATTACCACTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACgtagta                                                                                >  1:4186664/1‑93 (MQ=255)
         tAATGGATTACCACTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACgtagta                                                                                <  2:5798301/93‑1 (MQ=255)
          aaTGGATTACCACTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTaa                                                                               <  2:1051614/93‑1 (MQ=255)
          aaTGGATTACCACTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTaa                                                                               >  1:5642019/1‑93 (MQ=255)
           aTGGATTACCACTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAg                                                                              <  1:7935396/93‑1 (MQ=255)
            tGGATTACCACTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGc                                                                             >  1:3701419/1‑93 (MQ=255)
            tGGATTACCACTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGc                                                                             >  2:4568984/1‑93 (MQ=255)
             ggATTACCACTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGcc                                                                            >  1:2123337/1‑93 (MQ=255)
              gATTACCACTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCt                                                                           >  1:9115322/1‑93 (MQ=255)
              gATTACCACTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCt                                                                           >  2:6083611/1‑93 (MQ=255)
               aTTACCACTTAAATGAAGCGGAAGAAGGGGGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGCAAGCCTc                                                                          >  2:3152252/1‑93 (MQ=255)
               aTTACCACTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTc                                                                          >  1:2929274/1‑93 (MQ=255)
                 tACCACTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGt                                                                        <  2:8196423/93‑1 (MQ=255)
                   ccACTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTa                                                                      <  2:9611899/93‑1 (MQ=255)
                    cACTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTAt                                                                     <  1:10015857/93‑1 (MQ=255)
                    cACTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTAt                                                                     >  1:2229348/1‑93 (MQ=255)
                     aCTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATc                                                                    >  2:4870524/1‑93 (MQ=255)
                     aCTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATc                                                                    >  2:6238052/1‑93 (MQ=255)
                       ttAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGt                                                                  <  2:7885060/93‑1 (MQ=255)
                       ttAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGt                                                                  >  2:311556/1‑93 (MQ=255)
                        tAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGtt                                                                 >  2:6800373/1‑93 (MQ=255)
                          aaTGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGc                                                               <  2:5632160/93‑1 (MQ=255)
                           aTGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGcc                                                              >  1:1404193/1‑93 (MQ=255)
                           aTGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGcc                                                              >  1:730312/1‑93 (MQ=255)
                            tGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCa                                                             <  1:6586016/93‑1 (MQ=255)
                              aaGCGGAAGAACGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGa                                                           <  1:4021383/93‑1 (MQ=255)
                                gCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTTGCGAGAGGGCACGTAGTAAGCCTCGTTATCGTTGCCAGaaa                                                         >  2:6961042/1‑93 (MQ=255)
                                gCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGaaa                                                         <  1:2537563/93‑1 (MQ=255)
                                gCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGaaa                                                         <  1:3582427/93‑1 (MQ=255)
                                 cGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAAt                                                        <  2:9035263/93‑1 (MQ=255)
                                 cGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAAt                                                        <  2:3973778/93‑1 (MQ=255)
                                    aagaagGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAAtgtg                                                     <  2:4278358/93‑1 (MQ=255)
                                       aagGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGtt                                                  <  2:3377481/93‑1 (MQ=255)
                                        agGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGttt                                                 <  2:5337065/93‑1 (MQ=255)
                                           gCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTAc                                              <  1:9594862/93‑1 (MQ=255)
                                           gCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTACCAGAAATGTGGTTTTAc                                              <  2:4090984/93‑1 (MQ=255)
                                             gTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCg                                            >  1:5213543/1‑93 (MQ=255)
                                               gCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCa                                          <  1:1018243/93‑1 (MQ=255)
                                               gCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCa                                          <  1:506470/93‑1 (MQ=255)
                                                 cAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATa                                        <  2:4846519/93‑1 (MQ=255)
                                                 cAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATa                                        <  1:6914276/93‑1 (MQ=255)
                                                  aaaGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATac                                       >  2:3727615/1‑93 (MQ=255)
                                                  aaaGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATac                                       <  2:5736431/93‑1 (MQ=255)
                                                   aaGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATaca                                      >  1:6261732/1‑93 (MQ=255)
                                                    agagagCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCATAAATGTGGTTTTACCGCATacac                                     >  2:7226448/1‑93 (MQ=255)
                                                     gagagCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACAcc                                    >  1:5332775/1‑93 (MQ=255)
                                                     gagagCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACAcc                                    <  1:8161824/93‑1 (MQ=255)
                                                       gagCCTGGTGGTCCTGGGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCct                                  <  1:7935941/93‑1 (MQ=255)
                                                       gagCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACAACct                                  >  2:2796501/1‑93 (MQ=255)
                                                        agCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCctc                                 <  1:5466179/93‑1 (MQ=255)
                                                         gCCTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCctct                                <  2:5218636/93‑1 (MQ=255)
                                                          ccTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTa                               <  1:4506784/93‑1 (MQ=255)
                                                          ccTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTa                               >  2:8814373/1‑93 (MQ=255)
                                                          ccTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTa                               >  1:3170867/1‑93 (MQ=255)
                                                           cTGATGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTAc                              <  1:7722542/93‑1 (MQ=255)
                                                            tGATGGTGCTGCGTGACATGGTTGGCGAGAAGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACt                             <  2:8361904/93‑1 (MQ=255)
                                                              aTGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTgg                           >  2:8079067/1‑93 (MQ=255)
                                                              aTGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTgg                           >  2:7836367/1‑93 (MQ=255)
                                                              aTGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTgg                           >  1:9955202/1‑93 (MQ=255)
                                                              aTGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTgg                           >  2:1522202/1‑93 (MQ=255)
                                                              aTGGTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTgc                           >  1:287444/1‑92 (MQ=255)
                                                                ggTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCa                         >  2:7333676/1‑93 (MQ=255)
                                                                 gTGCTGCGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCAt                        >  1:2189008/1‑93 (MQ=255)
                                                                    ctgcGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGt                     <  1:73771/93‑1 (MQ=255)
                                                                     tgcGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTc                    >  1:937366/1‑93 (MQ=255)
                                                                      gcGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTCCGGCATACCCCCTCTACTGCCATTGTcc                   >  1:8267552/1‑93 (MQ=255)
                                                                      gcGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTcc                   <  2:5708517/93‑1 (MQ=255)
                                                                      gcGTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTcc                   >  1:7946799/1‑93 (MQ=255)
                                                                        gTGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCCGt                 >  1:8374950/1‑93 (MQ=255)
                                                                         tGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCCGTc                >  1:6029812/1‑93 (MQ=255)
                                                                         tGACATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCCGTc                >  2:3682270/1‑93 (MQ=255)
                                                                           aCATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCCGTCtt              >  1:7434178/1‑93 (MQ=255)
                                                                            cATGGTTGGCGAGAGGGTACGTCGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCCGTCTTg             <  1:8923973/93‑1 (MQ=255)
                                                                            cATGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCCGTCTTg             <  1:6119979/93‑1 (MQ=255)
                                                                             aTGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCCGTCTTGt            <  2:8781995/93‑1 (MQ=255)
                                                                              tGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCCGTCTTGTc           >  1:236469/1‑93 (MQ=255)
                                                                              tGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCCGTCTTGTc           <  2:9757159/93‑1 (MQ=255)
                                                                              tGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCCGTCTTGTc           >  2:7590200/1‑93 (MQ=255)
                                                                              tGGTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCCGTCTTGTc           <  2:431015/93‑1 (MQ=255)
                                                                               ggTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCCGTCTTGTCg          <  1:7016657/93‑1 (MQ=255)
                                                                               ggTTGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCCGTCTTGTCg          <  1:8458931/93‑1 (MQ=255)
                                                                               ggTGGGCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATGCACCCTCTACTGGCATTGTCCGTCTTGTCg          <  1:5887871/93‑1 (MQ=255)
                                                                                 ttGGCGAGAGGTTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCCGTCTTGTCggg        >  1:9738971/1‑93 (MQ=255)
                                                                                   ggCGAGAGGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCCGTCTTGTCGGGcc      >  1:4913958/1‑93 (MQ=255)
                                                                                       agagGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCCGTCTTGTCGGGCCTGGt  >  2:7129608/1‑93 (MQ=255)
                                                                                          |                                                                                         
TCCATAAGTTAATGGATTACCACTTAAATGAAGCGGAAGAAGGGCGTGCCAAAGAGAGCCTGATGGTGCTGCGTGACATGGTTGGCGAGAAGGTACGTAGTAAGCCTCGTTATCGTTGCCAGAAATGTGGTTTTACCGCATACACCCTCTACTGGCATTGTCCGTCTTGTCGGGCCTGGT  >  NC_000913/1339534‑1339713

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