Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 302,516 C→A 8.5% D20Y (GAT→TAT)  paoA ← PaoABC aldehyde oxidoreductase, 2Fe‑2S subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913302,5160CA 8.5% 288.4 / 2.8 104D20Y (GAT→TAT) paoAPaoABC aldehyde oxidoreductase, 2Fe‑2S subunit
Reads supporting (aligned to +/- strand):  ref base C (39/53);  new base A (9/2);  total (49/55)
Fisher's exact test for biased strand distribution p-value = 2.20e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 105 total aligned reads displayed.

CCGTTCACCTTCAGTGTCAGGGGCATTATCTCTGGCGCGGGTGTAGCTGCCGGAACGCTTGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAACCCGATTATCTTCGGGGTATTCGCCTTGGTTGCTCATACCAGGCCTCCGGTATTCTAAAGGGGAAAATAAGAGTGCGTTGCGAGAGTGCTGATGCTTTTTTAAACGTTAAGCATAGTCGGCGACAGCAG  >  NC_000913/302368‑302665
                                                                                                                                                    |                                                                                                                                                     
ccGTTCACCTTCAGTGTCAGGGGCATTATCTCTGGCGCGGGTGTAGCTGCCGGAACGCTTGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCg                                                                                                                                                      <  2:407718/150‑1 (MQ=255)
ccGTTCACCTTCAGTGTCAGGGGCATTATCTCTGGCGCGGGTGTAGCTGCCGGAACGCTTGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCg                                                                                                                                                      >  2:1103811/1‑150 (MQ=255)
 cGTTCACCTTCAGTGTCAGGGGCATTATCTCTGGCGCGGGTGTAGCTGCCGGAACGCTTGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCGGTTGCGGCGCGCCCTTTTAAACGATAGCGCCGGGGCAAACTTTAAACGt                                                                                                                                                     >  2:1035080/1‑150 (MQ=255)
  gTTCACCTTCAGTGTCAGGGGCATTATCTCTGGCGCGGGTGTAGCTGCCGGAACGCTTGCCGCCAGCGTAGAATTAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTg                                                                                                                                                    >  1:29318/1‑150 (MQ=255)
  gTTCACCTTCAGTGTCAGGGGCATTATCTCTGGCGCGGGTGTAGCTGCCGGAACGCTTGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTg                                                                                                                                                    <  1:1103811/150‑1 (MQ=255)
   ttCACCTTCAGTGGCAGGGGCATTATCTCTGGCGCGGGTGTAGCTGCCGGAACGCTTGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACTAAAATCGTGc                                                                                                                                                   <  2:355784/150‑1 (MQ=255)
     cACCTTCAGTGTCAGGGGCATTATCTCTGGCGCGGGTGTAGCTGCCGGAACGCTTGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCCCCGCTGTTTCGGCGCGCACTTTTATCAGAAAGCGACGGGTCAAACTTAAATCGTGCgg                                                                                                                                                 >  1:695990/1‑150 (MQ=255)
       ccTTCAGTGTCAGGGGCATTATCTCTGGCGCGGGGGAAGCTGCCGGAACGCTTGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCt                                                                                                                                               <  2:1511422/150‑1 (MQ=255)
       ccTTCAGTGTCAGGGGCATTAGCTCTGGCGCGGGGGTAGCTGCCGGAACGCTGGCCGCCAGCGTAGAATGAGGATCAAGCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTATTCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCt                                                                                                                                               <  2:1517971/150‑1 (MQ=255)
         ttCAGTGGCATGGGCATTATCGCTGGCGCGGGTGAAGCTGCCGGAACGCGTGCCGCCAGCGGAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCg                                                                                                                                             <  1:260312/150‑1 (MQ=255)
           cAGTGTCAGGGGCATTATCGCTGGCGCGGGTGTAGCTGCCGGAACGCTTGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGACTCGTg                                                                                                                                           <  1:839669/150‑1 (MQ=255)
            aGTGTCAGGGGCATTATCTCTGGCGCGGGTGTAGCTGCCGGAACGCTTGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCAGTTTGGGCGGGAACTATAATAACATCAAGACCGGGCAAACAAAAAAAGTGCGGcaccttc                                                                                                                                          >  2:1537666/1‑144 (MQ=255)
             gtgtCAGGGGCATTATCTCTGGCGCGGGTGTAGCTGCCGGAACGCTTGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCt                                                                                                                                         <  1:1411637/150‑1 (MQ=255)
              tgtCAGGGGCATTATCTCTGGCGCGGGTGTAGCTGCCGGAACGCTTGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCCCGGTTTAGGCCGTAACATTTATGATTTGGCGTCGGGTCGAACTTAAATTGGTCGGCTTGTGGtt                                                                                                                                        >  2:1038272/1‑150 (MQ=255)
               gtCAGGGGCATTATCTCTGGCGCGGGTGTAGCTGCCGGAACGCTTGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTc                                                                                                                                       >  2:101143/1‑150 (MQ=255)
               gtCAGGGGCATTATCTCTGGCGCGGGTGTAGCTGCCGGAACGCTTGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTc                                                                                                                                       >  2:1219450/1‑150 (MQ=255)
                tCAGGGGCATTATCTCTGGCGCGGGTGTAGCTGCCGGAACGCTTGCCGCCAGCGTATAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTcc                                                                                                                                      <  2:587252/150‑1 (MQ=255)
                    gggCATTATCTCTGGCGCGGGTGTAGCTGCCGGAACGCTTGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAAc                                                                                                                                  <  1:514495/150‑1 (MQ=255)
                      gCATTATCTCTGGCGCGGGTGTAGCTGCCGGAACGCTTGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCTGGTGCGGCGCTCACTGTAATCAGATCGCGACGGGTCAAACTTAAATCGTGGGGGTCGTGCTTCCCGAccc                                                                                                                                >  2:1140238/1‑150 (MQ=255)
                         ttATCTCTGGCGCGGGTGTAGCTGCCGGAACGCTTGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACATTAATCAGATCCCGACGGGGGAAAATTAAATACTGCGGCTCGGGCGTCCCAACCCGaa                                                                                                                             >  2:734870/1‑149 (MQ=255)
                           aTCTCTGGCGCGGGTGTAGCTGCCGGAACGCTTGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAACCCGATTa                                                                                                                           <  2:736169/150‑1 (MQ=255)
                            tctcTGGCGCGAGTGTAGCTGCAGGAACGCTTGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATAAGCTCGCGCCGGGTGAAACTTAAATCGTGCGGCTCGTTCTTACCCACCCGATTAt                                                                                                                          >  1:314792/1‑150 (MQ=255)
                             ctctGGCGCGGGTGTAGCTGCCGGAACGCTTGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCCCGGTTGGGGGGGGCCCTTTAAAAAAATAGCGCCGGGGGAAAAATAAAAAGTGCGGCGCGTGGTTACCAACCCCAatata                                                                                                                         >  2:369852/1‑149 (MQ=255)
                                 ggCGCGGGTGTAGCTGCCGGAACGCTTGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAACCCGATTATCTTCg                                                                                                                     <  2:552917/150‑1 (MQ=255)
                                    gcgTGTGTAGCTGCCGGAACGCTTTCCACCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCTCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGGCAAACTTAAATCGTGCGGCTCGTGCTTCCCATCCCGATTATCTTCgggg                                                                                                                  <  2:1210991/150‑1 (MQ=255)
                                       ggTGTAGCTGCCGGAACGCTTGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATGGCGCCGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCCACCCGATTATCTTCGGGGTAt                                                                                                               >  1:1185665/1‑150 (MQ=255)
                                       ggTGTAGCTGCCGGAACGCTTGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCGGGTTCGGCGCGCACATTAATAAGATAGCGGCGGGGCAAACATAAAAACTGCGGCGCGTGCTTCCCAACCCCAATATATTAGGGGTAt                                                                                                               >  2:402124/1‑150 (MQ=255)
                                         tgtAGCTGCCGGAACGCTTGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCGGTTGCGGCGCTCACGTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTAGTGCGTCCCCACCCGATTATCTTCGGGGTATTc                                                                                                             >  2:1203320/1‑150 (MQ=255)
                                          gtAGCTGCCGGAACGCTTGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAACCCGATTATCTTCGGGGTATTCg                                                                                                            <  2:166273/150‑1 (MQ=255)
                                             gCTGCCGGAACGCTTGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCAAACCCGATTATCTTCGGGGTATTCGCCt                                                                                                         <  2:543655/150‑1 (MQ=255)
                                               tGCCGGAACGCTTGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACGTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAACCCGATTATCTTCGGGGTATTCGCCGTg                                                                                                       >  2:190177/1‑150 (MQ=255)
                                                 ccGGAACGCTTGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTTCGGCCCTCCCATTAATAACAAAGCGACGGGGCAAACATAAATAGTGCGGCTAGTGCTTACCAAACCGAATATATTAGGGGTATTAACCTTTGt                                                                                                     >  2:832868/1‑150 (MQ=255)
                                                 ccGGAACGCTTGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTGCCTTTAATCAGCTAGCGCCCGGGCCAACTTAAATAGTGCGGCGCGTGGTTCCCCACCCGAAGATCATCGGGGGATTAGccattgt                                                                                                     >  2:320283/1‑145 (MQ=255)
                                                       cGCTTGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACATTAATCAGATAGCCACGGGTCAAACTTAAATCTTGCGGCTAGTTCGTCCCAACCCGATTATATTAGGGGTTAACGCCTTGGTTGcgaa                                                                                               >  2:819463/1‑147 (MQ=255)
                                                        gCTTGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAACCCGATTATCTTCGGGGTATTCGCCTTGGTTGCTCAt                                                                                              <  1:701339/150‑1 (MQ=255)
                                                          ttGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAACCCGATTATCTTCGGGGTATTCGCCTTGGTTGCTCATAc                                                                                            <  1:1119607/150‑1 (MQ=255)
                                                           tGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAACCCGATTATCTTCGGGGTATTCGCCTTGGTTGCTCATAcc                                                                                           <  1:101143/150‑1 (MQ=255)
                                                           tGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCGCCCTTTAATCAGATCGCGACGGGTCAAACTTAAATCGGGCGGCGCGTGCTTCCCAACCCGATTATCTTCGGGGGATTCGCCTTGGTTGCTCCTAcc                                                                                           >  1:1126081/1‑150 (MQ=255)
                                                             ccgccAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAACCCGATTATCTTCGGGGTATTCGCCTTGGTTGCTCATACCAg                                                                                         <  2:1691552/150‑1 (MQ=255)
                                                                ccAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCAGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAACCCGATTATCTTCGGGGTATTCGCCTTGGTTGCTCATACCAGGcc                                                                                      <  2:399067/150‑1 (MQ=255)
                                                                   gCGTAGAATGAGGATAAACCACGGCTGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCTCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAACCCGATTATCTTCGGGGTATTCGCCTTGGTTGCTCATACCAGGCCTcc                                                                                   <  1:1622707/150‑1 (MQ=255)
                                                                        gAATGAGGATAATCCACGGCGGTCGCCGCTGGTGCGGCGCTCACTTTAATCAGATCGCGACGGGGCAAACTTAAATCGTGCGGCTCGTGCTTCCCAACCCGATTATCTTCGGGGTATTCGCCTTGGTTGCTCATACCAGGCCTCCGGTAt                                                                              <  1:240253/150‑1 (MQ=255)
                                                                          aTGAGTATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAACCCGATTATCTTCGGGGTATTCGCCTTGGTTGCTCATACCAGGCCTCCGGTATTc                                                                            <  1:1038272/150‑1 (MQ=255)
                                                                            gAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAACCCGATTATCTTCGGGGTATTCGCCTTGGTTGCTCATACCAGGCCTCCGGTATTCTa                                                                          <  1:112920/150‑1 (MQ=255)
                                                                            gAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAACCCGATTATCTTCGGGGTATTCGCCTTGGTTGCTCATACCAGGCCTCCGGTATTCTa                                                                          <  2:371086/150‑1 (MQ=255)
                                                                              ggATAAACCACGGCGGTCGCCGCAGTTGCGGCGCTCCCTGTATTCAGGTCGCGACGGGTCAAGCTTAAATCGTGCGGCTCGTGCTTCCCGCCCCGATTATGTGCGGGGTATTCGCCGTGGGTGGGAAGACCAGGCCTCCGGGATTCGaat                                                                        >  2:273748/1‑149 (MQ=255)
                                                                               gATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAACCCGATTATCTTCGGGGTATTCGCCTTGGTTGCTCATACCAGGCCTCCGGTATTCTAAAg                                                                       <  2:1131716/150‑1 (MQ=255)
                                                                                  aaaCCCCGGGGGGGGCCGCTGGTGCGGCGGTCACTTTAATCAGATCGGGACGGGTCAAACTTAAATAGTGCGGCGCGTGCTTCCCAACCCGATTATCTTCGGGGTATTCGCCTTGGTTGCTCTTACCAGGCCTCCGGTTTTGTAAAgggg                                                                    <  2:1291473/150‑1 (MQ=255)
                                                                                  aaaCCACGGCGGTCGCCGCTGTTGGGGCGGGCACATTTATAAGATAGCGACGGGGCAAAATTTAATAGTGCGGCGCGTGGGTGCCAACCCGATTTTGTTAGGGGTTTTAGCCGTGGTTGGTAATAACAGGGCTCCGGGTTTGTAAAaggg                                                                    >  2:398123/1‑150 (MQ=255)
                                                                                  aaaCCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAACCCGATTATCTTCGGGGTATTCGCCTTGGTTGCTCATACCAGGCCTCCGGTATTCTAAAgggg                                                                    <  1:1702120/150‑1 (MQ=255)
                                                                                  aaaCCACGGCGGTCGCCGCTGTTGCGGCGCACAAGTTAATAAGATAGCGACGGGGCATACGTAAAAAGTGCGGCGCGAGCTTCCCACCCCGCTTTTCTTAGGGGTTTTTGGCGTGGTTGGGCATACCAGGGCGCCGGGTTTCTAAAgggg                                                                    >  2:321030/1‑150 (MQ=255)
                                                                                        cggcggTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAACCCGATTATCTTCGGGGTATTCGCCTTGGTTGCTCATACCAGGCCTCCGGTATTCTAAAGGGGAAAATa                                                              <  1:734870/150‑1 (MQ=255)
                                                                                              gcgccgcTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAACCCGATTATCTTCGGGGTATTCGCCTTGGTTGCTCATACCAGGCCTCCGGTATTCTAAAGGGGAAAATAAGAGTg                                                        <  1:384300/149‑1 (MQ=255)
                                                                                                 ccgcTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCCACCCGATGAGCTTCGGGGTATTCGACTTGGTTGCTCATACCAGGCCTCCCGTATTCTAATGGGGAAAATAATAGTGCGt                                                     >  1:1089525/1‑150 (MQ=255)
                                                                                                  cgcTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAACCCGATTATCTTCGGGGTATTCGCCTTGGTTGCTCATACCAGGCCTCCGGTATTCTAAAGGGGAAAATAAGAGTGCGtt                                                    >  1:1219661/1‑150 (MQ=255)
                                                                                                  cgcTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAACCCGATTATCTTCGGGGTATTCGCCTTGGTTGCTCATACCAGGCCTCCGGTATTCTAAAGGGGAAAATAAGAGTGCGtt                                                    <  2:1177803/150‑1 (MQ=255)
                                                                                                   gcTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAACCCGATTATCTTCGGGGTATTCGCCTTGGTTGCTCATACCAGGCCTCCGGTATTCTAAAGGGGAAAATAAGAGTGCGTTg                                                   >  1:54980/1‑150 (MQ=255)
                                                                                                   gcTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAACCCGATTATCTTCGGGGTATTCGCCTTGGTTGATCATACCAGGCATCAGGTATTCTAAAGGGGTAAATAAGAGTGCGTTg                                                   >  1:49634/1‑150 (MQ=255)
                                                                                                   gcTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGGCAAACTTAAATCGTGCGGCTCGTGCTTCCCAACCCGATTTTCTTCGGGGTATTCGCCTTGGTTGCTCATACCAGGCCTCCGGTATTCTAAAGGGGAAAATAAGAGTGCGTTg                                                   <  1:813839/150‑1 (MQ=255)
                                                                                                      gTTGCGGCGCTAACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAACCCGATTATCTTCGGGGTATTCGCCTTGGTTGCTCATACCAGGCCTCCGGTATTCTAAAGGGGAAAATAAGAGTGCGTTGCga                                                >  2:158708/1‑150 (MQ=255)
                                                                                                       ttGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAACCCGATTATCTTCGGGGTATTCGCCTTGGTTGCTCATACCAGGCCTCCGGTATTCTAAAGGGGAAAATAAGAGTGCGTTGCgag                                               <  1:984808/150‑1 (MQ=255)
                                                                                                       ttGCGGCGCTCAATTTAATCAGATCGCGACTGGTCAAACTTAAATAGTTCGGCTCGTGCGTCCCAACCCGATTATCTTCGGGGTATTCGCCTTGGTTGCTCATACCAGGCCTCCGGTATTCTAAAGGGGAAAATAAGAGTGCGTTGCgag                                               <  1:903058/150‑1 (MQ=255)
                                                                                                        tGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAACCCGATTATCTTCGGGGTATTCGCCTTGGTTGCTCATACCAGGCCTCCGGTATTCTAAAGGGGAAAATAAGAGTGCGTTGCgaga                                              >  2:1711714/1‑150 (MQ=255)
                                                                                                            gcgcTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAACCCGATTATCTTCGGGGTATTCGCCTTGGTTGCTCATACCAGGCCTCCGGTATTCTAAAGGGGAAAATAAGAGTGCGTTGCGAGAGTGc                                          <  1:832868/150‑1 (MQ=255)
                                                                                                             cgcTCACTTTAATCAGATCGCGACGGGTCAAACTTAATTCGTGCGGCTCGTGCTTCCCAACCCGATTATCTTCGGGGTTTTCGCCTTGGTTGCTCATACCAGGCCTCCGGTATTCTAAAGGGGAAAATAAGAGTGCGTTGCGAGAGTGCt                                         <  1:1203320/150‑1 (MQ=255)
                                                                                                             cgcTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAACCCGATTATCTTCGGGGTATTCGCCTTGGTTGCTCATACCAGGCCTCCGGTATTCTAAAGGGGAAAATAAGAGTGCGTTGCGAGAGTGCt                                         >  1:35665/1‑150 (MQ=255)
                                                                                                              gcTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAACCCGATTATCTTCGGGGTATTCGCCTTGGTTGCTCATACCAGGCCTCCGGTATTCTAAAGGGGAAAATAAGAGTGCGTTGCGAGAGTGCTg                                        <  2:104385/150‑1 (MQ=255)
                                                                                                               cTCACTTTAATCAGCTCGCGACGGGTCAAACTTACCTCGTGCGGCTCGTGCGACCCAACCCGATTATCTTCGGGGTATTCGCCTTGGTTGCTCATACCAGGCCTCCGGTATTCTAAAGGGGAATTTAAGAGTGCGTTGCGAGAGTGCTGa                                       <  2:1076176/150‑1 (MQ=255)
                                                                                                               cTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAACCCGATTATCTTCGGGGTATTCGCCTTGGTTGCTCATACCAGGCCTCCGGTATTCTAAAGGGGAAAATAAGAGTGCGTTGCGAGAGTGCTGa                                       >  2:910823/1‑150 (MQ=255)
                                                                                                                tCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAACCCGATTATCTTCGGGGTATTCGCCTTGGTTGCTCATACCAGGCCTCCGGTATTCTAAAGGGGAAAATAAGAGTGCGTTGCGAGAGTGCTGAt                                      <  1:419320/150‑1 (MQ=255)
                                                                                                                    tttAATCAGATCGCGACGGGGCAAACTTAAATCGTGCGGCTCGTGCTTCCCAACCCGATTATCTTCGGGGTATTCGCCTTGGTTGCTCATACCAGGCCTCCGGTATTCTAAAGGGGAAAATAAGAGTGCGTTGCGAGAGTGCTGATGCtt                                  <  1:881150/150‑1 (MQ=255)
                                                                                                                     ttAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAACCCGATTATCTTCGGGGTATTCGCCTTGGTTGCTCATACCAGGCCTCCGGTATTCTAAAGGGGAAAATAAGAGTGCGTTGCGAGAGTGCTGATGCttt                                 >  1:1567692/1‑150 (MQ=255)
                                                                                                                     ttAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAACCCGATTATCTTCGGGGTATTCGCCTTGGTTGCTCATACCAGGCCTCCGGTATTCTAAAGGGGAAAATAAGAGTGCGTTGCGAGAGTGCTGATGCttt                                 >  1:1078887/1‑150 (MQ=255)
                                                                                                                       aaTCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAACCCGATTATCTTCGGGGTATTCGCCTTGGTTGCTCATACCAGGCCTCCGGTATTCTAAAGGGGAAAATAAGAGTGCGTTGCGAGAGTGCTGATGCttttt                               >  2:1139055/1‑150 (MQ=255)
                                                                                                                       aaTCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAACCCGATTATCTTCGGGGTATTCGCCTTGGTTGCTCATACCAGGCCTCCGGTATTCTAAAGGGGAAAATAAGAGTGCGTTGCGAGAGTGCTGATGCttttt                               <  1:1125695/150‑1 (MQ=255)
                                                                                                                        aTCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAACCCGATTATCTTCGGGGTATTCGCCTTGGt                                                                                                     <  1:207301/79‑1 (MQ=255)
                                                                                                                        aTCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAACCCGATTATCTTCGGGGTATTCGCCTTGGt                                                                                                     >  2:207301/1‑79 (MQ=255)
                                                                                                                        aTCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAACCCGATTATCTTCGGGGTATTCGCCTTGGTTGCTCATACCAGGCCTCCGGTATTCTAAAGGGGAAAATAAGAGTGCGTTGCGAGAGTGCTGATGCtttttt                              >  1:259308/1‑150 (MQ=255)
                                                                                                                         tCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAACCCGATTATCTTCGGGGTATTCGCCTTGGTTGCTCATACCAGGCCTCCGGTATTCTAAAGGGGAAAATAAGAGTGCGTTGCGAGAGTGCTGATGCTTTTTTa                             <  1:397794/150‑1 (MQ=255)
                                                                                                                         tCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAACCCGATTATCTTCGGGGTATTCGCCTTGGTTGCTCATACCAGGCCTCCGGTATTCTAAAGGGGAAAATAAGAGTGCGTTGCGAGAGTGCTGATGCTTTTTTa                             <  2:1492859/150‑1 (MQ=255)
                                                                                                                         tCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAACCCGATTATCTTCGGGGTATTCGCCTTGGTTGCTCATACCAGGCCTCCGGTATTCTAAAGGGGAAAATAAGAGTGCGTTGCGAGAGTGCTGATGCTTTTTTa                             >  2:605275/1‑150 (MQ=255)
                                                                                                                         tCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAACCCGATTATCTTCGCGGTATTCGCCTTGGTTGCTCATACCAGGCCTCCGGTATTCTAAAGGGGAAAATAAGAGTGCGTTGCGAGAGTGCTGATGCTTTTTTa                             >  2:86162/1‑150 (MQ=255)
                                                                                                                            taTCGCGCCGGGTCAAACTTAACTCGTGCGGCTCGTGCTTCCCAACCCGATTATCTTAGGGGTATTCGCCTTGGTTGCTCATACCAGGCCTCCGGTATTCTAAAGGGGAAAATAAGAGTGCGTTGCGAGAGTGCTGATGCTTTTTTAAAc                          <  1:819463/149‑1 (MQ=255)
                                                                                                                            gATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAACCCGATTATCTTCGGGGTATTCGCCTTGGTTGCTCATACCAGGCCTCCGGTATTCTAAAGGGGAAAATAAGAGTGCGTTGCGAGAGTGCTGATGCTTTTTTAAAc                          >  2:1543965/1‑150 (MQ=255)
                                                                                                                            gATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAACCCGATTATCTTCGGGGTATTCGCCTTGGTTGCTCATACCAGGCCTCCGGTATTCTAAAGGGGAAAATAAGAGTGCGTTGCGAGAGTGCTGATGCTTTTTTAAAc                          <  2:1185665/150‑1 (MQ=255)
                                                                                                                             aTCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAACCCGATTATCTTCGGGGTATTCGCCTTGGTTGCTCATACCAGGCCTCCGGTATTCTAAAGGGGAAAATAAGAGTGCGTTGCGAGAGTGCTGATGCTTTTTTAAACg                         >  2:349114/1‑150 (MQ=255)
                                                                                                                                 cgACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAACCCGATTATCTTCGGGGTATTCGCCTTGGTTGCTCATACCAGGCCTCCGGTATTCTATAGGGGAAAATAAGAGTGCGTTGCGAGAGTGCTGATGCTTTTTTAAACGTTat                     >  1:1068730/1‑149 (MQ=255)
                                                                                                                                  gACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAACCCGATTATCTTCGGGGTATTCGCCTTGGTTGCTCATACCAGGCCTCCGGTATTCTAAAGGGGAAAATAAGAGTGCGTTGCGAGAGTGCTGATGCTTTTTTAAACGTTAAg                    <  1:788636/150‑1 (MQ=255)
                                                                                                                                   ccGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAACCCGATTATCTTCGGGGTATTCGCCTTGGTTGCTCATAACAGGCCTCCGGTATTCTAAAGGGTAAAATAAGAGTGCGTTGCGAGAGTGCTGATGCTTTTTTAAACGTTAAGc                   <  1:1553459/149‑1 (MQ=255)
                                                                                                                                   aCGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAACCCGATTATCTTCGGGGTATTCGCCTTGGTTGCTCATACCAGGCCTCCGGTATTCTAAAGGGGAAAATAAGAGTGCGTTGCGAGAGTGCTGATGCTTTTTTAAACGTTAAGc                   <  1:402124/150‑1 (MQ=255)
                                                                                                                                   aCGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAACCCGATTATCTTCGGGGTATTCGCCTTGGTTGCTCATACCAGGCCTCCGGTATTCTAAAGGGGAAAATAAGAGTGCGTTGCGAGAGTGCTGATGCTTTTTTAAACGTTAAGc                   >  1:577002/1‑150 (MQ=255)
                                                                                                                                     gggTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAACCCGATTATCTTCGGGGTATTCGCCTTGGTTGCTCATACCAGGCCTCCGGTATTCTAAAGGGGAAAATAAGAGTGCGTTGCGAGAGTGCTGATGCTTTTTTAAACGTTAAGCAt                 <  1:1016370/150‑1 (MQ=255)
                                                                                                                                       gTCAAACTTAAATCGTGCGGCTCTTGCTTCCCAACCCGATTATCTTCGGGGTTTTCGCCTTGTTTGCTCATTCCAGGCCTCCGGTATTCTAAAGGGGAAAATAAGAGTGCGTTGCGAGAGTGCTGATGCTTTTTTAAACGTTAAGCATAg               <  1:948758/150‑1 (MQ=255)
                                                                                                                                         cAAACTTAAATCGTGCGGCTCGTGCTTCCCAACCCGATTATCTTCGGGGTATTCGCCTTGGTTGCTCATACCAGGCCTCCGGTATTCTAAAGGGGAAAATAAGAGTGCGTTGCGAGAGTGCTGATGCTTTTTTAAACGTTAAGCATAGTc             >  1:13802/1‑150 (MQ=255)
                                                                                                                                         aaaaCTTAAATCGTGCGGCTCGTGCTTCCCAACCCGATTATCTTCGGGGTATTCGCCTTGGTTGCTCATACCAGGCCTCCGGTATTCTAAAGGGGAAAATAAGAGTGCGTTGCGAGAGTGCTGATGCTTTTTTAAACGTTAAGCATAGTc             <  1:1537666/149‑1 (MQ=255)
                                                                                                                                           aaCTTAAATCGTGCGGCTCGTGCTTCCCAACCCGATTATCTTCGGGGTATTCGCCTTGGTTGCTCATACCAGGCCTCCGGTATTCTAAAGGGGAAAATAAGAGTGCGTTGCGAGAGTGCTGATGCTTTTTTAAACGTTAAGCATAGTCgg           >  2:495214/1‑150 (MQ=255)
                                                                                                                                             cTTAAATCGTGCGGCTCGTGCTTCCCAACCCGATTATCTTCGGGGTATTCGCCTTGGTTGCTCATACCAGGCCTCCGGTATTCTAAAGGGGAAAATAAGAGTGCGTTGCGAGAGTGCTGATGCTTTTTTAAACGTTAAGCATAGTCGGCg         >  2:354717/1‑150 (MQ=255)
                                                                                                                                             cTTAAATCGTGCGGCTCGTGCTTCCCAACCCGATTATCTTCGGGGTATTCGCCTTGGTTGCTCATACCAGGCCTCCGGTATTCTAAAGGGGAAAATAAGAGTGCGTTGCGAGAGTGCTGATGCTTTTTTAAACGTTAAGCATAGTCGGCg         >  2:1473044/1‑150 (MQ=255)
                                                                                                                                                   tCGTGCGGCTCGTGCTTCCCAACCCGATTATCTTCGGGGTATTCGCCTTGGTTGCTCATACCAGGCCTCCGGTATTCTAAAGGGGAAAATAAGAGTGCGTTGCGAGAGTGCTGATGCTTTTTTAAACGTTAAGCATAGTCGGCGAcagca   <  1:398123/150‑1 (MQ=255)
                                                                                                                                                    cGTGCGGCTCGTGCTTCCCAACCCGATTATCTTCGGGGTATTCGCCTTGGTTGCTCATACCAGGCCTCCGGTATTCTAAAGGGGAAAATAAGAGTGCGTTGCGAGAGTGCTGATGCTTTTTTAAACGTTAAGCATAGTCGGCGAcagcag  <  2:1219661/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CCGTTCACCTTCAGTGTCAGGGGCATTATCTCTGGCGCGGGTGTAGCTGCCGGAACGCTTGCCGCCAGCGTAGAATGAGGATAAACCACGGCGGTCGCCGCTGTTGCGGCGCTCACTTTAATCAGATCGCGACGGGTCAAACTTAAATCGTGCGGCTCGTGCTTCCCAACCCGATTATCTTCGGGGTATTCGCCTTGGTTGCTCATACCAGGCCTCCGGTATTCTAAAGGGGAAAATAAGAGTGCGTTGCGAGAGTGCTGATGCTTTTTTAAACGTTAAGCATAGTCGGCGACAGCAG  >  NC_000913/302368‑302665

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