Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 1,408,984 G→C 5.5% G101R (GGC→CGC)  bamC → Outer membrane protein assembly factor BamC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1221,408,9840GC 5.5% 212.0 / 3.9 73G101R (GGC→CGC) bamCOuter membrane protein assembly factor BamC
Reads supporting (aligned to +/- strand):  ref base G (36/33);  new base C (2/2);  total (38/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.79e-01

AACGCCTTGACCCGCTGCTGCCGATTCCGCGTAACGTCGCTGACGATACCGTCAAGGGTGAGTTTGAAGTGCCTCGCCCGCAACCGTTGACGGGCGGTGCCGCCCAGGTCACCGACTACAGCCTGCAGCGCAGCGGCAGCAGCCGTTGGGTGCTGGCCCAGCATTCGCCAGCCGAAGTGTGGCCGGTTGCCCGGCAGTTCTTCGAGGACAACGGCTTCCGCATTGCTGAGGAGCGCCCGCAGACGGGTGAGTTCAATACCGCCT  >  ACB122/1408850‑1409113
                                                                                                                                      |                                                                                                                                 
aaCGCCTTGACCCGCTGCTGCCGATTCCGCGTAACGTCGCTGACGATACCGTCAAGGGTTAGTTTGAAGTGCCTCGCCCGCAACCGTTGACGGGCGGTGCCGCCCAGGTCACCGACTACAGCCTGCAGCGCAGCgg                                                                                                                                  <  8:73629/136‑1 (MQ=255)
 aCGCCTTGACCCGCTGCTGCCGATTCCGCGTAACGTCGCTGACGATACTGTCAAGGGTGAGTTTGAAGTGCCTCGCCCGCAACCGTTGACGGGCGGTGCCGCCCAGGTCACCGACTACAGCCTGCAGCGCAGCgg                                                                                                                                  <  4:27558/135‑1 (MQ=255)
         ccccGCGGCGGCCGAATCCCCGTAACGTCCCGGACGATACCGTCATGGGTGAGTTTGATGTGCCTCGCCCGCAACCGTTGACGGGCGGTGCCGCCCAGGTCACCGACTACAGCCTGCAGCGCAGCGGCAGCAGcc                                                                                                                          <  2:180103/134‑1 (MQ=255)
         aCCCGCGGCGGCCGATCCCGCGTAACGGCGCTGACGATACCGTCAAGGGTGAGTTTGAAGTGCCTCGCCCGCAACCGATGACGGGCGGTGCCGCCCAGGTCACCGACTACAGCCTGCAGCGCAGCGGCagcagc                                                                                                                           <  1:103491/134‑1 (MQ=255)
          cccGCTGCTGCCGATTCCGCGTAACGTCGCTGACGATACCGTCAAGGGTGAGTTTGAAGTGCCTCGCCCGCAACCGTTGACGGGCGGTGCCGCCCAGGTCCCCGACTACAGCCTGCAGCGCAGCGGCAGCAGCCGt                                                                                                                        >  8:144976/1‑136 (MQ=255)
          cccGCTGCTGCCGATTCCGCGTAACGTCGCTGACGATACCGTCAAGGGTGAGGTTGAAGTGCCTCGCCCGCAACCGTTGACGGGCGGTGCCGCCCAGGTCACCGACTACAGCCTGCAGCGCAGCgg                                                                                                                                  >  1:177926/1‑126 (MQ=255)
          cccGCTGCTGCCGATTCCGCGTAACGTCGCTGACGATACCGTCAAGGGTGAGGTTGAAGTGCCTCGCCCGCAACCGTTGACGGGCGGTGCCGCCCAGGTCACCGACTACAGCCTGCAGCGCAGCgg                                                                                                                                  <  2:177926/126‑1 (MQ=255)
               tgctgcCGATTCCGCGTAACGTCGCTGACGATACCGTCAAGGGTGAGTTTGAAGTGCCTCGCCCGCAACCGTTGACGGGCGGTGCCGCCCAGGTCACCGACTACAGCCTGCAGCGCAGCGGCAGCAGCCGTTggg                                                                                                                    <  6:237748/135‑1 (MQ=255)
                gctgcCGATTCCGCGTAACGTCGCTGACGATACCGTCAAGGGTGAGTTTGAAGTGCCTCGCCCGCAACCGTTGACGGGCGGTGCCGCCCAGGTCACCGACTACAGCCTGCAGCGCAGCGGCAGCAGCCGATGGGt                                                                                                                   >  1:248108/1‑135 (MQ=255)
                 ctgcCGATTCCGCGTAACGTCGCTGACGATACCGTCAAGGGTGAGTTTGAAGTGCCTCGCCCGCAACCGTTGACGGGCGGTGCCGCCCAGGTCACCGACTTCAGCCTGCAGCGCAGCGGCAGCCGCCggtgggtg                                                                                                                  >  6:270006/1‑135 (MQ=255)
                 ctgcCGATTCCGCGTAACGTCGCTGACGATACCGTCAAGGGTGAGTTTGAAGTGCCTCGCCCGCAACCGTTGACGGGCGGTGCCGCCCAGGTCACCGACTACAGCCTGCAGCGCAGCGGCAGCAGCCGTTGGGTg                                                                                                                  >  1:132820/1‑135 (MQ=255)
                 ctgcCGATTCCGCGTAACGTCGCTGACGATACCGTCAAGGGTGAGTTTGAAGTGCCTCGCCCGCAACCGTTGACGGGCGGTGCCGCCCAGGTCACCGACTACAGCCTGCAGCGCAGCGGCAGCAGCCGTTGGGTGc                                                                                                                 >  7:233658/1‑136 (MQ=255)
                  tgcCGATTCCGCGTAACGTCGCTGACGATACCGTCAAGGGTGAGTTTGAAGTGCCTCGCCCGCAACCGTTGACGGGCGGTGCCGCCCAGGTCACCGACTACAGCCTGCAGCGCAGCGGCAGCAGCCGTTGGGTGc                                                                                                                 >  4:213538/1‑135 (MQ=255)
                  tgcCGATTCCGCGTAACGTCGCTGACGATACCGTCAAGGGTGAGTTTGAAGTGCCTCGCCCGCAACCGTTGACGGGCGGTGCCGCCCAGGTCACCGACTACAGCCTGCAGCGCAGCGGCAGCAGCCGTTGGGTGCt                                                                                                                >  8:10915/1‑136 (MQ=255)
                   gcCGATTCCGCGTAACGTCGCTGACGATACCGTCAAGGGTGAGTTTGAAGTGCCTCGCCCGCAACCGTTGACGGGCGGTGCCGCCCAGGTCACCGACTACAGCCTGCAGCGCAGCGGCAGCAGCCGTTGGGTGCt                                                                                                                <  5:226416/135‑1 (MQ=255)
                    ccGATTCCGCGTAACGTCGCTGACGATACCGTCAAGGGTGAGTTTGAAGTGCCTCGCCCGCAACCGTTGACGGGCGGTGCCGCCCAGGTCACCGACTACAGCCTGCAGCGCAGCGGCAGCAGCCGTTGGGTGCTg                                                                                                               <  4:41916/135‑1 (MQ=255)
                    ccGATTCCGCGTAACGTCGCTGACGATACCGTCAAGGGTGAGTTTGAAGTGCCTCGCCCGCAACCGTTGACGGGCGGTGCCGCCCAGGTCACCGACTACAGCCTGCAGCGCAGCGGCAGCAGCCGTTGGGTGCTg                                                                                                               <  7:132150/135‑1 (MQ=255)
                       aTTCCGCGTAACGTCGCTGACGATACCGTCAAGGGTGAGTTTGAAGTGCCTCGCCCGCAACCGTTGACGGGCGGTGCCGCCCAGGTCACCGACTACAGCCTGCAGCGCAGCGGCAGCAGCCGTTGGGTGCTGGcc                                                                                                            >  1:226828/1‑135 (MQ=255)
                            gcgTAACGTCGCTGACGATACCGTCAAGGGTGAGTTTGAAGTGCCTCGCCCGCAACCTTTGACGGGCGGTGCCGCCAAGGTCACCGACTACAGCCTGCAGCGCAGCGGCAGCAGCCGTTGGGTGCTGGCCCAGCa                                                                                                       >  1:387242/1‑135 (MQ=255)
                              gTAACGTCGCTGACGATACCGTCAAGGGTGAGTTTGAAGTGCCTCGCCCGCAACCGTTGACGGGCGGTGCCGCCCAGGTCACCGACTACAGCCTGCAGCGCAGCGGCAGCAGCCGTTGGGTGCTGGCCCAGCAtt                                                                                                     >  4:303716/1‑135 (MQ=255)
                              gTAACGTCGCTGACGATACCGTCAAGGGTGAGTTTGAAGTGCCTCGCCCGCAACCGTTGACGGGCGGTGCCGCCCAGGTCACCGACTACAGCCTGCAGCGCAGCGGCAGCAGCCGTTGGGTGCTGGCCCAGCATTc                                                                                                    >  7:270351/1‑136 (MQ=255)
                                aaCGTCGCTGACGATACCGTCAAGGGTGAGTTTGAAGTGCCTCGCCCGCAACCGTTGACGGGCGGTGCCGCCCAGGTCACCGACTACAGCCTGCAGCGCAGCGGCAGCAGCCGTTGGGTGCTGGCCCAGCATTCg                                                                                                   >  6:119757/1‑135 (MQ=255)
                                  cGTCGCTGACGATACCGTCAAGGGTGAGTTTGAAGTGCCTCGCCCGCAAGCGTTGACGGGCGGTGCCGCCCAGGTCACCGACTACAGCCTGCAGCTCAGCGGCAGCAGCCGTTGGGTGCTGGCCCAGCATTCGcc                                                                                                 <  1:323144/135‑1 (MQ=255)
                                    tCGCTGACTATACCGTCTAGGGTGAGTTTGAAGTGCCTCGCCCGCATCCGTTGACGGTCGGTGCCGCCCAGGTCACCGACTACAGCCTGCAGCGCAGCGGCAGCAGCCGTTGGTTGCTGGCCCAGCATTCGCCAg                                                                                               <  2:26670/135‑1 (MQ=255)
                                    tCGCTGACGATACCGTCAAGGGTGAGTTTGAAGTGCCTCGCCCGCAACCGTTGACGGGCGGTGCCGCCCAGGTCACCGACTACAGCCTGCAGCGCAGCGGCAGCAGCCGTTGGGTGCTGGCCCAGCATTCGCCAg                                                                                               <  2:7645/135‑1 (MQ=255)
                                                      aGGGTGAGTTTGAAGTGCCTCGCCCGCAACCGTTGACGGGCGGTGCCGCCCAGGTCACCGACTACAGCCTGCAGCGCAGCGGCAGCAGCCGTTGGGTGCTGGCCCAGCATTCGCCAGCCGAAGTGTGGCCGGTTg                                                                             >  5:128019/1‑135 (MQ=255)
                                                       gggTGAGTTTGAAGTGCCTCGCCCGCAACCGTTCACGGGCGGTGCCGCCCAGGTCACCGACTACAGCCTGCAGCGCAGCGGCAGCAGCCGTTGGGTGCTGGCCCAGCATTCGCCAGCCGAAGTGTGGCCGGTTGc                                                                            >  3:52535/1‑135 (MQ=255)
                                                         gTGAGTTTGAAGTGCCTCGCCCGCAACCGTTGACGGGCGGTGCCGCCCAGGTCACCGACTACAGCCTGCAGCGCAGCGGCAGCAGCCGTTGGGTGCTGGCCCAGCATTCGCCAGCCGAAGTGTGGCCGGTTGCCCg                                                                         <  7:201354/136‑1 (MQ=255)
                                                             gTTTGAAGTGCCTCGCCCGCAACCGTTGACGGGCGGTGCCGCCCAGGTCACCGACTACAGCCTGCAGCGCAGCGGCAGCAGCCGTTGGGTGCTGGCCCAGCATTCGCCAGCCGAAGTGTGGCCGGTTGCCCGGc                                                                       <  4:323464/134‑1 (MQ=255)
                                                               ttGAAGTGCCTCGCCCGCAACCGTTGACGGGCGGTGCCGCCCAGGTCACCGACTACAGCCTGCAGCGCAGCGGCAGCAGCCGTTGGGTGCTGGCCCAGCATTCGCCAGCCGAAGTGTGGCCGGTTGCCCGGCAGtt                                                                   >  7:205031/1‑136 (MQ=255)
                                                               ttGAAGTGCCTCGCCCGCAACCGTTGACGGGCGGTGCCGCCCAGGTCACCGACTACAGCCTGCAGCGCAGCGGCAGCAGCCGTTGGGTGCTGGCCCAGCATTCGCCAGCCGAAGTGTGGCCGGTTGCCCGGCAGt                                                                    >  4:205991/1‑135 (MQ=255)
                                                               ttGAAGTGCCTCGCCCGCAACCGTTGACGGGCGGTGCCGCCCAGGTCACCGACTACAGCCTGCAGCGCAGCGGCAGCAGCCGTTGGGTGCTGGCCCAGCATCCGCCAGCCGAAGTGTGGCCGGTTGCCCGGCAGt                                                                    >  2:7743/1‑135 (MQ=255)
                                                                 gAAGTGCCTCGCCCGCAACCGTTGACGGGCGGTGCCGCCCAGGTCACCGACTACAGCCTGCAGCGCAGCGGCAGCAGCCGTTGGGTGCTGGCCCAGCATTCGCCAGCCGAAGTGTGGCCGGTTGCCCGGCAGttct                                                                 >  7:173289/1‑136 (MQ=255)
                                                                  aaGTGCCTCGCCCGCAACCGTTGACGGGCGGTGCCGCCCAGGTCACCGACTACAGCCTGCAGCGCAGCGGCAGCAGCCGTTGGGTGCTGGCCCAGCATTCGCCAGCCGAAGTGTGGCCGGTTGCCCGGCAGttc                                                                  <  5:372255/134‑1 (MQ=255)
                                                                          cGCCCGCAACCGTTGACGGGCGGTGTCGCCCAGGTCACCGACTACAGCCTGCAGCGCAGCGGCAGCAGCCGTTGGGTGCTGGCCCAGCATTCGGCAGCCGAAGTGTGGCCGGTTGCCCCGCAGTTCTTCGAGGAc                                                         >  2:242952/1‑135 (MQ=255)
                                                                          cGCCCGCAACCGTTGACGGGCGGTGCCGCCCAGGTCACCGACTACAGCCTGCAGCGCAGCGGCAGCAGCCGTTGGGTGCTGGCCCAGCATTCGCCAGCCGAAGTGTGGCCGGTTGCCCGGCAGTTCTTCGAGGa                                                          >  3:258471/1‑134 (MQ=255)
                                                                           gCCCGCAACCGTTGACGGGCGGTGCCGCCCAGGTCACCGACTACAGCCTGCAGCGCAGCGGCAGCAGCCGTTGGGTGCTGGCCCAGCATTCGCCAGCCGAAGTGTGGCCGGTTGCCCGGCAGTTCTTCGAGGACaa                                                       >  8:140334/1‑136 (MQ=255)
                                                                            cccGCAACCGTTGACGGGCGGTGCCGCCCAGGTCACCGACTACAGCCTGCAGCGCAGCGGCAGCAGCCGTTGGGTGCTGGCCCAGCATTCGCCAGCCGAAGTGTGGCCGGTTGCCCGGCAGTTCTTCGAGGAc                                                         <  4:127447/133‑1 (MQ=255)
                                                                             ccGCAACCGTTGACGGGCGGTGCCGCCCAGGTCACCGACTACAGCCTGCAGCGCAGCGGCAGCAGCCGTTGGGTGCTGGCCCAGCATTCGCCAGCCGAAGTGTGGCCGGTTGCCCGGCAGTTCTTCGAGGACAAc                                                      <  3:36709/135‑1 (MQ=255)
                                                                                  aCCGTTGACGGGCGGTGCCGCCCAGGTCACCGACTACAGCCTGCAGCGCAGCGGCAGCAGCCGTTGGGTGCTGGCCCAGCATTCGCCAGCCGAAGTGTGGCCGGTTGCCCGGCAGTTCTTCGAGGACAACGGCtt                                                 <  1:301843/135‑1 (MQ=255)
                                                                                  aCCGGTGACGGGCGGGGCCGCCCAGGTCACCGACTACAGCCTGCAGCGCAGCGGCAGCAGCCGTTGGGTGCTGGCCCAGCATTCGCCAGCCGAAGTGTGGCCGGGTGCCCGGCAGTTCTTCGAGGACAACGGCtt                                                 <  3:15044/135‑1 (MQ=255)
                                                                                    cGTTGACGGGCGGTGCCGCCCAGGTCACCGACTACAGCCTGCAGCGCAGCGGCAGCAGCCGTTGGGTGCTGGCCCAGCATTCGCCAGCCGAAGTGTGGCCGGTTGCCCGGCAGTTCTTCGAGGACAACGGCTTcc                                               <  6:128019/135‑1 (MQ=255)
                                                                                    cGTTGACGGGCGGTGCCGCCCAGGTCACCGACTACAGCCTGCAGCGCAGCGGCAGCAGCCGTTGGGTGCTGGCCCAGCATTCGCCAGCCGAAGTGTGGCCGGTTGCCCGGCAGTTCTTCGAGGACAACGGCTTcc                                               <  5:237411/135‑1 (MQ=255)
                                                                                    cGTTGACGGGCGGTGCCGCCCAGGTCACCGACTACAGCCTGCAGCGCAGCGGCAGCAGCCGTTGGGTGCTGGCCCAGCATTCGCCAGCCGAAGTGTGGCCGGTTGCCCGGCAGTTCTTCGAGGACAACGGCTTcc                                               <  2:5913/135‑1 (MQ=255)
                                                                                               ggTGCCGCCCAGGTCACCGACTACAGCCTGCAGCGCAGCGGCAGCAGCCGTTGGGTGCTGGCCCAGCATTCGCCAGCCGAAGTGTGGCCGGTTGCCCGGCAGTTCTTCGAGGACAACGGCTTCCGCATTGCTgag                                    <  2:226828/135‑1 (MQ=255)
                                                                                                gTGCCGCCCAGGTCACCGACTACAGCCTGCAGCGCAGCGGCAGCAGCCGTTGGGTGCTGGCCCAGCATTCGCCAGCCGAAGTGTGGCCGGTTGCCCGGCAGTTCTTCGAGGACAACGGCTTCCGCATTGCTga                                     >  1:288454/1‑133 (MQ=255)
                                                                                                gTGCCGCCCAGGTCACCGACTACAGCCTGCAGCGCAGCGGCAGCAGCCGTTCGGTGCTGGCCCAGCATTCGCCAGCCGAAGTGTGGCCGGTTGCCCCGCAGTTCTTCGAGGACACCGGCTTCCGCATTGCTga                                     >  4:90688/1‑133 (MQ=255)
                                                                                                 tGCCGCCCAGGTCACCGACTACAGCCTGCAGCGCAGCGGCAGCAGCCGTTGGGTGCTg                                                                                                               <  7:29296/58‑1 (MQ=255)
                                                                                                 tGCCGCCCAGGTCACCGACTACAGCCTGCAGCGCAGCGGCAGCAGCCGTTGGGTGCTg                                                                                                               >  8:29296/1‑58 (MQ=255)
                                                                                                 tGCCGCCCAGGTCACCGACTACAGCCTGCAGCGCAGCGGCAGCAGCCGTTGGGTGCTGGCCCAGCATTCGCCAGCCGAAGTGTGGCCGGTTGCCCGGCAGTTCTTCGAGGACTACGGCTTCCGCATTGCTgaggag                                 <  8:100891/136‑1 (MQ=255)
                                                                                                    cgccCAGGTCACCGACTACAGCCTGCAGCGCAGCGGCAGCAGCCGTTGGGTGCTGGCCCAGCATTCGCCAGCCGAAGTGTGGCCGGTTGCCCGGCAGTTCTTCGAGGACAACGGCTTCCGCATTGCTGAGGAgcg                               >  1:100245/1‑135 (MQ=255)
                                                                                                    cgccCAGGTCACCGACTACAGCCTGCAGCGCAGCGGCAGCAGCCGTTGGGTGCTGGCCCAGCATTCGCCAGCCGAAGTGTGGCCGGTTGCCCGGCAGTTCTTCGAGGACAACGGCTTCCGCATTGCTGAGGAgcg                               >  1:211141/1‑135 (MQ=255)
                                                                                                      cccAGGTCACCGACTACAGCCTGCAGCGCAGCGGCAGCAGCCGTTGGGTGCTGGCCCAGCATTCGCCAGCCGAAGTGTGGCCGGTTGCCCGGCAGTTCTTCGAGGACAACGGCTTCCGCATTGCTGAGGAGCGcc                             <  4:258471/135‑1 (MQ=255)
                                                                                                        cAGGTCACCGACTACAGCCTGCAGCGCAGCGGCAGCAGCCGTTGGGTGCTGGCCCAGCATTCGCCAGCCGAAGTGTGGCCGGTTGCCCGGCAGTTCTTCGAGGACAACGGCTTCCGCATTGCTGAGGAGCGCCCGc                          <  8:173289/136‑1 (MQ=255)
                                                                                                         aGGTCACCGACTACAGCCTGCAGCGCAGCGGCAGCAGCCGTTGGGTGCTGGCCCAGCATTCGCCAGCCGAAGTGTGGCCGGTTGCCCGGCAGTTCTTCGAGGACAATGGCTTCCGCATTGCTGAGGAGCGccc                            >  5:385600/1‑133 (MQ=255)
                                                                                                         aGGTCACCGACTACAGCCTGCAGCGCAGCGGCAGCAGCCGTTGGGTGCTGGCCCAGCATTCGCCAGCCGAAGTGTGGCCGGTTGCCCGGCAGTTCTTCGAGGACAACGGCTTCCGCATTGCTGAGGAGCGCCCGc                          <  6:395779/135‑1 (MQ=255)
                                                                                                              aCCGACTACAGCCTGCAGCGCAGCGGCAGCAGCCGTTGGGTGCTGGCCCAGCATTCGCCAGCCGAAGTGTGGCCGGTTGCCCGGCAGTTCTTCGAGGACAACGGCTTCCGCATTGCTGAGGAGCGCCCGCAGACg                     <  2:255272/135‑1 (MQ=255)
                                                                                                              aCCGACTACAGCCTGCAGCGCAGCGGCAGCAGCCGTTGGGTGCTGGCCCAGCATTCGCCAGCCGAAGTGTGGCCGGTTGCCCGGCAGTTCTTCGAGGACAACGGCTTCCGCATTGCTGAGGAGCGCCCGCAGACg                     <  1:195840/135‑1 (MQ=255)
                                                                                                              aCCGACTACAGCCTGCAGCGCAGCGGCAGCAGCCGTTGGGTGCTGGCCCAGCATTCGCCAGCCGAAGTGTGGCCGGTTGCCCGGCAGTTCTTCGAGGACAACGGCTTCCGCATTGCTGAGGAGCGCCCGCAGACg                     <  7:205729/135‑1 (MQ=255)
                                                                                                              aCCGACTACAGCCTGCAGCGCAGCGGCAGCAGCCGTTGGGTGCTGGCCCAGCATTCGCCAGCCGAAGTGTGGCCGGTTGCCCGGCAGTTCTTCGAGGACAACGGCTTCCGCATTGCTGAGGAGCGCCCGCAGACg                     >  1:255272/1‑135 (MQ=255)
                                                                                                              aCCGACTACAGCCTGCAGCGCAGCGGCAGCAGCCGTTGGGTGCTGGCCCAGCATTCGCCAGCCGAAGTGTGGCCGGTTGCCCGGCAGTTCTTCGAGGACAACGGCTTCCGCATTGCTGAGGAGCGCCCGCAGACg                     >  1:27166/1‑135 (MQ=255)
                                                                                                              aCCGACTACAGCCTGCAGCGCAGCGGCAGCAGCCGTTGGGTGCTGGCCCAGCATTCGCCAGCCGAAGTGTGGCCGGTTGCCCGGCAGTTCTTCGAGGACAACGGCTTCCGCATTGCTGAGGAGCGCCCGCAGACg                     <  4:393233/135‑1 (MQ=255)
                                                                                                              aCCGACTACAGCCTGCAGCGCAGCGGCAGCAGCCGTTGGGTGCTGGCCCAGCATTCGCCAGCCGAAGTGTGGCCGGTTGCCCGGCAGTTCTTCGAGGACAACGGCTTCCGCATTGCTGAGGAGCGCCCGCAGACg                     >  6:366047/1‑135 (MQ=255)
                                                                                                              aCCGACTACAGCCTGCAGCGCAGCGGCAGCAGCCGTTGGGTGCTGGCCCAGCATTCGCCAGCCGAAGTGTGGCCGGTTGCCCGGCAGTTCTTCGAGGACAACGGCTTCCACATTGCTGAGGAGCGCCCGCAGCCg                     >  2:140773/1‑135 (MQ=255)
                                                                                                              aCCGACTACAGCCTGCAGCGCAGCGGCAGCAGCCGGTGGGTGCTGGCCCAGCATTCGCCAGCCGAAGTGTGGCCGGTTGCCCGGCAGTTCTTCGAGGACAACGGCTTCCGCATTGCTGAGGAGCGCCCGCAGACg                     <  1:179339/135‑1 (MQ=255)
                                                                                                               ccGACTACAGCCTGCAGCGCAGCGGCAGCAGCCGTTGGGTGCTGGCCCAGCATTCGCCAGCCGAAGTGTGGCCGGTTGCCCGGCAGTTCTTCGAGGACAACGGCTTCCGCATTGCTGAGGAGCGCCCGCAGACg                     <  2:132820/134‑1 (MQ=255)
                                                                                                               ccGACTACAGCCTGCAGCGCAGCGGCAGCAGCCGTTGGGTGCTGGCCCAGCATTCGCCAGCCGAAGTGTGGCCGGTTGCCCGGCAGTTCTTCGAGGACAACGGCTTCCGCATTGCTGAGGAGCGCCCGCAGACg                     <  2:166640/134‑1 (MQ=255)
                                                                                                               ccGACTACAGCCTGCAGCGCAGCGGCAGCAGCCGTTGGGTGCTGGCCCAGCATTCGCCAGCCGAAGTGTGGCCGGTTGCCCGGCAGTTCTTCGAGGACAACGGCTTCCGCATTGCTGAGGAGCGCCCGCAGACg                     <  3:213538/134‑1 (MQ=255)
                                                                                                                 gACTACAGCCTGCAGCGCAGCGGCAGCAGCCGTTGGGTGCTGGCCCAGCATTCCCCAGCCGAAGTGTGGCCGGTTGCCCGGCAGTTCTTCGAGGACAACGGCTACCGCATTGCTGAGGCGCGCCCCCAGCCGGGt                  >  5:306618/1‑135 (MQ=255)
                                                                                                                        gCCTGCAGCGCAGCGGCAGCAGCCGTTGGGTGCTGGCCCAGCATTCGCCAGCCGAAGTGTGGCCGGTTGCCCGGCAGTTCTTCGAGGACAACGGCTTCCGCATTGCTGAGGAGCGCCCGCAGCCGGGTGAGTTCa           >  6:246642/1‑135 (MQ=255)
                                                                                                                        gCCTGCAGCGCAGCGGCAGCAGCCGTTGGGTGCTGGCCCAGCATTCGCCAGCCGAAGTGTGGCCGGTTGCCCGGCAGTTCTTCGAGGACAACGGCTTCCGCATTGCTGAGGAGCGCCCGCAGACGGGTGAGTTc            >  3:143883/1‑134 (MQ=255)
                                                                                                                        gCCTGCAGCGCAGCCGCAGCAGCCGTTGGGTGCTGGCCCAGCATTCGCCAGCCGAAGTGTGGCCGGTTGCCCGGCAGTTCTTCGAGGACAACGGCTTCCGCATTGCTGAGGAGCGCCCGCAGACGGGTGa                >  6:52900/1‑130 (MQ=255)
                                                                                                                        gCCTGCAGCGCAGCCGCAGCAGCCGTTGGGTGCTGGCCCAGCATTCGCCAGCCGAAGTGTGGCCGGTTGCCCGGCAGTTCTTCGAGGACAACGGCTTCCGCATTGCTGAGGAGCGCCCGCAGACGGGTGa                <  5:52900/130‑1 (MQ=255)
                                                                                                                        gCCTGCAGCGCAGCCGCAGCAGCCGTTGGGTGCTGGCCCAGCATTCGCCAGCCGAAGTGTGGCCGGTTGCCCGGCAGTTCTTCGAGGACAACGGCTTCCGCATTGCTGAGGAGCGCCCGCAGACGGGTGa                <  5:47523/130‑1 (MQ=255)
                                                                                                                        gCCTGCAGCGCAGCCGCAGCAGCCGTTGGGTGCTGGCCCAGCATTCGCCAGCCGAAGTGTGGCCGGTTGCCCGGCAGGTCTTCGAGGACAACGGCTTCCGCATTGCTGAGGAGCGCCCGCAGACGGGTGa                >  6:47523/1‑130 (MQ=255)
                                                                                                                            gcagcgcagcGGCAGCAGCCGTTGGGTGCTGGCCCAGCATTCGCCAGCCGAAGTGTGGCCGGTTGCCCGGCAGTTCTTCGAGGACAACGGCTTCCGCATTGCTGAGGAGCGCCCGCAGCCGGGTGAGTTCAATAc       >  2:378255/1‑135 (MQ=255)
                                                                                                                            gcagcgcagcGGCAGCAGCCGTTGGGTGCTGGCCCAGCATTCGACAGCCGAAGTGTGGCAGGTTGCCCGGCAGTTCTTCGAGGACAACTGCTTCCGCATTGCTGAGGAGCTCCCGCAGCCGGGTGAGTTCAATAc       >  6:164636/1‑135 (MQ=255)
                                                                                                                             cagcgcagcgGCAGCAGCCGTTGGGTGCTGGCCCAGCATTCGCCAGCCGAAGTGTGGCCGGTTGCCCGGCAGTTCTTCGAGGACAACGGCTTCCGCATTGCTGAGGAGCGCCCGCAGACGGGTGAGTTCAATAcc      <  2:273398/135‑1 (MQ=255)
                                                                                                                             cagcgcagcgGCAGCAGCCGTTGGGTGCTGGCCCAGCATTCGCCAGCCGAAGTGTGGCCGGTTGCCCGGCAGTTCTTCGAGGACAACGGCTTCCGCATTGCTGAGGAGCGCCCGCAGACGGGTGAGTTCAATAcc      >  3:42459/1‑135 (MQ=255)
                                                                                                                                cgcagcgGCAGCAGCCGTTGGGTGCTGGCCCAGCATTCGCCAg                                                                                               >  2:268570/1‑43 (MQ=255)
                                                                                                                                cgcagcgGCAGCAGCCGTTGGGTGCTGGCCCAGCATTCGCCAg                                                                                               <  1:268570/43‑1 (MQ=255)
                                                                                                                                cgcagcgGCAGCAGCCGTTGGGTGCTGGCCCAGCATTCGCCAGCCGAAGTGTGGCCGGTTGCCCGGCAGTTCTTCGAGGACATCGGCTTCCGCATTGCTGAGGAGCGCCCGCAGACGGGTGAGTTCAATACCGCCt  >  7:46176/1‑136 (MQ=255)
                                                                                                                                cgcagcgGCAGCAGCCGTTGGGTGCTGGCCCAGCATTCGCCAGCCGAAGTGTGGCCGGTTGCCCGGCAGTTCTTCGAGGACAACGGCTTCCGCATTGCTGAGGAGCGCCCGCAGACGGGTGAGTTCAATACCGcc   >  3:201943/1‑135 (MQ=255)
                                                                                                                                cgcagcgGCAGCAGCCGTTGGGTGCTGGCCCAGCATTCGCCAGCCGAAGTGTGGCCGGTTGCCCGGCAGTTCTTCGAGGACAACGGCTTCCGCATTGCTGAGGAGCGCCCGCAGACGGGTGAGTTCAATACCGcc   >  3:85069/1‑135 (MQ=255)
                                                                                                                                cgcagcgGCAGCAGCCGTTGGGTGCTGGCCCAGCATTCGCCAGCCGAAGTGTGGCCGGTTGCCCGGCAGTTCTTCGAGGACAACGGCTTCCGCATTGCTGAGGAGCGCCCGCAGACGGGTGAGTTCAATACCGcc   >  5:150378/1‑135 (MQ=255)
                                                                                                                                cgcagcgGCAGCAGCCGTTGGGTGCTGGCCCAGCATTCGCCAGCCGAAGTGTGGCCGGTTGCCCGGCAGTTCGTCGCGGACAACGGCTTCCGCATTGCTGCGGAGCGCCCGCAGCCGGGTGAGTTCAATACCGCCt  >  8:22062/1‑136 (MQ=255)
                                                                                                                                 gcagcgGCAGCAGCCGTTGGGTGCTGGCCCAGCATTCGCCAGCCGAAGTGTGGCCGGTTGCCCGGTAGTTCTTCGAGGACAACGGCTTCCGCATTGCTGAGGAGCGCCCGCAGACGGGTGAGTTCAATACCGCCt  >  3:229848/1‑135 (MQ=255)
                                                                                                                                      |                                                                                                                                 
AACGCCTTGACCCGCTGCTGCCGATTCCGCGTAACGTCGCTGACGATACCGTCAAGGGTGAGTTTGAAGTGCCTCGCCCGCAACCGTTGACGGGCGGTGCCGCCCAGGTCACCGACTACAGCCTGCAGCGCAGCGGCAGCAGCCGTTGGGTGCTGGCCCAGCATTCGCCAGCCGAAGTGTGGCCGGTTGCCCGGCAGTTCTTCGAGGACAACGGCTTCCGCATTGCTGAGGAGCGCCCGCAGACGGGTGAGTTCAATACCGCCT  >  ACB122/1408850‑1409113

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: ‑