Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 3,697,188 G→A 5.4% P624L (CCG→CTG)  paaZ_2 ← Bifunctional protein PaaZ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1223,697,1880GA 5.4% 245.4 / 3.0 74P624L (CCG→CTG) paaZ_2Bifunctional protein PaaZ
Reads supporting (aligned to +/- strand):  ref base G (31/39);  new base A (2/2);  total (33/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.52e-01

CCCAGCTGATTGGTTACTTCCACATCCCACGCCACCACCCCCTGCGGCTGCCCCAGCGGGCTGGTCTTGCCCTGGTCGATCTTGCGCTTGCAGGTCAGGCGCGCCTGGATGGTATCGCCAATGCCCACCGGGTTGATGAAGCGCAGGGTGTCCAGCCCATAGTTGGCCAATACCGGTCCGGCGCCGGGCGACACGAACAGCCCGGCCGCTGCCGACAGTACGAAGTAACCGTGGGCAATGCGTTTACCGAACTGCGAGGCCTTG  >  ACB122/3697059‑3697322
                                                                                                                                 |                                                                                                                                      
cccAGCTGATTGGTTACTTCCACATCGCACGCCACCACCCCCTGCGGCTGCCCCAGCGGGCGGGTCTTGCCCTGGTCGATCTTGCGCTTGCAGGTCAGGCGCGCCTGGATGGTATCGCCAATGCCCACCGGGtt                                                                                                                                    <  4:203836/134‑1 (MQ=255)
cccAGCTGATTGGTTACTTCCACATCCCACGCCACCACCCCCTGCGGCTGCCCCAGCGGGCTGGTCTTGCCCTGGTCGATCTTGCGCTTGCAGGTCAGGCGCGCCTGGATGGTATCGCCAATGCCCACCGGGtt                                                                                                                                    <  8:117287/134‑1 (MQ=255)
 ccAGCTGATTGGTTACTTCCACATCCCACGCCACCACCCCCTGCGGCTGCCCCAGCGGGCTGGTCTTGCCCTGGTCGATCTTGCGCTTGCAGGTCAGGCGCGCCTGGATGGTATCGCCAATGCCCACCGGGTtga                                                                                                                                  >  8:345743/1‑135 (MQ=255)
      tGATTGGTTACTTCCACATCCCACGCCACCACCCCCTGCGGCTGCCCCAGCGGGCTGGTCTTGCCCTGGTCGATCTTGCGCTTGCAGGTCAGGCGCGCCTGGATGGTATCGCCAATGCCCACCGGGTTGATGAAg                                                                                                                             <  4:102951/135‑1 (MQ=255)
              tACTTCCACCTCCCACGCCACCACCCCCTGCGGCTGCCCCAGCGGGCTGGTCTTGCCCTGGTCGATCTTGCGCTTGCAGGTCAGGCGCGCCTGGATGGTATCGCCAATGCCCACCGGGTTGATGAAGCGCAGGgt                                                                                                                     <  1:324271/135‑1 (MQ=255)
              tACTTCCACATCCCACGCCACCACCCCCTGCGGCTGCCCCAGCGGGCTGGTCTTGCCCTGGTCGATCTTGCGCTTGCAGGTCAGGCGCGCCTGGATGGTATCGCCAATGCCCACCGGGTTGATGAAGCGCAGGgt                                                                                                                     <  3:239222/135‑1 (MQ=255)
              tACTTCCACATCCCACGCCACCACCCCCTGCGGCTGCCCCAGCGGGCTGGTCTTGCCCTGGTCGATCTTGCGCTTGCAGGTCAGGCGCGCCTGGATGGTATCGCCAATGCCCACCGGGTTGATGAAGCGCAGGgt                                                                                                                     <  7:245969/135‑1 (MQ=255)
              tACTTCCACATCCCACGCCACCACCCCCTGCGGCTGCCCCAGCGGGCTGGTCTTGCCCTGGTCGATCTTGCGCTTGCAGGTCAGGCGCGCCTGGATGGTATCGCCAATGCCCACCGGGTTGATGAAGCGCAGGgt                                                                                                                     <  1:363783/135‑1 (MQ=255)
              tACTTCCACATCCCACGCCACCACCCCCTGCGGCTGCCCCAGCGGGCTGGTCTTGCCCTGGTCGATCTTGCGCTTGCAGGTCAGGCGCGCCTGGATGGTATCGCCAATGCCCACCGGGTTGATGAAGCGCAGGgt                                                                                                                     <  2:28252/135‑1 (MQ=255)
              tACTTCCACATCACACGCCACCACCACCTGCGGCAGCCCCAGCGGGCTGGTCTTGCCCTGGTCGAGCTTGCGCTTGCAGGTCAGGCGCGCCTGGATGGTATCGCCAATGCCCACCGGGTTGATGAAGCGCAGGgt                                                                                                                     <  2:299156/135‑1 (MQ=255)
                  tCCACATCCCACGCCACCACCCCCTGCGGCTGCCCCAGCGGGCTGGTCTTGCCCTGGTCGATCTTGCGCTTGCAGGTCAGGCGCGCCGGGATGGTATCGCCAATGCCCACCGGGTTGATGAAGCGCAGGGTGTcc                                                                                                                 >  3:13567/1‑135 (MQ=255)
                    cacaTCCCACGCCACCACCCCCTGCGGCTGCCCCAGCGGGCTGGTCTTGCCCTGGTCGATCTTGCGCTTGCAGGTCAGGCGCGCCTGGATGGTATCGCCAATGCCCACCGGGTTGATGAAGCGCAGGGTGTCCAg                                                                                                               <  3:303190/135‑1 (MQ=255)
                        tCCCACGCCACCACCCCCTGCGGCTGCCCCAGCGGGCTGGTCTTGCCCTGGTCGATCTTGCGCTTGCAGGTCAGGCGCGCCTGGATGGTATCGCCAATGCCCACCGGGTTGATGAAGCGCAGGGTGTCCAGCCCa                                                                                                           <  5:422976/135‑1 (MQ=255)
                        tCCCACGCCACCACCCCCTGCGGCTGCCCCAGCGGGCTGGTCTTGCCCTGGTCGATCTTGCGCTTGCAGGTCAGGCGCGCCTGGATGGTATCGCCAATGCCCACCGGGTTGATGAAGCGCAGGGTGTCCAGCCCa                                                                                                           <  4:141867/135‑1 (MQ=255)
                        tCCCACGCCACCACCCCCTGCGGCTGCCCCAGCGGGCTGGTCTTGCCCTGGTCGATCTTGCGCTTGCAGGTCAGGCGCGCCTGGATGGTATCGCCAATGCCCACCGGGTTGATGAAGCGCAGGGTGTCCAGCCCa                                                                                                           <  4:13567/135‑1 (MQ=255)
                        tCCCACGCCACCACCCCCTGCGGCTGCCCCAGCGGGCTGGTCTTGCCCTGGTCGATCTTGCGCTTGCAGGTCAGGCGCGCCTGGATGGTATCGCCAATGCCCACCGGGTTGATGAAGCGCAGGGTGTCCAGCCCa                                                                                                           <  7:158800/135‑1 (MQ=255)
                        tCCCACGCCACCACCCCCTGCGGCTGCCCCAGCGGGCTGGTCTTGCCCTGGTCGATCTTGCGCTTGCAGGTCAGGCGCGCCTGGATGGTATCGCCAATGCCCACCGGGTTGATGAAGCGCAGGGTGTCCAGCCCa                                                                                                           <  2:388896/135‑1 (MQ=255)
                        tCCCACGCCACCACCCCCTGCGGCTGCCCCAGCGGGCTGGTCTTGCCCTGGTCGATCTTGCGCTTGCAGGTCAGGCGCGCCTGGATGGTATCGCCAATGCCCACCGGGGTGATGAAGCGCAGGGTGTCCAGCCCa                                                                                                           <  7:118677/135‑1 (MQ=255)
                        tCCCACGCCACCACCCCCTGCGGCTGCCCCAGCGGGCTGGTCGTGCCCTGGTCGATCTTGCGCTTGCAGGTCAGGCGCGCCTGGATGGTATCGCCAATGCCCACCGGGTTGATGAAGCGCAGGGTGTCCAGCCCa                                                                                                           <  5:415051/135‑1 (MQ=255)
                        tCCCACGCCACCACCCCCTGCGGCGGCCCCAGCGGGCTGGTCTTGCCCTGGTCGATCTTGCGCTTGCAGGTCAGGCGCGCCTGGATGGTATCGCCAATGCCCACCGGGTTGATGAAGCGCAGGGTGTCCAGCCCa                                                                                                           <  3:311766/135‑1 (MQ=255)
                         cccACGCCACCACCCCCTGCGGCTGCCCCAGCGGGCTGGTCTTGCCCTGGTCGATCTTGCGCTTGCAGGTCAGGCGCGCCTGGATGGTATCGCCAATGCCCACCGGGTTGATGAAGCGCAGGGTGTCCAGCCCAt                                                                                                          <  1:365131/135‑1 (MQ=255)
                              gCCACCACCACCTGCGGCTGCCCCAGCGGGCTGGTCTTGCCCTGGGCGATCTTGCGCTTGCAGGTCTGGCGCGCCTGGATTGTATCGCCAATGCCCACCGGGTTGATGAAGCGCAGGGTGTCCAGCCCATAGTTg                                                                                                     <  1:153250/135‑1 (MQ=255)
                                      ccccTGCGGCTGCCCCAGCGGGCTGGTCTTGCCCTGGTCGATCTTGCGCTTGCAGGTCAGGCGCGCCTGGATGGTATCGCCAATGCCCACCGCGTTGATGAACCGCAGGGTGTCCAGCCCATAGTTGGCCAAtcc                                                                                             >  5:119323/1‑133 (MQ=255)
                                           gCGGCTGCCCCAGCGGGCTGGTCTTGCCCTGGTCGATCTTGCGCTTGCAGGTCAGGCGCGCCTGGATGGTATCGCCAATGCCCACCGGGTTGATGAAGCGCAGGGTGTCCAGCCCATAGTTGGCCAATACCGGTc                                                                                        >  3:59914/1‑135 (MQ=255)
                                             ggCTGCCCCCGCGGGCTGGTCTTGCCCTGGTCGCTCTTGCGCGTGCAGGTCCGGCGCGCCTGGATGGTATCGCCACTGCCCCCCGGGTTGCTGCAGCGCAGGGTGGACCGCCCATAGTTGGCCAATCCCGGTCCg                                                                                      >  6:106419/1‑135 (MQ=255)
                                             ggCTGCCCCAGCGGGCTGGTCTTGCCCTGGTCGATCTTGCGCTTGCAGGTCAGGCGCGCCTGGATGGTATCGCCAATGCCCACCGGGTTGATGAAGCGCAGGGTGTCCAGCCCATAGTTGGCCAATACCGGTCCg                                                                                      >  4:334547/1‑135 (MQ=255)
                                             ggCTGCCCCAGCGGGCTGGTCTTGCCCTGGTCGATCTTGCGCTTGCAGGTCAGGCGCGCCTGGATGGTATCGCCAATGCCCACCGGGTTGATGAAGCGCAGGGTGTCCAGCCCATAGTTGGCCAATACCGGTCCg                                                                                      >  6:192047/1‑135 (MQ=255)
                                             ggCTGCCCCAGCGGGCTGGTCTTGCCCTGGTCGATCTTGCGCTTGCAGGTCAGGCGCGCCTGGATGGTATCGCCAATGCCCACCGGGTTGATGAAGCGCAGGGTGTCAAGCCGATAGTTGGCCAATACAGGTCCg                                                                                      >  8:97477/1‑135 (MQ=255)
                                             ggCTGCCCCAGCGGGCTGGGCTTGCCCTGGTCGATCTTGCGCTTGCAGGTCCGGCGCGCCGGGATGGTATCGCCAATGCCCACCGGGGTGATGAAGCGCAGGGTGTCCAGCCCATAGGTGGCCAAGACCCGTCCg                                                                                      >  1:422475/1‑135 (MQ=255)
                                               cTGCCCCAGCGGGCTGGTCTTGCCCTGGTCGATCTTGCGCTTGCAGGTCAGGCGCGCCTGGATGGTATCGCCAATGCCCACCGGGTTGATGAAGCGCAGGGTGTCCAGCCCATAGTTGGCCAATACCGGTCCGgc                                                                                    >  1:111335/1‑135 (MQ=255)
                                                  ccccAGCGGGCTGGTCTTGCCCTGGTCGATCTTGCGCTTGCAGGTCAGGCGCGGCTGGATGGTATCGCCAATGCCCACCGGGTTGATGAAGCGCAGGGTGTCCAGCCCATAGTTGGCCAATACCGGTCCGGCGcc                                                                                 >  6:63552/1‑135 (MQ=255)
                                                  ccccAGCGGGCTGGTCTTGCCCTGGTCGATCTTGCGCTTGCAGGTCAGGCGCGCCTGGATGGTATCGCCAATGCCCACCGGGTTGATGAAGCGCAGGgt                                                                                                                     <  1:383572/99‑1 (MQ=255)
                                                  ccccAGCGGGCTGGTCTTGCCCTGGTCGATCTTGCGCTTGCAGGTCAGGCGCGCCTGGATGGTATCGCCAATGCCCACCGGGTTGATGAAGCGCAGGgt                                                                                                                     >  2:383572/1‑99 (MQ=255)
                                                   ccgaGCGGGCTGGTCTTGCCCTGGTCGATCTTGCGCTTGCAGGTCAGGCGCGCCTGGATGGTATCGCCAATGCCCACCGGGTTGATGAAGCGCAGGGTGTCGAGCCCTTAGTTGGCCAATACCGGTCCGGCGCCg                                                                                <  6:163129/132‑1 (MQ=255)
                                                     cAGCGGGCTGGTCTTGCCCTGGTCGATCTTGCGCTTGCAGGTCAGGCGCGCCTGGATGGTATCGCCAATGCCCACCGGGTTGATGAAGCGCAGGGTGTCCAGCCCATAGTTGGCCAATACCGGTCCGGCGCCggg                                                                              >  5:95397/1‑135 (MQ=255)
                                                     cAGCGGGCTGGTCTTGCCCTGGTCGATCTTGCGCTTGCAGGTCAGGCGCGCCTGGATGGTATCGCCAATGCCCACCGGGTTGATGAAGCGCAGGGTGTCCAGCCCATAGTTGGCCAATACCGGTCCGGCGCCggg                                                                              >  5:14596/1‑135 (MQ=255)
                                                     cAGCGGGCTGGTCTTGCCCTGGTCGATCTTGCGCTTGCAGGTCAGGCGCGCCTGGATGGTATCGCCAATGCCCACCGGGTTGATGAAGCGCAGGGTGTCCAGCCCATAGTTGGCCAATACCGGTCCGGCGCCggg                                                                              <  4:219699/135‑1 (MQ=255)
                                                     cAGCGGGCTGGTCTTGCCCTGGTCGATCTTGCGCTTGCAGGTCAGGCGCGCCTGGATGGTATCGCCAATGCCCACCGGGGTGATGAAGCGCAGGGTGTCCAGCCCATAGTTGGCCAAGACCGGTCCGGCGCCggg                                                                              >  4:216954/1‑135 (MQ=255)
                                                                tCTTGCCCTGGTCGATCTTGCGCTTGCAGGTCAGGCGCGCCTGGATGGTATCGCCAATGCCCACCGGGTTGATGAAGCGCAGGGTGTCCAGCCCATAGTTGGCCAATACCGGTCCGGCGCCGGGCGACACGAACa                                                                   >  2:148368/1‑135 (MQ=255)
                                                                tCTTCCCCTGGTCCATCTAGCGCTTGCAGGTCAGGCGCGCCTGGATGGTATCGCAAATTCCCACAGGGTTGATGAAGCGCAGGGTGACCAGCCCATAGTTGGCCAATACAGGTCCGGCGCCGGACGACACGAAcc                                                                   >  2:142906/1‑134 (MQ=255)
                                                                      ccTGGTCGATCTTGCGGTTGCAGGTCAGGCGCGCCTGTATGGTATCCCCAATGCCCACCGGGTTGATGAAGCGCAGGGTGTCCCCCCCATAGTTGGCCACTACCGGTCCGGCGCCGGGCGACACGAACAGCCCgg                                                             <  2:400437/135‑1 (MQ=255)
                                                                            cGATCTTGCGCTTGCAGGTCTGGCGCGCCTGGAGGGTATCGCCAATGCCCACCGGGTTGATGACGCGCTGGGTGTCCAGCCCATAGTTGGCCAATACCGGTCCGGCGCCGGGCGACACGAACAGCCCGGCCGCTg                                                       <  3:420958/135‑1 (MQ=255)
                                                                            cGATCTTGCGCTTGCAGGTCAGGCGCGCCTGGATGGTATCGCCAATGCCCACCGGGTTGATGAAGCGCAGGGTGTCCAGCCCATAGTTGGCCAATACCGGTCCGGCGCCGGGCGACACGAACAGCCCGGCCGCTg                                                       >  5:461855/1‑135 (MQ=255)
                                                                            cGATCTTGCGCTTGCAGGTCAGGCGCGCCTGGATGGTATCGCCAATGCCCACCAGGTTGATGAAGCGCAggg                                                                                                                      <  8:86851/72‑1 (MQ=255)
                                                                            cGATCTTGCGCTTGCAGGTCAGGCGCGCCTGGATGGTATCGCCAATGCCCACCAGGTTGATGAAGCGCAggg                                                                                                                      <  8:448913/72‑1 (MQ=255)
                                                                            cGATCTTGCGCTTGCAGGTCAGGCGCGCCTGGATGGTATCGCCAATGCCCACCAGGTTGATGAAGCGCAggg                                                                                                                      >  7:86851/1‑72 (MQ=255)
                                                                            cGATCTTGCGCTTGCAGGTCAGGCGCGCCTGGATGGTATCGCCAATGCCCACCAGGTTGATGAAGCGCAggg                                                                                                                      >  7:448913/1‑72 (MQ=255)
                                                                             gATCTTGCGCTTTCAGGTCAGGAGCGCCTGGATGGTATCGCCAATGCCCACCGGGTTGATGAAGCGCTCGGTGTCCAGCCCATAGTTGGCCAATACCGGTACGGCGCCGGGCGACACGAACAGCCCTGCCGCTGc                                                      >  5:55018/1‑135 (MQ=255)
                                                                             gATCTTGCGCTTGCAGGTCAGGCGCGCCTGGATGGTATCGCCAATGCCCACCGGGTTGATGAAGCGCAGGGTGTCCAGCCCATAGTTGGCCAATACCGGTCCGGCGCCGGGCGACACGAACAGCCCGGCCGCTGc                                                      >  8:141268/1‑135 (MQ=255)
                                                                             gATCTTGCGCTTGCAGGTCAGGCGCGCCTGGATGGTATCGCCAATGCCCACCGGGTTGATGAAGCGCAGGGTGTCCAGCCCATAGTTGGCCAATACCGGTCCGGCGCCGGGCGACACGAACAGCCCGGCCGCTGc                                                      <  7:345743/135‑1 (MQ=255)
                                                                             gATCTTGCGCTTGCAGCTCAGGCGCGCCTGGATGGTATCGCCAATGCCCACCGGCTTGATGAAGCGCAGGGTGTCCAGCCCATAGTTGGCCAATACCGCTCCGCCGCCGGGCGCCACGAACAGCCCGGCCGCTGc                                                      >  4:162816/1‑135 (MQ=255)
                                                                             gATCTTGCGCTGGCAGGGCAGGCGCGCCTGGATGGTATCGCCAATGCCCACCGGGTTGATGAAGCGCAGGGTGTCCAGCCCATAGTTGGCCAATACCGGTCCGGCGCCGGGCGACACGAACAGCCCGGCCGCTGc                                                      <  3:346772/135‑1 (MQ=255)
                                                                                  tGCGCTTGCAGGTCAGGCGCGCCTGGATGGTATCGCCAATGCCCACCGGGTTGATGAAGCGCAGGgt                                                                                                                     >  6:68933/1‑67 (MQ=255)
                                                                                  tGCGCTTGCAGGTCAGGCGCGCCTGGATGGTATCGCCAATGCCCACCGGGTTGATGAAGCGCAGGgt                                                                                                                     <  5:68933/67‑1 (MQ=255)
                                                                                    cgcTTGCAGGTCAGGCGCGCCTGGATGGGATCGCCAATGCCCACCGGGTTGATGAAGCGCAGGGTGTCCAGCCCATAGTTGGCCAATACCGGTCCGGCGCCGGGCGACACGAACAGCCCGGCCGCTGCCGAc                                                  <  7:105288/132‑1 (MQ=255)
                                                                                    cgcTTGCAGGTCAGGCGCGCCTGGAGGGTCTCGCCAATGCCCACCGGGTTGATGAAGCGCAGGGTGTCCAGCCCATAGTTGGCCAATACCGGTCCGGCGCCGGGCGACACGAACAGCCCGGCCGCTGCCGACa                                                 <  4:231906/133‑1 (MQ=255)
                                                                                    cgcTTGCAGGGCAGGCGCGCCTGGATGGTATCGCCAATGCCCACCGGGGTGATGAAGCGCAGGGTGTCCAGCCCATAGTTGGCCAATACCGGTCCGGCGCCGGGCGACACGAACAGCCCGGCCGCTGCCGACa                                                 <  5:368569/133‑1 (MQ=255)
                                                                                       ttGCCGGTCAGGCGCGCCTGGATGGTATCGCCAATGCCCACCGGGGTGATGAAGCGCAGGGTGTCCAGCCCATAGTTGGCCAATACCGGTCCGGCGCCGGGCGACACGAACAGCCCGGCCGCTGCCGACAGTACg                                            <  7:387931/135‑1 (MQ=255)
                                                                                       ttGCAGGGCAGGCGCGCCTGGATGGTATCGCCAATGCCCACCGGGTTGATGAAGCGCAGGGTGTCCAGCCCATAGTTGGCCAATACCGGTCCGGCGCCGGGCGACACGAACAGCCCGGCCGCTGCCGACACTACg                                            <  1:386019/135‑1 (MQ=255)
                                                                                           aGGTCAGGCGCGCCTGGATGGTATCGCCAATGCCCACCGGGTTGATGAAGCGCAGGGTGTCCAGCCCATAGTTGGCCAATACCGGTCCGGCGCCGGGCGACACGAACAGCCCGGCCGCTGCCGACAGTACGAAGt                                        <  2:111335/135‑1 (MQ=255)
                                                                                            ggTCAGGCGCGCCTGGATGGTATCGCCAATGCCCACCGGGTTGATGAAGCGCATGGTGTCCAGCCCATAGTTGGCCAATACCGGTCCGGCGCCGGGCGACACGAACAGCCCGGCCGCTGCCGACAGTCCGAAGTa                                       >  2:352853/1‑135 (MQ=255)
                                                                                            ggTCAGGCGCGCCTGGATGGTATCGCCAATGCCCACCGGGTTGATGAAGCGCAGGGTGTCCAGCCCATAGTTGGCCAATACCGGTCCGGCGCCGGGCGACACGAACAGCCCGGCCGCTGCCGACAGTACGAAGTa                                       >  1:67762/1‑135 (MQ=255)
                                                                                                    gcgcGTGGATGGTATCGCCAATGCCCACCGGGTTGATGAAGCGCAGGGTGTCCAGCCCATAGTTGGCCAATACCGGTCCGGCGCCGGGCGACACGAACAGCCCGGCCGCTGCCGACAGTACGAAGTAACCGTggg                               <  2:67762/135‑1 (MQ=255)
                                                                                                            aTGGTATCGCCAATGCCCACCGGGTTGATGAAGCGCAGGGTGTCCAGCCCATAGTTGGCCAATACCGGTCCGGCGCCGGGCGACACGAACAGCCCGGCCGCTGCCGACAGTACGAAGTAACCGTGGGCAATGCGt                       <  1:148368/135‑1 (MQ=255)
                                                                                                                tATCGCCAATGCCCACCGGGTTGATGAAGCGCAGGGTGTCCAGCCCATAGTTGGCCAATACCGGTCCGGCGCCGGGCGACACGAACAGCCCGGCCGCTGCCGACAGTACGAAGTAACCGTGGGCAATGCGTTTa                    <  4:59914/134‑1 (MQ=255)
                                                                                                                   cGCCAATGCCCACCGGGTTGATGAAGCGCAGGGTGTCCAGCCCATAGTTGGCCAATACCGGTCCg                                                                                      >  3:86609/1‑65 (MQ=255)
                                                                                                                   cGCCAATGCCCACCGGGTTGATGAAGCGCAGGGTGTCCAGCCCATAGTTGGCCAATACCGGTCCg                                                                                      <  4:86609/65‑1 (MQ=255)
                                                                                                                   cGCCAATGCCCACCGGGTTGATGAAGCGCAGGGTGTCCAGCCCATAGTTGGCCAATACCGGTCCGGCGCCGGGCGCCACGAACAGCCCGGCCGCTGCCGACAGTACGAAGTAACCGTGGGCAATGCGTTTACCGa                >  4:403632/1‑135 (MQ=255)
                                                                                                                   cGCCAATGCCCACCGGGTTGATGAAGCGCAGGGTGTCCAGCCCATAGTTGGCCAATACCGGTCCGGCGCCGGGCGACACGAACAGCCCGGCCGCTGCCGACAGTACGAAGTAACCGTGGGCAATGCGTTTACCgc                >  2:17720/1‑134 (MQ=255)
                                                                                                                   cGCCAATGCCCACCGGGTTGATGAAGCGCAGGGTGTCCAGCCCATAGTTGGCCAATACCGGTCCGGCGCCGCGCGCCACGAACAGCCCGTCCGCTGCAGACAGTACGAAGTAACCGTGGGCAATGCGTTTACCGa                >  4:163885/1‑135 (MQ=255)
                                                                                                                     ccAATGCCCACCGGGTTGATGAAGCGCAGGGTGTCCAGCCCATAGTTGGCCAATACCGGTCCGgcgc                                                                                  >  7:395173/1‑67 (MQ=255)
                                                                                                                     ccAATGCCCACCGGGTTGATGAAGCGCAGGGTGTCCAGCCCATAGTTGGCCAATACCGGTCCGgcgc                                                                                  <  8:395173/67‑1 (MQ=255)
                                                                                                                        aTGCCCACCGGGTTGATGAAGCGCAGGGTGTCCAGCCCATAGTTGGCCAATACCGGTCCGGCGCCGGGCGACACGAACAGCCCGGCCGCGGCCGACAGTACGAAGTACCCGGGGGCAATGCGTTTACCGAACTGc           >  1:230282/1‑135 (MQ=255)
                                                                                                                            ccACCGGGTTGATGAAGCGCAGGGTGTCCAGCCCATAGTTGGCCAATACCGGTCCGGCGCCGGGCGCCACGAACAGCCCGGCCGCTGCCGACAGTACGAAGTAACCGTTGGCAATGCGTTTACCGAACTGCGAgg       >  1:400382/1‑135 (MQ=255)
                                                                                                                                 gggTTGATGAAGCGCCGGGTGTACAGCCCATAGATGGCCAATCCCGGGCCGGCGCCGGGCGACACGAACAGCCCGGCCGCGGCCGACAGTACGAAGTAACCGTGGGCAATGCGTTTACCGAACTGCGAGGCCTTg  <  2:339190/135‑1 (MQ=255)
                                                                                                                                 gggTTGATGAAGCGCAGGGTGTCCAGCCCATAGTTGGCCAATACCGGTCCGGCGCCGGGCGACACGAACAGCCCGGCCGCTGCCGCCAGTACGAAGTAACCGTGGGCAATGCGTTTACCGAACTGCGAGGCCTTg  >  1:144722/1‑135 (MQ=255)
                                                                                                                                 gggTTGATGAAGCGCAGGGTGTCCAGCCCATAGTTGGCCAATACCGGTCCGGCGCCGGGCGACACGAACAGCCCGGCCGCTGCCGACAGTACGAAGTAACCGTGGGCAATGCGTTTACCGAACTGCGAGGCCTTg  >  5:302631/1‑135 (MQ=255)
                                                                                                                                 |                                                                                                                                      
CCCAGCTGATTGGTTACTTCCACATCCCACGCCACCACCCCCTGCGGCTGCCCCAGCGGGCTGGTCTTGCCCTGGTCGATCTTGCGCTTGCAGGTCAGGCGCGCCTGGATGGTATCGCCAATGCCCACCGGGTTGATGAAGCGCAGGGTGTCCAGCCCATAGTTGGCCAATACCGGTCCGGCGCCGGGCGACACGAACAGCCCGGCCGCTGCCGACAGTACGAAGTAACCGTGGGCAATGCGTTTACCGAACTGCGAGGCCTTG  >  ACB122/3697059‑3697322

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: