Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 3,773,715 C→A 5.8% T370K (ACA→AAA)  nikA → Nickel‑binding periplasmic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1223,773,7150CA 5.8% 202.2 / 4.0 70T370K (ACA→AAA) nikANickel‑binding periplasmic protein
Reads supporting (aligned to +/- strand):  ref base C (39/27);  new base A (2/2);  total (41/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.83e-01

TCAAGCCCTACCCTTACGACCCGGCGCGCGCCGCCCATGCCCTTGATGCTGCCGGGTGGACGCTGGCCCCAGGCGCGTCTGTGCGCAGTAAAGGCGGCGAACCCTTGGCCATCGAGCTGGTATTCCTCGGCACAAGTGCCTTGCAGAAGAGCCTGGCCGAAGTGCTGCAGGGCGAGCTGGCGAAAATCGGCATCCACATCGACCTGCGTGCGGTTGAAGAAGGGGCTTTGGTGCGGCGCCAACGCGACGGCGACTTCGGCATGA  >  ACB122/3773583‑3773846
                                                                                                                                    |                                                                                                                                   
tCAAGCGCTACCCTTACGACCCGGCGCGCGCCGCCCATGCCCTTGATGCTGCCGGGTGGACGCTGGCCCCAGTCGCGTCTGTCCGCAGTAAATGCGGCGAACCGTTGTCCATCGAGCTGGTATTCCTCGGCACaa                                                                                                                                   <  8:317350/135‑1 (MQ=255)
tCAAGCCCTACCCTTACGACCCGGCGCGCGCCGCCCATGCCCTTGATGCTGCCGGGTGGACGCTGGCCCCAGGCGCGTCTGTGCGCAGTAAAGGCGGCGAACCCTTGGCCATCGAGCTGGTATTCCTCGGCACaa                                                                                                                                   <  7:350924/135‑1 (MQ=255)
tCAAGCCCTACCCTTACGACCCGGCGCGCGCCGCCCATGCCCTTGATGCTGCCGGGTGGACGCTGGCCCCAGGCGCGTCTGTGCGCAGTAAAGGCGGCGAACCCTTGGCCATCGAGCTGGTATTCCTCGGCACaa                                                                                                                                   <  3:2918/135‑1 (MQ=255)
tCAAGCCCTACCCTTACGACCCGGCGCGCGCCGCCCATGCCCTTGATGCTGCCGGGTGGACGCTGGCCCCAGGCGCGGCTGTGCGCAGTAAAGGCGGCGAACCCTTGGCCATCGAGCTGGTATTCCTCGGCCCaa                                                                                                                                   >  8:416682/1‑135 (MQ=255)
  aaGCCCTACCCTTACGACCCGGCGCGCGCCGCCCATGCCCTTGATGCTGCCGGGTGGACGCTGGCCCCAGGCGCGGCTGTGCGCAGGAAAGGCGGCGAACCCTTGGCCATCGAGCTGGTTTTCCTCGGCCCaact                                                                                                                                 >  4:1540/1‑133 (MQ=255)
   aGCCCTACCCTTACGACCCGGCGCGCGCCGCCCATGCCCTTGATGCTGCCGGGTGGACGCTGGCCCCAGGCGCGTCTGTGCGCAGTAAAGGCGGCGAACCCTTGGCCATCGAGCTGGTATTCCTCGGCACAAGTg                                                                                                                                <  5:418577/135‑1 (MQ=255)
    gCCCTACCCTTACGACCCGGCGCGCGCCGCCCATGCCCTTGATGCTGCCGGGTGGACGCTGGCCCCAGGCGCGTCTGTGCGCAGTAAAGGCGGCGAACCCTTGGCCATCGAGCTGGTATTCCTCGGCACaa                                                                                                                                   <  3:354549/131‑1 (MQ=255)
    gCCCTACCCTTACGACCCGGCGCGCGCCGCCCATGCCCTTGATGCTGCCGGGTGGACGCTGGCCCCAGGCGCGTCTGTGCGCAGTAAAGGCGGCGAACCCTTGGCCATCGAGCTGGTATTCCTCGGCACaa                                                                                                                                   >  4:354549/1‑131 (MQ=255)
           ccTTACGACCCGGCGCGCGCCGCCCATGCCCTTGATGCTGCCGGGTGGACGCTGGCCCCAGGCGCGTGTGTTCGGAGTAAAAGCGGCGCACCCGTGGCCCTCGCGCTGGTATTCCTCGGCACAAGTGCCTTGCa                                                                                                                         >  5:23899/1‑134 (MQ=255)
           ccTTACGACCCGGCGCGCGCCGCCCATGCCCTTGATGCTGCCGGGTGGACGCTGGCCCCAGGCGCGTCTGTGCGGAGTAAAGGCGGCGAACCCTGGGCCCTCGAGCTGGTATTACTCGGCCCAGGTGCCTTGCa                                                                                                                         >  1:4283/1‑134 (MQ=255)
                  aCCCGGCGCGCGCCGCCCATGCCCTTGATGCTGCCGGGTGGACGCTGGCCCCAGGCGCGTCTGTGCGCAGTAAAGGCGGCGACCCCTTGGCCATCGAGCTGGTATTCCTCGGCACAAGTGCCTTGCAGAAGAGcc                                                                                                                 >  4:334147/1‑135 (MQ=255)
                  aCCCGGCGCGCGCCGCCCATGCCCTTGATGCTGCCGGGTGGACGCTGGCCCCAGGCGCGTCTGTGCGCAGTAAAGGCGGCGAACCCTTGGCCATCGAGCTGGTATTCCTCGGCACCAGTGCCTTGCAGCAGAGcc                                                                                                                 >  8:102029/1‑135 (MQ=255)
                  aCCCGGCGCGCGCCGCCCATGCCCTTGATGCTGCCGGGTGGACGCTGGCCCCAGGCGCGTCTGTGCGCAGTAAAGGCGGCGAACCCTTGGCCAGCGAGCTGGTATTCCTCGGCACAAGTGCCTTGCAGAAGAGcc                                                                                                                 >  8:368198/1‑135 (MQ=255)
                     cGGCGCGCGCCGCCCATGCCCTTGATGCTGCCGGGTGGACGCTGGCCCCAGGCGCGTCTGTGCGCAGTAAAGGCGGCGAACCCTTGGCCATCGAGCTGGTATTCCTCGGCACAAGTGCCTTGCAGAAGAGCCTgg                                                                                                              <  7:60156/135‑1 (MQ=255)
                     cGGCGCGCGCCGCCCATGCCCTTGATGCTGCCGGGTGGACGCTGGCCCCAGGCGCGTCTGTGCGCAGTAAAGGCGGCGAACCCTTGGCCATCGAGCTGGTATTCCTCGGCACAAGTGCCTTGCAGAAGAGCCTgg                                                                                                              <  8:235706/135‑1 (MQ=255)
                       gcgcgcgcCGCCCATGCCCTTGATGCTGCCGGGTGGACGCTGGCCCCCGGCGCGTCTGTGCGCAGTAAAGGCGGCGAACCCTTGGCCATCGAGCTGGTATTCCTCGGCACAAGTGCCTTGCAGAAGAGCCTGGc                                                                                                             >  1:311415/1‑134 (MQ=255)
                       gcgcgcgcCGCCCATGCCCTTGATGCTGCCGGGTGGACGCTGGCCCCAGGCGCGTCTGTGCGCAGTAAAGGCGGCGAACCCTTGGCCATCGAGCTGGTATTCCTCGGCACAAGTGCCTTGCAGAAGAGCCTGGcc                                                                                                            <  4:181386/135‑1 (MQ=255)
                       gcgcgcgcCGCCCATGCCCTTGATGCTGCCGGGTGGACGCTGGCCCCAGGCGCGGCTGTGCGCAGTAAAGGCGGCGCACCCTTGGCCCTCGCGCGGGGATTCCTCGGCCCAAGTGCCGTGCAGCAGCGCCTGGcc                                                                                                            >  2:396801/1‑135 (MQ=255)
                       gcgcgcgcCGCCCATGCCCTTGATGCTGCCGGGTGGACGCGGGCCCCAGGCGCGTCTGTGCGCAGTAAAGGCGGCGAACCCTTGGCCATCGAGCTGGTATTCCTCGGCACAAGGGCCTTGCAGAAGAGCCTGGcc                                                                                                            <  7:90791/135‑1 (MQ=255)
                        cgcgcgcCGCCCATGCCCTTGATGCTGCCGGGTGGACGCTGGCCCCAGGCGCGTCTGTGCGCAGTAAAGGCGGCGAACCCTTGGCCATCGAGCTGGTATTCCTCGGCACAAGTGCCTTGCAGAAGAGCCTGGcc                                                                                                            >  1:269220/1‑134 (MQ=255)
                                  ccATGCCCTTGATGCTGCCGGGGGGACGCTGGCCCCAGGCGCGTCTGTGCGCAGTTAAGGCGGCGAACACTTGGCCCTCGAGCTGGTATTCCTCGGGACACCTGCCTTGCAGAAGAGCCTGGCCGAAGTGCTGCa                                                                                                 <  1:295286/135‑1 (MQ=255)
                                               gctgcCGGGTGGACGCTGGCCCCAGGCGCGTCTGTGCGCAGTAAAGGCGGCGACCCCTTGGCCATCGAGCTGGTATTCCTCGGCACAAGTGCCTTGCAGAAGAGCCTTGCCGAAGTGCTGCAGGGCGCGCTGGc                                                                                     >  1:323787/1‑134 (MQ=255)
                                               gctgcCGGGTGGACGCTGGCCCCAGGCGCGTCTGTGCGCAGTAAAGGCGGCGAACCCTTGGCCATCGAGCTGGTATTCCTCGGCACAAGTGCCTTGCAGAAGAGCCTGGCCGAAGTGCTGCAGGGCGAGCTGGc                                                                                     >  8:74092/1‑134 (MQ=255)
                                                    cGGGTGGACGCTGGCCCCAGGCGCGTCTGTGCGCAGTAAAGGCGGCGAACCCTTGGCCATCGAGCTGGTATTCCTCGGCACAAGTGCCTTGCAGAAGAGCCTGGCCGAAGTGCTGCAGGGCGAGCTGGCGAAAAt                                                                               >  1:310864/1‑135 (MQ=255)
                                                    cGGGTGGACGCTGGCCCCAGGCGCGTCTGTGCGCAGTAAAGGCGGCGAACCCTTGGCCATCGAGCTGGTATTCCTCGGCACAAGTGCCTTGCAGAAGAGCCTGGCCGAAGTGCTGCAGGGCGAGCTGGCGAAAAt                                                                               >  6:157243/1‑135 (MQ=255)
                                                    cGGGTGGACGCTGGCCCCAGGCGCGTCTGTGCGCAGTAAAGGCGGCGAACCCTTGGCCATCGAGCTGGTATTCCTCGGCACAAGTGCCTTGCAGAAGAGCCTGGCCGAAGTGCTGCAGGGCGAGCTGGCGAAAAt                                                                               >  3:231470/1‑135 (MQ=255)
                                                     gggTGGACGCTGGCCCCAGGCGCGTCTGTGCGCAGTAAAGGCGGCGACCCCTTGGCCATCGAGCTGGGATTCCTCGGCACAAGTGCCTTGCAGAAGAGCCTGGCCGAAGTGCTGCAGGGCGAGCTGGCGAATATc                                                                              >  1:420607/1‑135 (MQ=255)
                                                     gggTGGACGCTGGCCCCAGGCGCGTCTGTGCGCAGTAAAGGCGGCGACCCCATGGCCATCGAGCTGGTATTACTCGGCACAAGTGCCTTGCAGAAGAGCCTGGCCGAAGTGCTGCAGGGCGATCTGGCGAATATc                                                                              >  1:229512/1‑135 (MQ=255)
                                                     gggTGGACGCTGGCCCCAGGCGCGTCTGTGCGCAGTAAAGGCGGCGAACCCTTGGCCATCGCGCTGGTATTCCTCGGCACAAGTGCCTTGCAGAAGAGCCTGGCCGAAGTGCTGCAGGGCGAGCTGGCGAAAATc                                                                              >  5:412484/1‑135 (MQ=255)
                                                     gggTGGACGCTGGCCCCAGGCGCGTCTGTGCGCAGTAAAGGCGGCGAACCCTTGGCCATCGAGCTGGTATTCCTCGGCACAAGTGCCTTGCAGAAGCGCCTGGCCGAAGTGCTGCAGGGCGAGCTGGCGAAAATc                                                                              >  5:133564/1‑135 (MQ=255)
                                                     gggTGGACGCTGGCCCCAGGCGCGTCTGTGCGCAGTAAAGGCGGCGAACCCTTGGCCATCGAGCGGGTATTCCTCGGCACAAGTGCCTTGCAGAAGCGCCTGGCCGAAGTGCTGCAGGGCGAGCTGGCGAAAATc                                                                              >  2:326738/1‑135 (MQ=255)
                                                       gTGGACGCTGGCCCCAGGCGCGTCTGTGCGCAGTAAAGGCGGCGAACCCTTGGCCATCGAGCTGGTATTCCTCGGCACAAGTGCCTTGCAGAAGAGCCTGGCCGAAGTGCTGCAGGGCGAGCTGGCGAAAATCgg                                                                            <  1:138799/135‑1 (MQ=255)
                                                       gTGGACGCGGGCCCCAGGCGCGTCTGTGCGCAGTAAAGGCGGCGAACCCTTGGCCATCGAGCTGGTATTCCTCGGCACAAGTGCCTTGCAGAAGAGCCTGGCCGAAGTGCTGCAGGGCGAGCTGGCGAAAATCgg                                                                            <  8:249438/135‑1 (MQ=255)
                                                              cTGGCCCCAGGCGCGTCTGTGCGCAGTAAAGGCGGCGAACCCTTGGCCATCGAGCTGGTATTCCTCGGCACAAGTGCCTTGCAGAAGAGCCTGGCCGAAGTGCTGCAGGGCGAGCTGGCGAAAATCGGCATCcac                                                                     >  4:316768/1‑135 (MQ=255)
                                                                ggCCCCAGGCGCGTCTGTGCGCAGTAAAGGCGGCGAACCCTTGGCCATCGAGCTGGTATTCCTCGGCACAAGTGCCTTGCAGAAGAGCCTGGCCGAAGTGCTGCAGGGCGAGCTGGCGAAAATCGGCATCCACAt                                                                   >  4:212156/1‑135 (MQ=255)
                                                                 gCCCCAGGCGCGTCGGTGCGCAGTAAAGGCGGCGAACCCTTGGCCAGCGAGCTGGTATTCCTCGGCACAAGTGCCTTGCAGAAGAGCCTGGCCGAAGTGCTGCAGGGCGAGCTGGCGAAAATCGGCATCCACATc                                                                  <  7:243427/135‑1 (MQ=255)
                                                                      aGGCGCGTCTGTGCGCAGTAAAGGCGGCGAACCCTTGGCCATCGAGCTGGTATTCCTCGGCACAAGTGCCTTGCAGAAGAGCCTGGCCGAAGTGCTGCAGGGCGAGCTGGCGAAAATCGGCATCCACATCGACCt                                                             <  3:417719/135‑1 (MQ=255)
                                                                        gcgcGTCTGTGCGCAGTAAAGGCGGCGAACCCTTGGCCATCGAGCTGGTATTCCTCGGCACAAGTGCCTTGCAGAAGAGCCTGGCCGAAGTGCTGCAGGGCGAGCTGGCGAAAATCGGCATCCACATCGACCtgc                                                           <  2:310225/135‑1 (MQ=255)
                                                                        gcgcGTCTGTGCGCAGTAAAGGCGGCGAACCCTTGGCCATCGAGCTGGTATTCCTCGGCACAAGTGCCTTGCAGAAGAGCCTGGCCGAAGTGCTGCAGGGCGAGCTGGCGAAAATCGGCATCCACATCGACCtgc                                                           <  6:243435/135‑1 (MQ=255)
                                                                          gcgTCTGTGCGCAGTAAAGGCGGCGACCCCTTGGCCATCGAGCTGGTATTCCTCGGCACAAGTGCCTTGCAGAAGAGCCTGGCCGAAGTGCTGCAGGGCGAGCTGGCGAAAATCGGCATCCACATCGACCtgagt                                                         >  6:237451/1‑135 (MQ=255)
                                                                          gcgTCTGTGCGCAGTAAAGGCGGCGAACCCTTGGCCATCGAGCTGGTATTCCTCGGCACAAGTGCCTTGCAGAAGAGCCTGGCCGAAGTGCTGCAGGGCGAGCTGGCGAAAATCGGCATCCACATCGACCtgcgt                                                         >  7:32354/1‑135 (MQ=255)
                                                                              cTGTGCGCAGTAAAGGCGGCGAACCCTTGGCCATCGAGCTGGTATTCCTCGGCACAAGTGCCTTGCAGAAGAGCCTGGCCGAAGTGCTGCAGGGCGAGCTGGCGAAAATCGGCATCCACATCGACCTGCGTGCgg                                                     <  2:229512/135‑1 (MQ=255)
                                                                              cTGTGCGCAGTAAAGGCGGCGAACCCTTGGCCATCGAGCTGGTATTACTCGGCACAAGTGCCTTGCAGAAGAGCCTGGCCGAAGTGCTGCAGGGCGAGCTGGCGAAAATCGGCATCCACATCGCCCTGCGTGCgg                                                     >  6:93155/1‑135 (MQ=255)
                                                                                gtgCGCAGTAAAGGCGGCGAACCCTTGGCCATCGAGCTGGTATTCCTCGGCACAAGTGCCTTGCAGAAGAGCCTGGCCGAAGTGCTGCAGGGCGAGCTGGCGAAAATCGGCATCCACATCGACCTGCGTGCGGt                                                    >  1:78553/1‑134 (MQ=255)
                                                                                  gcgcAGTATAGGCGGCGAACCCTTGGCCATCGAGCTGGTATTCCTCGGCACAAGTGCCTTGCAGAAGAGCCTGGCCGAAGTGCTGCAGGGCGAGCTGGCGAAAATCGGCATCCACATCGACCTGCGTGCGGTTga                                                 <  2:5547/135‑1 (MQ=255)
                                                                                    gcAGTAAAGGCGGCGAACCCTTGGCCATCGAGCTGGTATTCCTCGGCACAAGTGCCTTGCAGAAGAGCCTGGCCGAAGTGCTGCAGGGCGAGCTGGCGAAAATCGGCATCCACATCGACCTGCGTGCGGTTgaag                                               <  6:23899/135‑1 (MQ=255)
                                                                                    gcAGTAAAGGCGGCGAACCCTTGGCCATCGAGCTGGTATTCCTCGGCACAAGTGCCTTGCAGAAGAGCCTGGCCGAAGTGCTGCAGGGCGAGCTGGCGAAAATCGGCATCCACATCGACCTGCGTGCGGTTgaag                                               <  2:78553/135‑1 (MQ=255)
                                                                                    gcAGTAAAGGCGGCGAACCCTTGGCCATCGAGCTGGTATTCCTAGGCACAAGTGCCTTGCAGAAGAGACTGGCCGAAGTGCTGCAGGGCGAGCTGGCGAAAATCGGCATCCACATCGACCTGCGTGCGGTTgaa                                                >  8:384749/1‑134 (MQ=255)
                                                                                         aaaGGCGGCGAACCCTTGGCCATCGAGCTGGTATTCCTCGGCACAAGTGCCTTGCAGAAGAGCCTGGCCGAAGTGCTGCAGGGCGAGCTGGCGAAAATCGGCATCCACATCGACCTGCGTGCGGTTGAAGAAggg                                          >  2:316805/1‑135 (MQ=255)
                                                                                           aGGCGGCGAACCCTTGGCCATCGAGCTTGTATTCCTCGGCACAGGTGCCTTGCAGAAGAGCCTGGCCGAAGTGCTGCAGGGCGAGCTGGCGAAAATCGGCATCCACATCGACCTGCGTGCGGTTGAAGAAgggg                                         <  2:4283/134‑1 (MQ=255)
                                                                                              cggcgAACCCTTGGCCATCGAGCTGGTATTCCTCGGCAAAAGTGCCTTGCAGAAGAGCCTGGCCGAAGTGCTGCAGGGCGAGCTGGCGAAAATCGGCATCCACATCGACCTGCGTGCGGTTga                                                 <  7:75646/123‑1 (MQ=255)
                                                                                              cggcgAACCCTTGGCCATCGAGCTGGTATTCCTCGGCAAAAGTGCCTTGCAGAAGAGCCTGGCCGAAGTGCTGCAGGGCGAGCTGGCGAAAATCGGCATCCACATCGACCTGCGTGCGGTTga                                                 >  8:68311/1‑123 (MQ=255)
                                                                                              cggcgAACCCTTGGCCATCGAGCTGGTATTCCTCGGCAAAAGTGCCTTGCAGAAGAGCCTGGCCGAAGTGCTGCAGGGCGAGCTGGCGAAAATCGGCATCCACATCGACCTGCGTGCGGTTga                                                 <  7:68311/123‑1 (MQ=255)
                                                                                              cggcgAACCCTTGGCCATCGAGCTGGTATTCCTCGGCAAAAGTGCCTTGCAGAAGAGCCTGGCCGAAGTGCTGCAGGGCGAGCTGGCGAAAATCGGCATCCACATCGACCTGCGTGCGGTTga                                                 >  8:75646/1‑123 (MQ=255)
                                                                                               ggcgCACCCTTGGCCATCGAGCGGGTATGCCTCGGCACAAGTGCCTGGCAGAAGAGCCTGGCCGAAGTGCTGCAGGGCGAGCTGGCGAAAATCGGCATCCACATCGACCTGCGTTCGGTTGAAGAAGGGGCTTTg                                    <  7:87403/135‑1 (MQ=255)
                                                                                                gcgAACCCTTGGCCATCGAGCTGGTATTCCTCGGCACAAGTGCCTTGCAGAAGAGCCTGGCCGAAGTGCTGCAGGGCGAGCTGGCGAAAATCGGCATCCACATCGACCTGCGTGCGGTTGAAGAAGGGGCTTTgg                                   >  4:435664/1‑135 (MQ=255)
                                                                                                   aaCCCTTGGCCATCGAGCTGGTATTCCTCGGCACAAGTGCCTTGCAGAAGAGCCTGGCCGAAGTGCTGCAGGGCGAGCTGGCGAAAATCgg                                                                            <  1:87828/91‑1 (MQ=255)
                                                                                                   aaCCCTTGGCCATCGAGCTGGTATTCCTCGGCACAAGTGCCTTGCAGAAGAGCCTGGCCGAAGTGCTGCAGGGCGAGCTGGCGAAAATCgg                                                                            >  2:87828/1‑91 (MQ=255)
                                                                                                   aaCCCTTGGCCATCGAGCTGGTATTCCTCGGCACAAGTGCCTTGCAGAAGAGCCTGGCCGAAGTGCTGCAGGGCGAGCTGGCGAAAATCGGCATCCACATCGACCTGCGTGCGGTTGAAGAAGGGGCTTTGGTgc                                >  2:369522/1‑135 (MQ=255)
                                                                                                   aaCCCTTGGCCATCGAGCTGGTATTCCTCGGCACAAGTGCCTTGCAGAAGAGCCTGGCCGAAGTGCTGCAGGGCGAGCTGGCGAAAATCGGCATCCACATCGACCTGCGTGCGGTTGAAGAAGGGGCTTTGGTgc                                >  2:378640/1‑135 (MQ=255)
                                                                                                   aaCCCTTGGCCATCGAGCTGGTATTCCTCGGCACAAGTGCCTTGCAGAAGAGCCTGGCCGAAGTGCTGCAGGGCGAGCTGGCGAAAATCGGCATCCACATCGACCTGCGTGCGGTTGAAGAAGGGGCTTTGGTgc                                >  5:152444/1‑135 (MQ=255)
                                                                                                             cATCGAGCTGGTATTCCTCGGCACAAGTGCCTTGCAGAAGAGCCTGGCCGAAGTGCTGCAGGGCGCGCTGGCGAAAATCGGCATCCACATCGACCTGCGTGCGGTTGAAGAAGGGGCTTTGGTGCGGCGCCAAcg                      >  8:384735/1‑135 (MQ=255)
                                                                                                                  aGCTGGTATTCCTCGGCACAAGTGCCTTGCAGAAGAGCCTGGCCGAAGTGCTGCAGGGCGAGCTGGCGAAAATCGGCATCCACATCGACCTGCGTGCGGTTGAAGAAGGGGCTTTGGTGCGGCGCCAACGCGACg                 >  2:399241/1‑135 (MQ=255)
                                                                                                                   gCTGGTATTCCTCGGCACAAGTGCCTTGCAGAAGAGCCTGGCCGAAGTGCTGCAGGGCGAGCTGGCGAAAATCGGCATCCACATCGACCTGCGTGCGGTTGAAGAAGGGGCTTTGGTGCGGCGCCAACGCGACgg                >  6:417691/1‑135 (MQ=255)
                                                                                                                        tATTCCTCGGCACAAGTGCCTTGCAGAAGAGCCTGGCCGAAGTGCTGCAGGGCGAGCTGGCGAAAATCGGCATCCACATCGACCTGCGTGCGGTTGAAGAAGGGGCTTTGGTGCGGCGCCAACGCGACGGCGa             <  7:127700/133‑1 (MQ=255)
                                                                                                                          ttCCTCGGCACACGGGCCTGGCAGAAGAGCCTGGCCGAAGGGCTGCAGGGCGAGCTGGCGAAAATCGGCATCCACATCGCCCGGGGTGCGGTTGAAGAAGGGGCTTTGGTGCGGCGCCAACGCGACGGCGACTTc         <  8:263010/135‑1 (MQ=255)
                                                                                                                           tCCTGGGCACTTGTGCCTTGCTGACGTGCCTGGCAGAAGTGCTGCAGGGGGAGCTGGCGAAAATGGGAATCCTCCTCGTCCTGCGTGCGGTTGAAGATGGGGCTTTGGTGCGGCGCCAACGCGACGGCGACTTCg        <  5:212446/135‑1 (MQ=255)
                                                                                                                               cGGCACAAGTGCCTTGCAGAAGAGCCTGGCCGAAGTGCTGCAGGGCGAGCTGGCGAAAATCGGCATCCACATCGACCTGCGTGCGGTTGAAGAAGGGGCTTTGGTGCGGCGCCAACGCGACGGCGACTTCGGCa     >  2:335371/1‑134 (MQ=255)
                                                                                                                                ggCACAAGTGCCTTGCAGAAGAGCCTGGCCGAAGTGCTGCAGGGCGAGCTGGCGAAAATCGGCATCCACATCGACCTGCGTGCGGTTGAAGAAGGGGCTTTGGTGCGGCGCCAACGCGACGGCGACTTCGGCATg   >  5:248716/1‑135 (MQ=255)
                                                                                                                                 gCACAAGTGCCTTGCAGAAGAGCCTGGCCGAAGTGCTGCAGGGCGAGCTGGCGAAAATCGGCATCCACATCGACCTGCGTGCGGTTGAAGAAGGGGCTTTGGTGCGGCGCCAACGCGACGGCGACTTCGGCATGa  <  3:124964/135‑1 (MQ=255)
                                                                                                                                  cacaAGTGCCTTGCAGAAGAGCCTGGCCGAAGTGCTGCAGGGCGAGCTGGCGAAAATCGGCATCCACATCGACCTGCGTGCGGTTGAAGAAGGGGCTTTGGTGCGGCGCCAACGCGACGGCg              <  3:439958/122‑1 (MQ=255)
                                                                                                                                  cacaAGTGCCTTGCAGAAGAGCCTGGCCGAAGTGCTGCAGGGCGAGCTGGCGAAAATCGGCATCCACATCGACCTGCGTGCGGTTGAAGAAGGGGCTTTGGTGCGGCGCCAACGCGACGGCg              >  4:439958/1‑122 (MQ=255)
                                                                                                                                    |                                                                                                                                   
TCAAGCCCTACCCTTACGACCCGGCGCGCGCCGCCCATGCCCTTGATGCTGCCGGGTGGACGCTGGCCCCAGGCGCGTCTGTGCGCAGTAAAGGCGGCGAACCCTTGGCCATCGAGCTGGTATTCCTCGGCACAAGTGCCTTGCAGAAGAGCCTGGCCGAAGTGCTGCAGGGCGAGCTGGCGAAAATCGGCATCCACATCGACCTGCGTGCGGTTGAAGAAGGGGCTTTGGTGCGGCGCCAACGCGACGGCGACTTCGGCATGA  >  ACB122/3773583‑3773846

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: