Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,504,504 G→A 5.7% S263N (AGC→AAC)  insQ → IS609 transposase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,504,5040GA 5.7% 222.3 / 2.3 70S263N (AGC→AAC) insQIS609 transposase B
Reads supporting (aligned to +/- strand):  ref base G (36/30);  new base A (2/2);  total (38/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 = 9.73e-01

TACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAG  >  NC_000913/1504360‑1504652
                                                                                                                                                |                                                                                                                                                    
tACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTcagccag                                                                                                                                                  <  1:1239909/149‑1 (MQ=255)
 aCAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTcagccagc                                                                                                                                                 <  2:1023835/149‑1 (MQ=255)
            aTTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAAt                                                                                                                                      <  2:398323/149‑1 (MQ=255)
            aTTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAAt                                                                                                                                      <  2:1597680/149‑1 (MQ=255)
                  gTATCGCAAATATCCGCAGGGACGACCTTCACAAAGTCACACCGGCCGTCAGCAAAAACCACGCAATGATACTCATTGAGGAGTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCATCCGGGGCGCAATGTCCgg                                                                                                                                <  1:138096/149‑1 (MQ=255)
                        cAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCATTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCCGCAGCTGTTACGGTCAGCCAGCCGGGGCGCAATTTCCGGGCaaaa                                                                                                                          >  1:1056393/1‑149 (MQ=255)
                               ccGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCCGCAAAAACCACGCAATGATAGGCCTTGAGGATTTGAAGGGCAGCATCATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGtt                                                                                                                   <  1:404056/149‑1 (MQ=255)
                                        cTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAACCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGtt                                                                                                          >  1:1359893/1‑149 (MQ=255)
                                         tACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTc                                                                                                         <  1:828792/149‑1 (MQ=255)
                                            cTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGTTTTAAACCGTTCGAt                                                                                                      >  1:185573/1‑149 (MQ=255)
                                            cTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACcgatcgat                                                                                                      >  1:991784/1‑149 (MQ=255)
                                            cTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGAt                                                                                                      >  1:730883/1‑149 (MQ=255)
                                            cTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGAt                                                                                                      >  1:1195969/1‑149 (MQ=255)
                                             ttCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCCGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAGCCGTTCGATa                                                                                                     >  2:317342/1‑149 (MQ=255)
                                             ttCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATa                                                                                                     >  2:475622/1‑149 (MQ=255)
                                               cacaAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACt                                                                                                   >  2:813731/1‑149 (MQ=255)
                                                acaAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTg                                                                                                  <  2:101314/149‑1 (MQ=255)
                                                 caAAGTCACAACGGCCGTCAGCAAAAACCACGCATTGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTgg                                                                                                 <  2:721189/149‑1 (MQ=255)
                                                 caAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTgg                                                                                                 >  1:1594084/1‑149 (MQ=255)
                                                  aaaGTCACAACGGCCATCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGa                                                                                                >  2:793805/1‑149 (MQ=255)
                                                   aaGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGAt                                                                                               <  1:750698/149‑1 (MQ=255)
                                                        acaACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAggg                                                                                          <  2:1400172/149‑1 (MQ=255)
                                                         ccaCGGCCGTCCGCAAAAACCACGAAATGAGAGTCATGGAGGATGTGAAGGGCAGCAACATGTCAAAGTCAGCAGAGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGGTTAAACCGTTCGATACTGGATCAGGGc                                                                                         <  1:100802/147‑1 (MQ=255)
                                                          aaCGGCCGTCAGCATACCCCACGGAATGATAGTCATTGAGGATTTGAAGGGCAGGAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCt                                                                                        <  1:156262/149‑1 (MQ=255)
                                                          aaCGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCt                                                                                        <  1:392053/149‑1 (MQ=255)
                                                          aaCGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCt                                                                                        <  2:185573/149‑1 (MQ=255)
                                                               ccGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTAt                                                                                   >  1:513703/1‑149 (MQ=255)
                                                                cGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATg                                                                                  <  1:488363/149‑1 (MQ=255)
                                                                cGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATg                                                                                  <  1:596246/149‑1 (MQ=255)
                                                                cGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATg                                                                                  <  1:625591/149‑1 (MQ=255)
                                                                  tCAGCAAAACCCACGCAATGATAGTCATTGAGGAGTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGGCCGGTCAACCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGaa                                                                                <  2:1581819/149‑1 (MQ=255)
                                                                  tCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGaa                                                                                >  1:324574/1‑149 (MQ=255)
                                                                   cAGCAAAAACCTCGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGaaa                                                                               <  2:1122469/149‑1 (MQ=255)
                                                                   cAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAACCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGaaa                                                                               >  2:1445458/1‑149 (MQ=255)
                                                                    aGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAAt                                                                              >  2:579296/1‑149 (MQ=255)
                                                                      cAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGt                                                                                     <  2:232126/140‑1 (MQ=255)
                                                                      cAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGt                                                                                     >  1:232126/1‑140 (MQ=255)
                                                                       aaaaaCCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATgcg                                                                           >  1:1500278/1‑149 (MQ=255)
                                                                           aCCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGcgccgc                                                                       <  2:513703/149‑1 (MQ=255)
                                                                                    tGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTGCGGGCAAAATCAGGGTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCg                                                              <  2:257081/149‑1 (MQ=255)
                                                                                      aTAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCATAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTCCGta                                                            >  1:1174610/1‑149 (MQ=255)
                                                                                        aGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGtata                                                          >  1:192221/1‑149 (MQ=255)
                                                                                          tCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCATTGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAag                                                        >  2:390942/1‑149 (MQ=255)
                                                                                          tCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAag                                                        <  1:1100695/149‑1 (MQ=255)
                                                                                          tCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAag                                                        <  1:491999/149‑1 (MQ=255)
                                                                                                aGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGctc                                                  >  1:749290/1‑149 (MQ=255)
                                                                                                 ggATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGctct                                                 >  2:391145/1‑149 (MQ=255)
                                                                                                    tttGAAGGTCATCAACATGTCAAAGTCCGCAGCGGGTACGGTCAGCAAGACGGGGCGCAATGTCCGGGCAAAATCAGCTTTAAACCGTTCGATCCTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCaggc                                              >  1:1135261/1‑149 (MQ=255)
                                                                                                       taatgtgaGCAACATGTCAAAGTCAGCAGCGGGTACGTTCAGCCAGACGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGTGCTGGTATTAAATGCGCCGCCAGCTTGCGTTTAAGCAGCTCtggcgtg                                           <  2:479739/142‑1 (MQ=255)
                                                                                                        aaGGTCAGCAACATGTCAAAGTCCGCAGCGGGTACGGTCAGCCAGCCGGGGAGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTAGGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCtggcgtgg                                          <  2:1065618/149‑1 (MQ=255)
                                                                                                         aGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTAAGCCAGCCGGGTCGAAATGTCCTGGCAAAATCAGTTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGATCtggcgtggc                                         <  2:1056393/149‑1 (MQ=255)
                                                                                                            tCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGt                                      >  2:249335/1‑149 (MQ=255)
                                                                                                            tCAGCAACATGTCAAAGTCAGCAGCGGGAACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGt                                      <  2:1344048/149‑1 (MQ=255)
                                                                                                             cAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTc                                     >  2:1306466/1‑149 (MQ=255)
                                                                                                             cAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTc                                     >  1:1133207/1‑149 (MQ=255)
                                                                                                                 aaCATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTACACCGTTCGCTACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGt                                 >  2:1156938/1‑149 (MQ=255)
                                                                                                                     tGTCAATGTAAGCAGCGTGTACGTTCAGCCAGCCGGGGCGAAATGTCCGGGAAAAATCTGGTGTAAACCGTTCGATACTGGATCAGGGGGGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGtgctt                             <  2:1019237/149‑1 (MQ=255)
                                                                                                                     tGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGtgctt                             >  1:28850/1‑149 (MQ=255)
                                                                                                                      ggcAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGtgcttg                            <  1:249335/147‑1 (MQ=255)
                                                                                                                      gTCAAACTCAGCAGCGGGTACGGGCAGCCAGCCGGGGCGCAATGTCCGGGCAAACTCAGGTTTAAACCGTGCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGtgcttg                            <  2:1007818/149‑1 (MQ=255)
                                                                                                                         aaaGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCAGGTATGATATGCGCCGCCAGCTTGCGTATAAGCCGCTCTGGCGTGGCGGTCAGGTGCTTGCTg                         >  2:1060801/1‑149 (MQ=255)
                                                                                                                                 cagcGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGccc                 >  2:203276/1‑148 (MQ=255)
                                                                                                                                        tACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGtata          >  1:533052/1‑149 (MQ=255)
                                                                                                                                           ggTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACaa       >  2:289659/1‑149 (MQ=255)
                                                                                                                                           ggTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACaa       >  1:473322/1‑149 (MQ=255)
                                                                                                                                           ggTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACaa       >  1:1307882/1‑149 (MQ=255)
                                                                                                                                            gTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTTGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACaac      >  1:690621/1‑148 (MQ=255)
                                                                                                                                            gTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACaac      >  1:1383768/1‑148 (MQ=255)
                                                                                                                                            gTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTAAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAg      >  2:1135246/1‑149 (MQ=255)
                                                                                                                                             tCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGc     <  2:324574/149‑1 (MQ=255)
                                                                                                                                             tCAGCCAGCCGGGGCGCAATGTCCCGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGACAGGTGCTTGCTGTTCCGCCAGCGTATACAAGc     >  1:727037/1‑149 (MQ=255)
                                                                                                                                             tCAACCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGc     <  2:1359893/149‑1 (MQ=255)
                                                                                                                                                gccagccGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAg  >  2:940687/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
TACAGCGACTGCATTCCTGTATCGCAAATATCCGCAGGGACTACCTTCACAAAGTCACAACGGCCGTCAGCAAAAACCACGCAATGATAGTCATTGAGGATTTGAAGGTCAGCAACATGTCAAAGTCAGCAGCGGGTACGGTCAGCCAGCCGGGGCGCAATGTCCGGGCAAAATCAGGTTTAAACCGTTCGATACTGGATCAGGGCTGGTATGAAATGCGCCGCCAGCTTGCGTATAAGCAGCTCTGGCGTGGCGGTCAGGTGCTTGCTGTTCCGCCAGCGTATACAAGCCAG  >  NC_000913/1504360‑1504652

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