Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,504,667 T→G 21.9% C317W (TGT→TGG)  insQ → IS609 transposase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,504,6670TG21.9% 87.3 / 3.4 47C317W (TGT→TGG) insQIS609 transposase B
Reads supporting (aligned to +/- strand):  ref base T (23/13);  new base G (3/8);  total (26/21)
Fisher's exact test for biased strand distribution p-value = 4.33e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGG‑TCAGCCAGCCGGGGCGCAATGTCCGGGCA‑AAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAG‑CGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGG  >  NC_000913/1504436‑1504901
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                          
ccACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGG‑TCAGCCAGCCGGGGCGCAATGTCCGGGCA‑AAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAG‑CGTTGCGCGTACTGTGGTCATa                                                                                                                                                                                                                                     <  1:436368/239‑1 (MQ=255)
   cGCAATGATAGTCATGGAGGATTTGAAGGGCGGCAACATGGCAAAGTCAGAAGCGGGTACGG‑TCAGCCAGCCGGGGCGCAATGTCCGGGCA‑AAATCAGGTGTAAACCGTTCGATACTGGATCAGGGCTGGTAGGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAG‑CGTTGCGCGTACTGTGGTCATACAg                                                                                                                                                                                                                                  <  2:1688504/239‑1 (MQ=255)
   cGCAATGATAGACATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGG‑TCAGCCAGCCGGGGCGCAATGTCCGGGCA‑AAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAG‑CGTTGCGCGTACTGTGGTCATACAg                                                                                                                                                                                                                                  <  2:167840/239‑1 (MQ=255)
   cGCAATGAGAGTCATTAAGGATTAGAAGGTCAGGAACATGTCAAAGTCAGCAGCGGGTACGG‑TCAGCCAGCCGGGGCGCAATGTCCGGGCA‑AAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAG‑CGTTGCGCGTACTGTGGTCATACAg                                                                                                                                                                                                                                  <  2:798640/239‑1 (MQ=255)
                  tGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGG‑TCAGCCAGCCGGGGCGCAATGTCCGGGCA‑AAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACCAGCCAG‑CGTTGCGCGTACTTTGGGCATACAGCGAAAGAGAACCGc                                                                                                                                                                                                                    >  2:879037/1‑238 (MQ=255)
                         ttGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGG‑TCAGCCAGCCGGGGCGCAATGTCCGGGCA‑AAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAG‑CGTTGCGCGTACTGTGGGCCTACAGCGAAAGAGAACCGCCTGTcaca                                                                                                                                                                                                            >  1:984881/1‑239 (MQ=255)
                               gTCAGCAACATGTCAAAGTCAGCAGCGGGTACGG‑TCAGCCAGCCGGGGCGCAATGTCCGGGCA‑AAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAG‑CGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTaa                                                                                                                                                                                                      >  1:891234/1‑239 (MQ=255)
                                       cATGTCAAAGTCAGCAGCGGGTACGG‑TCAGCCAGCCGGGGCGCAATGTCCGGGCA‑AAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGGTCCGCCAGCGTATACAAGCCAG‑CGTTGTGCGTACTGGGGTCATACAGCGAAAGAGAAACGCCTGTCACAAAAGAAATTCAGAt                                                                                                                                                                                              >  1:912330/1‑239 (MQ=255)
                                         tGTCAAAGTCAGCAGCGGGTACGG‑TCAGCCAGCCGGGGCGCAATGTCCGGGCA‑AAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGTGGCGTGCTGTTCCGCCCGCGTATACAAGCCAG‑CGTTGCGCGTACTGGGGTCATCAAGCAAAAGGGAAACGCCTGTCCCCAAGAAAAATCAAATg                                                                                                                                                                                             >  2:721603/1‑238 (MQ=255)
                                            cAAAGTCAGCAGCGGGTACGG‑TCAGCCAGCCGGGGCGCAATGTCCGGGCA‑AAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCTTGGCGGTCAGGGGCTTGCTGTTCCGCCAGCGTATACAAGCCAG‑GTTTGCGCGTACTGGGGTCATACAGCGAAAGGGAACCGCCTGTCCCAAAGTAAATTCAGATGCCAg                                                                                                                                                                                         >  1:735576/1‑239 (MQ=255)
                                             aaaGGCAGCGGCGGGTAGGG‑TCAGCCGGCCGGGGGGAAAGGGCGGGGCA‑AAAGCGGGTTTAAACGGTTGGATACGGGATCGGGGCGGGGATGAAGTGCGCCGCCAGCTGGCGTATAAGCAGCTCGGGGGGGGGGGTGGGGTGCTTGCTGTTCCGCCAGCATATACAAGCCGG‑TGTTGCGAGTATTGTGTTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAgg                                                                                                                                                                                        <  2:1562383/239‑1 (MQ=255)
                                               aGTCAGGAGGGGGTACGG‑TCAGCCCGCCGGGGCGCAGTGTCGGGGAA‑AATACGGGTTTAACCGGTGCGATTGGGGATCGGGGCGGGTAGGAAATGCGCCGCCGGTTGGCGTATAAGCAGCTGTGGCGGGGGGGTCGGGTGTTGGCTGTCCCGCCAGGGTATACAAGCCGG‑CGTGGCGCGTATTGGGGTCATACAGCAAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTa                                                                                                                                                                                      <  2:472553/239‑1 (MQ=255)
                                                   cagcgggggTACGGCACAGCGCGCGGGGGCGCATTG‑CCGGGCACAAAGCGGTTTAAAACCTTTCGAGACGGGACGAGGGGTGGTATGAAGTGGGCCGCCGGCTTGCGTATAAGAGGCTGTGGCTTGGGGGTCGGGTGTTTGCTTTTCCGCCAGCGTATACAAGCCAG‑TGTTGCGCGTACTGTGGTAATAAAGCGAAAGAAAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTAtg                                                                                                                                                                                    <  2:657711/232‑1 (MQ=255)
                                                                   acgccgcccGGGGAGCAATGGCCGGACA‑AATCCGGGTTTAACTCGTTCGATGGTGGAGTGGGGTGGTTATGAAGTCCGCCGCAGGCTTTCTTTAAAGAACCTCTGGCGTGGCGGTAGGGGGCTGGCTGTCCCGCCAGGGTAAACAAGCAAG‑CTTTGGGCGATTTGTGGTCACACAGGAAAAAAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACg                                                                                                                                                                   <  2:1075939/232‑1 (MQ=255)
                                                                     gccagccGGGGCGCAATGTCCGGGCA‑AAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAG‑CGTTGCGCGTACTGTGGGCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGAAGCCCGGTATGTGGATATACAGCGAACGcc                                                                                                                                                                 >  2:1491885/1‑239 (MQ=255)
                                                                          ccGGGGCGCAATGTCCGGGCA‑AAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAG‑CGTTGCGCGTACTGTGGTCATACAGCGTAAGAAAACCGCCTGGCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCCAACGCCGATGt                                                                                                                                                            >  1:841948/1‑239 (MQ=255)
                                                                              ggCGCAATGTCCGGGCA‑AAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAG‑CGTTGCGCGTACTGTGGTCATACAGCGGAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAAt                                                                                                                                                        >  1:864979/1‑239 (MQ=255)
                                                                                cgcAATGTCCGGGCA‑AAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCTGCGTATACAAGCCAG‑CGTTGCGCGTACTGTGGGCATACAGCGAAAGAGAACCGCCTGTTACAAAGTAAATTCAGATGCCAGGTATGTGGATCTACAGCGAACGCCGATGTAAATg                                                                                                                                                       >  1:171491/1‑238 (MQ=255)
                                                                                   aaTGTCCGGGCA‑AAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCGGTTCCGCCAGCGTATACAAGCCCG‑CGTTGCGCGGACGGTGGTCCTACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGgcg                                                                                                                                                    >  2:698552/1‑238 (MQ=255)
                                                                                         cGGGCA‑AAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGGGCTTGCTGTTCCGCCAGCGTATACAAGACAG‑CGTTGCGCGTACTGTGGTTATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCCGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCa                                                                                                                                              >  1:9325/1‑238 (MQ=255)
                                                                                                                 ttCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAG‑CGTTGCGCGTACTGTGGTCATa                                                                                                                                                                                                                                     <  1:398693/128‑1 (MQ=255)
                                                                                                                 ttCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAG‑CGTTGCGCGTACTGTGGTCATa                                                                                                                                                                                                                                     >  2:398693/1‑128 (MQ=255)
                                                                                                                 ttCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAG‑CGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGGACGcc                                                                                                                       >  1:391900/1‑238 (MQ=255)
                                                                                                                 ttCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGGGCTTGCTGTTCCGCCAGCGTATACAAGCCAG‑CGTTGCGCGTACTGTGGTCATACAGCGGAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAAAATTTTAGCGGCGGGGCACGcc                                                                                                                       >  2:845539/1‑238 (MQ=255)
                                                                                                                         tGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAG‑CGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGGACCCCgggccttc                                                                                                               >  1:167746/1‑231 (MQ=255)
                                                                                                                           gATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAG‑CGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGGGGATATACAGCGAACCCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCGCGCCGCTCTTGc                                                                                                               >  2:780420/1‑236 (MQ=255)
                                                                                                                                gggCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGCGCTTGCTGTTCCGCCAACGTATACCAACCAG‑CGTTGCGCGTACTGTGGTCATACAGCGCAAGAGAACCGCCTGTCCCAAAGGAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGGTGTAAAAGGCGCTCCGAACATTTTAGCGGCCGGGCACGCCGCTCttcgcctcggag                                                                                                       >  2:995808/1‑229 (MQ=255)
                                                                                                                                        agaaaaTGCGCCGCCGGTTGGGGTATAAGCAGCTCGGGGGGGGGGGTCGGGTGGTGGCTGTCCCGCAGGCGAAAAAAAGCCAG‑GGTTGCGCGGACGGGGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTg                                                                                               <  2:107513/236‑1 (MQ=255)
                                                                                                                                                cgccgcCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAG‑CGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGGACGCCGCTCTTTGCTCTGGTGAGGCGGgtcagtcaa                                                                                        >  2:166238/1‑237 (MQ=255)
                                                                                                                                                     ccAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAG‑CGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGGTCTTGCCCGTGGAGAGGTGGGGGAGACAAgccgc                                                                                   >  1:1722442/1‑238 (MQ=255)
                                                                                                                                                                  aaGCAGCTCTGGCGTGGCGGTCGGGTGCTTGCTGTTCCGCCAGCTTATAAAAGCCAG‑GGTTGCGCGTACTGGGGTCATCCAGCAAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGa                                                                     <  1:1201933/239‑1 (MQ=255)
                                                                                                                                                                   agcggcTCGGGCGGGGCGTTAGGGTTCTGGCTGTCCCGCAGGCGTATAAAAGCCGG‑CGTGGGCCGTACGGGGGTCAAACGGCGAAAAAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGaa                                                                    <  1:1336824/239‑1 (MQ=255)
                                                                                                                                                                             gcgtggcgGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAG‑CGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATg                                                                                                                                                             >  1:854833/1‑140 (MQ=255)
                                                                                                                                                                                  gcgGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAG‑CGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGGCCGCCGGTCTTTTCTGGGGGGGGATGGTGTGGTCAGGGCGCCCGCTGGAGCAGGGAcccccccccaaagtt                                                     >  1:1108176/1‑230 (MQ=255)
                                                                                                                                                                                          gtgtTTGCTGTCCCGCCAGCGTATACAAGCCAG‑AGTTGCGCGTACTGGGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGa                                              <  1:857263/238‑1 (MQ=255)
                                                                                                                                                                                          gTGCTGGCGGTCCCGCCGGCGTAAACAAGCCGG‑CGTGGCGCTTACTGTGGCCATACAGCAAAAAAAAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGa                                              <  1:962587/238‑1 (MQ=255)
                                                                                                                                                                                           tgcttg‑tgTTCCCCCAGTGTTTACAAGCCCGATGTTGCGGGACCTGGGGCCATAGACAGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCg                                               <  1:177260/236‑1 (MQ=255)
                                                                                                                                                                                               tcgttttCCGCCAGCGTATACAAGCAAG‑CGTTGCGCGTATTGTGGTCAAACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCt                                        <  1:1717341/234‑1 (MQ=255)
                                                                                                                                                                                                  tGTCCCGCCGGCGTATAAAAGCCCG‑GTTTGCGGGTATTGTGGTAATACAGCAAAAGAAAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGaa                                     <  1:1718510/239‑1 (MQ=255)
                                                                                                                                                                                                  tGTCCCGCCAGCGTATCCAAGCGAG‑GGTTGCGCGTATTGGGGTAATACAGCAAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGaa                                     <  1:1111438/239‑1 (MQ=255)
                                                                                                                                                                                                                aaaCAAGCCGG‑GGTTGCGCGTCCTGTGTACAAACGGGGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGACCCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATcc                       <  2:1720600/237‑1 (MQ=255)
                                                                                                                                                                                                                  aCAAGCCAG‑CGTTGCGCGTACTGTGGTCATAAAGCAAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCAc                     <  1:879037/239‑1 (MQ=255)
                                                                                                                                                                                                                    aaGCCAG‑CGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTca                                                                                                                                                                                                              >  1:918180/1‑52 (MQ=255)
                                                                                                                                                                                                                       ccAG‑CGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGGTCTTTCCTTTGGAGAGATGGGGCAGGCACGCCGCCCGCTGGAGCAGGAACCCCCCCAAATAATTCAAGCGGCCGCCTGAAAGCAGCAGGGCTCCCCGTCCt                >  1:674935/1‑239 (MQ=255)
                                                                                                                                                                                                                       ccAG‑CGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCCGTCTTGCCTGTGGGGAGATGGGGCCGTCAAGCCGCCCGCTTAAGCAGGAACCCCCCCAAATGATTTCGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCt                >  1:173481/1‑239 (MQ=255)
                                                                                                                                                                                                                            cGTTGCGCTTCCTGGGGACTTACAGGAAAGGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCa             <  1:721603/238‑1 (MQ=255)
                                                                                                                                                                                                                                  gcgaattgtgGTCATACAGAGAAAGAGACCCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTg      <  2:167746/233‑1 (MQ=255)
                                                                                                                                                                                                                                      aCTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGGTCTTGCCTGTGGAGAGATGGGGCAGGCAGGCCCCCCGTTGCAGCAGGAACCCACCGAAATGATTCAGGGGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTggggg  >  1:1498591/1‑239 (MQ=255)
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                          
CCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGG‑TCAGCCAGCCGGGGCGCAATGTCCGGGCA‑AAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAG‑CGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGG  >  NC_000913/1504436‑1504901

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑