Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 1,209,298 G→T 5.7% G238G (GGC→GGA)  rhlB_1 ← ATP‑dependent RNA helicase RhlB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1221,209,2980GT 5.7% 218.1 / 3.8 70G238G (GGC→GGA) rhlB_1ATP‑dependent RNA helicase RhlB
Reads supporting (aligned to +/- strand):  ref base G (43/23);  new base T (2/2);  total (45/25)
Fisher's exact test for biased strand distribution p-value = 6.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.69e-01

ACCAGCGTGGTGCGGTGTTCACGGGCAAGGGTGGCCAAGCGGTCATACACCAGGCCCCAGACATCAATGGCCATCACCGCCCCCAGCGGTACCGGCGGCACTTCGATGGCCAGGTCACGCTTGCGGGCATGGCCGATATCGACAATGGCGCAGGGACGCCCGCTACCCACCAGAAACTGTGCCACCCGTTCGACCGGGCGTTGGGTAGCGGACAGGCCGATACGCCGCAGTGGCTTTTTGCACAAGGCTTGCAGGCGTTCCAAG  >  ACB122/1209167‑1209430
                                                                                                                                   |                                                                                                                                    
aCCAGCGTGGTGCGGTGTTCACGGGCAAGGGTGGCCAAGCGGTCATACACCAGGCCCCAGACATCAATGGCCATCACCGCCCCCAGCGGTACCGGCGGCACTTCGATGGCCAGGTCACGCTTGCGGGCATGGCCg                                                                                                                                   <  5:340361/135‑1 (MQ=255)
  cAGCGTGGTGCGGTGTTCACGGGCAAGGGTGGCCAAGCGGTCATACACCAGGCCCCAGACATCAATGGCCATCACCGCCCCCAGCGGTACCGGCGGCACTTCGATGGCCAGGTCACGCTTGCGGGCATGGCCGat                                                                                                                                 >  6:292791/1‑135 (MQ=255)
  cAGCGTGGTGCGGTGTTCACGGGCAAGGGTGGCCAAGCGGTCATACACCAGGCCCCAGACATCAATGGCCATCACCGCCCCCAGCGGTACCGGCGGCACTTCGATGGCCAGGTCACGCTTGCGGGCATGGCCGat                                                                                                                                 <  3:90868/135‑1 (MQ=255)
       tggtgCGGTGTTCACGGGCAAGGGTGGCCAAGCGGTCATACACCAGGCCCCAGACATCAATGGCCATCACCGCCCCCAGCGGTACCGGCGGCACTTCGATGGCCAGGTCACGCTTGCGGGCATGGCCGATATCg                                                                                                                             >  7:321792/1‑134 (MQ=255)
       tggtgCGGTGTTCACGGGCAAGGGTGGCCAAGCGGTCATACACCAGGCCCCAGACATCAATGGCCATCACCGCCCCCAGCGGTACCGGCGGCACTTCGATGGCCAGGTCACGCTTGCGGGCATGGCCGATATCGa                                                                                                                            <  6:88581/135‑1 (MQ=255)
             ggTGTTCACGGGCAAGGGTGGCCAAGCGGTCATACACCAGGCCCCAGACATCAATGGCCATCACCGCCCCCAGCGGTACCGGCGGCACTTCGATGGCCAGGTCACGCTTGCGGGCATGGCCGATATCGACAATgg                                                                                                                      >  7:333825/1‑135 (MQ=255)
               tgtTCACGGGCAAGGGTGGCCAAGCGGTCATACACCAGGCCCCAGCCATCAATGGCCATCACCGCCCCCAGCGGTACCGGCGGCACTTCGATGGCCAGGTCACGCTTGCGGGCATGGCCGATATCGACAATGgcg                                                                                                                    >  1:290242/1‑135 (MQ=255)
                    aCGGGCAAGGGTGGCCAAGCGGTCATACACCAGGCCCCAGACATCAATGGCCATCACCGCCCCCAGCGGTACCGGCGGCACTTCGATGGCCAGGTCACGCTTGCGGGCATGGCCGATATCGACAATGGCGCAggg                                                                                                               >  6:192718/1‑135 (MQ=255)
                       ggCAAGGGTGTCCAAGCGCTCATACACCAGGCCCCAGACATCAATGGCCATCACCGCCCCCAGCGGTACCGGCGGCACTTCGATGGCCAGGTCACGCTTGCGGGCATGGCCGATATCGACAATGGCGCAGGGACg                                                                                                            >  6:278168/1‑135 (MQ=255)
                       ggCAAGGGTGGCCAAGCGGTCATACACCAGGCCCCAGACATCAATGGCCATCACCGCCCCCAGCGGTACCGGCGGCACTTCGATGGCCAGGTCACGCTTGCGGGCATGGCCGATATCGACAATGGCGCAGGGACg                                                                                                            >  8:252237/1‑135 (MQ=255)
                           aGGGTGGCCAAGCGGTCATACACCAGGCCCCAGACATCAATGGCCATCACCGCCCCCAGCGGTACCGGCGGCACTTCGATGGCCAGGTCACGCTTGCGGGCATGGCCGATATCGACAATGGCGCAGGGACGCCCg                                                                                                        <  6:349499/135‑1 (MQ=255)
                              gTGGCCAAGCGGTCATACACCAGGCCCCAGACATCAATGGCCATCACCGCCCCCAGCGGTACCGGCGGCACTTCGATGGCCAGGTCACGCTTGCGGGCATGGCCGATATCGACAATGGCGCAGGGACGCCCGCTa                                                                                                     >  1:199839/1‑135 (MQ=255)
                              gTGGCCAAGCGGTCATACACCAGGCCCCAGACATCAATGGCCATCACCGCCCCCAGCGGTACCGGCGGCACTTCGATGGCCAGGTCACGCTTGCGGGCATGGCCGATATCGACAATGGCGCAGGGACGCCCGCTa                                                                                                     >  5:33562/1‑135 (MQ=255)
                                    aaGCGGTCATACACCAGGCCCCAGACATCAATGGCCATCACCGCCCCCAGCGGTACCGGCGGCACTTCGATGGCCAGGTCACGCTTGCGGGCATGGCCGATATCGACAATGGCGCAGGGACGCCCGCTAcccg                                                                                                 >  5:284322/1‑132 (MQ=255)
                                    aaGCGGTCATACACCAGGCCCCAGACATCAATGGCCATCACCGCCCCCAGCGGTACCGGCGGCACTTCGATGGCCAGGTCACGCTTGCGGGCATGGCCGATATCGACAATGGCGCAGGGACGCCCGCTAcccc                                                                                                 >  2:307628/1‑132 (MQ=255)
                                    aaGCGGTCATACACCAGGCCCCAGACATCAATGGCCATCACCGCCCCCAGCGGTACCGGCGGCACTTCGATGGCCAGGTCACGCTTGCGGGCATGGCCGATATCGACAATGGCGCAGGGACGCCCGCTACcca                                                                                                 >  2:302645/1‑133 (MQ=255)
                                    aaGCGGTCATACACCAGGCCCCAGACATCAATGGCCATCACCGCCCCCAGCGGTACCGGCGGCACTTCGATGGCCAGGTCACGCTTGCGGGCATGGCCGATATCGACAATGGCGCAGGGACGCCCGCTACcca                                                                                                 >  2:307624/1‑133 (MQ=255)
                                     aGCGGTCATACACCAGGCCCCAGACATCAATGGCCATCACCGCCCCCAGCGGTACCGGCGGCACTTCGATGGCCAGGTCACGCTTGCGGGCATGGCCGATATCGACAATGGCGCAGGGACGCCCGCTACcca                                                                                                 >  7:14235/1‑132 (MQ=255)
                                         gTCAGACACCAGGCCCCAGACATCAATGGCCATCACCGCCCCCAGCGGTCCCGGCGGCACTTCGATGGCCAGGTCACGCGTGCGGGCATGGCCGATATCGACAATGGCGCAGGGACGCCCGCGGCCCACCAGacc                                                                                          >  5:47518/1‑133 (MQ=255)
                                          tCATACACCAGGCCCCAGACATCAATGGCCATCACCGCCCCCAGCGGTACCGGCGGCACTTCGATGGCCAGGTCACGCTTGCGGGCATGGCCGATATCGACAATGGCGCAGGGACGCCCGCTACCCACCAGAAAc                                                                                         >  1:85861/1‑135 (MQ=255)
                                               cacCAGGCCCCAGACATCAATGGTCATCACCGCCCCCAGCGGTACCGGCGGCACTTCGATGGCCAGGTCACGCTTGCGGGCATGGCCGATATCGACAATGGCGCAGGGACGCCCGCTACCCACCAGAAACtgtg                                                                                     >  1:257206/1‑134 (MQ=255)
                                                       cccAGACATCAATGGCCATCACCGCCCCCAGCGGTACCGGCGGCACTTCGATGGCTAGGTCACGCTTGCGGGCATGGCCGATATCGACAATGGCGCAGGGACGCCCGCt                                                                                                      >  3:33564/1‑109 (MQ=255)
                                                       cccAGACATCAATGGCCATCACCGCCCCCAGCGGTACCGGCGGCACTTCGATGGCTAGGTCACGCTTGCGGGCATGGCCGATATCGACAATGGCGCAGGGACGCCCGCt                                                                                                      <  4:33564/109‑1 (MQ=255)
                                                       cccAGACATCAATGGCCATCACCGCCCCCAGCGGTACCGGCGGCACTTCGATGGCCAGGTCACGCTTGCGGTCATGGCCGATATCGACAATGGCGCAGGGACGCCCGCTACCCACCAGAAACTGTGCCCCCCGtt                                                                            >  6:113545/1‑135 (MQ=255)
                                                       cccAGACATCAATGGCCATCACCGCCCCCAGCGGTACCGGCGGCACTTCGATGGCCAGGTCACGCTTGCGGGCATGTCCGATATCGACAATGGCGCAGGGACGCCCGCTACCCACCAGAAACTGTGCCACCCGtt                                                                            >  7:309807/1‑135 (MQ=255)
                                                       cccAGACATCAATGGCCATCACCGCCCCCAGCGGTACCGGCGGCACTTCGATGGCCAGGTCACGCTTGCGGGCATGTCCGATATCGACAATGGCGCAGGGACGCCCGCTACCCACCAGAAACTGTGCCACCCGtt                                                                            >  7:314983/1‑135 (MQ=255)
                                                       cccAGACATCAATGGCCATCACCGCCCCCAGCGGTACCGGCGGCACTTCGATGGCCAGGTCACGCTTGCGGGCATGGCCGATATCGACAATGGCGCAGGGACGCCCGCTACCCACCAGAAACTGTGCCACCCGtt                                                                            >  3:166314/1‑135 (MQ=255)
                                                       cccAGACATCAATGGCCATCACCGCCCCCAGCGGTACCGGCGGCACTTCGATGGCCAGGTCACGCTTGCGGGCATGGCCGATATCGACAATGGCGCAGGGACGCCCGCTACCCACCAGAAACTGTGCCACCCGtt                                                                            >  7:188712/1‑135 (MQ=255)
                                                       cccAGACATCAATGGCCATCACCGCCCCCAGCGGTACCGGCGGCACTTCGATGGCCAGGTCACGCTTGCGGGCATGGCCGATATCGACAATGGCGCAGGGACGCCCGCTACCCACCAGAAACTGTGCCACCCGtt                                                                            >  5:71524/1‑135 (MQ=255)
                                                       cccAGACATCAATGGCCATCACCGCCCCCAGCGGTACCGGCGGCACTTCGATGGCCAGGTCACGCTTGCGGGCATGGCCGATATCGACAATGGCGCAGGGACGCCCGCTACCCACCAGAAACTGTGCCACCCGtt                                                                            >  2:260051/1‑135 (MQ=255)
                                                       cccAGACATCAATGGCCATCACCGCCCCCAGCGGTACCGGCGGCACTTCGATGGCCAGGTCACGCTTGCGGGCATGGCCAATATCGACAATGGCGCAGGGACGc                                                                                                           <  8:60057/104‑1 (MQ=255)
                                                       cccAGACATCAATGGCCATCACCGCCCCCAGCGGTACCGGCGGCACTTCGATGGCCAGGTCACGCTTGCGGGCATGGCCAATATCGACAATGGCGCAGGGACGc                                                                                                           >  7:60057/1‑104 (MQ=255)
                                                       cccAGACATCAATGGCCATCACCG‑CCCCAGCGGTACCGGCGGCACTTCGATGGCCAGGTCACGCTTGCGGGCATGGCCGATATCGACAATGGCGCAGGGACGCCCGCTACCCACCAGAAACTGTGCCACCCGtt                                                                            >  3:285126/1‑134 (MQ=255)
                                                       cccAGACATCAATGGACATCAACGCCCCCAGCTGTACCGGCGGCACTTCGATGGCCAGCACACGCTTGCGGGCATGGCCGATATCGACAATGGCGAAGGGACGCCCGCTCACCACCAGACACTGTTCCACCCgaa                                                                            >  5:165895/1‑133 (MQ=255)
                                                             cATCAATGGCCATCACCGCCCCCAGCGGTACCGGCGGCACTTCGATGGCCAGGTCACGCTTGCGGGCATGGCCGATATCGACAATGGCGCAGGGACGCCCGCTACCCACCAGAAACTGTGCCACCCGTTCGACCg                                                                      >  7:330599/1‑135 (MQ=255)
                                                              aTCAATGGCCATCACCGCCCCCAGCGTTACCGGCGGCACTTCGATGGCCAGGTCACGCTTGCGGGCATGGCCGATATCGACAATGGCGCAGGGACGCCCGCTACCCACCAGAAACTGTGCCACCCGTTCGACCgg                                                                     <  4:95756/135‑1 (MQ=255)
                                                                cAATGGCCATCACCGCCCCCAGCGGTACCGGCGGCACTTCGATGGCCAGGTCACGCTTCCTGGCATGGCCGATATCGACAATGGCGCAGGGACGCCCGCTACCCACCAGAAACTGTGCCACCCGTTCGACCGGGc                                                                   <  3:180241/135‑1 (MQ=255)
                                                                    ggCCATCACCGCCCCCAGCGGTACCGGCGGCACTTCGATGGCCAGGTCACGCTTGCGGGCATGGCCGATATCGACAATGGCGCAGGGACGCCCGCTACCCACCAGAAACTGTGCCACCCGTTCGACCGGGCGTTg                                                               >  3:300075/1‑135 (MQ=255)
                                                                      ccATCACCGCCCCCAGCGGTACCGGCGGCACTTCGATGGCCAGGTCACGCTTGCGGGCATGTCCGATATCGACAATGGCGCAGGGACGCCCGCTACCCACCAGAAACTGTGCCACCCGTTCGACCGGGCGTTgg                                                              <  8:314983/134‑1 (MQ=255)
                                                                      ccATCACCGCCCCCAGCGGTACCGGCGGCACTTCGATGGCCAGGTCACGCTTGCGGGCATGTCCGATATCGACAATGGCGCAGGGACGCCCGCTACCCACCAGAAACTGTGCCACCCGTTCGACCGGGCGTTgg                                                              <  8:309807/134‑1 (MQ=255)
                                                                          ccccGCCCCCCGCGGTCCCGGCGGCCCTGCGAGGCCCAGGGCAGGCTGGCGGGCAGGGCCTATAGCGACAAGGGCGCCGGGACGCCCGCTACCCACCAGAAACTGTGCCACCCGTGCGACCGGGCGGTGGGTAGc                                                         <  5:302457/133‑1 (MQ=255)
                                                                          ccccGCCCCCAGCGGTACCGGCGGCACTTCGAGGGCCAGGTCACGCTTGCGGGCATGGCCGATATCGACAATGGCGCAGGGACGCCCGCTACCCACCAGAAACTGTGCCACCCGTTCGACCGGGCGTTGGGTAGc                                                         <  5:381999/133‑1 (MQ=255)
                                                                          cACCGCCCCCAGCGGTACCGGCGGCACTTCGATGGCCAGGTCACGCTTGCGTGCATGGCCGATATCGACAATGGCGCAGGGACGCCCGCTACCCACCAGAAACTGTGCCACCCGTTCGACCGGGCGTTGGGTAGc                                                         <  1:260051/135‑1 (MQ=255)
                                                                          cACCGCCCCCAGCGGTACCGGCGGCACTTCGATGGCCAGGTCACGCTTGCGGGCATGGCCGATATCGACAATGGCGCAGGGAGGCCCGCTACCCACCAGAAACTGTGCCACCCGTTCGACCGGGCGTTGGGTAGc                                                         <  6:33562/135‑1 (MQ=255)
                                                                          cACCGCCCCCAGCGGTACCGGCGGCACTTCGATGGCCAGGTCACGCTTGCGGGCATGGCCGATATCGACAATGGCGCAGGGACGCCCGCTACCCACCAGAAACTGTGCCACCCGTTCGACCGGGCGTTGGGTAGc                                                         >  8:304216/1‑135 (MQ=255)
                                                                           aCCGCCCCCAGCGGTACCGGCGGCACTTCGATGGCCAGGTCACGCTTGCGGGCATGGCCGATATCGACAATGGCGCAGGGACGCCCGCTACCCACCAGACACAGTGCCACCCGGTCGCCCGGGCGTTGGGTAGCg                                                        >  5:184278/1‑135 (MQ=255)
                                                                           aCCGCCCCCAGCGGTACCGGCGGCACTTCGATGGCCAGGTCACGCTTGCGGGCATGGCCGATATCGACAATGGCGCAGGGACGCCCGCTACCCACCAGAAACTGTGCCACCCGTTCGACCGGGCGTTGGGTAGc                                                         >  2:102623/1‑134 (MQ=255)
                                                                           aCCGCCCCCAGCGGTACCGGCGGCACTTCGATGGCCAGGTCACGCTTGCGGGCATGGCCGATATCGACAATGGCGCAGGGACGCCCGCTACCCACCAGAAACTGTGCCACCCGTTCGACCGGGCGTTGGGTAGCg                                                        >  2:251404/1‑135 (MQ=255)
                                                                           aCCGCCCCCAGCGGTACCGGCGGCACTTCGATGGCCAGGTCACGCTTGCGGGCATGGCCGATATCGACAATGGCGCAGGGACGCCCGCTACCCACCAGAAACTGTGCCACCCGTTCGACCGGGCGTTGGGTAGCg                                                        >  3:67278/1‑135 (MQ=255)
                                                                           aCCGCCCCCAGCGGTACCGGCGGCACTTCGATGGCCAGGTCACGCTTGCGGGCATGGCCGATATCGACAATGGCGCAGGGACGCCCGCTACCCACCACAAAATGTGCCACCCGTTGGACCGGGCGGTGGGGAGCg                                                        >  5:242107/1‑135 (MQ=255)
                                                                                  ccAGCGGTACCGGCGGCACTTCGATGGCCAGGTCACGCTTGCGGGCATGGCCGATATCGACAATGGCGCAGGGACGCCCGCTACCCACCAGAAACTGTGCCACCCGTTCGACCGGGCGTTGGGTAGCGGACAGGc                                                 >  7:175350/1‑135 (MQ=255)
                                                                                    aggggTACCGGGGTCACTTGGATGGCCAGGTCACGCTTGCGGGCATGGCCGATATCTACAATGGCGCAGGGACGCCCGCTACCCACCAGAAACTGTGCCACCCGTTCGACCGGGCGTTGGGTAGCGGACAGGCCg                                               <  5:329455/132‑1 (MQ=255)
                                                                                       ggTACCGGCGGCACTTCGATGGCCAGGTCACGCTTGCGGGCATGGCCGATATCGACAATGGCGCAGGGACGCCCGCTACCCACCAGAAACTGTGCCACCCGTTCGACCGGGCGTTGGGTAGCGGACAGGCCGATa                                            <  7:91820/135‑1 (MQ=255)
                                                                                       ggTACCGGCGGCACTTCGATGGCCAGGTCACGCTTGCGGGCATGGCCGATATCGACAATGGCGCAGGGACGCCCGCTACCCACCAGAAACTGTGCCACCCGTTCGACCGGGCGTTGGGTAGCGGACAGGCCGATa                                            >  8:91820/1‑135 (MQ=255)
                                                                                       ggTACCGGCGGCACTTCGATGGCCAGGTCACGCTTGCGGGCATGGCCGATATCGACAATGGCGCAGGGACGCCCGCTACCCACCAGAAACTGTGCCACCCGTTCGACCGGGCGTTGGGTAGCGGACAGGCCGATa                                            >  8:169837/1‑135 (MQ=255)
                                                                                       ggTACCGGCGGCACTTCGATGGCCAGGTCACGCTTGCGGGCATGGCCGATATCGACAATGGCGCAGGGACGCCCGCTACCCACCAGAAACTGTGCCACCCGTTCGACCGGGCGTTGGGTAGCGGACAGGCCGATa                                            >  6:17813/1‑135 (MQ=255)
                                                                                       ggTACCGGCGGCAATTCGATGGCCAGGTCACGCTTGCGGGCATGGCCGATATCGACAATGGCGCAGGGCCGCCCGCTACCCCCCAGCAACTGTGCCACCCGTTCGACCGGGCTGTGGGGAGCGGACAGGCCGATa                                            >  5:373322/1‑135 (MQ=255)
                                                                                                      tCGATGGCCAGGTCACGCTTGCGGGCATGGCCGATATCGACAATGGCGCAGGGACGCCCGCTACCCACCAGAAACTGTGCCACCCGTTCGACCGGGCGTTGGGTAGCGGACAGGCCGATACGCCGCAGTGGCttt                             >  3:313666/1‑135 (MQ=255)
                                                                                                      tCGATGGCCAGGTCACGCTTGCGGGCATGGCCGATATCGACAATGGCGCAGGGACGCCCGCTACCCACCAGAAACTGTGCCACCCGTTCGACCGGGCGTTGGGTAGCGGACAGGCCGATACGCCGCAGTGGCttt                             >  5:350407/1‑135 (MQ=255)
                                                                                                        gATGGCCAGGTCACGCTTGCGGGCATGGCCGATATCGACAATGGCGCAGGGACGCCCGCTACCCACCAGAAACTGTGCCACCCGTTCGACCGTGCGTTGGGTAGCGGACAGGCCGATACGCCGCAGTGGCttttt                           <  3:326921/135‑1 (MQ=255)
                                                                                                         aTGGCCAGGTCACGCTTGCGGGCATGGCCGATATCGACAATGGCGCAGGGACGCCCGCTACCCACCAGAAACTGTGCCACCCGTTCGACCGGGCGTTGGGTAGCGGACAGGCCGATACGCCGCAGTGGCttttt                           <  1:236753/134‑1 (MQ=255)
                                                                                                          tGGCCAGGTCACGCTTGCGGGCATGGCCGATATCGACAATGGCGCAGGGACGCCCGCTACCCACCAGAAACTGTGCCACCCGTTCGACCGGGCGTTGGGTAGCGGACAGGCCGATACGCCGCAGTGGCTTTTTGc                         <  5:292791/135‑1 (MQ=255)
                                                                                                          tGGCCAGGTCACGCTTGCGGGCATGGCCGATATCGACAATGGCGCAGGGACGCCCGCTACCCACCAGAAACTGTGCCACCCGTTCGACCGGGCGTTGGGTAGCGGACAGGCCGATACGCCGCAGTGGCTTTTTGc                         <  2:199839/135‑1 (MQ=255)
                                                                                                          tGGCCAGGTCACGCTTGCGGGCATGGCCGATATCGACAATGGCGCAGGGACGCCCGCTACCCACCAGAAACTGTGCCACCCGTTCGACCGGGCGTTGGGTAGCGGACAGGCCGATACGCCGCAGTGGCTTTTTGc                         <  1:251404/135‑1 (MQ=255)
                                                                                                                     cGCTTGCGGGCATGGCCGATATCGACAATGGCGCAGGGACGCCCGCTACCCACCAGAAACTGTGCCACCCGTTCGACCGGGCGTTGGGTAGCGGACAGGCCGATACGCCGCAGTGGCTTTTTGCACAAGGCTTGc              <  8:297543/135‑1 (MQ=255)
                                                                                                                     cGCTTGCGGGCATGGCCGATATCGACAATGGCGCAGGGACGCCCGCTACCCACCAGAAACTGTGCCACCCGTTCGACCGGGCGTTGGGTAGCGGACAGGCCGATACGCCGCAGTGGCTTTTTGCACAAGGCTTGc              <  2:53770/135‑1 (MQ=255)
                                                                                                                            gggCATGGCCGATATCGACAATGGCGCAGGGACGCCCGCTACCCACCAGAAACTGTGCCACCCGTTCGACCGGGCGTTGGGTAGCGGACAGGCCGATACGCCGCAGTGGCTTTTTGCACAAGGCTTGCAGGCGtt       <  5:17813/135‑1 (MQ=255)
                                                                                                                              gCATGGCCGATATCGACAATGGCGCAGGGACGCCCGCTACCCACCAGAAACTGTGCCACCCGTTCGACCGGGCGTTGGGTAGCGGACAGGCCGATACGCCGCAGTGGCTTTTTGCACAAGGCTTGCAGGCGTTc      >  3:310667/1‑134 (MQ=255)
                                                                                                                              gCATGGCCGATATCGACAATGGCGCAGGGACGCCCGCTACCCACCAGAAACTGTGCCACCCGTTCGACAGGGCGTTGGGGAGCGGACAGGCCGATACGCCGCAGTGGCTTGTTGCACAAGCCGTGCAGGCGtta      >  5:338498/1‑133 (MQ=255)
                                                                                                                                 tGGCCGATATCGACAATGGCGCAGGGACGCCCGCTACCCACCAGAAACTGTGCCACCCGTTCGACCGGGCGTTGGGTAGCGGACAGGCCGATACGCCGCAGTGGCTTTTTGCACAAGGCTTGCAGGCGTTCCAAg  <  1:102623/135‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                    
ACCAGCGTGGTGCGGTGTTCACGGGCAAGGGTGGCCAAGCGGTCATACACCAGGCCCCAGACATCAATGGCCATCACCGCCCCCAGCGGTACCGGCGGCACTTCGATGGCCAGGTCACGCTTGCGGGCATGGCCGATATCGACAATGGCGCAGGGACGCCCGCTACCCACCAGAAACTGTGCCACCCGTTCGACCGGGCGTTGGGTAGCGGACAGGCCGATACGCCGCAGTGGCTTTTTGCACAAGGCTTGCAGGCGTTCCAAG  >  ACB122/1209167‑1209430

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