Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,376,291 T→C 13.1% intergenic (+26/‑26) efp → / → ecnA polyproline‑specific translation elongation factor EF‑P/entericidin A membrane lipoprotein, antidote entericidin B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,376,2910TC13.1% 244.9 / 23.9 84intergenic (+26/‑26)efp/ecnApolyproline‑specific translation elongation factor EF‑P/entericidin A membrane lipoprotein, antidote entericidin B
Reads supporting (aligned to +/- strand):  ref base T (36/37);  new base C (5/6);  total (41/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.89e-01

GGTGGCAAACCGGCTACCCTGTCTACTGGCGCTGTGGTTAAAGTTCCGCTGTTTGTACAAATCGGCGAAGTCATCAAAGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTTTGCTTGCCAGCACGCTGCTCACGGGCTGTAACACCGCTCGCGGTTTCGGCGAAGACATCAAACATCTCG  >  NC_000913/4376146‑4376413
                                                                                                                                                 |                                                                                                                          
ggtggCAAACCGGCTACCCTGTCTACTGGCGCTGTGGTTAAAGTTCCGCTGTTTGTACAAATCGGCGAAGTCATCAAAGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCAGCCTGCAGGCCGc                                                                                                                          <  1:131756/148‑1 (MQ=255)
                 ccTGTCTACTGGCGCTGTGGTTAAAGTTCCGCTGTTTGTACAAATCGGCGAAGTCATCAAAGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGc                                                                                                                          >  2:756418/1‑131 (MQ=255)
                 ccTGTCTACTGGCGCTGTGGTTAAAGTTCCGCTGTTTGTACAAATCGGCGAAGTCATCAAAGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGc                                                                                                                          <  1:756418/131‑1 (MQ=255)
                   tGTCTACTGGCGCTGTGGTTAAAGTTCCGCTGTTTGTACAAATCGGCGAAGTCATCAAAGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTc                                                                                                       >  1:483803/1‑148 (MQ=255)
                    gTCTACTGGCGCTGTGGTTAAAGTTCCGCTGTTTGTACAAATCGGCGAAGTCATCAAAGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCa                                                                                                      >  1:603046/1‑148 (MQ=255)
                       tACTGGCGCTGTGGTTAAAGTTCCGCTGTTTGTACAAATCGGCGAAGTCATCAAAGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCagag                                                                                                   <  2:603046/148‑1 (MQ=255)
                           ggCGCTGTGGTTAAAGTTCCGCTGTTTGTACAAATCGGCGAAGTCATCAAAGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGc                                                                                                                          >  4:1106662/1‑121 (MQ=255)
                           ggCGCTGTGGTTAAAGTTCCGCTGTTTGTACAAATCGGCGAAGTCATCAAAGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGc                                                                                                                          <  3:1106662/121‑1 (MQ=255)
                                      tAAAGTTCCGCTGTTTGTACAAATCGGCGAAGTCATCAAAGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTAt                                                                                                              >  4:1095972/1‑122 (MQ=255)
                                      tAAAGTTCCGCTGTTTGTACAAATCGGCGAAGTCATCAAAGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTAt                                                                                                              <  3:1095972/122‑1 (MQ=255)
                                        aaGTTCCGCTGTTTGTACAAATCGGCGAAGTCATCAAAGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCAGCCTGCAGGCCGCa                                                                                                                         >  3:653593/1‑109 (MQ=255)
                                        aaGTTCCGCTGTTTGTACAAATCGGCGAAGTCATCAAAGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCAGCCTGCAGGCCGCa                                                                                                                         <  4:653593/109‑1 (MQ=255)
                                            tCCGCTGTTTGTACAAATCGGCGAAGTCATCAAAGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGtt                                                                              <  2:483803/148‑1 (MQ=255)
                                                 tGTTTGTACAAATCGGCGAAGTCATCAAAGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCAGCCTGCAGGCCGCa                                                                                                                         <  1:314268/100‑1 (MQ=255)
                                                 tGTTTGTACAAATCGGCGAAGTCATCAAAGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCAGCCTGCAGGCCGCa                                                                                                                         >  2:314268/1‑100 (MQ=255)
                                                        aCAAATCGGCGAAGTCATCAAAGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCa                                                                                                           >  1:300676/1‑107 (MQ=255)
                                                        aCAAATCGGCGAAGTCATCAAAGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCa                                                                                                           <  2:300676/107‑1 (MQ=255)
                                                                    aGTCATCAAAGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTTTGCTTGCCAGCACgctg                                                      <  3:364094/148‑1 (MQ=255)
                                                                     gTCATCAAAGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGtt                                                                              >  4:754335/1‑123 (MQ=255)
                                                                     gTCATCAAAGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGtt                                                                              <  3:754335/123‑1 (MQ=255)
                                                                       catcaAAGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCa                                                                                                           >  2:448688/1‑92 (MQ=255)
                                                                       catcaAAGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCa                                                                                                           <  1:448688/92‑1 (MQ=255)
                                                                          caAAGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCa                                                                                                                   <  1:771229/81‑1 (MQ=255)
                                                                          caAAGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCa                                                                                                                   >  2:771229/1‑81 (MQ=255)
                                                                          caAAGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAgg                                                                                                         >  4:423586/1‑91 (MQ=255)
                                                                          caAAGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAgg                                                                                                         <  3:423586/91‑1 (MQ=255)
                                                                          caAAGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGtttt                                                                      <  3:460029/126‑1 (MQ=255)
                                                                          caAAGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGtttt                                                                      >  4:460029/1‑126 (MQ=255)
                                                                          caAAGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTTTGCTTGCCAGCACGCTGCTCACg                                                <  3:493305/148‑1 (MQ=255)
                                                                          caAAGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTTTGCTTGCCAGCACGCTGCTCACg                                                >  4:493305/1‑148 (MQ=255)
                                                                          caAAGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCAGCCTGCAGGCCGCa                                                                                                                         <  4:95047/75‑1 (MQ=255)
                                                                          caAAGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCAGCCTGCAGGCCGCa                                                                                                                         >  3:95047/1‑75 (MQ=255)
                                                                           aaaGTGGATAGCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAg                                                                                                          >  1:592831/1‑89 (MQ=255)
                                                                           aaaGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAg                                                                                                          <  2:592831/89‑1 (MQ=255)
                                                                               tGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTc                                                                                                            >  4:384115/1‑83 (MQ=255)
                                                                               tGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTc                                                                                                            <  3:384115/83‑1 (MQ=255)
                                                                                 gATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTTTGCTTGCCAGCACGCTGCTCACGGGCTGTa                                         >  1:497380/1‑148 (MQ=255)
                                                                                  aTACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTTTGCTTGCCAGCACGCTGCTCACg                                                <  4:808657/140‑1 (MQ=255)
                                                                                  aTACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTTTGCTTGCCAGCACGCTGCTCACg                                                >  3:808657/1‑140 (MQ=255)
                                                                                       cGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGc                                                                                       <  1:731392/96‑1 (MQ=255)
                                                                                       cGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGc                                                                                       <  3:955713/96‑1 (MQ=255)
                                                                                       cGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGc                                                                                       >  4:955713/1‑96 (MQ=255)
                                                                                       cGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGc                                                                                       >  2:731392/1‑96 (MQ=255)
                                                                                       cGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGttt                                                                       >  3:1061058/1‑112 (MQ=255)
                                                                                       cGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGttt                                                                       <  4:1061058/112‑1 (MQ=255)
                                                                                       cGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTttgc                                                                    >  4:421161/1‑115 (MQ=255)
                                                                                       cGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTttgc                                                                    <  3:421161/115‑1 (MQ=255)
                                                                                             ggTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTttgct                                                                   <  4:360157/110‑1 (MQ=255)
                                                                                             ggTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTttgct                                                                   >  3:360157/1‑110 (MQ=255)
                                                                                               tGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCa                                                                                                           <  4:443203/68‑1 (MQ=255)
                                                                                               tGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCa                                                                                                           >  3:443203/1‑68 (MQ=255)
                                                                                               tGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCa                                                                                                           <  4:314843/68‑1 (MQ=255)
                                                                                               tGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCa                                                                                                           >  3:314843/1‑68 (MQ=255)
                                                                                                  aTACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTc                                                                                                            <  1:129285/64‑1 (MQ=255)
                                                                                                  aTACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTc                                                                                                            >  2:129285/1‑64 (MQ=255)
                                                                                                  aTACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAgatga                                                                                               <  3:546583/77‑1 (MQ=255)
                                                                                                  aTACGTCTCTCGCGTGAAGTAATGCGGTTGTGGGGCGGCCTGGAGGCTGCACCATCACTTATTCAGGTCAGAgatga                                                                                               >  4:546583/1‑77 (MQ=255)
                                                                                                    aCGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCag                                                                                                     >  4:485071/1‑69 (MQ=255)
                                                                                                    aCGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCag                                                                                                     <  3:485071/69‑1 (MQ=255)
                                                                                                    aCGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTAt                                                                                  <  3:1131011/88‑1 (MQ=255)
                                                                                                    aCGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTAt                                                                                  >  4:1131011/1‑88 (MQ=255)
                                                                                                     cGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGaaa                                                                                          <  1:818666/79‑1 (MQ=255)
                                                                                                     cGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGaaa                                                                                          >  2:818666/1‑79 (MQ=255)
                                                                                                       tctctcGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTc                                                                             <  3:715923/90‑1 (MQ=255)
                                                                                                       tctctcGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTc                                                                             >  4:715923/1‑90 (MQ=255)
                                                                                                       tctctcGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTTTGCTTGCCAGCACgctgc                                                     >  3:590872/1‑114 (MQ=255)
                                                                                                       tctctcGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTTTGCTTGCCAGCACgctgc                                                     <  4:590872/114‑1 (MQ=255)
                                                                                                       tctctcGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTTTGCTTGCCAGCACGCTGCTc                                                   <  1:488937/116‑1 (MQ=255)
                                                                                                       tctctcGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTTTGCTTGCCAGCACGCTGCTc                                                   >  2:488937/1‑116 (MQ=255)
                                                                                                          ctcGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGc                                                                                       <  2:78010/77‑1 (MQ=255)
                                                                                                          ctcGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGc                                                                                       >  1:78010/1‑77 (MQ=255)
                                                                                                              cgTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGtttt                                                                      >  3:509961/1‑90 (MQ=255)
                                                                                                              cgTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGtttt                                                                      <  4:509961/90‑1 (MQ=255)
                                                                                                                tGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTTTGCTTGCCAGCACGCTGCTCACGGGCTGTAACACCGCTCGCGGTTTCGGCGAAGACATCaa          <  2:497380/148‑1 (MQ=255)
                                                                                                                 gAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTTTGCTTGCCAGCACGCTGCTCACGGGCTg                                           >  3:643466/1‑114 (MQ=255)
                                                                                                                 gAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTTTGCTTGCCAGCACGCTGCTCACGGGCTg                                           <  4:643466/114‑1 (MQ=255)
                                                                                                                                   tgcagCCTGCAGGCCGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTTTGCTTGCCAGCACGCTGCTCACGGGCTGTAACACCGCTCGCGGTTTCGGCGAAGACATCaa          <  1:1118324/125‑1 (MQ=255)
                                                                                                                                   tgCAGCCTGCAGGCCGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTTTGCTTGCCAGCACGCTGCTCACGGGCTGTAACACCGCTCGCGGTTTCGGCGAAGACATCaa          >  2:1118324/1‑129 (MQ=255)
                                                                                                                                    gcagCCTGCAGGCCGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGt                                                                         >  3:455007/4‑65 (MQ=255)
                                                                                                                                    gcagCCTGCAGGCCGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGt                                                                         <  4:455007/62‑1 (MQ=255)
                                                                                                                                        ccTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTTTGCTTGCCAGCACGCTGCTCACGGGCt                                            <  1:798279/90‑1 (MQ=255)
                                                                                                                                        ccTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTTTGCTTGCCAGCACGCTGCTCACGGGCt                                            >  2:798279/1‑90 (MQ=255)
                                                                                                                                                cTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTTTGCTTGCCAGCACGCTGCTCACGGGCTGTAACACCGCTCGCGGTTTCGGCGAAGACATCAAACATCTCg  <  1:649906/124‑1 (MQ=255)
                                                                                                                                                cTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTTTGCTTGCCAGCACGCTGCTCACGGGCTGTAACACCGCTCGCGGTTTCGGCGAAGACATCAAACATCTCg  >  2:649906/1‑124 (MQ=255)
                                                                                                                                                 tGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTTTGCTTGCCAGCAc                                                          <  4:551849/67‑1 (MQ=255)
                                                                                                                                                 tGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTTTGCTTGCCAGCAc                                                          >  3:551849/1‑67 (MQ=255)
                                                                                                                                                 |                                                                                                                          
GGTGGCAAACCGGCTACCCTGTCTACTGGCGCTGTGGTTAAAGTTCCGCTGTTTGTACAAATCGGCGAAGTCATCAAAGTGGATACCCGCTCTGGTGAATACGTCTCTCGCGTGAAGTAATGCGGTTGTGGTGCGGCCTGCAGGCTGCACCATCACTTATTCAGGTCAGAGATGATGAAACGCCTTATCGTTCTTGTTTTGCTTGCCAGCACGCTGCTCACGGGCTGTAACACCGCTCGCGGTTTCGGCGAAGACATCAAACATCTCG  >  NC_000913/4376146‑4376413

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: