Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 120,975 Δ1 bp 100% intergenic (+272/+100) 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.90.8% 248.2 / 19.3 65intergenic (+272/+100)polA/engBDNA polymerase I/putative GTP‑binding protein EngB
Reads supporting (aligned to +/- strand):  ref base T (3/3);  new base . (39/20);  total (42/23)
Fisher's exact test for biased strand distribution p-value = 6.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.

CTCAGAGTAACTGAATGGCTGGTGAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCAGTTCCATC  >  ACB122/120841‑121105
                                                                                                                                      |                                                                                                                                  
cTCAGAGTAACTGAATGGCTGGTGAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGCTATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACGCCGAACTACCCCCTCCCCATACGAAGTCCGGGGt                                                                                                                                    >  7:381359/1‑135 (MQ=255)
           tGAATGGCTGGTGAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCTGACATGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGttttttttgccc                                                                                                                         >  5:56565/1‑131 (MQ=255)
           tGAATGGCTGGTGAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGttttttttgccc                                                                                                                         >  2:4215/1‑131 (MQ=255)
                  cTGGTGAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAAttt                                                                                                                 >  3:300771/1‑135 (MQ=255)
                        aaGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCg                                                                                                             >  6:243524/1‑133 (MQ=255)
                        aaGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCACTAGCGGTCCTGACTTGGACACCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTGTTGAGCGGAATTTGCCg                                                                                                             >  1:363478/1‑133 (MQ=255)
                              ttCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAATTCCCGGGTTTTTTTTT‑CGAGGAATTTGACGGGAGCcgg                                                                                                     >  5:83177/1‑135 (MQ=255)
                                  aTGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGttt                                                                                                  <  4:300771/134‑1 (MQ=255)
                                        cTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTcgc                                                                                            >  7:213849/1‑134 (MQ=255)
                                        cTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTcgc                                                                                            >  7:340064/1‑134 (MQ=255)
                                        cTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTcgc                                                                                            >  7:61563/1‑134 (MQ=255)
                                         tATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGa                                                                                          <  8:30563/135‑1 (MQ=255)
                                            gttgtGTTACGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCGATAGGAAGTCCGGTG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCTAAt                                                                                        <  5:129367/134‑1 (MQ=255)
                                            gttgtGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAAt                                                                                        <  4:256852/134‑1 (MQ=255)
                                            gttgtGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAAt                                                                                        <  5:37378/134‑1 (MQ=255)
                                              tgtgTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGc                                                                                      <  5:243524/134‑1 (MQ=255)
                                              tgtgTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGc                                                                                      <  6:109800/134‑1 (MQ=255)
                                               gtgtTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGGCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCg                                                                                     <  4:221389/134‑1 (MQ=255)
                                                tgtTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGAt                                                                                   >  5:266394/1‑135 (MQ=255)
                                                  ttAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGt                                                                                 >  8:279044/1‑135 (MQ=255)
                                                   tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTT‑CCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGt                                                                                 >  3:12789/1‑134 (MQ=255)
                                                   tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGt                                                                                 >  2:125092/1‑134 (MQ=255)
                                                   tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGt                                                                                 >  5:137378/1‑134 (MQ=255)
                                                   tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGt                                                                                 >  5:200265/1‑134 (MQ=255)
                                                   tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGAGAATGCTATGt                                                                                 >  3:162549/1‑134 (MQ=255)
                                                   tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCACGGAATTTGACGGGAGCGGGTTTTCTCGCGAATGCGAGGt                                                                                 >  6:315116/1‑134 (MQ=255)
                                                     aGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGa                                                                              >  4:78074/1‑135 (MQ=255)
                                                     aGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGa                                                                              >  7:43098/1‑135 (MQ=255)
                                                      gtgtTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGa                                                                              <  3:265527/134‑1 (MQ=255)
                                                        gtTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAAt                                                                            >  3:322539/1‑134 (MQ=255)
                                                        gtTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAAt                                                                            >  1:270532/1‑134 (MQ=255)
                                                          tGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATaca                                                                         >  2:144500/1‑135 (MQ=255)
                                                             cAGATATCTGGCCCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGCTGTCGAATgcacag                                                                      >  2:82167/1‑133 (MQ=255)
                                                             cAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACAccg                                                                      <  6:102925/135‑1 (MQ=255)
                                                                        cccccGCCCTAGCGGCCCGGACTTGGCCCCCGAACTTCCCCGTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGTATTTTGCCGGGAGCGGGTTTTCTCGCGAATGCTATGCCGAATACACCGCCGCATTcgcg                                                           <  1:109860/134‑1 (MQ=255)
                                                                        aCCCCGCCCTCGCGGCCCGGACTTGGACCCCGAACTTCCCCCTCCCCATCCGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTcgcg                                                           <  1:181962/135‑1 (MQ=255)
                                                                        aCCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTcgcg                                                           <  6:200265/135‑1 (MQ=255)
                                                                         ccccGCCCTAGCGGTCCGGACTTGGAACCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTT‑CCCGGAATTTGGCGGGAGCGGGTTTTATCGCCAATGCGATGTAGAATACCCCGCCTCATTCGCgg                                                          >  4:228168/1‑135 (MQ=255)
                                                                             gCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTaa                                                      <  4:322539/135‑1 (MQ=255)
                                                                             ccccAAGCGGTCCGGACTTGCACCCCTAACTTCCCCCTCCCCATCCGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTaa                                                      <  2:270532/134‑1 (MQ=255)
                                                                               ccTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAAc                                                    >  1:367645/1‑135 (MQ=255)
                                                                               ccTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAAc                                                    >  6:49728/1‑135 (MQ=255)
                                                                               ccTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGGGTTCTCGCGACTGCGATGTCGAATACACCGCCGCATTCGCGGGTaacc                                                    >  7:191903/1‑133 (MQ=255)
                                                                               ccTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑‑TTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAAcc                                                   >  8:287609/1‑135 (MQ=255)
                                                                                 tAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAcccc                                                  >  2:128287/1‑135 (MQ=255)
                                                                                    cGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCg                                                 <  3:62297/133‑1 (MQ=255)
                                                                                     gggccGGCCTTGGACCCCGAACTTCCCCCTCCCCATACGACGCCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATCCACCGCCGCATTCGCGGGTAAACCCGGGt                                              <  2:337809/132‑1 (MQ=38)
                                                                                       tCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCa                                            >  3:190193/1‑135 (MQ=255)
                                                                                        ccGGACTTGGACCCCGAACTTCCCCCTCCCCATACGCAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAg                                           <  4:270512/135‑1 (MQ=255)
                                                                                         cGGACTTGGACCCCGCACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTGTTTTTGACCGGAATTTGCCGGGAGCGGGTTTTCGCGCGAAAGCGAGGTCGCATTCACCGCCGCCTGAGCGGGGGAACCCggggccgg                                          >  5:306230/1‑130 (MQ=255)
                                                                                         cGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑‑TTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGgc                                         <  7:287609/135‑1 (MQ=255)
                                                                                            aCTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGGATTTGACGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAATCCCGGGTCAGgcgc                                       >  3:280095/1‑135 (MQ=255)
                                                                                            aCTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGgcg                                        >  2:245116/1‑134 (MQ=255)
                                                                                                   ccccGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCa                                 >  3:266694/1‑134 (MQ=255)
                                                                                                    cccGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑GTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTACACCCGGCTCAGGCGCAAACCAc                                >  6:15570/1‑134 (MQ=255)
                                                                                                         aCTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGCATGCGATGTCGAATACACCCCCCCATTCGCGGGTAAACCCGGGTCCGGCGCAAACCACTAAACg                          >  1:369836/1‑135 (MQ=255)
                                                                                                         aCTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCCGGCGCAAACCACTAAACg                          >  3:155046/1‑135 (MQ=255)
                                                                                                         aCTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCCAACCACTAAACg                          >  3:162521/1‑135 (MQ=255)
                                                                                                         aCTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACg                          >  6:195210/1‑135 (MQ=255)
                                                                                                         aCTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACg                          >  8:128299/1‑135 (MQ=255)
                                                                                                         aCTTCCCCCTCCCCATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTAGCGGATGCGATGACGTATACACCCCCGCATTCGACGGTAAACCCGGGGAATCCCCAAACCACAAAACg                          >  2:374569/1‑135 (MQ=255)
                                                                                                         aCTTCCCCCTCCCCATACGAAGTCCGGGG‑TGTTTTTTGCGCGGAATTTGCCGGGAGCGTGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGAATTCGCGGGTACACCCGGGTCAGGCGCAAACCACTAAACg                          >  4:19460/1‑135 (MQ=255)
                                                                                                               cccGCCCCATACGACG‑GCGGGGTTTTTTTTTCCCGGTAATTTGCGGGTCGGGGTTTTTCTCGCGATTGCTTTGGCGAATCCACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGc                    <  8:340064/135‑1 (MQ=255)
                                                                                                                    cccATACGAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTc               <  7:279044/135‑1 (MQ=255)
                                                                                                                           gAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCAGt        >  1:38333/1‑135 (MQ=255)
                                                                                                                           gAAGTCCGGGG‑TTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCAGt        >  7:186185/1‑135 (MQ=255)
                                                                                                                                  ggggTTTTTTTTTGCCCGTATTTTGCCGGGAGGGGTTTTTCTCGCGAATGCGATTGCCAAAACACCGCCTCATTCGGGGTAAAACCGGGGTCAGGCGAAAACCACAACACGACCGGCTTGTCTTCCAGTtccatc  <  3:363245/135‑1 (MQ=255)
                                                                                                                                  ggggTTTTTTTTTGCCCGTAATTTGCCGGGAGGGGTTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCAGTtccatc  <  6:137378/135‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                  
CTCAGAGTAACTGAATGGCTGGTGAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCAGTTCCATC  >  ACB122/120841‑121105

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