Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,529,926 G→A 5.7% A51A (GCG→GCA)  hslO → heat shock protein Hsp33

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,529,9260GA 5.7% 232.0 / 2.9 70A51A (GCG→GCA) hslOheat shock protein Hsp33
Reads supporting (aligned to +/- strand):  ref base G (37/29);  new base A (2/2);  total (39/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

AACATGACCAATTACATCGCTATCTGTTTGAAAACTTTGCCGTGCGCGGCGAACTGGTAACCGTTTCGGAAACCCTGCAACAGATCCTTGAGAACCACGATTATCCGCAGCCCGTTAAAAACGTGCTGGCAGAACTGCTGGTTGCGACCAGCCTGTTAACCGCTACGCTGAAGTTTGATGGTGATATCACCGTACAGCTGCAGGGCGACGGTCCGATGAATCTGGCGGTTATTAACGGTAACAATAACCAGCAGATGCGCGGTGTGGCGCGCGTGCAGGGCGAAATTCCAGAAA  >  NC_000913/3529781‑3530074
                                                                                                                                                 |                                                                                                                                                    
aaCATGACCAATTACATCGCTATCTGTTTGAAAACTTTGCCGTGCGCGGCGAACTGGTAACCGTTTCGGAAACCCTGCAACAGATCCTTGAGAACCACGATTATCCGCAGCCCGTTAAAAACGTGCTGGCAGAACTGCTGGTTGCGAcc                                                                                                                                                   >  2:439914/1‑149 (MQ=255)
   aTGACCAATTACATCGCTATCTGTTTGAAAACTTTGCCGTGCGCGGCGAACTGGTAACCGTTTCGGAAACCCTGCAACAGATCCTTGAGAACCACGATTATCCGCAGCCCGTTAAAAACGTGCTGGCAGAACTGCTGGTTGCGACCAGc                                                                                                                                                >  2:143426/1‑149 (MQ=255)
    tGACCAATTACATCGCTATCTGTTTGAAAACTTTGCCGTGCGCGGCGAACTGGTAACCGTTTCGGAAACCCTGCAACAGATCCTTGAGAACCACGATTATCCGCAGCCCGTTAAAAACGTGCTGGCAGAACTGCTGGTTGCGACCAGcc                                                                                                                                               >  1:74196/1‑149 (MQ=255)
     gACCAATTACATCGCTATCTGTTTGAAAACTTTGCCGTGCGCGGCGAACTGGTAACCGTTTCGGAAACCCTGCAACAGATCCTTGAGAACCACGATTATCCGCAGCCCGTTAAAAACGTGCTGGCAGAACTGCTGGTCGCAACCAGcag                                                                                                                                              <  2:289853/149‑3 (MQ=255)
       ccAATTACATCGCTATCTGTTTGAAAACTTTGCCGTGCGCGGCGAACTGGTAACCGTTTCGGAAACCCTGCAACAGATCCTTGAGAACCACGATTATCCGCAGCCCGTTAAAAACGTGCTGGCAGAACTGCTGGTTGCGACCAGCCTGt                                                                                                                                            <  2:443845/149‑1 (MQ=255)
           ttACATCGCTATCTGTTTGAAAACTTTGCCGTGCGCGGCGAACTGGTAACCGTTTCGGAAACCCTGCAACAGATCCTTGAGAACCACGATTATCCGCAGCCCGTTAAAAACGTGCTGGCAGAACTGCTGGTTGCGACCAGCCTGTTAAc                                                                                                                                        >  3:429031/1‑149 (MQ=255)
             aCATCGCTATCTGTTTGAAAACTTTGCCGTGCGCGGCGAACTGGTAACCGTTTCGGAAACCCTGCAACAGATCCTTGAGAACCACGATTATCCGCAGCCCGTTAAAAACGTGCTGGCAGAACTGCTGGTTGCGACCAGCCTGTTAACCg                                                                                                                                      >  2:493319/1‑149 (MQ=255)
             aCATCGCTATCTGTTTGAAAACTTTGCCGTGCGCGGCGAACTGGTAACCGTTTCGGAAACCCTGCAACAGATCCTTGAGAACCACGATTATCCGCAGCCCGTTAAAAACGTGCTGGCAGAACTGCTGGTTGCGACCAGCCTGTTAACCg                                                                                                                                      >  2:325454/1‑149 (MQ=255)
               aTCGCTATCTGTTTGAAAACTTTGCCGTGCGCGGCGAACTGGCAACCGTTTCGGAAACCCTGCAACAGATCCTTGAGAACCACGATTATCCGCAGCCCGTTAAAAACGTGCTGGCAGAACTGCTGGTTGCGACCAGCCTGTTAACCGCt                                                                                                                                    <  1:477386/149‑1 (MQ=255)
                 cGCTATCTGTTTGAAAACTTTGCCGTGCGCGGCGAACTGGTAACCGTTTCGGAAACCCTGCAACAGATCCTTGAGAACCACGATTATCCGCAGCCCGTTAAAAACGTGCTGGCAGAACTGCTGGTTGCGACCAGCCTGTTAACCGCTAc                                                                                                                                  <  2:74196/149‑1 (MQ=255)
                    tATCTGTTTGAAAACTTTGCCGTGCGCGGCGAACTGGTAACCGTTTCGGAAACCCTGCAACAGATCCTTGAGAACCACGATTATCCGCAGCCCGTTAAAAACGTGCTGGCAGAACTGCTGGTTGCGACCAGCCTGTTAACCGCTACGCt                                                                                                                               >  1:194644/1‑149 (MQ=255)
                        tGTTTGAAAACTTTGCCGTGCGCGGCGAACTGGTAACCGTTTCGGAAACCCTGCAACAGATCCTTGAGAACCACGATTATCCGCAGCCCGTTAAAAACGTGCTGGCAGAACTGCTGGTTGCGACCAGCCTGTTAACCGCTACGCTGAAg                                                                                                                           >  1:460465/1‑149 (MQ=255)
                         gTTTGAAAACTTTGCCGTGCGCGGCGAACTGGTAACCGTTTCGGAAACCCTGCAACAGATCCTTGAGAACCACGATTATCCGCAGCCCGTTAAAAACGTGCTGGCAGAACTGCTGGTTGCGACCAGCCTGTTAACCGCTACGCTGAAGt                                                                                                                          >  4:292506/1‑149 (MQ=255)
                          tttGAAAACTTTGCCGTGCGCGGCGAACTGGTAACCGTTTCGGAAACCCTGCAACAGATCCTTGAGAACCACGATTATCCGCAGCCCGTTAAAAACGTGCTGGCAGAACTGCTGGTTGCGACCAGCCTGTTAACCGCTACGCTGAAGtt                                                                                                                         >  4:100180/1‑149 (MQ=255)
                          tttGAAAACTTTGCCGTGCGCGGCGAACTGGTAACCGTTTCGGAAACCCTGCAACAGATCCTTGAGAACCACGATTATCCGCAGCCCGTTAAAAACGTGCTGGCAGAACTGCTGGTTGCGACCAGCCTGTTAACCGCTACGCTGAAGtt                                                                                                                         >  1:368495/1‑149 (MQ=255)
                          tttGAAAACTTTGCCGTGCGCGGCGAACTGGTAACCGTTTCGGAAACCCTGCAACAGATCCTTGAGAACCACGATTATCCGCAGCCCGTTAAAAACGTGCTGGCAGAACTGCTGGTTGCGACCAGCCTGTTAACCGCTACGCTGAAGtt                                                                                                                         <  3:292506/149‑1 (MQ=255)
                              aaaaCTTTGCCGTGCGCGGCGAACTGGTAACCGTTTCGGAAACCCTGCAACAGATCCTTGAGAACCACGATTATCCGCAGCCCGTTAAAAACGTGCTGGCAGAACTGCTGGTTGCGACCAGCCTGTTAACCGCTACGCTGAAGTTTGAt                                                                                                                     >  3:379325/1‑149 (MQ=255)
                                  cTTTGCCGTGCGCGGCGAACTGGTAACCGTTTCGGAAACCCTGCAACAGATCCTTGAGAACCACGATTATCCGCAGCCCGTTAAAAACGTGCTGGCAGAACTGCTGGTTGCGACCAGCCTGTTAACCGCTACGCTGAAGTTTGATGGTg                                                                                                                 >  2:527663/1‑149 (MQ=255)
                                    ttGCCGTGCGCGGCGAACTGGTAACCGTTTCGGAAACCCTGCAACAGATCCTTGAGAACCACGATTATCCGCAGCCCGTTAAAAACGTGCTGGCAGAACTGCTGGTTGCGACCAGCCTGTTAACCGCTACGCTGAAGTTTGATGGTGat                                                                                                               <  3:100180/149‑1 (MQ=255)
                                         gTGCGCGGCGAACTGGTAACCGTTTCGGAAACCCTGCAACAGATCCTTGAGAACCACGATTATCCGCAGCCCGTTAAAAACGTGCTGGCAGAACTGCTGGTTGCGACCAGCCTGTTAACCGCTACGCTGAAGTTTGATGGTGATATCAc                                                                                                          <  2:368495/149‑1 (MQ=255)
                                                  gAACTGGTAACCGTTTCGGAAACCCTGCAACAGATCCTTGAGAACCACGATTATCCGCAGCCCGTTAAAAACGTGCTGGCAGAACTGCTGGTTGCGACCAGCCTGTTAACCGCTACGCTGAAGTTTGATGGTGATATCACCGTACAGCt                                                                                                 <  3:29977/149‑1 (MQ=255)
                                                        gTAACCGTTTCGGAAACCCTGCAACAGATCCTTGAGAACCACGATTATCCGCAGCCCGTTAAAAACGTGCTGGCAGAACTGCTGGTTGCGACCAGCCTGTTAACCGCTACGCTGAAGTTTGATGGTGATATCACCGTACAGCTGCAggg                                                                                           >  2:59445/1‑149 (MQ=255)
                                                          aaCCGTTTCGGAAACCCTGCAACAGATCCTTGAGAACCACGATTATCCGCAGCCCGTTAAAAACGTGCTGGCAGAACTGCTGGTTGCAACCAGCCTGTTAACCGCTACGCTGAAGTTTGATGGTGATATCACCGTACAGCTGCAgg                                                                                            <  4:270188/146‑1 (MQ=255)
                                                          aaCCGTTTCGGAAACCCTGCAACAGATCCTTGAGAACCACGATTATCCGCAGCCCGTTAAAAACGTGCTGGCAGAACTGCTGGTTGCAACCAGCCTGTTAACCGCTACGCTGAAGTTTGATGGTGATATCACCGTACAGCTGCAgg                                                                                            >  3:270188/1‑146 (MQ=255)
                                                             cGTTTCGGAAACCCTGCAACAGATCCTTGAGAACCACGATTATCCGCAGCCCGTTAAAAACGTGCTGGCAGAACTGCTGGTTGCGACCAGCCTGTTAACCGCTACGCTGAAGTTTGATGGTGATATCACCGTACAGCTGCAGGGCGACg                                                                                      >  3:558556/1‑149 (MQ=255)
                                                               tttCGGAAACCCTGCAACAGATCCTTGAGAACCACGATTATCCGCAGCCCGTTAAAAACGTGCTGGCAGAACTGCTGGTTGCGACCAGCCTGTTAACCGCTACGCTGAAGTTTGATGGTGATATCACCGTACAGCTGCAGGGCGACGGt                                                                                    >  3:255546/1‑149 (MQ=255)
                                                               tttCGGAAACCCTGCAACAGATCCTTGAGAACCACGATTATCCGCAGCCCGTTAAAAACGTGCTGGCAGAACTGCTGGTTGCGACCAGCCTGTTAACCGCTACGCTGAAGTTTGATGGTGATATCACCGTACAGCTGCAGGGCGACGGt                                                                                    >  3:456497/1‑149 (MQ=255)
                                                                   ggAAACCCTGCAACAGATCCTTGAGAACCACGATTATCCGCAGCCCGTTAAAAACGTGCTGGCAGAACTGCTGGTTGCGACCAGCCTGTTAACCGCTACGCTGAAGTTTGATGGTGATATCACCGTACAGCTGCAGGGCGACGGTCCGa                                                                                <  1:418177/149‑1 (MQ=255)
                                                                   ggAAACCCTGCAACAGATCCTTGAGAACCACGATTATCCGCAGCCCGTTAAAAACGTGCTGGCAGAACTGCTGGTTGCGACCAGCCTGTTAACCGCTACGCTGAAGTTTGATGGTGATATCACCGTACAGCTGCAGGGCGACGGTCCGa                                                                                >  2:250566/1‑149 (MQ=255)
                                                                    gAAACCCTGCAACAGATCCTTGAGAACCACGATTATCCGCAGCCCGTTAAAAACGTGCTGGCAGAACTGCTGGTTGCGACCAGc                                                                                                                                                >  2:481054/1‑84 (MQ=255)
                                                                    gAAACCCTGCAACAGATCCTTGAGAACCACGATTATCCGCAGCCCGTTAAAAACGTGCTGGCAGAACTGCTGGTTGCGACCAGc                                                                                                                                                <  1:481054/84‑1 (MQ=255)
                                                                    gAAACCCTGCAACAGATCCTTGAGAACCACGATTATCCGCAGCCCGTTAAAAACGTGCTGGCAGAACTGCTGGTTGCGACCAGCCTGTTAACCGCTACGCTGAAGTTTGATGGTGATATCACCGTACAGCTGCAGGGCGACGGTCCGAt                                                                               <  4:439936/149‑1 (MQ=255)
                                                                    gAAACCCTGCAACAGATCCTTGAGAACCACGATTATCCGCAGCCCGTTAAAAACGTGCTGGCAGAACTGCTGGTTGCGACCAGCCTGTTAACCGCTACGCTGAAGTTTGATGGTGATATCACCGTACAGCTGCAGGGCGACGGTCCGAt                                                                               >  1:239985/1‑149 (MQ=255)
                                                                    gAAACCCTGCAACAGATCCTTGAGAACCACGATTATCCGCAGCCCGTTAAAAACGTGCTGGCAGAACTGCTGGTTGCGACCAGCCTGTTAACCGCTACGCTGAAGTTTGATGGTGATATCACCGTACAGCTGCAGGGCGACGGTCCGAt                                                                               >  3:423862/1‑149 (MQ=255)
                                                                    gAAACCCAGCAACAGCTCCTTGAGCACCACGATTATCCGCAGCCCGTTAAAAACGTGCTGTCAGAACTGCTGCTTGAGACCAGCCTGTTAACAACAACGCTGCAGTTGGATGGTGATATCACCGTACAGCTGCAGGGCTCCGGTCCGAt                                                                               >  2:105416/1‑149 (MQ=255)
                                                                     aaaCCCTGCAACAGATCCTTGAGAACCACGATTATCCGCAGCCCGTTAAAAACGTGCTGGCAGAACTGCTGGTTGCGACCAGCCTGTTAACCGCTACGCTGAAGTTTGATGGTGATATCACCGTACAGCTGCAGGGCGACGGTCCGATg                                                                              >  1:131659/1‑149 (MQ=255)
                                                                       aCCCTGCAACAGATCCTTGAGAACCACGATTATCCGCAGCCCGTTAAAAACGTGCTGGCAGAACTGCTGGTTGCGACCAGCCTGTTAACCGCTACGCTGAAGTTTGATGGTGATATCACCGTACAGCTGCAGGGCGACGGTCCGATGaa                                                                            <  4:416251/149‑1 (MQ=255)
                                                                        cccTGCAACAGATCCTTGAGAACCACGATTATCCGCAGCCCGTTAAAAACGTGCTGGCAGAACTGCTGGTTGCGACCAGCCTGTTAACCGCTACGCTGAAGTTTGATGGTGATATCACCGTACAGCTGCAGGGCGACGGTCCGATGAAt                                                                           >  1:399299/1‑149 (MQ=255)
                                                                              aaCAGATCCTTGAGAACCACGATTATCCGCAGCCCGTTAAAAACGTGCTGGCAGAACTGCTGGTTGCGACCAGCCTGTTAACCGCTACGCTGAAGTTTGATGGTGATATCACCGTACAGCTGCAGGGCGACGGTCCGATGAATCt                                                                         >  3:169900/1‑145 (MQ=255)
                                                                              aaCAGATCCTTGAGAACCACGATTATCCGCAGCCCGTTAAAAACGTGCTGGCAGAACTGCTGGTTGCGACCAGCCTGTTAACCGCTACGCTGAAGTTTGATGGTGATATCACCGTACAGCTGCAGGGCGACGGTCCGATGAATCt                                                                         <  4:169900/145‑1 (MQ=255)
                                                                              aaCAGATCCTTGAGAACCACGATTATCCGCAGCCCGTTAAAAACGTGCTGGCAGAACTGCTGGTTGCGACCAGCCTGTTAACCGCTACGCTGAAGTTTGATGGTGATATCACCGTACAGCTGCAGGGCGACGGTCCGATGAATCTGGCg                                                                     >  2:305371/1‑149 (MQ=255)
                                                                                      cTTGAGAACCACGATTATCCGCAGCCCGTTAAAAACGTGCTGGCAGAACTGCTGGTTGCGACCAGCCTGTTAACCGCTACGCTGAAGTTTGATGGTGATATCACCGTACAGCTGCAGGGCGACGGTCCGATGAATCTGGCGGTTATTaa                                                             <  2:399299/149‑1 (MQ=255)
                                                                                      cTTGAGAACCACGATTATCCGCAGCCCGTTAAAAACGTGCTGGCAGAACTGCTGGTTGCGACCAGCCTGTTAACCGCTACGCTGAAGTTTGATGGTGATATCACCGTACAGCTGCAGGGCGACGGTCCGATGAATCTGGCGGTTATTaa                                                             <  4:429031/149‑1 (MQ=255)
                                                                                         gagaACCACGATTATCCGCAGCCCGTTAAAAACGTGCTGGCAGAACTGCTGGTTGCGACCAGCCTGTTAACCGCTACGCTGAAGTTTGATGGTGATATCACCGTACAGCTGCAGGGCGACGGTCCGATGAATCTGGCGGTTATTAACgg                                                          <  1:439914/149‑1 (MQ=255)
                                                                                             aCCACGATTATCCGCAGCCCGTTAAAAACGTGCTGGCAGAACTGCAGGTTGCGACCAGCCTGTTAACCGCTTCGCTGAAGTTTGATGGTGATATCACCGTTCAGCTGCAGGGCGACGGTCCGATGAATCTGGCGGTTATTAACGGTAAc                                                      >  2:22311/1‑149 (MQ=255)
                                                                                               cACGATTATCCGCAGCCCGTTAAAAACGTGCTGGCAGAACTGCTGGTTGCGACCAGCCTGTTAACCGCTACGCTGAAGTTTGATGGTGATATCACCGTACAGCTGCAGGGCGACGGTCCGATGAATCTGGCGGTTATTAACGGTAACaa                                                    <  2:376420/149‑1 (MQ=255)
                                                                                                aCGATTATCCGCAGCCCGTTAAAAACGTGCTGGCAGAACTTCTGGTTGCGACCAGCCTGTTAACCGCTACGCTGAAGTTTGATGGTGATATCACCGTACAGCTGCAGGGCGACGGTCCGATGAATCTGGCGGTTATTAACGGTAACAAt                                                   >  1:469100/1‑149 (MQ=255)
                                                                                                aCGATTATCCGCAGCCCGTTAAAAACGTGCTGGCAGAACTGCTGGTTGCGACCAGCCTGTTAACCGCTACGCTGAAGTTTGATGGTGATATCACCGTACAGCTGCAGGGCGACGGTCCGATGAATCTGGCGGTTATTAACGGTAACAAt                                                   <  4:498129/149‑1 (MQ=255)
                                                                                                aCGATTACCCGCAGCCCGTTAAAAACGTGCTGGCAGAACTGCTGGTTGTGACCAGCCTGTTAACCGCTACGCTGAATTTTGATGGTGATATCACCGTACAGCTGCAGGGCGACGGTCCGATGAATCTGGCGGTTATTAACGGTAACAAt                                                   <  4:29706/149‑1 (MQ=255)
                                                                                                  gATTATCCGCAGCCCGTTAAAAACGTGCTGGCAGAACTGCTGGTTGCGACCAGCCTGTTAACCGCTACGCTGAAGTTTGATGGTGATATCACCGTACAGCTGCAGGGCGACGGTCCGATGAATCTGGCGGTTATTAACGGTAACAATaa                                                 >  2:263803/1‑149 (MQ=255)
                                                                                                  gATTATCCGCAGCCCGTTAAAAACGTGCTGGCAGAACTGCTGGTTGCGACCAGCCTGTTAACCGCTACGCTGAAGTTTGATGGTGATATCACCGTACAGCTGCAGGGCGACGGTCCGATGAATCTGGCGGTTATTAACGGTAACAATaa                                                 <  3:368453/149‑1 (MQ=255)
                                                                                                       tCCGCAGCCCGTTAAAAACGTGCTGGCAGAACTGCTGGTTGCGACCAGCCTGTTAACCGCTACGCTGAAGTTTGATGGTGATATCACCGTACAGCTGCAGGGCGACGGTCCGATGAATCTGGCGGTTATTAACGGTAACAATAACcagc                                            <  1:250566/149‑1 (MQ=255)
                                                                                                       tCCGCAGCCCGTTAAAAACGTGCTGGCAGAACTGCTGGTTGCGACCAGCCTGTTAACCGCTACGCTGAAGTTTGATGGTGATATCACCGTACAGCTGCAGGGCGACGGTCCGATGAATCTGGCGGTTATTAACGGTAACAATAACcagc                                            >  1:543061/1‑149 (MQ=255)
                                                                                                       tCCGCAGCCCGTTAAAAACGTGCTGGCAGAACTGCTGGTTGCGACCAGCCTGTTAACCGCTACGCTGAAGTTTGATGGTGATATCACCGTACAGCTGCAGGGCGACGGTCCGATGAATCTGGCGGTTATTAACGGTAACAATAACcagc                                            >  4:464443/1‑149 (MQ=255)
                                                                                                              cccGTTAAAAACGTGCTGGCAGAACTGCTGGTTGCGACCAGCCTGTTAACCGCTACGCTGAAGTTTGATGGTGATATCACCGTACAGCTGCAGGGCGACGGTCCGATGAATCTGGCGGTTATTAACGGTAACAATAACCAGCAGATgcg                                     >  3:247896/1‑149 (MQ=255)
                                                                                                                  ttAAAAACGTGCTGGCAGAACTGCTGGTTGCGACCAGCCTGTTAACCGCTACGCTGAAGTTTGATGGTGATATCACCGTACAGCTGCAGGGCGACGGTCCGATGAATCTGGCGGTTATTAACGGTAACAATAACCAGCAGATGCGCGgt                                 <  1:59445/149‑1 (MQ=255)
                                                                                                                      aaaCGTGCTGGCAGAACTGCTGGTTGCGACCAGCCTGTTAACCGCTACGCTGAAGTTTGATGGTGATATCACCGTACAGCTGCAGGGCGACGGTCCGATGAAt                                                                           >  2:200259/1‑103 (MQ=255)
                                                                                                                      aaaCGTGCTGGCAGAACTGCTGGTTGCGACCAGCCTGTTAACCGCTACGCTGAAGTTTGATGGTGATATCACCGTACAGCTGCAGGGCGACGGTCCGATGAAt                                                                           <  1:200259/103‑1 (MQ=255)
                                                                                                                         cGTGCTGGCAGAACTGCTGGTTGCGACCAGCCTGTTAACCGCTACGCTGAAGTTTGATGGTGATATCACCGTACAGCTGCAGGGCGACGGTCCGATGAATCTGGCGGTTATTAACGGTAACAATAACCAGCAGATGCGCGGTGTGgcgc                          <  4:423862/149‑1 (MQ=255)
                                                                                                                         cGTGCTGGCAGAACTGCTGGTTGCGACCAGCCTGTTAACCGCTACGCTGAAGTTTGATGGTGATATCACCATACAACTGCAGGGCGACGGTCCGATGAATCTGGCGGTTATTAACGGTAACAATAACCAGCAGATGCGCGGTGTGgcgc                          <  2:239985/149‑1 (MQ=255)
                                                                                                                             cTGGCAGAACTGCTGGTTGCGACCAGCCTGTTAACCGCTACGCTGAAGTTTGATGGTGATATCACCGTACAGCTGCAGGGCGACGGTCCGATGAATCTGGCGGTTATTAACGGTAACAATAACCAGCAGATGCGCGGTGTGGCGCGCGt                      <  1:258020/149‑1 (MQ=255)
                                                                                                                                  aGAACTGCTGGTTGCGACCAGCCTGTTAACCGCTACGCTGAAGTTTGATGGTGATATCACCGTACAGCTGCAGGGCGACGGTCCGATGAATCTGGCGGTTATTAACGGTAACAATAACCAGCAGATGCGCGGTGTGGTGCGCGTGCAgg                 <  1:304327/149‑1 (MQ=255)
                                                                                                                                  aGAACTGCTGGTTGCGACCAGCCTGTTAACCGCTACGCTGAAGTTTGATGGTGATATCACCGTACAGCTGCAGGGCGACGGTCCGATGAATCTGGCGGTTATTAACGGTAACAATAACCAGCAGATGCGCGGTGTGGCGCGCGTGCAgg                 >  2:394348/1‑149 (MQ=255)
                                                                                                                                  aGAACTGCTGGTTGCGACCAGCCTGTTAACCGCTACGCTGAAGTTTGATGGTGATATCACCGTACAGCTGCAGGGCGACGGTCCGATGAATCTGGCGGTTATTAACGGTAACAATAACCAGCAGATGCGCGGTGTGGCGCGCGTGCAgg                 >  4:151967/1‑149 (MQ=255)
                                                                                                                                  aGAACTGCTGGTTGCGACCAGCCTGTTAACCGCTACGCTGAAGTTTGATGGTGATATCACCGTACAGCTGCAGGGCGACGGTCCGATGAATCTGGCGGTTATTAACGGTAACAATAACCAGCAGATGCGCGGTGTGGCGCGCGTGCAgg                 >  2:254686/1‑149 (MQ=255)
                                                                                                                                      aggctgGTCGCAACCAGCCTGTTAACCGCTACGCTGAAGTTTGATGGTGATATCACCGTACAGCTGCAGGGCGACGGTCCGATGAATCTGGCGGTTATTAACGGTAACAATAACCAGCAGATGCGCGGTGTGGCGCGCGTGCAGGGCGa             >  3:305778/3‑149 (MQ=255)
                                                                                                                                        gctgGTTGCGACCAGCCTGTTAACCGCTACGCTGAAGTTTGATGGTGATATCACCGTACAGCTGCAGGGCGACGGTCCGATGAATCTGGCGGTTATTAACGGTAACAATAACCAGCAGATGCGCGGTGTGGCGCGCGTGCAGGGCGaaa           <  1:493319/149‑1 (MQ=255)
                                                                                                                                        gctgGTTGCGACCAGCCTGTTAACCGCTACGCTGAAGTTTGATGGTGATATCACCGTACAGCTGCAGGGCGACGGTCCGATGAATCTGGCGGTTATTAACGGTAACAATAACCAGCAGATGCGCGGTGTGGCGCGCGTGCAGGGCGaaa           <  1:325454/149‑1 (MQ=255)
                                                                                                                                          tgGTTGCGACCAGCCTGTTAACCGCTACGCGGAAGTTTGTTGGTGATATCACCGTACGGCTGCAGGGCGATGGTCCGATGAATCTGGGGGTTATTAACGGTAACAATAACCAGCAGATGCGCGGTGTGGCGCGCGTGCAGGACGAAAtt         >  2:381448/1‑149 (MQ=255)
                                                                                                                                              tGCGACCAGCCTGTTAACCGCTACGCTGAAGTTTGATGGTGATATCACCGTACAGCTGCAGGGCGACGGTCCGATGAATCTGGCGGTTATTAACGGTAACAATAACCAGCAGATGCGCGGTGTGGCGCGCGTGCAGGGCGAAATTCCAg     <  4:558556/149‑1 (MQ=255)
                                                                                                                                                 gACCAGCCTGTTAACCGCTACGCTGAAGTTTGATGGTGATATCACCGTACAGCTGCAGGGCGACGGTCCGATGAATCTGGCGGTTATTAACGGTAACAATAACCAGCAGATGCGCGGTGTGGCGCGCGTGCAGGGCGAAATTCCAGaaa  >  2:364066/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
AACATGACCAATTACATCGCTATCTGTTTGAAAACTTTGCCGTGCGCGGCGAACTGGTAACCGTTTCGGAAACCCTGCAACAGATCCTTGAGAACCACGATTATCCGCAGCCCGTTAAAAACGTGCTGGCAGAACTGCTGGTTGCGACCAGCCTGTTAACCGCTACGCTGAAGTTTGATGGTGATATCACCGTACAGCTGCAGGGCGACGGTCCGATGAATCTGGCGGTTATTAACGGTAACAATAACCAGCAGATGCGCGGTGTGGCGCGCGTGCAGGGCGAAATTCCAGAAA  >  NC_000913/3529781‑3530074

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑