Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 1,443,010 T→C 6.0% T17A (ACC→GCC)  prlF ← Antitoxin PrlF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1221,443,0100TC 6.0% 213.7 / 3.9 82T17A (ACC→GCC) prlFAntitoxin PrlF
Reads supporting (aligned to +/- strand):  ref base T (28/47);  new base C (2/3);  total (32/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.55e-01

CGTCCTCGGTTGCCTCTGCGCGGGTGAGTATCACCCCACCGTCGGGGCGAATCGAGTAATGGATCTTGTCACGTTTTTTCAAGCCCAAGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGA  >  ACB122/1442893‑1443143
                                                                                                                     |                                                                                                                                     
cGTCCTCGGTTGCCTCTGCGCGGGTGAGTATCACCCCACCGTCGGGGCGAATCGAGTAATGGATCTTGTCACGTTTTTTCAAGCCCAACGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAg                                                                                                                     >  8:213419/1‑136 (MQ=255)
     tCGGTTGCCTCTGCGCGGGTGAGTATCACCCCACCGTCGGGGCGAATCGAGTAATGGATCTTGTCACGTTTTTTCAAGCCCACGGCGCGCCGCCCCGTGTCGGGAACAGGGGTTGGATTGCGCTCGGTCAGCGTg                                                                                                                 >  6:256651/1‑135 (MQ=255)
     tCGGTTGCCTCTGCGCGGGTGAGTATCACCCCACCGTCGGGGCGAATCGAGTAATGGATCTTGTCACGTTTTTTCAAGCCCAAGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGt                                                                                                                  >  5:284554/1‑134 (MQ=255)
     tCGGTTGCCTCTGCGCGGGTGAGTATCACCCCACCGTCGGGGCGAATCGAGTAATGGATCTTGTCACGTTTTTTCAAGCCCAAGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTgg                                                                                                                >  3:315500/1‑136 (MQ=255)
     tCGGTTGCCTCTGCGCGGGTGAGTATCACCCCACCGTCGGGGCGAATCGAGTAATGGATCTTGTCACGTTTTTTCAAGCCCAAGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTg                                                                                                                 >  2:296605/1‑135 (MQ=255)
     tCGGTTGCCTCTGCGCGGGTGAGTATCACCCCACCGTCGGGGCGAATCGAGTAATGGATCTTGTCACGTTTTTTCAAGCCCAAGGCGCGCCGCACCGTTTCGGGAACAGTGGCTTGATAGCGATCGGTCAGCGTgg                                                                                                                >  7:14704/1‑136 (MQ=255)
     tCGGTTGCCTCTGCGCGGGTGAGTATCACCCCACCGTCGGGGCGAATCGAGTAATGGATCTTGTCACGTTTTTTCAAGCCCAAGGCGCGCCGCACCGTTTCGGGAACAGTGGCTTGATAGCGATCGGTCAGCGTg                                                                                                                 >  5:394107/1‑135 (MQ=255)
     tCGGTTGCCTCTGCGCGGGTGAGTATCACCCCACCGTCGGGGCGAATCGAGTAATGGATCTTGTCACGTTTTTTCAAGACCAAGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTgg                                                                                                                >  3:288236/1‑136 (MQ=255)
       ggTTGCCTCTGCGCGGGTGAGTATCACCCCACCGTCGGGGCGAATCGAGTAATGGATCTTGTCACGTTTTTTCAAGCCCAAGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGAc                                                                                                              >  8:91880/1‑136 (MQ=255)
              caggcgcGGGGGAGTATCCCCCCACGGTGGGGGCTAATCGAGTAATGGATCTTGTCAGGTTTTTTCACGCCCCAGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACg                                                                                                        <  5:252083/132‑1 (MQ=255)
               ctGCGCGGGTGAGTATCACCCCACCGTCGGGGCGAATCGAGTAATGGATCTTGTCACGTTTTTTCAAGCCCAAGGCGCTCCGCACCGTTTCGGGAACAGTGGGTTGATAGAGATCGGTCAGCGTGGACTCGACGt                                                                                                       >  1:231300/1‑135 (MQ=255)
               ctGCGCGGGTGAGTATCACCCCACCGTCGGGGCGAATCGAGTAATGGATCTTGTCACGTTTTTTCAAGCCCAAGGCGCGCCGCACCGTTTCGGGATCAGTGGTTTGATAGCGATCGGTCAGCGGGGACTCGACGt                                                                                                       >  6:365435/1‑135 (MQ=255)
               ctGCGCGGGTGAGTATCACCCCACCGTCGGGGCGAATCGAGTAATGGATCTTGTCACGTTTTTTCAAGCCCAAGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGt                                                                                                       >  2:338036/1‑135 (MQ=255)
                   gcgGGGGAGTATCACCCCACCGTCGGGGCGAATCGAGTAATGGATCTTGTCACGTTTTTTCAAGCCCAAGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTc                                                                                                      <  2:294890/132‑1 (MQ=255)
                    cgGGTGAGTATCACCCCACCGTCGGGGCGAATCGAGTAATGGATCTTGTGACGTTTTTTCAAGCCCAAGGCGCGCCGCACCGGTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAgg                                                                                                 >  3:247837/1‑136 (MQ=255)
                    cgGGTGAGTATCACCCCACCGTCGGGGCGAATCGAGTAATGGATCTTGTCACGTTTTTTCAAGCCCAAGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAgg                                                                                                 <  7:75556/136‑1 (MQ=255)
                    cgGGTGAGTATCACCCCACCGTCGGGGCGAATCGAGTAATGGATCTTGTCACGTTTTTTCAAGCCCAAGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAgg                                                                                                 <  7:91880/136‑1 (MQ=255)
                     gggTGAGTATCACCCCACCGTCGGGGCGAATCGAGTAATGGATCTTGTCACGTTTTTTCAAGCCCAAGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAgg                                                                                                 <  6:266522/135‑1 (MQ=255)
                     gggTGAGTATCACCCCACCGTCGGGGCGAATCGAGTAATGGATCTTGTCACGTTTTTTCAAGCCCAAGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAgg                                                                                                 <  1:131252/135‑1 (MQ=255)
                      ggTGAGTATCACCCCACCGTCGGGGCGAATCGAGTAATGGATCTTGTCACGTTTTTTCAAGCCCAAGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTg                                                                                               <  3:134702/136‑1 (MQ=255)
                         caGTAGCACCCCGCCGGCGGGGCGAATCGAGTAATTGATCTTGTCACGTTTTTTCAAGCCCAAGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGTACTCGACGTCGAAGGGGGCa                                                                                            <  8:19937/135‑1 (MQ=255)
                          aGTATCACCCCACCGTCGGGGCGAATCGAGTAATGGATCTTGTCACGTTTTTTCAAGCCCAAGGCGCGGCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCa                                                                                            <  5:284684/135‑1 (MQ=255)
                          aGTATCACCCCACCGTCGGGGCGAATCGAGTAATGGATCTTGTCACGTTTTTTCAAGCCCAAGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCa                                                                                            <  1:296605/135‑1 (MQ=255)
                          aGTATCACCCCACCGTCGGGGCGAATCGAGTAATGGATCTTGTCACGTTTTTTCAAGCCCAAGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCa                                                                                            >  6:90268/1‑135 (MQ=255)
                          aGTATCACCCCACCGTCGGGGCGAATCGAGTAATGGATCTTGTCACGTTTTTTCAAGCCCAAGGCGCGCCGCACAGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCa                                                                                            <  2:371025/135‑1 (MQ=255)
                          aGTATCACCCCACCGTCGGGGCGAATCGAGTAATGGATCTTGTCACGTTTTTGCAAGCCCAAGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCa                                                                                            >  1:221037/1‑135 (MQ=255)
                          aGTATCACCCCACCGTCGGGGCGAATCGAGTAATGGATCTTGTCACGTTGTTTCAAGCCCAAGGCGCGCCGGACCGTTTCGGGAACAGTTGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCa                                                                                            <  4:325726/135‑1 (MQ=255)
                            tATCACCCCACCGTCGGGGCGAATCGAGTAATGGATCTTGTCACGTTTTTTCCAGCCCAAGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGc                                                                                          <  2:51619/135‑1 (MQ=255)
                               cACCCCACCGTCGGGGCGAATCGAGTAATGGATCTTGTCACGGTTTTTGAAGCCCAAGGCGCGCCGCCCCGGTTCGGGCACCGGGGGTTGCTAGCGATAGGGGCGCGTGGACTCGACCTCGAAGGTGgccgccgg                                                                                       >  6:252083/1‑133 (MQ=255)
                                 ccccACCGTCGGGGCGAATCGAGTAATGGATCTTGTCACGTTTTTTCAAGCCCAAGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGa                                                                                     <  5:90268/135‑1 (MQ=255)
                                  cccTCCGTGGGGCCGAATCGAGTAATGGATCTTGTCACGTTTTTTCTAGCCCTAGGCGCGCGGCTCCGTTTCTGGAAGAGTGGTTTGTTAGGGAGCGGTCAGCGTGGACTCGACGTCGATGGTGGCAGCCATTAg                                                                                    <  2:168042/135‑1 (MQ=255)
                                  cccACCGTCGGGGCGAATCGAGTAATGGATCTTGTCACGTTTTTTCAAGCCCAAGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAg                                                                                    <  1:52419/135‑1 (MQ=255)
                                    cACCGTCGGGGCGAATCGAGTAATGGATCTTGTCACGTTTTTTCAAGCCCAAGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTaa                                                                                 >  4:146851/1‑136 (MQ=255)
                                     aCCGTCGGGGCGAATCGAGTAATGGATCTTGTCACGTTTTTTCAAGCCCAAGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAAt                                                                                <  7:213419/136‑1 (MQ=255)
                                       cGTCGGGGCGAATCGAGTAATGGATCTTGTCACGTTTTTTCAAGCCCCAGGCGCGCCGCCCCGTTTCGGGAACAGTGGTGTGATAGCGATCGGTCCGCGTGGCCTCGACGTCGAAGGTGGCAGCCATGAGTAATg                                                                               >  1:155092/1‑135 (MQ=255)
                                       cGTCGGGGCGAATCGAGTAATGGATCTTGTCACGTTTTTTCAAGCCCAAGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGt                                                                              <  3:75843/136‑1 (MQ=255)
                                         tCGGGGCGAATCGAGTAATGGATCTTGTCACGTTTTTTCAAGCCCAAGGCGCGCCGCTCCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGCCTCGACGTCGAAGGTGGCAGCCATGAGTAATGTcc                                                                            >  3:294568/1‑136 (MQ=255)
                                         tCGGGGCGAATCGAGTAATGGATCTTGTCACGTTTTTTCAAGCCCAAGGCGCGCCGCACCGTTTCGGGAACAGTGGttt                                                                                                                                     >  5:325255/1‑79 (MQ=255)
                                         tCGGGGCGAATCGAGTAATGGATCTTGTCACGTTTTTTCAAGCCCAAGGCGCGCCGCACCGTTTCGGGAACAGTGGttt                                                                                                                                     <  6:325255/79‑1 (MQ=255)
                                         tCGGGGCGAATCGAGTAATGGATCTTGTCACGTTTTTTCAAGCCCAAGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTc                                                                             >  6:134318/1‑135 (MQ=255)
                                            gggCGAATCGAGTAATGGATCTTGTCACGTTTTTTCAAGCCCAAGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCtgt                                                                         >  8:213984/1‑136 (MQ=255)
                                            gggCGAATCGAGTAATGGATCTTGTCACGTTTTTTCAAGCCCAAGCCGCGCCGCACCGTTTCGGGAACAGAGGTTTGATAGCGATCGGTCAGCGTAGACTCTACGTCCAAGGTGGCAGCCATGACTAATGTCgtgt                                                                         >  7:148557/1‑136 (MQ=255)
                                             ggCGAATCGAGTAATGGATCTTGTCACGTTTTTTCAAGCCCAAGGCGCGCCGCACCGTTTCGGGAACAGTGGCTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCtgtg                                                                        <  8:14704/136‑1 (MQ=255)
                                             ggCGAAGCGAGTAATGGATCTTGTCACGTTTTTTCATGCCCATGGCGCGCCGCACCGTTTCGGGAACATGGGTGTGATAGCGATCGGTCAGCGTGGACTCGACGTCGACGTGGGCAGCCATGAGTAATGTCCtgtg                                                                        <  8:199267/136‑1 (MQ=255)
                                              gCGAATCGAGTAATGGATCTTGTCACGTTTTTTCAAGCCCAAGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCtgtg                                                                        <  6:284554/135‑1 (MQ=255)
                                              gCGAATCGAGTAATGGATCTTGTCACGTTTTTTCAAGCCCAAGGCGCGCCGCACCGTTTCGGGAACAGTGGCTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCtgtg                                                                        <  6:394107/135‑1 (MQ=255)
                                              gCGAATCGAGTAATGGATCTTGTCACGTTTTTTCAAGCCCAAGGCGCGCCCCACCGTTTCCGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTgg                                                                       >  8:333680/1‑136 (MQ=255)
                                               cGAATCGAGTAATGGATCTTGTCACGTTTTTTCAAGCCCAAGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATc                                                                                                                           >  4:220623/1‑83 (MQ=255)
                                               cGAATCGAGTAATGGATCTTGTCACGTTTTTTCAAGCCCAAGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATc                                                                                                                           <  3:220623/83‑1 (MQ=255)
                                               cGAATCGAGTAATGGATCTTGTCACGTTTTTTCAAGCCCAAGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTgg                                                                       >  1:387805/1‑135 (MQ=255)
                                                      aGTAATGGATCTTGTCACGTTTTTTCAAGCCCAAGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCa                                                                                            <  4:2697/107‑1 (MQ=255)
                                                      aGTAATGGATCTTGTCACGTTTTTTCAAGCCCAAGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCa                                                                                            >  3:2697/1‑107 (MQ=255)
                                                       gTAATGGATCTTGTCACGTTTTTTCAAGCCCAAGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCg                                                              <  8:175229/136‑1 (MQ=255)
                                                       gTAATGGATCTTGTCACGTTTTTTCAAGCCCAAGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCg                                                              <  4:226299/136‑1 (MQ=255)
                                                       gTAATGGATCTTGTCACGTTTTTTCAAGCCCAAGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGGCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCg                                                              <  4:144854/136‑1 (MQ=255)
                                                        tAATGGATCTTGTCACGTTTTTTCAAGCCCAAGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCg                                                              <  1:105537/135‑1 (MQ=255)
                                                        tAATGGATCTTGTCACGTTTTTTCAAGCCCAAGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCg                                                              <  1:33144/135‑1 (MQ=255)
                                                        tAATGGATCTTGTCACGTTTTTTCAAGCCCAAGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCg                                                              <  5:65388/135‑1 (MQ=255)
                                                        tAATGGATCTTGTCACGTTTTTTCAAGCCCAAGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCg                                                              <  6:1694/135‑1 (MQ=255)
                                                              tgcTTGTCAGGTTTTTTACCGCCCCCGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACa                                                         <  5:256651/132‑1 (MQ=255)
                                                               tCTTGTCACGTTTTTTCAAGCCCAAGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAAc                                                       <  6:211341/135‑1 (MQ=255)
                                                                cTTGTCACGTTTTTTCCAGCCCAAGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAAc                                                       <  1:336956/134‑1 (MQ=255)
                                                                   gTCACGTTTTTTCAAGCCCAAGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCaa                                                   >  5:268733/1‑135 (MQ=255)
                                                                   gTCACGTTTTTTCAAGCCCAAGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCaa                                                   <  8:300253/135‑1 (MQ=255)
                                                                    tCACGTTTTTTCAAGCCCAATGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGTTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATTAGTAATGTCCTGTGGGTGACAGCGATACAAc                                                       <  6:204048/130‑1 (MQ=255)
                                                                      aCGTTTTTTCAAGCCCAAGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCa                                                                                            <  3:36479/91‑1 (MQ=255)
                                                                      aCGTTTTTTCAAGCCCAAGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCa                                                                                            >  4:36479/1‑91 (MQ=255)
                                                                               cAAGCCCAAGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGatt                                       <  4:288236/135‑1 (MQ=255)
                                                                                    ccAAGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGattgcatt                                  <  5:336008/135‑1 (MQ=255)
                                                                                        ggCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCt                              <  6:190248/135‑1 (MQ=255)
                                                                                        ggCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCt                              <  1:311583/135‑1 (MQ=255)
                                                                                         gcgcgcCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTg                            >  7:217769/1‑136 (MQ=255)
                                                                                          cgcgcCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGattgcat                                   >  6:213125/1‑128 (MQ=255)
                                                                                          cgcgcCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGattgcat                                   <  5:213125/128‑1 (MQ=255)
                                                                                             gccgcACCGTTTCGGGAACAGTGGTTTGATAGCGAGCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCaa                        <  8:295698/136‑1 (MQ=255)
                                                                                              cggcACCGTTTCGGGAACAGTGGTTTGATAGCGAGCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCaa                        <  6:127354/133‑1 (MQ=255)
                                                                                              ccgcACCGTTTCGGGAACAGTGGTGTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCaa                        <  1:95744/135‑1 (MQ=255)
                                                                                               cgcACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATg                      >  4:271664/1‑136 (MQ=255)
                                                                                               cgcACCGTTTCGGGAACAGTGGCTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAAt                       <  1:338036/135‑1 (MQ=255)
                                                                                                  aCCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCATGGTAATGCGATTGCATTGCCTTGTCAATGCAg                   <  8:297152/136‑1 (MQ=255)
                                                                                                           ggAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTa           <  6:268733/135‑1 (MQ=255)
                                                                                                            gAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGTAGGTGCATTAcc         <  3:236036/136‑1 (MQ=255)
                                                                                                            gAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTAcc         <  8:309457/136‑1 (MQ=255)
                                                                                                              aCAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCa        >  5:127055/1‑135 (MQ=255)
                                                                                                                     tttGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGtga  >  6:135843/1‑134 (MQ=255)
                                                                                                                     |                                                                                                                                     
CGTCCTCGGTTGCCTCTGCGCGGGTGAGTATCACCCCACCGTCGGGGCGAATCGAGTAATGGATCTTGTCACGTTTTTTCAAGCCCAAGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGA  >  ACB122/1442893‑1443143

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: