Predicted mutation
evidence seq id position mutation annotation gene description
RA ACB122 120,983 (T)9→8 intergenic (+280/+92) polA → / ← engB DNA polymerase I/putative GTP‑binding protein EngB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB122120,9750T.95.0% 334.5 / 18.6 83intergenic (+272/+100)polA/engBDNA polymerase I/putative GTP‑binding protein EngB
Reads supporting (aligned to +/- strand):  ref base T (1/3);  new base . (47/29);  total (49/34)
Fisher's exact test for biased strand distribution p-value = 2.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACTGAATGGCTGGTGAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCAGTTCCATCCAG  >  ACB122/120849‑121108
                                                                                                                              |                                                                                                                                     
aaCTGAATGGCTGGTGAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGttttttttgcc                                                                                                                             >  3:197461/1‑134 (MQ=255)
     aaTGGCTGGTGAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGaa                                                                                                                       >  3:10292/1‑137 (MQ=255)
     aaTGGCTGGTGAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGaa                                                                                                                       >  3:94647/1‑137 (MQ=255)
     aaTGGCTGGTGAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGaa                                                                                                                       >  2:197927/1‑137 (MQ=255)
        ggCTGGTGAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAAttt                                                                                                                    <  1:265541/137‑1 (MQ=255)
          cTGGTGAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCCCCCTAGCGGTCCGGACTTGGACCCCGAACTTCACCCTCCCCATACGAAGTCCGGGG‑TTATTTTTGCCCGGAATTTGc                                                                                                                  >  6:303489/1‑137 (MQ=255)
                aaGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAg                                                                                                            >  5:314607/1‑137 (MQ=255)
                aaGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAg                                                                                                            >  5:24015/1‑137 (MQ=255)
                     cTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGt                                                                                                       <  8:270060/137‑1 (MQ=255)
                       tCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGttt                                                                                                     <  7:158700/137‑1 (MQ=255)
                       tCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGttt                                                                                                     <  6:350442/137‑1 (MQ=255)
                            gCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTcg                                                                                                <  4:94647/137‑1 (MQ=255)
                               ccTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGa                                                                                             <  8:7510/137‑1 (MQ=255)
                                 tATGTTGTGTTAAGTGTTGGCAGCAATCGGGCCCCCGCCCTAGCGGTCCGGACTTGGCCCCCGAACTTCCCCCTCCCCATACGAAGGCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAAt                                                                                           <  3:27412/137‑1 (MQ=255)
                                 tATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAAt                                                                                           <  8:29702/137‑1 (MQ=255)
                                  aTGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCACCATACGAAGTCCGGGG‑TTTTTTTTGCACGGAATTTGCCGGGAGCGGGTTTTATCGCGAATg                                                                                          >  8:230252/1‑137 (MQ=255)
                                    gttgtGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCg                                                                                        <  6:24015/137‑1 (MQ=255)
                                        tgtTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGt                                                                                    <  3:86024/137‑1 (MQ=255)
                                        tgtTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGt                                                                                    <  5:387653/137‑1 (MQ=255)
                                           tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGa                                                                                 >  8:398368/1‑137 (MQ=255)
                                           tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGa                                                                                 >  5:205543/1‑137 (MQ=255)
                                           tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGa                                                                                 >  5:154868/1‑137 (MQ=255)
                                           tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGa                                                                                 >  6:414339/1‑137 (MQ=255)
                                           tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGa                                                                                 >  4:335043/1‑137 (MQ=255)
                                           tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGa                                                                                 >  4:137204/1‑137 (MQ=255)
                                           tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGa                                                                                 >  7:440530/1‑137 (MQ=255)
                                           tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGa                                                                                 >  8:178947/1‑137 (MQ=255)
                                           tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGa                                                                                 >  1:207514/1‑137 (MQ=255)
                                           tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGa                                                                                 >  2:342174/1‑137 (MQ=255)
                                           tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTAGa                                                                                 >  7:373851/1‑137 (MQ=255)
                                           tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCACCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGa                                                                                 >  7:434547/1‑137 (MQ=255)
                                           tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCACCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGACGGGAGCGGGTTTTCTCGCGAATGCGATGTCGa                                                                                 >  6:277610/1‑137 (MQ=255)
                                           tAAGTGTTGCCAGAATTCTGCACCCCGCCCAAGCCGTCCGGACTTCGACCCCGAACTTCCCCCTCCCCATTAGATGTCCGGCG‑TTTTTTTTGCCCGGAATATCCCGGGACAGCGTTTTCTCGCGAATGCGATGTCGa                                                                                 >  8:365252/1‑137 (MQ=255)
                                             aGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAAt                                                                               >  1:230580/1‑137 (MQ=255)
                                             aGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAAt                                                                               >  5:213512/1‑137 (MQ=255)
                                                gtTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATaca                                                                            >  3:71132/1‑137 (MQ=255)
                                                  tGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACAcc                                                                          >  8:324977/1‑137 (MQ=255)
                                                   ggCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACAccg                                                                         <  5:76488/137‑1 (MQ=255)
                                                     cAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACAccgcc                                                                       >  1:82102/1‑137 (MQ=255)
                                                         tatCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCAt                                                                   <  7:398368/137‑1 (MQ=255)
                                                           tCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑GTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTc                                                                 >  7:176769/1‑137 (MQ=255)
                                                            cTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTcg                                                                >  8:159379/1‑137 (MQ=255)
                                                              ggACCCCTCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTcgcg                                                              <  7:188417/137‑1 (MQ=255)
                                                              ggACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTcgcg                                                              >  4:142301/1‑137 (MQ=255)
                                                              ggACCCCCCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTcgcg                                                              <  2:209214/137‑1 (MQ=255)
                                                                aCCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCggg                                                            >  4:403699/1‑137 (MQ=255)
                                                                  cccGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTa                                                          >  7:316199/1‑137 (MQ=255)
                                                                       ccTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAAccc                                                     >  7:354133/1‑137 (MQ=255)
                                                                       ccTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAAccc                                                     >  7:110593/1‑137 (MQ=255)
                                                                       ccTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAAccc                                                     >  6:195230/1‑137 (MQ=255)
                                                                       ccTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAAccc                                                     >  4:345120/1‑137 (MQ=255)
                                                                       ccTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAAccc                                                     >  1:257372/1‑137 (MQ=255)
                                                                       ccTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCAGTCGCGGGTAAAccc                                                     >  1:408175/1‑137 (MQ=255)
                                                                        cTAGCGGTCCGGACTTGCACCCCGAACATCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGTGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCg                                                    <  1:176090/137‑1 (MQ=39)
                                                                         tAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCgg                                                   >  4:371802/1‑137 (MQ=255)
                                                                                ccGGACTTGGACCCCGAACTTCCCCTTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGgc                                            <  1:16183/137‑1 (MQ=255)
                                                                                ccGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGgc                                            <  7:178947/137‑1 (MQ=255)
                                                                                ccGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGgc                                            <  4:10292/137‑1 (MQ=255)
                                                                                ccGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGgc                                            <  3:44587/137‑1 (MQ=255)
                                                                                ccGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGgc                                            <  5:85190/137‑1 (MQ=255)
                                                                                    aCTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCaa                                        >  2:417898/1‑137 (MQ=255)
                                                                                            cccGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTaa                                >  2:282990/1‑137 (MQ=255)
                                                                                            cccGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTaa                                >  2:228252/1‑137 (MQ=255)
                                                                                            cccGAACTTCCCCCTCCCCATACGAAGT‑CGGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTaa                                >  5:404401/1‑137 (MQ=255)
                                                                                                 aCTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGAc                           >  7:104872/1‑137 (MQ=255)
                                                                                                 aCTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGAc                           >  1:229572/1‑137 (MQ=255)
                                                                                                 aCTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGAc                           >  1:300121/1‑137 (MQ=255)
                                                                                                 aCTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACAGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGAc                           >  8:367296/1‑137 (MQ=255)
                                                                                                  cTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGAcc                          <  8:221054/137‑1 (MQ=255)
                                                                                                    tCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCgg                        <  6:205543/137‑1 (MQ=255)
                                                                                                     ccccccccccATACGAAGTCCGGGG‑TTTTTTTTGCCCGTAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACAAAACGACCGGc                       <  2:396878/137‑1 (MQ=39)
                                                                                                     cccccTCCCCTTACTAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCGGGGTGGGGGTTTTCTCGCGATTGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCGGGGTCTGGCGCAAACCACTACACGACCGGc                       <  6:154868/137‑1 (MQ=255)
                                                                                                     cccccTCCCCATCCGAAGTCC‑GGGTTTTTTTTTGCCCGTAATTTGCCGGGCGGGGGTTTTCTCGGTAATGCGATGTCTAAACCACCGCCGCATTCGCGGGTAAACCGGGGTCAGGCGCAAACCACTACACGACCGGc                       <  1:107992/137‑1 (MQ=255)
                                                                                                     cccccTCCCAATACGAAGTCC‑GGGTTTTTTTTTGCCCGGATTTTGCCGGTAGCGGTTTTTCTCGCGATTGGTATGTCGAATACACGGCCGCATTCGCGGTTAAACCCGGGTCAGGCGCAAACCAATAAACGACCGGc                       <  2:106664/137‑1 (MQ=255)
                                                                                                          cccccTTACGAAGT‑CGGGGGTTTTTTTTGCCCGTATTTTGCCGGGAGCGGGTTTTCTCGCGATTGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCGGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTc                  <  6:3115/136‑1 (MQ=255)
                                                                                                          cccccATACGAAGTCCGGGG‑TTTTTTTTGCCCGGATTTTGCCGGGTGCGGTTTTTCTCGCGAATGCGATGTCCAATACCCCGCGGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTACACGACCGGCTTGTc                  <  4:108220/136‑1 (MQ=255)
                                                                                                           ccccATACTACGTCCGGGGTTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATTTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTc                  <  8:220958/137‑1 (MQ=255)
                                                                                                               aTACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCTAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCa             <  3:137204/137‑1 (MQ=255)
                                                                                                                  cGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCAGtt          >  5:436298/1‑137 (MQ=255)
                                                                                                                   gAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCAGTtc         >  2:83973/1‑137 (MQ=255)
                                                                                                                   cgaattcGGGGGTTTTTTTTGCCCGGATTTTGCCGGGAGGGGGTTTTCTCGCGAATGCGATGTCTAATACACCGCCGCATTCGGGGTTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCAGtt          <  8:104872/131‑1 (MQ=255)
                                                                                                                    aaGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCAGTtcc        <  6:141516/137‑1 (MQ=255)
                                                                                                                       tCCGGGG‑TTTTTTTTGCCCGGATTTTGCCGGGAGCGGGTTTTCTCGCGAATGCTATGTCGAATACACCGCGGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCAGTtccatc     <  6:314607/137‑1 (MQ=255)
                                                                                                                         ccggggTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCAGTtccatcc    >  3:211928/3‑137 (MQ=255)
                                                                                                                          ggggTTTTTTTTT‑ACCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGAGATGTCGAATAAACAGCCGCATTAGCGGGTAAACCCGGGTCATGCGCAAACCAATAAACGACCGGCTTGTCTTCCAGTTCCATACAg  >  8:170070/1‑137 (MQ=255)
                                                                                                                          cggggTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCAGTtccatcca   >  7:443415/2‑137 (MQ=255)
                                                                                                                          cggggTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCAGTtccatcca   >  7:268232/2‑137 (MQ=255)
                                                                                                                          cggggTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGACGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCACACCACTAAACGACAGGCTTGTCTTCCAGTtccatcca   >  4:129759/2‑137 (MQ=255)
                                                                                                                              |                                                                                                                                     
AACTGAATGGCTGGTGAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCAGTTCCATCCAG  >  ACB122/120849‑121108

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