Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 4,609,122 G→A 5.2% D1293N (GAC→AAC)  CJ019_04355 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1224,609,1220GA 5.2% 252.1 / 2.8 77D1293N (GAC→AAC) CJ019_04355hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (54/19);  new base A (2/2);  total (56/21)
Fisher's exact test for biased strand distribution p-value = 2.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGACTATCCAGAGCCGCGTCGGCGGGTGCTGTTTCACTGGCAAGGGTTTCGCTTGCTTCAAGAAGTGCAGAGTGGTTTGGCTAGCTTGTATGTGTATGCCACTGTGGAGAGCTATGACCCGTTGGCGCGGGTTGACGGTAAACCTGGCAGCGAAGCAATTCACTATTTTCATACCACCTTGGTTGGGCAACCGGAGCAGCTGACTGATTCTGATGGCAACACGATCTGGCGTAGTGACCATCACGGTTGGGGAAAGATAATTAA  >  ACB122/4608989‑4609252
                                                                                                                                     |                                                                                                                                  
ggACTATCCAGAGCCGCGGCGGCGGGTGCTGTTTCACTGGCAAGGGTTTCGCTTGCTTCAAGAAGTGCAGAGTGGTTTGGCTAGCTTGTATGTGTATGCCACTGTGGAGAGCTATGACCCGTTGGCGCGGGTTGAc                                                                                                                                  <  3:10488/136‑1 (MQ=255)
ggACTATCCAGAGCCGCGGCGGCGGGTGCTGTTTCACTGGCAAGGGTATCGCTTGCTTCAAGAAGTGCAGAGTGGTTTGGCTAGCTTGTATGTGTATGCCACTGTGGAGAGCTATGACCCGTTGGCGCGGGTTGAc                                                                                                                                  <  8:373931/136‑1 (MQ=255)
gcaCTATCAAGAGCCGCGTCGGCGGGAGCTGTTTCACTGGCAAGCGATTCGCATCCATCAAGAAGTGCAGACAGGTTTGGCTAGCATGTATGTGTATGCCACTGTGGAGAGCTATGACCCCTATGCCCGGGTTgcc                                                                                                                                  >  8:8487/3‑134 (MQ=17)
  aCTATCCAGAGCCGCGTCGGCGGGTGCTGTTTCACTGGCAAGGGTTTCGCTTGCTTCAAGAAGTGCAGAGTGGTTTGGCTAGCGTGTATGTGTATGCCACGGTGGAGAGCTATGACCCGTTGGCGCGGGTTGACgg                                                                                                                                >  3:60161/1‑136 (MQ=255)
        cAGAGCCGCGTCGGCGGGTGCTGTTTCACTGGCAAGGGTTTAGCTTGCTGCAAGAATTGCAGAGTGGTTTGGGTACCTTGTTTGTGTATTCCACGGTTGAGAGCGATGACACGTTTGGGCGGGTTGCAGGTaaa                                                                                                                            >  6:28997/1‑134 (MQ=17)
              cgcgTCGGCGGGTGCTGTTTCACTGGCAAGGGTTTCGCTTGCTTCAAGAAGTGCAGAGTGGTTTGGCTAGCTTGTATGTGTATGCCACTGTGGAGCGCTATGACCCGTTGGCGCGGGTTGACGGTAAACCTGGCa                                                                                                                     >  5:189582/1‑135 (MQ=255)
              cgcgTCGGCGGGTGCTGTTTCACTGGCAAGGGTTTCGCTTGCTTCAAGAAGTGCAGAGTGGTTTGGCTAGCTTGTATGTGTATGCCACTGTGGAGAGCTATGACCCGTTGGGGCGGGTTGACGGTAAACCTGGCa                                                                                                                     >  6:161225/1‑135 (MQ=255)
              cgcgTCGGCGGGTGCTGTTTCACTGGCAAGGGTTTCGCTTGCTTCAAGAAGTGCAGAGTGGTTTGGCTAGCTTGTATGTGTATGCCACTGTGGAGAGCTATGACCCGTTGGCGCGGGTTGACGGTAAACCTGGCa                                                                                                                     >  2:36267/1‑135 (MQ=255)
              cgcgTCGGCGGGTGCTGTTTCACTGGCAAGGGTTTCGCTTGCTTCAAGAAGTGCAGAGTGGTTTGGCTAGCTTGTATGTGTATGCCACTGTGGAGAGCTATGACCCGTTGGCGCGGGTTGACGGTAAACCTGGCAg                                                                                                                    >  4:160474/1‑136 (MQ=255)
              cgcgTCGGCGGGTGCTGTTTAACTGGCAAGGGTTTCGCTTGCATCAAGAAGTGCAGAGTGGTTTGGCTAGCGTGTATGTGTATGCCACTGTGGCTAGCTATGACCCGTTGGCGCGGGGTGACGGTAAACCTGGCa                                                                                                                     >  1:284452/1‑135 (MQ=21)
                 gTCGGCGGGTGCTGTTTCACTGGCAAGGGTTTCGCTTGCTTCAAGAAGTGCAGAGTGGTTTGGCTAGCTTGTATGTGTATGCCACTGTGGAGAGCTATGACCCGTTGGCGCGGGTTGACGGTAAACCTGGCAGCg                                                                                                                  <  6:181814/135‑1 (MQ=255)
                   cggcggGTGCTGTTTCCCTGGCAAGGGTTTCGCTTGCTTCAAGAAGTGCAGAGTGGTTTGGCTAGCTTGTATGTGTATGCCACTTTGGAGAGCTATGACCCGTTGGCGCGGGTTGACGGTAAACCTGGCAGCGAAg                                                                                                               >  3:331343/1‑136 (MQ=255)
                   cggcggGTGCTGTTTCACTGGCAAGGGTTTCGCTTGCTTCAAGAAGTGCAGAGTGGTTTGGGTAGCTTGTATGTGTATTCCAATGTGGAGAGCTATGACCCGTTGGAGCGGGTTGCAGGTAAAACTGGCAGCgcat                                                                                                               >  4:215785/1‑133 (MQ=18)
                   cggcggGTGCTGTTTCACTGGCAAGGGTTTCGCTTGCTTCAAGAAGTGCAGAGTGGTTTGGGGAGCTTGTATGTGTATTGCACGGTGGCGAGCGATTCCCCGGTGGCGCGGGTTGACGGTTTACCAGGCCGCCaa                                                                                                                >  2:375792/1‑135 (MQ=17)
                   cggcggGTGCTGTTTCACTGGCAAGGGTTTCGCTTGCTTCAAGAAGTGCAGAGTGGTTTGGCTAGCTTGTATGTGTATGCCACTGTGGAGAGCTATGACCCGTTGGCGCGGGTTGACGGTAAACCTGGCAGCGAAg                                                                                                               >  8:92815/1‑136 (MQ=255)
                   cggcggGTGCTGTTTCACTGGCAAGGGTTTCGCTTGCTTCAAGAAGTGCAGAGTGGTTTGGCTAGCTTGTATGTGTATGCCACTGTGGAGAGCTATGACCCGTTGGCGCGGGTTGACGGTAAACCTGGCAGCGAAg                                                                                                               >  8:27842/1‑136 (MQ=255)
                   cggcggGTGCTGTTTCACTGGCAAGGGTTTCGCTTGCTTCAAGAAGTGCAGAGTGGTTTGGCTAGCTTGTATGTGTATGCCACTGTGGAGAGCTATGACCCGTTGGCGCGGGTTGACGGTAAACCTGGCAGCGAAg                                                                                                               >  8:313541/1‑136 (MQ=255)
                                   aCTGGCAAGGGTTTCGCTTGCTTCAAGAAGTGCAGAGTGGTTTGGCTAGCTTGTATGTGTATGCCACTGTGGAGAGCTATGACCCGTTGGCGCGGGTTGACGGTAAACCTGGCAGCGAAGCAATTCACTAtttt                                                                                                 <  1:36267/134‑1 (MQ=255)
                                    cTGGCAAGGGTTTCGCTTGCTTCAAGAAGTGCAGAGTGGTTTGGCTAGCTTGTATGTGTATGCCACTGTGGAGAGCTATGACCCGTTGGCGCGGGTTGACGGTAAACCTGGCAGCGAAGCAATTCACTATTTTCa                                                                                               <  2:242083/135‑1 (MQ=255)
                                          aGGGTTTCGCTTGCTTCAAGAAGTGCAGAGTGGTTTGGCTAGCTTGTATGTGTATGCCACTGTGGAGAGCTATGACCCGTTGGCGCGGGTTGACGGTAAACCTGGCAGCGAAGCAATTCACTATTTTCATaccac                                                                                         >  2:54404/1‑135 (MQ=255)
                                                         tCAAGAAGTGCAGAGTGGTTTGGCTAGCTTGTATGTGTATGCCACTGTGGAGAGCTATGACCCGTTGGCGCGGGTTGACGGTAAACCTGGCAGCGAAGCaa                                                                                                            <  5:388768/101‑1 (MQ=255)
                                                         tCAAGAAGTGCAGAGTGGTTTGGCTAGCTTGTATGTGTATGCCACTGTGGAGAGCTATGACCCGTTGGCGCGGGTTGACGGTAAACCTGGCAGCGAAGCaa                                                                                                            >  6:388768/1‑101 (MQ=255)
                                                         tCAAGAAGTGCAGAGTGGTTTGGCTAGCTTGTATGTGTATGCCACTGTGGAGAGCTATGACCCGTTGGCGCGGGTTGACGGTAAACCTGGCAGCGAAGCAATTCACTATTTTCATACCACCTTGGTTGGGCAAcc                                                                          >  1:323395/1‑135 (MQ=255)
                                                                    agagTGGTTTGGCTAGCTTGTATGTGTATGCCACTGTGGAGAGCTATGACCCGTTGGCGCGGGTTGACGGTAAACCTGGCAGCGCAGCAATTCACTATTTTCATACCACCTTGGTTGGGCAACCGGAGCAGctgcc                                                              >  3:277688/1‑134 (MQ=255)
                                                                    agagTGGTTTGGCTAGCTTGTATGTGTATGCCACTGTGGAGAGCTATGACCCGTTGGCGCGGGTTGACGGTAAACCTGGCAGCGAAGCAATTCACTATTTTCATACCACCTTGGTTGGGCAACCGGAGCAGctgac                                                              >  4:328612/1‑136 (MQ=255)
                                                                    agagTGGTTTGGCTAGCTTGTATGTGTATGCCACTGTGGAGAGCTATGACCCGTTGGCGCGGGTTGACGGTAAACCTGGCAGCGAAGCAATTCACTATTTTCATACCACCTTGGTTGGGCAACCGGAGCAGctga                                                               >  1:76311/1‑135 (MQ=255)
                                                                    agagTGGTTTGGCTAGCTTGTATGTGTATGCCACTGTGGAGAGCTATGACCCGTTGGCGCGGGTTAACGGTAAACCTGGCAGCGAAGCAATTCACTATTTTCATACCACCTTGGTTGGGCAACCGGAGCAGctgac                                                              >  8:24356/1‑136 (MQ=255)
                                                                    agagTGGTTTGGCTAGCTTGTATGTGTATGCCACTGTGGAGAGCTATGACCCGTTGGCGCGGGTTAACGGTAAACCTGGCAGCGAAGCAATTCACTATTTTCATACCACCTTGGTTGGGCAACCGGAGCAGctgac                                                              >  8:18314/1‑136 (MQ=255)
                                                                              ggCTAGCTTGTATGTGTATGCCACTGTGGAGCGCTATGACCCGTTGGCGCGGGTTGACGGTAAACCTGGCAGCGAAGCCATTCACTATTTTCATACCACCTTGGTTGGGCAACCGGAGCAGCTGACTGATTCTGa                                                     >  6:326380/1‑135 (MQ=255)
                                                                                 tAGCTTGTATGTGTATGCCACTGTGGAGAGCTATGACCCGTTGGCGCGGGTTGACGGTAAACCTGGCAGCGAAGCAATTCACTATTTTCATACCACCTTGGTTGGGCAACCGGAGCAGCTGACTGATTCTGATgg                                                  <  5:28997/135‑1 (MQ=255)
                                                                                    cTTGTATGTGTATGCCACTGTGGAGAGCTATGACCCGTTGGCGCGGGTTGACGGTAAACCTGGCAGCGAAGCAATTCACTATTTTCATACCACCTTGGTTGGGCAACCGGAGCAGCTGACTGATTCTGATGGCa                                                <  6:245102/134‑1 (MQ=255)
                                                                                       gtatgtgtatgCCACTGTGGAGAGCTATGACCCGTTGGCGCGGGTTGACGGTAAACCTGGCAGCGAAGCAATTCACTATTTTCATACCACCTTGGTTGGGCAACCGGAGCAGCTGACTGATTCTGATGGCAacac                                            >  7:125013/1‑135 (MQ=255)
                                                                                            tgtatgCCACTGTGGAGAGCTATGACCCGTTGGCGCGCGTTGACGGTAAACCTGGCAGCGAAGCAATTCACTATTTTCATACCACCTTGGTTGGGCAACCGGAGCAGCTGACTGATTCTGATGGCAACACGATCt                                       >  2:327489/1‑135 (MQ=255)
                                                                                                tgCCACTGTGGAGAGCTATGACCCGTTGGCGCGGGTTGACGGTAAACCTGGCAGCGAAGCAATTCACTAtttt                                                                                                 <  1:256507/73‑1 (MQ=255)
                                                                                                tgCCACTGTGGAGAGCTATGACCCGTTGGCGCGGGTTGACGGTAAACCTGGCAGCGAAGCAATTCACTAtttt                                                                                                 >  2:256507/1‑73 (MQ=255)
                                                                                                tgCCACTGTGGAGAGCTATGACCCGTTGGCGCGGGTTAACGGTAAACCTGGCAGCGAAGCAATTCACTATTTTCATACCACCTTGGTTGGGCAACCGGAGCAGCTGACTGATTCTGATGGCAACACGATCTGGCGt                                  <  7:18314/136‑1 (MQ=255)
                                                                                                tgCCACTGTGGAGAGCTATGACCCGTTGGCGCGGGTTAACGGTAAACCTGGCAGCGAAGCAATTCACTATTTTCATACCACCTTGGTTGGGCAACCGGAGCAGCTGACTGATTCTGATGGCAACACGATCTGGCGt                                  <  7:24356/136‑1 (MQ=255)
                                                                                                 gCCACTGTGGAGAGCTATGACCCGTTGGCGCGGGTTGACGGTAGACCTGGCAGCGAAGCAATTCACTATTTTCATACCACCTTGGTTGGGCAACCGGAGCAGCTGACTGATTCTGATGGCAACACGATCTGGCg                                   >  3:197064/1‑134 (MQ=255)
                                                                                                 gCCACTGTGGAGAGCTATGACCCGTTGGCGCGGGTTGACGGTAAACCTGGCAGCGAAGCAATTCACTATTTTCATACCACCTTGGTTGGGCAACCGGAGCAGCTGACTGAt                                                          <  2:233652/111‑1 (MQ=255)
                                                                                                 gCCACTGTGGAGAGCTATGACCCGTTGGCGCGGGTTGACGGTAAACCTGGCAGCGAAGCAATTCACTATTTTCATACCACCTTGGTTGGGCAACCGGAGCAGCTGACTGAt                                                          >  1:233652/1‑111 (MQ=255)
                                                                                                 gCCACTGTGGAGAGCTATGACCCGTTGGCGCGGGTTGACGGTAAACCTGGCAGCGAAGCAATTCACTATTTTCATACCACCTTGGTTGGGCAACCGGAGCAGCTGACTGATTCTGATGGCa                                                >  7:309265/1‑121 (MQ=255)
                                                                                                 gCCACTGTGGAGAGCTATGACCCGTTGGCGCGGGTTGACGGTAAACCTGGCAGCGAAGCAATTCACTATTTTCATACCACCTTGGTTGGGCAACCGGAGCAGCTGACTGATTCTGATGGCa                                                <  8:309265/121‑1 (MQ=255)
                                                                                                 gCCACTGTGGAGAGCTATGACCCGTTGGCGCGGGTTGACGGTAAACCTGGCAGCGAAGCAATTCACTATTTTCATACCACCTTGGTTGGGCAACCGGAGCAGCTGACTGATTCTGATGGCAACACGATCTGGc                                    >  5:52340/1‑133 (MQ=255)
                                                                                                 gCCACTGTGGAGAGCTATGACCCGTTGGCGCGGGTTGACGGTAAACCTGGCAGCGAAGCAATTCACTATTTTCATACCACCTTGGTTGGGCAACCGGAGCAGCTGACTGATTCTGATGGCAACACGATCTGGc                                    >  1:51207/1‑133 (MQ=255)
                                                                                                 gCCACTGTGGAGAGCTATGACCCGTTGGCGCGGGTTGACGGTAAACCTGGCAGCGAAGCAATTCACTATTTTCATACCACCTTGGTTGGGCAACCGGAGCAGCTGACTGATTCTGATGGCAACACGATCTGGCg                                   >  7:81561/1‑134 (MQ=255)
                                                                                                 gCCACTGTGGAGAGCTATGACCCGTTGGCGCGGGTTGACGGTAAACCTGGCAGCGAAGCAATTCACTATTTTCATACCACCTTGGTTGGGCAACCGGAGCAGCTGACTGATTCTGATGGCAACACGATCTGGCg                                   >  7:281746/1‑134 (MQ=255)
                                                                                                 gCCACTGTGGAGAGCTATGACCCGTTGGCGCGGGTTGACGGTAAACCTGGCAGCGAAGCAATTCACTATTTTCATACCACCTTGGTTGGGCAACCGGAGCAGCTGACTGATTCTGATGGCAACACGATCTGGCg                                   >  4:51873/1‑134 (MQ=255)
                                                                                                  ccACTGGGGAGAGCTCTGACCCGGTGGCGCGGGTTGACGGTAAACCTGGCAGCGAAGCAATTCACTATTTTCATACCACCTTGGTTGGGCAACCGGAGCAGCTGACTGATTCTGATGGCAACACGATCTGGCGTa                                 <  4:60161/135‑1 (MQ=255)
                                                                                                     cTGTGGAGAGCTATGACCCGTTGGCGCGGGTTGACGGTAAACCTGGCAGCGAAGCAATTCACTATTTTCATACCACCTTGGTTGGGCAACCGGAGCAGCTGCCTGATTCTGATGGCAACACGATCTGGCGTAGTGa                             >  3:271911/1‑136 (MQ=255)
                                                                                                        tgGAGAGCTATGACCCGTTGGCGCGGGTTGACGGTAAACCTGGCAGCGAAGCAATTCACTATTTTCATACCACCTTGGTTGGGCAACCGGAGCAGCTGACTGATTCTGATGGCAACACGATCTGGCGTAGTGAcc                           >  1:317973/1‑135 (MQ=255)
                                                                                                         ggAGAGCTATGACCCGTTGGCGCGGGTTGACGGTAAACCTGGCAGCGAAGCAATTCACTATTTTCATACCACCTTGGTTGGGCAACCGGAGCAGCTGACTGATTCTGATGGCAACACGATCTGGCGTAGTGACCa                          >  1:338472/1‑135 (MQ=255)
                                                                                                         ggAGAGCTATGACCCGTTGGCGCGGGTTGACGGTAAACCTGGCAGCGAAGCAATTCACTATTTTCATACCACCTTGGTTGGGCAACCGGAGCAGCTGACTGATTCTGATGGCAACACGATCTGGCGTAGTGACCAt                         >  7:261909/1‑136 (MQ=255)
                                                                                                           agagCTATGACCCGTTGGCGCGGGTTGACGGTAACCCTGGCAGCGAAGCAATTCACTATTTTCATACCACCTTGGTTGGGCAACCGGAGCAGCTGACTGATTCTGATGGCAACACGATCTGGCGTAGTGACCATCa                       >  8:24136/1‑136 (MQ=255)
                                                                                                           agagCTATGACCCGTTGGCGCGGGTTGACGGTAAACCTGGCAGCGAAGCAATTCACTATTTTCATACCACCTTGGTTGGGCAACCGGAGCAGCTGACTGATTCTGATGGCAACACGATCTGGCGTAGTGACCATCa                       >  4:64311/1‑136 (MQ=255)
                                                                                                           agagCTATGACCCGTTGGCGCGGGTTGACGGTAAACCTGGCAGCGAAGCAATTCACTATTTTCATACCACCTTGGTTGGGCAACCGGAGCAGCTGACTGATTCTGATGGCAACACGATCTGGCGTAGTGACCATCa                       >  3:55969/1‑136 (MQ=255)
                                                                                                             agCTATGACCCGTTGGCGCGGGTTGACGGTAAACCTGGCAGCGAAGCAATTCACTATTTTCATACCACCTTGGTTGGGCAACCGGAGCAGCTGACTGATTCTGATGGCAACACGATCTGGCGTAGTGACCATCAc                      >  1:265607/1‑135 (MQ=255)
                                                                                                              gCTATGACCCGTTGGCGCGGGTTGACGGTAAACCTGGCAGCGAAGCAATTCACTATTTTCATACCACCTTGGTTGGGCAACCGGAGCAGCTGACTGATTCTGATGGCAACACGATCTGGCGTAGTGACCATCACgg                    <  8:385219/136‑1 (MQ=255)
                                                                                                              gCTATGACCCGTTGGCGCGGGTTGACGGTAAACCTGGCAGCGAAGCAATTCACTATTTTCATACCACCTTGGTTGGGCAACCGGAGCAGCTGACTGATTCTGATGGCAACACGATCTGGCGTAGTGACCATCACg                     >  6:47532/1‑135 (MQ=255)
                                                                                                                tATGACCCGTTGGCGCGGGTTGACGGTAAACCTGGCAGCGAAGCAATTCACTATTTTCATACCACCTTGGTTGGGCAACCGGAGCAGCTGACTGATTCTGATGGCAACACGATCTGGCGTAGTGACCATCACGGt                   >  5:158515/1‑135 (MQ=255)
                                                                                                                  tGACCCGTTGGCGCGGGTTGACGGTAAACCTGGCAGCGAAGCAATTCACTATTTTCATACCACCTTGGTTGGGCAACCGGAGCAGCTGACTGATTCTGATGGCAACACGATCTGGCGTAGTGACCATCACGGTTg                 >  6:29077/1‑135 (MQ=255)
                                                                                                                   gACCCGTTGGCGCGGGTTGACGGTAAACCTGGCAGCGAAGCAATTCACTATTTTCATACCACCTTGGTTGGGCAACCGGAGCAGCTGACTGATTCTGATGGCAACACGATCTGGCGTAGTGACCATCACGGTTggg               >  7:265343/1‑136 (MQ=255)
                                                                                                                    aCCCGTTGGCGCGGGTTGACGGTAAACCTGGCAGCGAAGCAATTCACTATTTTCATACCACCTTGGTTGGGCAACCGGAGCAGCTGACTGATTCTGATGGCAACACGATCTGGCGTAGTGACCATCACGGTTgggg              <  3:64311/136‑1 (MQ=255)
                                                                                                                     cccGTTGGCGCGGGTTGACGGTAAACCTGGCAGCGAAGCAATTCACTATTTTCATACCACCTTGGTTGGGCAACCGGAGCAGCTGACTGATTCTGATGGCAACACGATCTGGCGTAGTGACCATCACGGTTgggg              <  2:152805/135‑1 (MQ=255)
                                                                                                                     cccGTTGGCGCGGGTTGACGGTAAACCTGGCAGCGAAGCAATTCACTATTTTCATACCACCTTGGTTGGGCAACCGGAGCAGCTGACTGATTCTGATGGCAACACGATCTGGCGTAGTGACCATCACGGTTgggg              >  7:81702/1‑135 (MQ=255)
                                                                                                                     cccGTTGGCGCGGGTTGACGGTAAACCTGGCAGCGAAGCAATTCACTATTTTCATACCACCTTGGTTGGGCAACCGGAGCAGCTGACTGATTCTGATGGCAACACGATCTGGCGTAGTGACCATCACGGTTggg               >  1:129779/1‑134 (MQ=255)
                                                                                                                     cccGTTGGCGCGGGTTGACGGTAAACCTGGCAGCGAAGCAATTCACTATTTTCATACCACCTTGGTTGGGCAACCGGAGCAGCTGACTGATTCTGATGGCAACACGATCTGGCGTAGTGACCATCACGGTTggg               >  2:354039/1‑134 (MQ=255)
                                                                                                                      ccGTTGGCGCGGGTTGACGGTAAACCTGGCAGCGAAGCAATTCACTATTTTCATACCACCTTGGTTGGGCAACCGGAGCAGCTGACTGATTCTGATGGCAACACGATCTGGCGTAGTGACCATCACGGTTGGGGaa            <  7:313541/136‑1 (MQ=255)
                                                                                                                       cggtGGTGCGGGTTGACGGTAAACCTGGCAGCGAAGCAATTCACTATTTTCATACCACCTTGGTTGGGCAACCGGAGCAGCTGACTGATTCTGATGGCAACACGATCTGGCGTAGTGACCATCACGGTTGGGGaa            <  5:161225/132‑1 (MQ=255)
                                                                                                                       cGTTGGCGCGGGTTGACGGTAAACCTGGCAGCGAAGCAATTCACTATTTTCATACCACCTTGGTTGGGCAACCGGAGCAGCTGACTGATTCTGATGGCa                                                <  8:175257/99‑1 (MQ=255)
                                                                                                                       cGTTGGCGCGGGTTGACGGTAAACCTGGCAGCGAAGCAATTCACTATTTTCATACCACCTTGGTTGGGCAACCGGAGCAGCTGACTGATTCTGATGGCa                                                >  7:175257/1‑99 (MQ=255)
                                                                                                                            gcgcGGGTTGACGGTAAACCTGGCAGCGAAGCAATTCACTATTTTCATACCACCTTGGTTGGGCAACCGGAGCAGCTGACTGATTCTGATGGCAACACGATCTGGCGTAGTGACCATCACGGTTGGGGAAAGAtaa      >  7:267984/1‑136 (MQ=255)
                                                                                                                            gcgcGGGTTGACGGTAAACCTGGCAGCGAAGCAATTCACTATTTTCATACCACCTTGGTTGGGCAACCGGAGCAGCTGACTGATTCTGATGGCAACACGATCTGGCGTAGTGACCATCACGGTTGGGGAAAGAtaa      >  7:237864/1‑136 (MQ=255)
                                                                                                                            gcgcGGGTTGACGGTAAACCTGGCAGCGAAGCAATTCACTATTTTCATACCACCTTGGTTGGGCAACCGGAGCAGCTGACTGATTCTGATGGCAACACGATCTGGCGTAGTGACCATCACGGTTGGGGAAAGAta       >  5:379897/1‑135 (MQ=255)
                                                                                                                            gcgcGGGTTGACGGTAAACCTGGCAGCGAAGCAATTCACTATTTTCATACCACCTTGGTTGGGCAACCGGAGCAGCTGACTGATTCTGATGGCAACACGATCTGGCGTAGTGACCATCACGGTTGGGGAAAGAta       >  5:178983/1‑135 (MQ=255)
                                                                                                                            gcgcGGGTTGACGGTAAACCTGGCAGCGAAGCAATTCACTATTTTCATACCACCTTGGTTGGGCAACCGGAGCAGCTGACTGATTCTGATGGCAACACGATCTGGCGTAGTGACCATCACGGTTGGGGAAAGAta       >  1:57389/1‑135 (MQ=255)
                                                                                                                            gcgcGGGTTGACGGTAAACCTGGCAGCGAAGCAATTCACTATTTTCATACCACCTTGGTTGGGCAACCGGAGCAGCTGACTGATTCTGATGGCAACACGATCTGGCGTAGTGACCATCACGGTTGGGGAAAGAta       >  2:279731/1‑135 (MQ=255)
                                                                                                                             cgcgGGTTGACGGTAAACCTGGCAGCGAAGCAATTCACTATTTTCATACCACCTTGGTTGGGCAACCGGAGCAGCTGACTGATTCTGATGGCAACACGATCTGGCGTAGTGACCATCACGGTTGGGGAAAGAtaa      <  6:379897/135‑1 (MQ=255)
                                                                                                                                  gTTGACGGTAAACCTGGCAGCGAAGCAATTCACTATTTTCATACCACCTTGGTTGGGCAACCGGAGCAGCTGACTGATTCTGATGGCAACACGATCTGGCGTAGTGACCATCACGGTTGGGGAAAGAtaattaa  >  5:371760/1‑134 (MQ=255)
                                                                                                                                     |                                                                                                                                  
GGACTATCCAGAGCCGCGTCGGCGGGTGCTGTTTCACTGGCAAGGGTTTCGCTTGCTTCAAGAAGTGCAGAGTGGTTTGGCTAGCTTGTATGTGTATGCCACTGTGGAGAGCTATGACCCGTTGGCGCGGGTTGACGGTAAACCTGGCAGCGAAGCAATTCACTATTTTCATACCACCTTGGTTGGGCAACCGGAGCAGCTGACTGATTCTGATGGCAACACGATCTGGCGTAGTGACCATCACGGTTGGGGAAAGATAATTAA  >  ACB122/4608989‑4609252

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: