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→TGGinsQIS609 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

CCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGG  >  NC_000913/1504436‑1504901
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                          
ccACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAGCGTTGCGCGTACTGTGGTCATa                                                                                                                                                                                                                                     <  1:436368/239‑1 (MQ=255)
   cGCAATGATAGTCATGGAGGATTTGAAGGGCGGCAACATGGCAAAGTCAGAAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTGTAAACCGTTCGATACTGGATCAGGGCTGGTAGGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAGCGTTGCGCGTACTGTGGTCATACAg                                                                                                                                                                                                                                  <  2:1688504/239‑1 (MQ=255)
   cGCAATGATAGACATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAGCGTTGCGCGTACTGTGGTCATACAg                                                                                                                                                                                                                                  <  2:167840/239‑1 (MQ=255)
   cGCAATGAGAGTCATTAAGGATTAGAAGGTCAGGAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAGCGTTGCGCGTACTGTGGTCATACAg                                                                                                                                                                                                                                  <  2:798640/239‑1 (MQ=255)
                  tGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACCAGCCAGCGTTGCGCGTACTTTGGGCATACAGCGAAAGAGAACCGc                                                                                                                                                                                                                    >  2:879037/1‑238 (MQ=255)
                         ttGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAGCGTTGCGCGTACTGTGGGCCTACAGCGAAAGAGAACCGCCTGTcaca                                                                                                                                                                                                            >  1:984881/1‑239 (MQ=255)
                               gTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTaa                                                                                                                                                                                                      >  1:891234/1‑239 (MQ=255)
                                       cATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGGTCCGCCAGCGTATACAAGCCAGCGTTGTGCGTACTGGGGTCATACAGCGAAAGAGAAACGCCTGTCACAAAAGAAATTCAGAt                                                                                                                                                                                              >  1:912330/1‑239 (MQ=255)
                                         tGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGTGGCGTGCTGTTCCGCCCGCGTATACAAGCCAGCGTTGCGCGTACTGGGGTCATCAAGCAAAAGGGAAACGCCTGTCCCCAAGAAAAATCAAATg                                                                                                                                                                                             >  2:721603/1‑238 (MQ=255)
                                            cAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCTTGGCGGTCAGGGGCTTGCTGTTCCGCCAGCGTATACAAGCCAGGTTTGCGCGTACTGGGGTCATACAGCGAAAGGGAACCGCCTGTCCCAAAGTAAATTCAGATGCCAg                                                                                                                                                                                         >  1:735576/1‑239 (MQ=255)
                                             aaaGGCAGCGGCGGGTAGGGTCAGCCGGCCGGGGGGAAAGGGCGGGGCAAAAGCGGGTTTAAACGGTTGGATACGGGATCGGGGCGGGGATGAAGTGCGCCGCCAGCTGGCGTATAAGCAGCTCGGGGGGGGGGGTGGGGTGCTTGCTGTTCCGCCAGCATATACAAGCCGGTGTTGCGAGTATTGTGTTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAgg                                                                                                                                                                                        <  2:1562383/239‑1 (MQ=255)
                                               aGTCAGGAGGGGGTACGGTCAGCCCGCCGGGGCGCAGTGTCGGGGAAAATACGGGTTTAACCGGTGCGATTGGGGATCGGGGCGGGTAGGAAATGCGCCGCCGGTTGGCGTATAAGCAGCTGTGGCGGGGGGGTCGGGTGTTGGCTGTCCCGCCAGGGTATACAAGCCGGCGTGGCGCGTATTGGGGTCATACAGCAAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTa                                                                                                                                                                                      <  2:472553/239‑1 (MQ=255)
                                                   cagcgggggTACGGCACAGCGCGCGGGGGCGCATTGCCGGGCACAAAGCGGTTTAAAACCTTTCGAGACGGGACGAGGGGTGGTATGAAGTGGGCCGCCGGCTTGCGTATAAGAGGCTGTGGCTTGGGGGTCGGGTGTTTGCTTTTCCGCCAGCGTATACAAGCCAGTGTTGCGCGTACTGTGGTAATAAAGCGAAAGAAAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTAtg                                                                                                                                                                                    <  2:657711/232‑1 (MQ=255)
                                                                   acgccgcccGGGGAGCAATGGCCGGACAAATCCGGGTTTAACTCGTTCGATGGTGGAGTGGGGTGGTTATGAAGTCCGCCGCAGGCTTTCTTTAAAGAACCTCTGGCGTGGCGGTAGGGGGCTGGCTGTCCCGCCAGGGTAAACAAGCAAGCTTTGGGCGATTTGTGGTCACACAGGAAAAAAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACg                                                                                                                                                                   <  2:1075939/232‑1 (MQ=255)
                                                                     gccagccGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAGCGTTGCGCGTACTGTGGGCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGAAGCCCGGTATGTGGATATACAGCGAACGcc                                                                                                                                                                 >  2:1491885/1‑239 (MQ=255)
                                                                          ccGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAGCGTTGCGCGTACTGTGGTCATACAGCGTAAGAAAACCGCCTGGCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCCAACGCCGATGt                                                                                                                                                            >  1:841948/1‑239 (MQ=255)
                                                                              ggCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAGCGTTGCGCGTACTGTGGTCATACAGCGGAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAAt                                                                                                                                                        >  1:864979/1‑239 (MQ=255)
                                                                                cgcAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCTGCGTATACAAGCCAGCGTTGCGCGTACTGTGGGCATACAGCGAAAGAGAACCGCCTGTTACAAAGTAAATTCAGATGCCAGGTATGTGGATCTACAGCGAACGCCGATGTAAATg                                                                                                                                                       >  1:171491/1‑238 (MQ=255)
                                                                                   aaTGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCGGTTCCGCCAGCGTATACAAGCCCGCGTTGCGCGGACGGTGGTCCTACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGgcg                                                                                                                                                    >  2:698552/1‑238 (MQ=255)
                                                                                         cGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGGGCTTGCTGTTCCGCCAGCGTATACAAGACAGCGTTGCGCGTACTGTGGTTATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCCGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCa                                                                                                                                              >  1:9325/1‑238 (MQ=255)
                                                                                                                 ttCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAGCGTTGCGCGTACTGTGGTCATa                                                                                                                                                                                                                                     <  1:398693/128‑1 (MQ=255)
                                                                                                                 ttCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAGCGTTGCGCGTACTGTGGTCATa                                                                                                                                                                                                                                     >  2:398693/1‑128 (MQ=255)
                                                                                                                 ttCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGGACGcc                                                                                                                       >  1:391900/1‑238 (MQ=255)
                                                                                                                 ttCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGGGCTTGCTGTTCCGCCAGCGTATACAAGCCAGCGTTGCGCGTACTGTGGTCATACAGCGGAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAAAATTTTAGCGGCGGGGCACGcc                                                                                                                       >  2:845539/1‑238 (MQ=255)
                                                                                                                         tGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGGACCCCgggccttc                                                                                                               >  1:167746/1‑231 (MQ=255)
                                                                                                                           gATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGGGGATATACAGCGAACCCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCGCGCCGCTCTTGc                                                                                                               >  2:780420/1‑236 (MQ=255)
                                                                                                                                gggCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGCGCTTGCTGTTCCGCCAACGTATACCAACCAGCGTTGCGCGTACTGTGGTCATACAGCGCAAGAGAACCGCCTGTCCCAAAGGAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGGTGTAAAAGGCGCTCCGAACATTTTAGCGGCCGGGCACGCCGCTCttcgcctcggag                                                                                                       >  2:995808/1‑229 (MQ=255)
                                                                                                                                        agaaaaTGCGCCGCCGGTTGGGGTATAAGCAGCTCGGGGGGGGGGGTCGGGTGGTGGCTGTCCCGCAGGCGAAAAAAAGCCAGGGTTGCGCGGACGGGGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTg                                                                                               <  2:107513/236‑1 (MQ=255)
                                                                                                                                                cgccgcCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGGACGCCGCTCTTTGCTCTGGTGAGGCGGgtcagtcaa                                                                                        >  2:166238/1‑237 (MQ=255)
                                                                                                                                                     ccAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGGTCTTGCCCGTGGAGAGGTGGGGGAGACAAgccgc                                                                                   >  1:1722442/1‑238 (MQ=255)
                                                                                                                                                                  aaGCAGCTCTGGCGTGGCGGTCGGGTGCTTGCTGTTCCGCCAGCTTATAAAAGCCAGGGTTGCGCGTACTGGGGTCATCCAGCAAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGa                                                                     <  1:1201933/239‑1 (MQ=255)
                                                                                                                                                                   agcggcTCGGGCGGGGCGTTAGGGTTCTGGCTGTCCCGCAGGCGTATAAAAGCCGGCGTGGGCCGTACGGGGGTCAAACGGCGAAAAAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGaa                                                                    <  1:1336824/239‑1 (MQ=255)
                                                                                                                                                                             gcgtggcgGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATg                                                                                                                                                             >  1:854833/1‑140 (MQ=255)
                                                                                                                                                                                  gcgGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGGCCGCCGGTCTTTTCTGGGGGGGGATGGTGTGGTCAGGGCGCCCGCTGGAGCAGGGAcccccccccaaagtt                                                     >  1:1108176/1‑230 (MQ=255)
                                                                                                                                                                                          gtgtTTGCTGTCCCGCCAGCGTATACAAGCCAGAGTTGCGCGTACTGGGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGa                                              <  1:857263/238‑1 (MQ=255)
                                                                                                                                                                                          gTGCTGGCGGTCCCGCCGGCGTAAACAAGCCGGCGTGGCGCTTACTGTGGCCATACAGCAAAAAAAAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGa                                              <  1:962587/238‑1 (MQ=255)
                                                                                                                                                                                           tgcttgtgTTCCCCCAGTGTTTACAAGCCCGATGTTGCGGGACCTGGGGCCATAGACAGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCg                                               <  1:177260/236‑1 (MQ=255)
                                                                                                                                                                                               tcgttttCCGCCAGCGTATACAAGCAAGCGTTGCGCGTATTGTGGTCAAACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCt                                        <  1:1717341/234‑1 (MQ=255)
                                                                                                                                                                                                  tGTCCCGCCGGCGTATAAAAGCCCGGTTTGCGGGTATTGTGGTAATACAGCAAAAGAAAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGaa                                     <  1:1718510/239‑1 (MQ=255)
                                                                                                                                                                                                  tGTCCCGCCAGCGTATCCAAGCGAGGGTTGCGCGTATTGGGGTAATACAGCAAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGaa                                     <  1:1111438/239‑1 (MQ=255)
                                                                                                                                                                                                                aaaCAAGCCGGGGTTGCGCGTCCTGTGTACAAACGGGGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGACCCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATcc                       <  2:1720600/237‑1 (MQ=255)
                                                                                                                                                                                                                  aCAAGCCAGCGTTGCGCGTACTGTGGTCATAAAGCAAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCAc                     <  1:879037/239‑1 (MQ=255)
                                                                                                                                                                                                                    aaGCCAGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTca                                                                                                                                                                                                              >  1:918180/1‑52 (MQ=255)
                                                                                                                                                                                                                       ccAGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGGTCTTTCCTTTGGAGAGATGGGGCAGGCACGCCGCCCGCTGGAGCAGGAACCCCCCCAAATAATTCAAGCGGCCGCCTGAAAGCAGCAGGGCTCCCCGTCCt                >  1:674935/1‑239 (MQ=255)
                                                                                                                                                                                                                       ccAGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCCGTCTTGCCTGTGGGGAGATGGGGCCGTCAAGCCGCCCGCTTAAGCAGGAACCCCCCCAAATGATTTCGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCt                >  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)
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                          
CCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAGCGTTGCGCGTACTGTGGTCATACAGCGAAAGAGAACCGCCTGTCACAAAGTAAATTCAGATGCCAGGTATGTGGATATACAGCGAACGCCGATGTAAATGGCGCTCGCAACATTTTAGCGGCGGGGCACGCCGTTCTTGCCTGTGGAGAGATGGTGCAGTCAGGCCGCCCGTTGAAGCAGGAACCCACCGAAATGATTCAGGCGACAGCCTGAACGTAGCAGGGATCCACGTCCTTCAGGGCGTGGAGG  >  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: