Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 31,613 G→A 5.8% R244Q (CGG→CAG)  oprP → Porin P

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB12231,6130GA 5.8% 223.6 / 2.0 69R244Q (CGG→CAG) oprPPorin P
Reads supporting (aligned to +/- strand):  ref base G (37/28);  new base A (2/2);  total (39/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTCGGATGGCGACAGCACCAAGCGCTACAACGTACGTGGCGTGTTCGCACCGCTGCACAGTGATGGCAACGTCCTGCACGTCGGTGCTCAGTACGCTTACCGCGACCTCAAAAACAGTGCGGTCGATACCCGGATTCGTTCGCGTCTGGGCGTCCGCGGTGTAGATACCAATGGCGGCGGCGATGCCGGCACCAACGGCAACCGCATGCTGTTTGGCGGTGCTGCTGCACAAGATGGCCTGTGGAAGGACGACTCGGTC  >  ACB122/31482‑31740
                                                                                                                                   |                                                                                                                               
cTCGGATGGCGACAGCACCAAGCGCTACAACGTACGTGGCGTGTTCGCACCGCTGCACAGTGATGGCAACGTCCTGCACGTCGGTGCTCAGTACGCTTACCGCGACCTCAAAAACAGTGCGGTCGATACCCGGAtt                                                                                                                             >  8:378036/1‑136 (MQ=255)
    gATGGCGACAGCACCAAGCGCTACAACGTACGTGGCGTGTTCGCACCGCTGCACAGTGATGGCAACGTCCTGCACGTCGGTGCTCAGTACGCTTACCGCGACCTCAAAAACAGTGCGGTCGATACCCGGAttc                                                                                                                            <  5:138149/133‑1 (MQ=255)
               cACCAAGCGCTACAACGTACGTGGCGTGTTCGCACCGCTGCACAGTGATGGCAACGTCCTGCACGTCGGTGCTCAGTACGCTTACCGCGACCTCAAAAACAGTGCGGTCGATACCCGGATTCGTTCGCGTCTggg                                                                                                               <  2:10041/135‑1 (MQ=255)
                  cAAGCGCTACAACGTACGTGGCGTGTTCGCACCGCTGCACAGTGATGGCAACGTCCTGCACGTCGGTGCTCAGTACGCTTACCGCGACCTCAAAAACAGTGCGGTCGATACCCGGATTCGTTCGCGTCTGGGCGTc                                                                                                           <  4:365442/136‑1 (MQ=255)
                  cAAGCGCTACAACGTACGTGGCGTGTTCGCACCGCTGCACAGTGATGGCAACGTCCTGCACGTCGGTGCTCAGTACGCTTACCGCGACCTCAAAAACAGTGCGGTCGATACCCGGATTCGTTCGCGTCTGGGCGTc                                                                                                           <  7:194635/136‑1 (MQ=255)
                   aaGCGCTACAACGTACGTGGCGTGTTCGCACCGCTGCACAGTGATGGCAACGTCCTGCACGTCGGTGCTCAGTACGCTTACCGCGACCTCAAAAACAGTGCGGTCGATACCCGGATTCGTTCGCGTCTGGGCGTcc                                                                                                          <  4:115832/136‑1 (MQ=255)
                    aGCGCTACAACGTACGTGGCGTGTTCGCACCGCTGCACAGTGATGGCAACGTCCTGCACGTCGGTGCTCAGTACGCTTACCGCGACCTCAAAAACAGTGCGGTCGATACCCGGATTCGTTCGCGTCTGGGCGTcc                                                                                                          <  6:147278/135‑1 (MQ=255)
                      cgcTACAACGTACGTGGCGTGTTCGCACCGCTGCACAGTGATGGCAACGTCCTGCACGTCGGTGCTCAGTACGCTTACCGCGACCTCAAAAACAGTGCGGTCGATACCCGGATTCGTTCGCGTCTGGGCGTCcgcg                                                                                                       >  3:340023/1‑136 (MQ=255)
                       gcTACAACGTACGTGGCGTGTTCGCACCGCTGCACAGTGATGGCAACGTCCTGCACGTCGGTGCTCAGTACGCTTACCGCGACCTCAAAAACAGTGCGGTCGATACCCGGATTCGTTCGCGTGTGGGCGTCcgcg                                                                                                       >  5:273590/1‑135 (MQ=255)
                       gcTACAACGTACGTGGCGTGTTCGCACCGCTGCACAGTGATGGCAACGTCCTGCACGTCGGTGCTCAGTACGCTTACCGCGACCTCAAAAACAGTGCGGTCGATACCCGGATTCGTTCGCGTCTGGGCGTCcgcg                                                                                                       >  6:129215/1‑135 (MQ=255)
                         tACAACGTACGTGGCGTGTTCGCACCGCTGCACAGTGATGGCAACGTCCTGCACGTCGGTGCTCAGTACGCTTACCGCGACCTCAAAAACAGTGCGGTCGATACCCGGATTCGTTCGCGTCTGGGCGTCCGCGgt                                                                                                     <  5:258515/135‑1 (MQ=255)
                            aaCGTACGTGGCGTGTTCGCACCGCTGCACAGTGATGGCAACGTCCTGCACGTCGGTGCTCAGTACGCTTACCGCGACCTCAAAAACAGTGCGGTCGATACCCGGATTCGTTCGCGTCTGGGCGTCCGCGGTGTa                                                                                                  >  1:292777/1‑135 (MQ=255)
                                   gtGGCGTGTTCGCACCGCTGCACAGTGATGGCAACGTCCTGCACGTCGGTGCTCAGTACGCTTACCGCGACCTCAAAAACAGTGCGGTCGATACCCGGATTCGTTCGCGTCTGGGCGTCCGCGGTGTAGATACCaa                                                                                          >  4:181973/1‑136 (MQ=255)
                                      gCGTGTTCGCACCGCTGCACAGTGATGGCAACGTCCTGCACGTCGGTGCTCAGTACGCTTACCGCGACCTCAAAAACAGTGCGGTCGATACCCGGATTCGTTCGCGTCTGGGCGTCCGCGGTGTAGATACCAATgg                                                                                       >  7:91017/1‑136 (MQ=255)
                                           ttCGCACCGCTGCACAGTGATGGCAACGTCCTGCACGTCGGTGCTCAGTACGCTTACCGCGACCTCAAAAACAGTGCGGTCGATACCCGGATTCGTTCGCGGCTGGGCGTCCGCGGTGTAGATACCAATggcggc                                                                                   <  7:171543/135‑1 (MQ=255)
                                             cGCACCGCTGCACAGTGATGGCAACGTCCTGCACGTCGGTGCTCAGTACGCTTACCGCGACCTCAAAAACAGTGCGGTCGATACCCGGATTCGTTCGCGTCTGGGCGTCCGCGGTGTAGATACCAATggcggcggc                                                                                >  7:236471/1‑136 (MQ=255)
                                              gCACCGCTGCACAGTGATGGCAACGTCCTGCACGTCGGTGCTCAGTACGCTTACCGCGCCCTCAAAAACAGTGCGGTCGATACCCGGATTCGTTCGCGGCTGGGCGTCCGCGGTGGAGATACCAATggccgcggc                                                                                >  6:221420/1‑135 (MQ=255)
                                                   gCTGCACAGTGATGGCAACGTCCTGCACGTCGGTGCTCAGTACGCTTACCGCGACCTCAAAAACAGTGCGGTCGATACCCGGATTCGTTCGCGTCTGGGCGTc                                                                                                           >  4:36223/1‑103 (MQ=255)
                                                   gCTGCACAGTGATGGCAACGTCCTGCACGTCGGTGCTCAGTACGCTTACCGCGACCTCAAAAACAGTGCGGTCGATACCCGGATTCGTTCGCGTCTGGGCGTc                                                                                                           <  3:36223/103‑1 (MQ=255)
                                                               tGGCAACGTCCTGCACGTCGGTGCTCAGTACGCTTACCGCGACCTCAAAAACAGTGCGGTCGATACCCGGATTCGTTCGCGTCTGGGCGTcc                                                                                                          >  4:365337/1‑92 (MQ=255)
                                                               tGGCAACGTCCTGCACGTCGGTGCTCAGTACGCTTACCGCGACCTCAAAAACAGTGCGGTCGATACCCGGATTCGTTCGCGTCTGGGCGTcc                                                                                                          <  3:365337/92‑1 (MQ=255)
                                                               tGGCAACGTCCTGCACGTCGGTGCTCAGTACGCTTACCGCGACCTCAAAAACAGTGCGGTCGATACCCGGATTCGTTCGCGTCTGGGCGTCCGCGGTGTAGATACCAATGGCGGCGGCGATGCCGGCACCAACgg                                                               >  2:47795/1‑135 (MQ=255)
                                                               tGGCAACGTCCTGCACGTCGGTGCTCAGTACGCTTACCGCGACCTCAAAAACAGTGCGGTCGATACCCGGATTCGTTCGCGTCTGGGCGTCCGCGGTGTAGATACCAATGGCGGCGGCGATGCCGGCACCAACgg                                                               >  7:293612/1‑135 (MQ=255)
                                                               tGGCAACGTCCTGCACGTCGGTGCTCAGTACGCTTACCGCGACCTCAAAAACAGTGCGGTCGATACCCGGATTCGTTCGCGTCTGGGCGTCCGCGGTGTAGATACCAATGGCGGCGGCGATGCCGGCACCAACg                                                                >  2:224322/1‑134 (MQ=255)
                                                                    aCGTCCTGCACGTCGGTGCTCAGTACGCTTACCGCGACCTCAAAAACAGTGCGGTCGATACCCGGATTCGTTCGCGTCTGGGCGTCCGCGGTGTAGATACCAATTGCGGCGGCGATGCCGGCACCAACGGCAACCg                                                         >  8:197257/1‑136 (MQ=255)
                                                                    aCGTCCTGCACGTCGGTGCTCAGTACGCTTACCGCGACCTCAAAAACAGTGCGGTCGATACCCGGATTCGTTCGCGTCTGGGCGTCCGCGGTGTAGATACCAATGGCGGCGGCGATGCCGGCACCAACGGCAAcc                                                          >  2:127053/1‑135 (MQ=255)
                                                                     cGTCCTGCACGTCGGTGCTCAGTACGCTTACCGCGACCTCAAAAACAGTGCGGTCGATACCCGGATTCGTTCGCGTCTGGGCGTCCGCGGTGTAGATACCAATGGCGGCGGCGATGCCGGCACCAACGGCAACCGc                                                        <  3:117093/136‑1 (MQ=255)
                                                                        ccTGCACGTCGGTGCTCAGTACGCTTACCGCGACCTCAAAAACAGTGCGGTCGATACCCGGATTCGTTCGCGTCTGGGCGTCCGCGGTGTAGATACCAATGGCGGCGGCGATGCCGGCACCAACGGCAACCGCATg                                                     >  3:176041/1‑136 (MQ=255)
                                                                         cTGCACGTCGGTGCTCAGTACGCTTACCGCGACCTCAAAAACAGTGCGGTCGATACCCGGATTCGTTCGCGTCTGGGCGTcc                                                                                                          >  8:208534/1‑82 (MQ=255)
                                                                         cTGCACGTCGGTGCTCAGTACGCTTACCGCGACCTCAAAAACAGTGCGGTCGATACCCGGATTCGTTCGCGTCTGGGCGTcc                                                                                                          <  7:208534/82‑1 (MQ=255)
                                                                         cTGCACGTCGGTGCTCAGTACGCTTACCGCGACCTCAAAAACAGTGCGGTCGATACCCGGATTCGTTCGCGTCTGGGCGTCCGCGGTGTAGATACCAATGGCGGCGGCGATGCCGGCACCAACGGCAACCGCATg                                                     >  6:226324/1‑135 (MQ=255)
                                                                         cTGCACGTCGGTGCTCAGTACGCTTACCGCGACCTCAAAAACAGTGCGGTCGATACCCAGATTCGTTCGCGTCTGGGCGTCCGCGGTGTAGATACCAATGGCGGCGGCGATGCCGGCACCAACGGCAACCTCATGc                                                    >  4:106053/1‑136 (MQ=255)
                                                                               gTCGGTGCTCAGTACGCTTACCGCGACCTCAAAAACAGTGCGGTCGATACCCGGATTCGTTCGCGTCTGGGCGTCCGCGGTGTAGATACCAATGGCGGCGGCGATGCCGGCACCAACGGCAACCGCATGCTGTTTg                                              >  3:150910/1‑136 (MQ=255)
                                                                                    tGCTCAGTACGCTTACCGCGACCTCAAAAACAGTGCGGTCGATACCCGGATTCGTTCGCGTCTGGGCGTCCGCGGTGTAGATACCAATGGCGGCGGCGATGCCGGCACCAACGGCAACCGCATGCTGTTTGGCgg                                          >  5:155372/1‑135 (MQ=255)
                                                                                    tGCTCAGTACGCTTACCGCGACCTCAAAAACAGTGCGGTCGATACCCGGATTCGTTCGCGTCTGGGCGTCCGCGGTGTAGATACCAATGGCGGCGGCGATGCCGGCACCAACGGCAACCGCATGCTGTTTGGCGGt                                         >  7:115049/1‑136 (MQ=255)
                                                                                    tGCTCAGTACGCTTACCGCGACCTCAAAAACAGTGCGGTCGATACCCGGATTCGTTCGCGTCTGGGCGTCCGCGGTGTAGATACCAATGGCGGCGGCGATGCCGGCACCAACGGCAACCGCATGCTGTTTGGCGGt                                         >  3:27233/1‑136 (MQ=255)
                                                                                       tCAGTACGCTTACCGCGACCTCAAAAACAGTGCGGTCGATACCCGGATTCGTTCGCGTCTGTGCGTCCGCGGTGTAGATACCAATGGCGGCGGCGATGCCGGCACCAACGGCAACCGCATGCTGTTTGGCGGtgct                                      >  4:97196/1‑136 (MQ=255)
                                                                                       tCAGTACGCTTACCGCGACCTCAAAAACAGTGCGGTCGATACCCGGATTCGTTCGCGTCTGGGCGTCCGCGGTGTAGATACCAATGGCGGCGGCGATGCCGGCACCAACGGCAACCGCATGCTGTTTGGCGGtg                                        <  1:31760/134‑1 (MQ=255)
                                                                                        cAGTACGCTTACCGCGACCTCAAAAACTGTGCGGTCGATACCCGGATTCGGTCGCGGCTGGGCGTCCGCGGTGTAGATACCAATGGCGGCGGCGATGCCGGCACCAACGGCAACCGCATGCTGTTTGGCGGtgctg                                     <  3:5767/136‑1 (MQ=255)
                                                                                             cGCTTACCGCGACCTCAAAAACAGTGCGGTCGATACCCGGATTCGTTCGCGTCTGGGCGTCCGCGGTGTAGATACCAATGGCGGCGGCGATGCCGGCACCAACGGCAACCGCATGCTGTTTGGCGGTGCTGCTGca                                >  3:289037/1‑136 (MQ=255)
                                                                                              gCTTACCGCGACCTCAAAAACAGTGCGGTCGATAGCCGGATTCGTTCGCGTCTGGGCGTCCGCGGTGTAGATACCAATGGCGGCGGCGATGCCGGCACCAACGGCAACCGCATGCTGTTTGGCGGTGCTGCTGcac                               <  8:3207/136‑1 (MQ=255)
                                                                                              gCTTACCGCGACCTCAAAAACAGTGCGGTCGATACCCGGATTCGTTCGCGTCTGGGCGt                                                                                                            <  3:291200/59‑1 (MQ=255)
                                                                                              gCTTACCGCGACCTCAAAAACAGTGCGGTCGATACCCGGATTCGTTCGCGTCTGGGCGt                                                                                                            >  4:291200/1‑59 (MQ=255)
                                                                                               cTTACCGCGACCTCAAAAACAGTGCGGTCGATACCCGGATTCGTTCGCGTCTGGGCGTCCGCGGTGTAGATACCAATGGCGGCGGCGATGCCGGCACCAACGGCAACCGCATGCTGTTTGGCGGTGCTGCTGcaca                              <  4:27233/136‑1 (MQ=255)
                                                                                               cTTACCGCGACCTCAAAAACAGTGCGGTCGATACCCGGATTCGTTCGCGTCTGGGCGTCCGCGGTGTAGATACCAATGGCGGCGGCGATGCCGGCACCAACGGCAACCGCATGCTGTTTGGCGGTGCTGCTGcaca                              <  3:280956/136‑1 (MQ=255)
                                                                                                   ccGCGACCTCAAAAACAGTGCGGTCGATACCCGGATTCGTTCGCGTCTGGGCGTCCGCGGTGTAGATACCAATGGCGGCGGCGATGCCGGCACCAACGGCAACCGCATGCTGTTTGGCGGTGCTGCTGCACACGAt                          <  4:340023/136‑1 (MQ=255)
                                                                                                   ccGCGACCTCAAAAACAGTGCGGTCGATACCCGGATTCGTTCGCGTCTGGGCGTCCGCGGTGTAGATACCAATGGCGGCGGCGATGCCGGCACCAACGGCAACCGCATGCTGTTTGGCGGTGCTGCTGCACAAGAt                          >  7:228799/1‑136 (MQ=255)
                                                                                                    cgcgACCTCAAACCCAGTGCGGTCGATCCCCGGCTTCGGTCGCGGCTGGGCGGCCGCGGTGTAGATACCAATGGCGGCGGCGATGCCGGCACCAACGGCAACCGCATGCTGTTTGGCGGTGCTGCTGCACAAGAt                          <  5:221420/135‑1 (MQ=255)
                                                                                                    cgcgACCTCAAAAACAGTGCGGTCGATACCCGGATTCGTTCGCGTCTGGGCGTCCGCGGTGTAGATACCAATGGCGGCGGCGATGCCGGCACCAACGGCAACCGCATGCTGTTTGGCGGTGCTGCTGCACAAGAt                          <  5:129215/135‑1 (MQ=255)
                                                                                                     gcgACCTCAAAAACAGTGCGGTCGATACCCGGATTCGTTCGCGTCTGGGCGTCCGCGGTGTAGATACCAATGGCGGCGGCGATGCCGGCACCaa                                                                  >  5:64754/1‑94 (MQ=255)
                                                                                                     gcgACCTCAAAAACAGTGCGGTCGATACCCGGATTCGTTCGCGTCTGGGCGTCCGCGGGGTAGATACCAATGGCGGCGGCGATGCCGGCACCaa                                                                  <  6:64754/94‑1 (MQ=255)
                                                                                                      cgACCTCAAAAACAGTGCGGTCGATACCCGGATTCGTTCGCGTCTGGGCGTCCGCGGTGTAGATACCAATGGCGGCGGCGATGCCGGCACCAACGGCAACCGCATGCTGTTTGGCGGTGCTGCTGCACAAGATgg                        >  5:114739/1‑135 (MQ=255)
                                                                                                      cgACCTCAAAAACAGTGCGGTCGATACCCGGATTCGTTCGCGTCTGGGCGTCCGCGGTGTAGATACCAATGGCGGCGGCGATGCCGGCACCAACGGCAACCGCATGCTGTTTGGCGGTGCTGCTGCACAAGATGGc                       >  8:55836/1‑136 (MQ=255)
                                                                                                      cgACCTCAAAAACAGTGCGGTCGATACCCGGATTCGTTCGCGTCTGGGCGTCCGCGGTGTAGATACCAATGGCGGCGGCGATGCCGGCACCAACGGCAACCGCATGCTCTTTGGCGGTGCTGCTGAACAAGATgt                        >  5:245782/1‑134 (MQ=255)
                                                                                                        aCCTCAAAAACAGTGCGGTCGATACCCGGATTCGTTCGCGTCTGGGCGTCCGCGGTGTAGATACCAATGGCGGCGGCGATGCCGGCACCAACg                                                                <  1:103164/93‑1 (MQ=255)
                                                                                                        aCCTCAAAAACAGTGCGGTCGATACCCGGATTCGTTCGCGTCTGGGCGTCCGCGGTGTAGATACCAATGGCGGCGGCGATGCCGGCACCAACg                                                                >  2:103164/1‑93 (MQ=255)
                                                                                                        aCCTCAAAAACAGTGCGGTCGATACCCGGATTCGTTCGCGTCTGGGCGTCCGCGGTGTAGATACCAATGGCGGCGGCGATGCCGGCACCAACGGCAACCGCATGCTg                                                  >  3:151405/1‑107 (MQ=255)
                                                                                                        aCCTCAAAAACAGTGCGGTCGATACCCGGATTCGTTCGCGTCTGGGCGTCCGCGGTGTAGATACCAATGGCGGCGGCGATGCCGGCACCAACGGCAACCGCATGCTg                                                  <  4:151405/107‑1 (MQ=255)
                                                                                                        aCCTCAAAAACAGTGCGGTCGATACCCGGATTCGTTCGCGTCTGGGCGTCCGCGGTGTAGATACCAATGGCGGCGGCGATGCCGGCACCAACGGCAACCGCATGCTGTTTGGCGGTGCTGCTGCACAAGATGGCCt                     >  7:116683/1‑136 (MQ=255)
                                                                                                        aCCTCAAAAACAGTGCGGTCGATACCCAGATTCGTTCGCGTCTGGGCGTCCGCGGTGTAGATACCAATGGCGGCGGCGATGCCGGCACCAACGGCAACCGCATGCTGTTTGGCGGTGCTGCTGCACAAGATGGCCt                     >  3:326463/1‑136 (MQ=255)
                                                                                                              aaaaCAGTGCGGTCGATACCCAGATTCGTTCGCGTCTGGGCGTCCGCGGTGTAGATACCAATGGGGGCGGCGATGCCGGCACCAACGGCAACCGCTTGCTGTTTGGCGGTGCTGCTGCACAAGATGGCCTGTGGaa               <  3:106053/136‑1 (MQ=255)
                                                                                                                  cAGTGGGGTCGATCCCCGGATTCGTTCGCGGCTGGGCGTCCGCGGTGTAGATACCATTGGGGGGGGCGATGCGGGAACCAACGGCAACCGCTGGCTTTTTGGCGGTGCTGCTGCACAAGATGGGCTGTTTAAGgac           <  4:176041/136‑1 (MQ=255)
                                                                                                                   aGTGCGGTCGATACCCGGATTCGTGCGCGTCTGGGCGTCCGCGGTGTAGATACCAATGGCGGCGGCGATGCCGGCACCAACGGCAACCGCATGCTGTTTGGCGGTGCTGCTGCACAAGATGGCCTGTGGAAGgacg          <  8:321448/136‑1 (MQ=255)
                                                                                                                        ggTCGATACCCGGATTCGTTCGCGTCTGGGCGTCCGCGGTGTAGATACCAATGGCGGCGGCGATGCCGGCACCAACGGCAACCGCATGCTGTTTGGCGGGGCTGCTGCACAAGATGGCCAGTGGAAGGACGCCTCg     >  3:269762/1‑136 (MQ=255)
                                                                                                                        ggTCGATACCCGGATTCGTTCGCGTCTGGGCGTCCGCGGTGTAGATACCAATGGCGGCGGCGATGCCGGCACCAACGGCAACCGCAGGCTGGTTGGCGGTGCTGCTGCACAAGATGGCCTGTGGAAGGACGACTCg     >  8:76644/1‑136 (MQ=255)
                                                                                                                          tCGATACCCGGATTCGTTCGCGTCTGGGCGTCCGCGGTGTAGATACCAATGGCGGCGGCGATGCCGGCACCAACGGCAACCGCATGCTGTTTGGCGGTGCTGCTGCACAAGATGGCCTGTGGAAGGACGACTCGGt   <  8:116683/136‑1 (MQ=255)
                                                                                                                          tCGATACCCGGATTCGTTCGCGTCTGGGCGTCCGCGGTGTAGATACCAATGGCGGCGGCGATGCCGGCACCAACGGCAACCGCATGCTGTTTGGCGGTGCTGCTGCACAAGATGGCCTGTGGAAGGACGACTCGGt   >  3:311004/1‑136 (MQ=255)
                                                                                                                           cGATACCCGGATTCGTTCGCGTCTGGGCGTCCGCGGTGTAGATACCAATGGCGGCGGCGATGCCGGCACCAACGGCAACCGCATGCTGTTTGGCGGTGCTGCTGCACAAGATGGCCTGTGGAAGGACGACTCGGTc  <  8:293612/136‑1 (MQ=255)
                                                                                                                           cGATACCCAGATTCGTTCGCGTCTGGGCGTCCGCGGTGTAGATACCAATGGCGGCGGCGATGCCGGCACCAACGGCAACCGCATGCTGTTTGGCGGTGCTGCTGCACAAGATGGCCTGTGGAAGGACGACTCGGTc  <  4:326463/136‑1 (MQ=255)
                                                                                                                                   |                                                                                                                               
CTCGGATGGCGACAGCACCAAGCGCTACAACGTACGTGGCGTGTTCGCACCGCTGCACAGTGATGGCAACGTCCTGCACGTCGGTGCTCAGTACGCTTACCGCGACCTCAAAAACAGTGCGGTCGATACCCGGATTCGTTCGCGTCTGGGCGTCCGCGGTGTAGATACCAATGGCGGCGGCGATGCCGGCACCAACGGCAACCGCATGCTGTTTGGCGGTGCTGCTGCACAAGATGGCCTGTGGAAGGACGACTCGGTC  >  ACB122/31482‑31740

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: