Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 4,247,323 T→C 5.3% M1V (ATG→GTG) † rocE ← Amino‑acid permease RocE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1224,247,3230TC 5.3% 246.0 / 2.8 75M1V (ATG→GTG) †rocEAmino‑acid permease RocE
Reads supporting (aligned to +/- strand):  ref base T (26/45);  new base C (2/2);  total (28/47)
Fisher's exact test for biased strand distribution p-value = 6.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

GCCCTGATTGATGGTTACGCCCGAGCCCATGAACAGGCCGGTTCCAATAACCCCACCCAGCGACAGCATGAAGATGTGCCGGCTTTTCAGGGCGCGGGTGAGCTGGATTTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTT  >  ACB122/4247195‑4247438
                                                                                                                                |                                                                                                                   
gcccTGATTGATGGTTACGCCCGAGCCCATGAACAGGCCGGTTCCAATAACCCCACCCAGCGACAGCATGAAGATGTGCCGGCTTTTCAGGGCGCGGGTGAGCTGGATTTTTTCAACGCTGTCAGTCATGGCacac                                                                                                              <  3:36017/136‑1 (MQ=255)
  ccTGATTGATGGTTACGCCCGAGCCCATGAACAGGCCGGTTCCAATAACCCCACCCAGCGACAGCATGAAGATGTGCCGGCTTTTCAGGGCGCGGGTGAGCTGGATTTTTTCAACGCTGTCAGTCATGGCACACCt                                                                                                            <  3:156751/136‑1 (MQ=255)
      attgatGGTTACGCCCGAGCCCATGAACAGGCCGGTTCCAATAACCCCACCCAGCGACAGCATTAAGATGTGCCGGCTTTTCAGGTCGCGTGTTAGCTGGATTTTTTCAACGCTGTCAGTCATGGCACACCTCGcc                                                                                                        <  4:10104/136‑1 (MQ=255)
           tGGTTACGCCCGAGCCCGTGAACAGGCCGGTTCCAATAACCCCACCCAGCGACAGCATGAAGATGTGCCGGCTTTTCAGGGCGCGGGTGAGCTGGATTTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGc                                                                                                   <  4:203838/136‑1 (MQ=255)
           tGGTTACGCCCGAGCCCATGAACAGGCCGGTTCCAATAACCCCACCCAGCGACAGCATGAAGATGTGCCGGCTTTTCAGGGCGCGGGTGAGCTGGATTTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGc                                                                                                   <  4:136879/136‑1 (MQ=255)
           tGGTTACGCCCGAGCCCATGAACAGGCCGGTTCCAATAACCCCACCCAGCGACAGCATGAAGATGTGCCGGCTTTTCAGGGCGCGGGTGAGCTGGATTTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGc                                                                                                   <  3:68349/136‑1 (MQ=255)
            ggTTACGCCCGACCCCATGCACAGGCCGGTTCCAATAACCCCACCCAGCGACAGCATGAAGATGTGCCGGCTTATCAGGGCGCGGGTGAGCTGGAATTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGc                                                                                                   <  2:299441/135‑1 (MQ=255)
             gTTACGCCCGAGCCCATGAACAGGCCGGTTCCAATAACCCCACCCAGCGACAGCATGAAGATGTGCCGGCTTTTCAGGGCGCGGGTGAGCTGGATTTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGCgg                                                                                                 <  7:317775/136‑1 (MQ=255)
              ttACGCCCGAGCCCATGAACAGGCCGGTTCCAATAACCCCACCCAGCGACAGCATGAAGATGTGCCGGCTTTTCAGGGCGCGGGTGAGCTGGATTTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGCgg                                                                                                 <  2:278739/135‑1 (MQ=255)
              ttACGCCCGAGCCCATGAACAGGCCGGTTCCAATAACCCCACCCAGCGACAGCATGAAGATGTGCCGGCTTTTCAGGGCGCGGGTGAGCTGGATTTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGCgg                                                                                                 <  6:213919/135‑1 (MQ=255)
              ttACGCCCGAGCCCATGAACAGGCCGGTTCCAATAACCCCACCCAGCGACAGCATGAAGATGTGCCGGCTTTTCAGGGCGCGGGTGAGCTGGATTTTTTCAACGCTGTCAGTCACGGCACACCTCGCCATTGCgg                                                                                                 >  5:252875/1‑135 (MQ=255)
              ttACGCCCGAGCCCATGAACAGGCCGGTTCCAATAACCCCACCCAGCGACAGCATGAAGATGTGCCGGCTTTTCAGGGCGCGGGTGAGCTGGATTTTTTCAACGCTGTCAGTCACGGCACACCTCGCCATTGCgg                                                                                                 >  5:257686/1‑135 (MQ=255)
               tACGCCCGAGCCCATGAACAGGCCGGTTCCAATAACCCCACCCAGCGACAGCATGAAGATGTGCCGGCTTTTCAGGGCGCGGGTGAGCTGGATTTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGc                                                                                                >  5:179710/1‑135 (MQ=255)
                aCGCCCGAGCCCATGAACAGGCCGGTTCCAATAACCCCACGCAGCGACAGCATGAAGATGTGCCGGCTTTTCAGGGAGCGGGTGAGCTGGATTTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGcc                                                                                               <  3:142018/135‑1 (MQ=255)
                  gCCCGAGCCCATGAACAGGCCGGTTCCAATAACCCCACCCAGCGACAGCATGAAGATGTGCCGGCTTTTCAGGGCGCGGGTGAGCTGGATTTTTTCAACGCTGTCAGTCACGGCACACCTCGCCATTGCGGCCgg                                                                                             <  6:252875/135‑1 (MQ=255)
                  gCCCGAGCCCATGAACAGGCCGGTTCCAATAACCCCACCCAGCGACAGCATGAAGATGTGCCGGCTTTTCAGGGCGCGGGTGAGCTGGATTTTTTCAACGCTGTCAGTCACGGCACACCTCGCCATTGCGGCCgg                                                                                             <  6:257686/135‑1 (MQ=255)
                   cccGAGCCCATGAACAGGCCGGTTCCAATAACCCCACCCAGCGACAGCATGAAGATGTGCCGGCTTTTCAGGGCGCGTGTGAGCTGGATTTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTg                                                                                           >  7:59383/1‑136 (MQ=255)
                       aGCCCATGAAGAGGCCGGTTCCAATAACCCCACCCAGCGACATCATGAAGATGTGCCGGCTTTTCAGGGCGCTGGTTAGCTGGATTTTTTCTACGCTGTCAGTCATTGCTCTCCTAGCCATTGCGGCCGGTGCGTg                                                                                       <  7:354516/136‑1 (MQ=255)
                       aGCCCATGAACAGGCCGGTTCCAATAACCCCACCCAGCGACAGCATGAAGATGTGCCGGCTTTTCAGGGCGCGGGTGAGCTGGATTTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGt                                                                                        <  5:175408/135‑1 (MQ=255)
                        gCCCATGAACAGGCCGGTTCCAATAACCCCACCCAGCGACAGCATGAAGATGTGCCGGCTTTTCAGGGCGCGGGTGAGCTGGATTTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTg                                                                                       <  1:390228/135‑1 (MQ=255)
                          ccATGAACAGGCCGGTTCCAATAACCCCACCCAGCGACAGCATGAAGATGTGCCGGCTTTTCAGGGCGCGGGTGAGCTGGATTTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTgg                                                                                      >  4:16848/1‑134 (MQ=255)
                          ccATGAACAGGCCGGTTCCAATAACCCCACCCAGCGACAGCATGAAGATGTGCCGGCTTTTCAGGGCGCGGGTGAGCTGGATTTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTg                                                                                       >  2:28398/1‑133 (MQ=255)
                             tGAACAGGCCGGTTCCAATCACCCCACCCAGCGACAGCATGAAGATGTGCCGGCTTTTCAGGGCGCGGGTGAGCTGGATTTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAggg                                                                                 <  4:34846/136‑1 (MQ=255)
                               aaCAGGCCGGTTCCAATAACCCCACCCAGCGACAGCATGAAGATGTGCCGGCTTTTCAGGGCGCGGGTGAGCTGGATTTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGt                                                                                >  1:29265/1‑135 (MQ=255)
                                       ggTTCCAATAACCCCACCCAGCGACAGCATGAAGATGTGCCGGCTTTTCAGGGCGCGGGTGAGCTGGATTTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGa                                                                            <  6:128750/131‑1 (MQ=255)
                                       ggTTCCAATAACCCCACCCAGCGACAGCATGAAGATGTGCCGGCTTTTCAGGGCGCGGGTGAGCTGGATTTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGa                                                                            >  5:128750/1‑131 (MQ=255)
                                       ggTTCCAATAACCCCACCCAGCGACAGCATGAAGATGTGCCGGCTTTTCAGGGCGCGGGTGAGCTGGATTTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGa                                                                        >  1:314309/1‑135 (MQ=255)
                                       ggTTCCAATAACCCCACCCAGCGACAGCATGAAGATGTGCCGGCTTTTCAGGGCGCGGGTGAGCTGGATTTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAg                                                                       >  3:197188/1‑136 (MQ=255)
                                       ggTTCCAATAACCCCACCCAGCGACAGCATGAAGATGTGCCGGCTTTTCAGGGCGCGGGTGAGCTGGATTTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAg                                                                       >  3:190777/1‑136 (MQ=255)
                                         ttCCAATAACCCCACCCAGCGACAGCATGAAGATGTGCCGGCTTTTCAGGGCGCGGGTGAGCTGGATTTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAggg                                                                     <  3:248579/136‑1 (MQ=255)
                                          tCCAATTACCCCACCAAGCGTAAGCATTAAGATGTGCCGGCTTTTGAGGGCGCGGGTGAGCTGGATTTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGTCCGGTGCGTGGGAGGGTTCGAATGAggg                                                                     <  2:189292/135‑1 (MQ=255)
                                          tCCAATAACCCCACCCAGCGACAGCATGAAGATGTGCCGGCTTTTCAGGGAGCGGGTGAGCTGGATTTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAggg                                                                     <  6:333532/135‑1 (MQ=255)
                                           ccAATAACCCCACCCAGCGACAGCATGAAGATGTGCCGGCTTTTCAGGGCGCGGGTGAGCTGGATTTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAgggg                                                                    <  8:131731/135‑1 (MQ=255)
                                                 aCCCCACCCAGCGACAGCATGAAGATGTGCCGGCTTTTCAGGGCGCGGGTGAGCTGGATTTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAgcg                                                             <  4:186827/136‑1 (MQ=255)
                                                     cacccaGCGACAGCATGAAGATGTGCCGGCTTTTCAGGGCGCGGGTGAGCTGGATTTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTAc                                                         >  8:268229/1‑136 (MQ=255)
                                                      acccaGCGTCAGAATGAAGATGTGCAGGCTTTTCAGGGCGCGGGTGAGCGGGATTTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGAGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTAc                                                         >  1:189292/1‑135 (MQ=255)
                                                      acccaGCGACAGCATGAAGATGTGCCGGCTTTTCAGGGCGCGTGTGAGCTGGATTTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACg                                                        <  8:59383/136‑1 (MQ=255)
                                                      acccaGCGACAGCATGAAGATGTGCCGGCTTTTCAGGGCGCGGGTGAGCTGGATTTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACg                                                        <  4:233832/136‑1 (MQ=255)
                                                      acccaGCGACAGCATGAAGATGTGCCGGCTTTTCAGGGCGCGGGTGAGCTGGATTTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACg                                                        <  8:43639/136‑1 (MQ=255)
                                                      acccaGCGACAGCATGAAGATGTGCCGGCTTTTCAGGGCGCGGGTGAGCTGGATTTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACg                                                        <  3:16848/136‑1 (MQ=255)
                                                      acccaGCGACAGCATGAAGATGTGCCGGCTTTTCAGGGCGCGGGTGAGCTGGATTTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACg                                                        <  3:317975/136‑1 (MQ=255)
                                                       cccaGCGACAGCATGAAGATGTGCCGGCTTTTCAGGGCGCGGGTGAGCTGGATTTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACgg                                                       <  7:268229/136‑1 (MQ=255)
                                                          aGCGACAGCATGAAGATGTGCCGGCTTTTCAGGGCGCGGGTGAGCTGGATTTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCa                                                     >  6:270910/1‑135 (MQ=255)
                                                          aGCGACAGCATGAAGATGTGCCGGCTTTTCAGGGCGCGGGTGAGCTGGATTTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCa                                                     >  1:386387/1‑135 (MQ=255)
                                                          aGCGACAGCATGAAGATGTGCCGGCTTTTCAGGGCGCGGGTGAGCTGGATTTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCAt                                                    >  7:129051/1‑136 (MQ=255)
                                                             gACAGCATGAAGATGTGCCGGCTTTTCAGGGCGCGGGTGAGCTGGATTTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGt                                                  >  4:88242/1‑135 (MQ=255)
                                                                     gAAGATGTGCCGGCTTTTCAGGGCGCGGGTGAGCTGGATTTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGa                                                                                    <  7:210890/93‑1 (MQ=255)
                                                                     gAAGATGTGCCGGCTTTTCAGGGCGCGGGTGAGCTGGATTTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGa                                                                                    >  8:210890/1‑93 (MQ=255)
                                                                     gAAGATGTGCCGGCTTTTCAGGGCGCGGGTGAGCTGGATTTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTg                                                                 <  5:109180/112‑1 (MQ=255)
                                                                     gAAGATGTGCCGGCTTTTCAGGGAGCGGGTGAGCTGGATTTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTg                                                                 >  6:109180/1‑112 (MQ=255)
                                                                       aGATGTGCCGGCTTTTCAGGGCGCGGGTGAGCTGGATTTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCtt                                        <  2:324041/135‑1 (MQ=255)
                                                                         aTGTGCCGGCTTTTCAGTGCGCGGGTGAGCTGGATTTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGtc                                     >  3:239077/1‑136 (MQ=255)
                                                                            tgCCGGCTTTTCTGGGCGCGGGTGAGCTGGATTTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCTGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTcc                                  <  8:338450/136‑1 (MQ=255)
                                                                            tgCCGGCTTTTCAGGGCGCGGGTGAGCTGGATTTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTcc                                  <  8:37809/136‑1 (MQ=255)
                                                                              ccGGCTTTTCAGGGCGCGGGTGAGCTTGATTTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCaa                                <  3:267869/136‑1 (MQ=255)
                                                                              ccGGCTTTTCAGGGCGCGGGTGAGCTGGATTTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCaa                                <  4:177856/136‑1 (MQ=255)
                                                                              ccGGCTCTTCAGGGCGCGGGGGAGCTGGATTTTATCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCaa                                <  8:338421/136‑1 (MQ=255)
                                                                               cGGCTTTTCAGGGCGCGGGTGAGCTGGATTTTTTCAACGCTGTCAGTCTTGGATGTCGTCGCCATTGCGGCCGGTGCGTGGGTGGGTTTGAATGAGGGGCTGAGGGTTACGGCATGTCATTATTCTTGTCTCCat                                <  6:151983/135‑2 (MQ=255)
                                                                               cGGCTTTTCAGGGCGCGGGTGAGCTGGATTTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAggg                                                                     <  8:332657/98‑1 (MQ=255)
                                                                               cGGCTTTTCAGGGCGCGGGTGAGCTGGATTTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAggg                                                                     >  7:332657/1‑98 (MQ=255)
                                                                               cGGCTTTTCAGGGCGCGGGTGAGCTGGATTTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCaa                                >  3:65711/1‑135 (MQ=255)
                                                                               cGGCTTTTCAGGGCGCGGGTGAGCTGGATTTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCaa                                <  5:170835/135‑1 (MQ=255)
                                                                               cGGCTTTTCAGGGCGCGGGTGAGCTGGATTTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCaa                                <  2:27611/135‑1 (MQ=255)
                                                                               cGGCTTTTCAGGGCGCGGGTGAGCTGGATTTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCa                                 >  1:27611/1‑134 (MQ=255)
                                                                               cGGCTTTTCAGGGCGCGGGTGAGCTGGATTTTTTCAACGCTGACAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCAGGTCATTATTCTTGTCTCCaa                                >  3:11705/1‑135 (MQ=255)
                                                                               aggCTTTTCAGGGCGGGTGTGAGCTGGATTTTTTCAACGCTTGCAGTCATGGCACACCTCGCCATTGCGGGCGGTGGGTGGGTGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCaa                                <  2:382133/134‑1 (MQ=255)
                                                                                ggCTTTTCAGGGCGCGGGTGAGCTGGATTTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATa                              <  8:348776/136‑1 (MQ=255)
                                                                                ggCTTTTCAGGGCGCGGGTGAGCTGGATTTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCAAAt                               >  3:331669/1‑135 (MQ=255)
                                                                                    tttCAGGGCGCGGGTGAGCTGGATTTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCaa                                >  4:197034/1‑130 (MQ=255)
                                                                                    tttCAGGGCGCGGGTGAGCTGGATTTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCaa                                <  3:197034/130‑1 (MQ=255)
                                                                                               gggTGAGCTGGATTTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCg                >  1:65581/1‑135 (MQ=255)
                                                                                                      cTGGATTTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCaa         <  2:386387/135‑1 (MQ=255)
                                                                                                            tttttCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACg                                                        <  3:332423/82‑1 (MQ=255)
                                                                                                            tttttCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACg                                                        >  4:332423/1‑82 (MQ=255)
                                                                                                              tttCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGGGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGtt  <  5:173206/134‑1 (MQ=255)
                                                                                                                                |                                                                                                                   
GCCCTGATTGATGGTTACGCCCGAGCCCATGAACAGGCCGGTTCCAATAACCCCACCCAGCGACAGCATGAAGATGTGCCGGCTTTTCAGGGCGCGGGTGAGCTGGATTTTTTCAACGCTGTCAGTCATGGCACACCTCGCCATTGCGGCCGGTGCGTGGGAGGGTTCGAATGAGGGGCTGAGCGCTACGGCATGTCATTATTCTTGTCTCCAATAAGCCTCATGGACCGTGTGCAACACGGTT  >  ACB122/4247195‑4247438

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