Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB143 3,013,661 G→A 5.3% A387T (GCG→ACG)  bcsC → Cellulose synthase operon protein C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1433,013,6610GA 5.3% 247.1 / 3.2 75A387T (GCG→ACG) bcsCCellulose synthase operon protein C
Reads supporting (aligned to +/- strand):  ref base G (41/30);  new base A (2/2);  total (43/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.37e-01

AGGACCTGGAGTCGTCCAGCCAGTATTGGGCCCTGCTGAAGAAAGCCGAGCGGGCCCTGGAGGTGAAAGACGTGCGCGGTGCGCGCGCACTTTACCAGCAGGCTCGCCAGCAACGGCCAAAGGATGAATACGCGTTGCTGGGCCTGGCAGATGTCTCGCTGGCCGAGGGCGATGTGGCGGCTGCCGAGCGGCAGTACCTGCAGGTACGGCGTATGGCGCCGGATAACGAAAGTGCCGTGCGTGGCTTGATGCGCGTTTACCAGGC  >  ACB143/3013530‑3013794
                                                                                                                                   |                                                                                                                                     
aGGACCTGGAGTCGTCCAGCCAGTATTGGGCCCTGCTGAAGAAAGCCGAGCGGGCCCTGGAGGTGAAAGACGTGCGCGGTGCGCGCGCACTTTACCAGCAGGCTCGCCAGCAACGGCCAAAGGATGAATACGCGtt                                                                                                                                   <  8:76193/136‑1 (MQ=255)
aGGACCTGGAGTCGTCCAGCCAGTATTGGGCCCTGCTGAAGAAAGCCGAGCGGGCCCTGGAGGTGAAAGACGTGCGCGGTGCGCGCGCACTTTACCAGCAGGCTCGCCAGCAACGGCCAAAGGATGAATACGCGtt                                                                                                                                   <  2:366565/136‑1 (MQ=255)
     cTGGAGTCGTCCAGCCAGTATTGGGCCCTGCTGAAGAAAGCCGAGCGGGCCCTGGAGGTGAAAGACGTGCGCGGTGCGCGCGCACTTTACCAGCAGGCTCGCCAGCAACGGCCAAAGGATGAATACGCGTTGCTgg                                                                                                                              >  8:174980/1‑136 (MQ=255)
         aGTCGTCCAGCCAGTATTGGGCCCTGCTGAAGAAAGCCGAGCGGGCCCTGGAGGTGAAAGACGTGCGCGGTGCGCGCGCACTTTACCAGCAGGCTCGCCAGCAACGGCCAAAGGATGAATACACGTTGCTGGGCCt                                                                                                                          >  4:72706/1‑136 (MQ=255)
         aGTCGTCCAGCCAGTATTGGGCCCTGCTGAAGAAAGCCGAGCGGGCCCTGGAGGTGAAAGACGTGCGCGGTGCGCGCGCACTTTACCAGCAGGCTCGCCAGCAACGGCCAAAGGATGAATAAGCGTTGCTGGGCCt                                                                                                                          >  1:273177/1‑136 (MQ=255)
         aGTCGTCCAGCCAGTATTGGGCCCTGCTGAAGAAAGCCGAGCGGGCCCTGGAGGTGAAAGACGTGCGCGGTGCGCGCGCACTTTACCAGCAGGCTAGCCAGCAACGGCCAAAGGATGAATAAACGTTTCTGGGCCt                                                                                                                          >  4:67082/1‑136 (MQ=255)
             gtcCAGCCAGTATTGGGCCCTGCTGAAGAAAGCCGAGCGGGCCCTGGAGGTGAAAGACGTGCGCGGTGCGCGCGCACTTTACCAGCAGGCTCGCCAGCAACGGCCAAAGGATGAATACGCGTTGCTGGGCCTGGCa                                                                                                                      <  1:349425/136‑1 (MQ=255)
                cagccagTATTGGGCCCTGCTGAAGAAAGCCGAGCGGGCCCTGGAGGTGAAAGACGTGCGCGGTGCGCGCGCACTTTACCAGCAGGCTCGCCAGCAACGGCCAAAGGATGAATACGCGTTGCTGGGCCTGGCAGAt                                                                                                                   >  4:374175/1‑136 (MQ=255)
                cagccagTATTGGGCCCTGCTGAAGAAAGCCGAGCGGGCCCTGGAGGTGAAAGACGTGCGCGGTGCGCGCGCACTTTACCAGCAGGCACGCCCGCAACGGCCAAAGGATGAATACGCGTTGCTGGGCCGGGCAGAt                                                                                                                   >  2:32047/1‑136 (MQ=255)
                cagccagTATTGGGCCCTGCTGAAGAAAGCCGAGCGGGCCCTGGAGGTGAAAGACGTGCGCGGTGCGCGCGCACATTTCCAGCCGGCTCGCCAGCAACGGCCAAACGATGAATACGCGTTGCAGGGCCGGGCAGAt                                                                                                                   >  6:120460/1‑136 (MQ=255)
                   ccagTATTGGGCCCTGCTGAAGAAAGCCGAGCGGGCCCTGGAGGTGAAAGACGTGCGCGGTGCGCGCGCACTTTACCAGCAGGCTCGCCAGCAACGGCCAAAGGATGAATACGCGTTGCTGGGCCTGGCAGATGtc                                                                                                                <  7:249601/136‑1 (MQ=255)
                     agTATTGGGCCCTGCTGAAGAAAGCCGAGCGGGCCCTGGAGGTGAAAGACGTGCGCGGTGCGCGCGCACTTTACCAGCAGGCTCGCCAGCAACGGCCAAAGGATGAATACGCGTTGCTGGGCCTGGCAGATGtctc                                                                                                              >  1:314190/1‑136 (MQ=255)
                          tGGGCCCTGCTGAACAAAGCCGAGCGGGCCCTGGAGGTGAAAGACGTGCGCGGTGCGCGCGCACTTTACCAGCAGGCTCCCCACCAACGCGCAAAGCATGAATAAGCGTTGCTGGGCCAGGCAAATGTCTCGCtcg                                                                                                         >  4:343305/1‑134 (MQ=255)
                               ccTGCTGAAGAAAGCCGAGCGGGCCCTGGAGGTGAAAGACGTGCGCGGTGCGCGCGCACTTTTCCAGCAGGCTCGCCAGCAACCGCCAAAGGATGAATACGCGTTGCTGGGCCAGGCAGATGTATCCCTGGCCGAg                                                                                                    >  4:181427/1‑136 (MQ=255)
                               ccTGCTGAAGAAAGCCGAGCGGGCCCTGGAGGTGAAAGACGTGCGCGGTGCGCGCGCACTTTACCAGCAGGCTCGCCAGCAACGGCCAAAGGATGAATACGCGTTGCTGGGCCTGGCAGATGTCTCGCTGGCCGAg                                                                                                    >  7:281052/1‑136 (MQ=255)
                               ccTGCTGAAGAAAGCCGAGCGGGCCCTGGAGGTGAAAGACGTGCGCGGTGCGCGCGCACTTTACCAGCAGGCTCGCCAGCAACGGCCAAAGGATGAATACGCGTTGCTGGGCCTGGCAGATGTCTCGCTGGCCGAg                                                                                                    >  2:323522/1‑136 (MQ=255)
                               ccTGCTGAAGAAAGCCGAGCGGGCCCTGGAGGTGAAAGACGTGCGCGGTGCGCGCGCACGTTACCAGCAGGCTCGCCAGCAACGGCCAAAGGATGAATACGCGTTGCTGGGCCTGGCAGATGTCTCGCTGGCCGAg                                                                                                    >  3:221915/1‑136 (MQ=255)
                                 tgctgAAGAAAGCCGAGCGGGCCCTGGAGGTCAAAGACGTGCGCGGTGCGCGCGCACTTTACCAGCAGGCTCGCCAGCAACCGCCAAACGATGAATACGCGTTGCTGGGCCTGGCAGATGTCTCGCTGGCCGAcgg                                                                                                  >  5:334640/1‑136 (MQ=255)
                                    tgAAGAAAGCCGAGCGGGCCCTGGAGGTGAAAGACGTGCGCGGTGCGCGCGCCCTTTACCAGCAGGCTCGCCAGCAACGGCCAAAGGATGAATACGCGTTGCTGGGCCTGGCAGATGTCTCGGTGGCCGAGGGCGa                                                                                               >  1:395866/1‑136 (MQ=255)
                                    tgAAGAAAGCCGAGCGGGCCCTGGAGGTGAAAGACGTGCGCGGTGCGCGCGCACTTTACCAGCAGGCTCGCCAGCAACGGCCAAAGGATGAATACGCGTTGCTGGGCCTGGCAGATGTCTCGCTGGCCGAGGGCGa                                                                                               >  1:22905/1‑136 (MQ=255)
                                    tgAAGAAAGCCGAGCGGGCCCTGGAGGTGAAAGACGTGCGCGGTGCGCGCGCACTTTACCAGCAGGCTCGCCAGCAACGGCCAAAGGATGAATACGCGTTGCTGGGCCTGGCAGATGTCTCGCTGGCCGAGGGCGa                                                                                               >  7:67401/1‑136 (MQ=255)
                                    tgAAGAAAGCCGAGCGGGCCCTGGAGGTGAAAGACGTGCGCGGTGCGCGCGCACTTTACCAGCAGGCTCGCCAGCAACGGCCAAAGGATGAATACGCGTTGCTGGGCCTGGCAGATGTCTCGCTGGCCGAGGGCGa                                                                                               >  8:173067/1‑136 (MQ=255)
                                    tgAAGAAAGCCGAGCGGGCCCTGGAGGTGAAAGACGTGCGCGGTGCGCGCGCACTTTACCAGCAGGCTCGCCAGCAACGGCCAAAGGATGAATACGCGTTGCTGGGCCTGGCAGATGTCTCGCTGGCCGAGGGCGa                                                                                               >  1:269870/1‑136 (MQ=255)
                                    tgAAGAAAGCCGAGCGGGCCCTGGAGGTGAAAGACGTGCGCGGTGCGCGCGCACTTTACCAGCAGGCTCCCCAGCAACGGCCAAAGGATGAATACGCGTTGCTGGGCCTGGCAGATGTCTCGCTGGCCGAGGGCGa                                                                                               >  4:221478/1‑136 (MQ=255)
                                    tgAAGAAAGCCGAGCGGGCCCTCGAGGTGAAAGACGTGCGCGGTGCGCGCGCACTTTACCAGCAGGCTCCCCAGCAACGCCCAAAGGATGAATACGCGTTGCAGGGCCTGGCAGATGTATCGCTGGCCGCGGGCGa                                                                                               >  6:310580/1‑136 (MQ=255)
                                      aagaaAGCCGAGCGGGCCCTGGAGGTGAAAGACGTGCGCGGTGCGCGCGCACTTTACCAGCAGGCTCGCCAGCAACGGCCAAAGGATGAATACGCGTTGCTGGGCCTGGCAGATGTCTCGCTGGCCGAGGGCGAtg                                                                                             >  8:40589/1‑136 (MQ=255)
                                      aagaaAGCCGAGCGGGCCCTGGAGGTGAAAGACGTGCGCGGTGCGCGCGCACTTTACCAGCAGGCTCGCCAGCAACGGCCAAAGGATGAATACGCGTTGCTGGGCCTGGCAGATGTCTCGCTGGCCGAGGGCGAtg                                                                                             >  4:328682/1‑136 (MQ=255)
                                               gAGCGGGCCCTGGAGGTGAAAGACGTGCGCGGTGCGCGCGCACTTTACCAGCAGGCTCGCCAGCAACGGCCAAAGGATGAATACGCGTTGCTGGGCCTGGCAGATGTCTCGCTGGCCGAGGGCGATGTGGCGGCTg                                                                                    <  4:54218/136‑1 (MQ=255)
                                                aGCGGGCCCTGGAGGTGAAAGACGTGCGCGGTGCGCGCGCACTTTACCAGCAGGCTCGCCAGCAACGGCCAAAGGATGAATACGCGTTGCTGGGCCTGGCAGATGTCTCGCTGGCCGAGGGCGATGTGGCGGCTGc                                                                                   <  5:79820/136‑1 (MQ=255)
                                                 gCGGGCCCTGGAGGTGAAAGACGTGCGCGGTGCGCGCGCACTTTACCAGCAGGCTCGCCAGCAACGGCCAAAGGATGAATACGCGTTGCTGGGCCTGGCAGATGTCTCGCTGGCCGAGGGCGATGTGGCGGCTGcc                                                                                  <  6:104857/136‑1 (MQ=255)
                                                 gCGGGCCCTGGAGGTGAAAGACGTGCGCGGTGCGCGCGCACTTTACCAGCAGGCACGCCAGCAACGGCCAAAGGATGAATACGCGTTGCTGGGCCTGGCAGATGTCTCGCTGGCCGAGGGCGATGTGGCGGCTGcc                                                                                  <  2:379920/136‑1 (MQ=255)
                                                  cGGGCCCTGGAGGTGAAAGACGTGCGCGGTGCGCGCGCACTTTACCAGCAGGCTCGCCAGCAACGGCCAAAGGATGAATACGCGTTGCTGGGCCTGGCAGATGTCTCGCTGGCCGAGGGCGATGTGGCGGCTGCCg                                                                                 <  7:289274/136‑1 (MQ=255)
                                                   gggCCCTGGAGGTGAAAGACGTGCGCGGTGCGCGCGCACTTTACCAGCAGGCTCGCCAGCAACGGCCAAAGGATGAATACGCGTTGCTGGGCCTGGCAGATGTCTCGCTGGCCGAGGGCGATGTGGCGGCTGcc                                                                                  >  1:316488/1‑134 (MQ=255)
                                                    ggCCCTGGAGGTGAAAGACGTGCGCGGTGCGCGCGCACTTTACCAGCAGGCTCGCCAGCAACGGCCAAAGGATGAATACGCGTTGCTGGGCCTGGCAGATGTCTCGCTGGCCGAGGGCGATGTGGCGGCTGCCGAg                                                                               <  2:233993/136‑1 (MQ=255)
                                                      cccTGGAGGTGAAAGACGTGCGCGGGGCGCGCGCACTTTACCAGCAGGCTCGCCAGCAACGGCCAAAGGATGAATACGCGTTGCTGGGCCTGGCAGATGTCTCGCTGGCCGAGGGCGATGTGGCGGCTGCCGAGCg                                                                             <  3:148503/136‑1 (MQ=255)
                                                        cTGGAGGTGAAAGACGTGCGCGGTGCGCGCGCACTTTACCAGCAGGCTCGCCAGCAACGGCCAAAGGATGAATACGCGTTGCTGGGCCTGGCAGATGTCTCGCTGGCCGAGGGCGATGTGGCGGCTGCCGAGCGGc                                                                           >  6:12964/1‑136 (MQ=255)
                                                           gAGGTGAAAGACGTGCGCGGTGCGCGCGCACTTTACCAGCAGGCTCGCCAGCAACGGCCAAAGGATGAATACGCGTTGCTGGGCCTGGCAGATGTCTCGCTGGCCGAGGGCGATGTGGCGGCTGCCGAGCGGCa                                                                          <  7:214318/134‑1 (MQ=255)
                                                             ggTGAAAGACGTGCGCGGTGCGCGCGCACTTTACCAGCAGGCTCGCCAGCAACGGCCAAAGGATGAATACGCGTTGCTGGGCCTGGCAGATGTCTCGCTGGCCGAGGGCGATGTGGCGGCTGCCGAGCGGCAGTAc                                                                      <  4:94981/136‑1 (MQ=255)
                                                                  aaGACGTGCGCGGTGCGCGCGCACTTTACCAGCAGGCTCGCCAGCAACGGCCAAAGGATGAATACGCGTTGCTGGGCCTGGCAGATGTCTCGCTGGCCGAGGGCGATGTGGCGGCTGCCGAGCGGCAGTACCTg                                                                   >  6:264082/1‑134 (MQ=255)
                                                                  aaGACGTGCGCGGTGCGCGCGCACTTTACCAGCAGGCTCGCCAGCAACGGCCAAAGGATGAATACGCGTTGCTGGGCCTGGCAGATGTCTCGCTGGCCGAGGGCGATGTGGCGGCTGCCGAGCGGCAGTACCTg                                                                   >  3:154106/1‑134 (MQ=255)
                                                                   aGACGTGCGCGGTGCGCGCGCACTTTACCAGCAGGCTCGCCAGCAACGGCCAAAGGATGAATACGCGTTGCTGGGCCTGGCAGATGTCTCGCTGGCCGAGGGCGATGTGGCGGCTGCCGAGCGGCAGTACCTGCAg                                                                <  3:161181/136‑1 (MQ=255)
                                                                    gACGTGCGCGGTGCGCGCGCACTTTACCAGCAGGCTCGCCAGCAACGGCCAAAGGATGAATACGCGTTGCTGGGCCTGGCAGATGTCTCGTTGGCCGAGGGCGATGTGGCGGCTGCCGAGCGGCAGTACCTGCAgg                                                               >  1:290792/1‑136 (MQ=255)
                                                                    gACGTGCGCGGTGCGCGCGCACTTTACCAGCAGGCTCGCCAGCAACGGCCAAAGGATGAATACGCGTTGCTGGGCCTGGCAGATGTCTCGCTGGCCGAGGGCGATGTGGCGGCTGCCGAGCGGCAGTACCTGCAgg                                                               <  1:291938/136‑1 (MQ=255)
                                                                     aCGTGCGCGGTGCGCGCGCACTTTACCAGCAGGCTCGCCAGCAACGGCCAAAGGATGAATACGCGTTGCTGGGCCTGGCAGATGTCTCGCTGGCCGAGGGCGATGTGGCGGCTGCCGAGCGGCAGTACCTGCAGGt                                                              <  8:329914/136‑1 (MQ=255)
                                                                          cgcggtgcgcgCGCACTTTACCAGCAGGCTCGCCAGCAACGGCCAAAGGATGAATACGCGTTGCTGGGCCTGGCAGATGTCTCGCTGGCCGAGGGCGATGTGGCGGCTGCCGAGCGGCAGTACCTGCAGGTACGGc                                                         <  2:143468/136‑1 (MQ=255)
                                                                           gcggtgcgcgCGCACTTTTCCAGCAGGCTCGCCAGCAACGGCCAAAGGCTGAATACGCGTTGCTGGGCCTGGCAGCTGTCTCGCTGGCCGAGGGCGATGTGGCGGCTGCCGAGCGGCAGTTCCTGCAGGTACGGCg                                                        >  5:212023/1‑136 (MQ=255)
                                                                           gcggtgcgcgCGCACTTTACCAGCAGGCTCGCCAGCACCGGCCAAAGGATGAATACGCGTTGCTGGGCCTGGCAGATGTCTCGCTGGCCGAGGGCGATGTGGCGGCTGCCGAGCGGCAGTACCTGCAGGTACGGCg                                                        >  1:233106/1‑136 (MQ=255)
                                                                           gcggtgcgcgCGCACTTTACCAGCAGGCTCGCCAGCAACGGCCAAAGGATGAATACGCGTTGCTg                                                                                                                               >  1:319195/1‑65 (MQ=255)
                                                                           gcggtgcgcgCGCACTTTACCAGCAGGCTCGCCAGCAACGGCCAAAGGATGAATACGCGTTGCTg                                                                                                                               <  2:319195/65‑1 (MQ=255)
                                                                               tgcgcgGGCACTTTACCAGCAGGCTCGCCAGCAACGGCCAAAGGATGAATACGCGTTGCTGGGCCTGGCAGATGTCTCGCTGGCCGAGGGCGATGTGGCGGCTGCCGAGCGGCAGTACCTGCAGGTACGGCGTATg                                                    <  5:256024/136‑1 (MQ=255)
                                                                               tgcgcgCGCACTTTACCAGCAGGCTCGCCAGCAACGGCCAAAGGATGAATACGCGTTGCTGGGCCTGGCAGATGTCTCGCTGGCCGAGGGCGATGTGGCGGCTGCCGAGCGGCAGTACCTGCAGGTACGGCGTATg                                                    <  7:254683/136‑1 (MQ=255)
                                                                               tgcgcgCGCACTTTACCAGCAGGCTCGCCAGCAACGGCCAAAGGATGAATACGCGTTGCTGGGCCTGGCAGATGTCTCGCTGGCCGAGGGCGATGTGGCGGCTGCCGAGCGGCAGTACCTGCAGGTACGGCGTATg                                                    <  6:188654/136‑1 (MQ=255)
                                                                               tgcgcgCGCACTTTACCAGCAGGCGCGCCAGCAACGGCCAAAGGATGAATACGCGTTGCTGGGCCTGGCAGATGTCTCGCTGGCCGAGGGCGATGTGGCGGCTGCCGAGCGGCAGTACCTGCAGGTACGGCGTATg                                                    <  1:295554/136‑1 (MQ=255)
                                                                                 cgcgcgcACTTTACCAGCAGGCTCGCCAGCAACGGCCAAAGGATGAATACACGTTGCTGGGCCTGGCAGATGTCTCGCTGGCCGAGGGCGATGTGGCGGCTGCCGAGCGGCAGTACCTGCAGGTACGGCGTATgg                                                   <  3:67082/135‑1 (MQ=255)
                                                                                 cgcgcgcACTTTACCAGCAGGCTCGCCAGCAACGGCCAAAGGATGAATACACGTTGCTGGGCCTGGCAGATGTCTCGCTGGCCGAGGGCGATGTGGCGGCTGCCGAGCGGCAGTACCTGCAGGTACGGCGTATgg                                                   <  3:72706/135‑1 (MQ=255)
                                                                                  gcgcgcACTTTAACAGCAGGCTCGCCCGCAACGGCCAAAGGATGAATACGCGTTGCTGGGCCTGGCAGATGTCTCGCTGGCCGAGGGCGATGTGGCGGCTGCCGAGCCGCAGTAACTGCAGGTACGGCGTATTgcg                                                 >  2:238109/1‑136 (MQ=255)
                                                                                    gcgcACTTTACCAGCAGGCTCGCCAGCAACGGCCAAAGGATGAATACGCGTTGCTGGGCCTGGCAGATGTCTCGCTGGCCGAGGGCGATGTGGCGGCTGCCGAGCGGCAGTACCTGCAGGTACGGCGTATGGCGcc                                               >  8:384504/1‑136 (MQ=255)
                                                                                     cgcACTTTACCAGCAGGCTCGCCAGCAACGGCCAAAGGATGAATACGCGTTGCTGGGCCTGGCAGATGTCTCGCTGGCCGAGGGCGATGTGGCGGCTGCCGAGCGGCAGTACCTGCAGGTACGGCGTATGGCGCCg                                              <  7:144289/136‑1 (MQ=255)
                                                                                             aCCAGCAGGCTCGCCAGCAACGGCCAAAGGATGAATACGCGTTGCTGGGCCTGGCAGATGTCTCGCTGGCCGAGGGCGATGTGGCGGCTGCCGAGCGGCAGTACCTGCAGGTACGGCGTATGGCGCCGGATAACg                                       <  1:218778/135‑1 (MQ=255)
                                                                                              ccAGCAGGCTCGCCAGCAACGGCCAAAGGATGAATACGCGTTGCTGGGCCTGGCAGATGTCTCGCTGGCCGAGGGCGATGTGGCGGCTGCCGAGCGGCAGTACCTGCAGGTACGGCGTATGGCGCCGGATAACGa                                      >  3:317166/1‑135 (MQ=255)
                                                                                               cagcagGCTCGCCAGCAACGGCCAAAGGATGAATACGCGTTGCTGGGCCTGGCAGATGTCTCGCTGGCCGAGGGCGATGTGGCGGCTGCCGAGCGGCAGTACCTGCAGGTACGGCGTATGGCGCCGGATAACGaaa                                    <  5:308691/136‑1 (MQ=255)
                                                                                               cagcagGCTCGACAGCAACGGCCAAAGGATGAATACGCGTTGCTGGGCCTGGCAGATGTCTCGCTGGCCGAGGGCGATGTGGCGGCTGCCGAGCGGCAGTACCTGCAGGTACGGCGTATGGCGCCGGATAACGaaa                                    <  4:186024/136‑1 (MQ=255)
                                                                                                     gCTCGCCAGCAACGGCCAAAGGATGAATACGCGTTGCTGGGCCTGGCAGATGTCTCGCTGGCCGAGGGCGATGTGGCGGCTGCCGAGCGGCAGTACCTGCAGGTACGGCGTATGGCGCCGGATAACGAAAGTGCcg                              >  1:256044/1‑136 (MQ=255)
                                                                                                     gCTCGCCAGCAACGGCCAAAGGATGAATACGCGTTGCTGGGCCTGGCAGATGTCTCGCTGGCCGAGGGCGATGTGGCGGCTGCCGAGCGGCAGTACCTGCAGGTACGGCGTATGGCGCCGGATAACGAAAGTGCcg                              <  4:299090/136‑1 (MQ=255)
                                                                                                       tCGCCAGCAACGGCCAAAGGATGAATACGCGTTGCTGGGCCTGGCAGATGTCTCGCTGGCCGAGGGCGATGTGGCGGCTGCCGAGCGGCAGTACCTGCAGGTACGGCGTATGGCGCCGGATAACGAAAGTGCcgtg                            >  7:127990/1‑136 (MQ=255)
                                                                                                         gCCAGCAACGGCCAAAGGATGAATACGCGTTGCTGGGCCTGGCAGATGTCTCGCTGGCCGAGGGCGATGTGGCGGCTGCCGAGCGGCAGTACCTGCAGGTACGGCGTATGGCGCCGGATAACGAAAGTGCcgtgcg                          >  5:260506/1‑136 (MQ=255)
                                                                                                         gCCAGCAACGGCCAAAGGATGAATACGCGTTGCTGGGCCTGGCAGATGTCTCGCTGGCCGAGGGCGATGTGGCGGCTGCCGAGCGGCAGTACCTGCAGGTACGGCGTATGGCGCCGGATAACGAAAGTGCcgtgcg                          >  7:152251/1‑136 (MQ=255)
                                                                                                         gCCAGCAACGGCCAAAGGATGAATACGCGTTGCTGGGCCTGGCAGATGTCTCGCTGGCCGAGGGCGATGTGGCGGCTGCCGAGCGGCAGTACCTGCAGGTACGGCGTATGGCGCCGGATAACGAAAGTGCcgtgcg                          >  1:281555/1‑136 (MQ=255)
                                                                                                               aaCGGCCAAAGGATGAATACGCGTTGCTGGGCCTGGCAGATGTCTCGCTGGCCGAGGGCGATGTGGCGGCTGCCGAGCGGCAGTACCTGCAGGTACGGCGTATGGCGCCGGATAACGAAAGTGCCGTGCGTGGCtt                    <  5:251488/136‑1 (MQ=255)
                                                                                                                    ccAAAGGATGAATACGCGTTGCTGGGCCTGGCAGATGTCTCGCTGGCCGAGGGCGATGTGGCGGCTGCCGAGCGGCAGTACCTGCAGGTACGGCGTATGGCGCCGGATAACGAAAGTGCCGTGCGTGGCTTGATgc               <  2:273177/136‑1 (MQ=255)
                                                                                                                      aaaGGATGAATACGCGTTGCTGGGCCTGGCAGATGTCTCGCTGGCCGAGGGCGATGTGGCGGCTGCCGAGCGGCAGTACCTGCAGGTACGGCGTATGGCGCCGGATAACGAAAGTGCcg                              <  4:308588/119‑1 (MQ=255)
                                                                                                                      aaaGGATGAATACGCGTTGCTGGGCCTGGCAGATGTCTCGCTGGCCGAGGGCGATGTGGCGGCTGCCGAGCGGCAGTACCTGCAGGTACGGCGTATGGCGCCGGATAACGAAAGTGCcg                              >  3:308588/1‑119 (MQ=255)
                                                                                                                                 aCGCGTTGCTGGGCCTGGCAGATGTCTCGCTGGCCGAGGGCGATGTGGCGGCTGCCGAGCGGCAGTACCTGCAGGTACGGCGTATGGCGCCGGATAACGAAAGTGCCGTGCGTGGCTTGATGCGCGTTTACCAGgc  >  1:80337/1‑136 (MQ=255)
                                                                                                                                 aCGCGTTGCTGGGCCTGGCAGATGTCTCGCTGGCCGAGGGCGATGTGGCGGCTGCCGAGCGGCAGTACCTGCAGGTACGGCGTATGGCGCCGGATAACGAAAGTGCCGTGCGTGGCTTGATGCGCGTTTACCAGgc  >  3:331714/1‑136 (MQ=255)
                                                                                                                                 aCGCGTTGCTGGGCCTGGCAGATGTCTCGCTGGCCGAGGGCGATGTGGCGGCTGCCGAGCGGCAGTACCTGCAGGTACGGCGTATGGCGCCGGATAACGAAAGTGCCGTGCGTGGCTTGATGCGCGTTTACCAGgc  >  3:237930/1‑136 (MQ=255)
                                                                                                                                   |                                                                                                                                     
AGGACCTGGAGTCGTCCAGCCAGTATTGGGCCCTGCTGAAGAAAGCCGAGCGGGCCCTGGAGGTGAAAGACGTGCGCGGTGCGCGCGCACTTTACCAGCAGGCTCGCCAGCAACGGCCAAAGGATGAATACGCGTTGCTGGGCCTGGCAGATGTCTCGCTGGCCGAGGGCGATGTGGCGGCTGCCGAGCGGCAGTACCTGCAGGTACGGCGTATGGCGCCGGATAACGAAAGTGCCGTGCGTGGCTTGATGCGCGTTTACCAGGC  >  ACB143/3013530‑3013794

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: