Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,187,658 T→G L770R (CTC→CGC)  rpoC → RNA polymerase, beta prime subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,187,6580TG95.1% 117.1 / ‑2.5 41L770R (CTC→CGC) rpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base G (23/16);  minor base A (1/1);  total (24/17)
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.14e-01

TGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGACTCGTCGTCTGGTTGACGTGGCGCAGGACCTGGTGGTTACCGAAGACGATTGTGGTACCCATGAAGGTATCATGATGACTCCGGTTATCGAGGGTGGTGACGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGAAGCCGGGTACTGCTGATATCCTCGTTCCGCGCAACACGCTGCTGCACGAACAGTGGTG  >  NC_000913/4187364‑4187955
                                                                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                                                                         
tggtAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTTGTCAGCTTGCTGGTATTCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACAGCGAACTTCCGTGAAGGTCTGAACGTACGCCAGTa                                                                                                                                                                                                                                                                                                     >  2:75277/1‑301 (MQ=255)
         cGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACGCCAGTACTTCAtct                                                                                                                                                                                                                                                                                             >  2:171869/1‑300 (MQ=255)
                         caaAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTTTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGCTGGCACCATCATCGAAACGCCACTCACCCCGAACTTCCGTGAAGTTCTGAACGTACGCCCGCCCTTCCTCTacactctcgttcttcat                                                                                                                                                                                                                                                                            >  1:264302/1‑288 (MQ=255)
                            aagacgtgtgctcttccgaactgcgAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACGCCAGTACTTCATCTCCACCCACGGTGCTCGTaaa                                                                                                                                                                                                                                                                         <  2:197348/279‑1 (MQ=255)
                              ttATCGATATCTGGGCTGCGGCGAACGAGCGTGTATCCAAAGCGATGATGGATAACCTGCAAGCTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCATGTTTCCTTCAACAGCATCTACATGATGGCCAACTCCGGTGCGCGTGTTTCTGCGGCACAGATTCGTCAGCTTGCTAGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACGCCAGTACTTCATCTCCACCCACGGTGCTCGTaaa                                                                                                                                                                                                                                                                         <  2:254997/299‑1 (MQ=255)
                                        cTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGTAGCCGGCTGGCTCCATCATCGAACCGCCAATTACCGCGAACTTCCGTGAAGGTCTGAACGTACGCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGAt                                                                                                                                                                                                                                                             >  2:202306/1‑301 (MQ=255)
                                         tGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAATCGTGACGGTCAGGAAGAGAAGTAGGTTTCCTTCAATAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGACCCATCATCGCAACGCCAACCACCGCCACCCCCCGTGACGGTCTGAACGTACGCCAGGACCTCATCCCTACCTACGGTGCACGCAAAGGCTTGACGGATa                                                                                                                                                                                                                                                            >  1:70640/1‑301 (MQ=255)
                                          gggCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACGCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATa                                                                                                                                                                                                                                                            <  2:85477/300‑1 (MQ=255)
                                                  gcgAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGACGGTCTGAACGTACGCCAGTACTTCATCTCCACCCACGTTGCTCGTaaaagatcggaagagcggattgtcg                                                                                                                                                                                                                                                   >  1:197342/1‑279 (MQ=255)
                                                      aCGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACGCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGGCTGGCGGATACCGCACTGAAAAc                                                                                                                                                                                                                                               <  2:62280/301‑1 (MQ=255)
                                                        gATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGACGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACGCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACtt                                                                                                                                                                                                                                             >  2:33580/1‑300 (MQ=255)
                                                         aTCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACGCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACtt                                                                                                                                                                                                                                             >  2:173367/1‑299 (MQ=255)
                                                              gtATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACACCAGTACTTAATCTCAACCCACGTAGCCGGTCAAGGTCTGGCGGATACCGCACTGAAAACTGcaacc                                                                                                                                                                                                                                        >  1:41702/1‑296 (MQ=255)
                                                                                                  aCCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTACAACAGCATCCACATGATGGCCGACTCCGGTGCGCGTGGTCCTGCGGCCCAGATTCGTCAGCTTGCTGGTCTGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACGCCCGTGAAGGTCTGAACGTACGCCAGTACTTCATCTCCACCCTCGGTGCTCGTAAAGGTCTGGCGGTGACCGCTCTGAAAACTGCGAACTCCGGTTCCCTgcacattcgccctgtttaagtggcgt                                                                                                                                                                                                   >  1:239954/1‑276 (MQ=255)
                                                                                                                         ggAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACGCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGACTCGTCGTCTGTTTGACGTGGCGCAGGACCTGGTGGTTACCGAAGag                                                                                                                                                                            >  2:150560/1‑300 (MQ=255)
                                                                                                                            cccgtcagcactcatttcctatcacagcTACATGATGGCCGAGTCCGGTGCGCGAGGGTATGCGGCACAGATTCGTCCACTTGCTGGTATGTGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACGCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGACTCGTCGTCTGGTTGACGTGGCGCAGGACCTGGTGGTTACCGAAGACGAt                                                                                                                                                                         <  1:211021/274‑1 (MQ=255)
                                                                                                                                gaGCAGGTTGCCTTCAACAGCATCTACATGATGGCCGGCTCCGGTGCGCGTGGTTCTGCGGGACAGATTCGTCAGCGTGCTGGTATCCGTGGTCTGCTGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACGTCCGTGAAGGTCTGAACGTACGCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGACTCGTCGTCTGGTTGACGTGGCGCAGGACCTGGTGGTTACCGAAGACGATtg                                                                                                                                                                       <  1:22791/298‑1 (MQ=255)
                                                                                                                                       tttCCTTCAACAGCATCTCCATGATGGCCGGCTCCGGTGCGCTTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACACCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGACTCGTCGTCTGGTTGACGTGGCACAGGACCTGGTGGTTACCGAAGACGATTGTGGTACCCa                                                                                                                                                              <  2:41703/301‑1 (MQ=255)
                                                                                                                                                    cATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACGCAAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGACTCGTCGTCTGGTTGACGTGGCGCAGGCCCTGGTGGTTACCGAAGACGacagatcggaagagaacacgcctgaa                                                                                                                                                 >  2:211027/1‑276 (MQ=255)
                                                                                                                                                      tCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACGCCAGTACTTCATCTCCACCCACGGTGCTCGTGAAGGTCTGGCGGATACCGCACTGACAACTGCGAACTCCGGTTACCTGACTCGTCGTCTGGTTGACGTTGCGCCGGACCCGGTGGTTACCGAAGACGATTGTGGTAACCATTAAGTGAACatgat                                                                                                                                               >  1:125071/1‑301 (MQ=255)
                                                                                                                                                               tgGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACGCCAGTACTTCATCTCCACCCACGGTGCGCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGACTCGTCGTCTGGTTGACGTGGCGCAGGACCGGTGGGTTACGGAAGACGATTGTGGTACTCATGAAGGTATAAGGATGAACCCggg                                                                                                                                      >  1:66199/1‑300 (MQ=255)
                                                                                                                                                                          ggTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACGCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGACTCGTCGTCTGGTTGACGTGGCGCAGGACCTGGTGGTTACCGAAGACGATTGTGGTACCCATGAAGGTATCATGATGACTCCGGTTATCGAGggtg                                                                                                                           >  2:95717/1‑301 (MQ=255)
                                                                                                                                                                           tgtggggTGGGTCTGCGGCACAGAGTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGGACGCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGACTCGTCGTCTGGTTGACGTGGCGCAGGACCTGGTGGTTACCGAAGACGATTGTGGTACCCATGAAGGTATCATGATGACTCCGGTTATCGAGggtgg                                                                                                                          <  1:179621/295‑1 (MQ=255)
                                                                                                                                                                            tGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACGCCAGTACTTCATCTACACCCACGGTGCTCGTAAAGGTTTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGCCTCGACGCCTGGTTGACGTGGCGCATGACCGGGTGGTTCACGAAGACGATTGTGGTACCCCAGATGTAATCATGATGATCCCAGTATTCGCGggtggt                                                                                                                         >  1:183385/1‑301 (MQ=255)
                                                                                                                                                                                    gTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACGCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGACTCGTCGTCTGGTTGACGTGGCGCAGGACCTGGTGGTTACCGAAGACGATTGTGGTACCCATGAAGGTATCATGATGACTCCGGTTATCGAGGGTGGTGACGTTa                                                                                                                  <  2:125074/300‑1 (MQ=255)
                                                                                                                                                                                        tGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCGCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGGACGCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGAGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGACTCGTCGTCTGGTTGACGTGGCGCAGGACCTGGTGGTTACCGAAGACGATTGTGGTACCCATGAAGGTATCATGATGACTCCGGTTATCGAGGGTGGTGACGTTAAagag                                                                                                             <  2:264310/301‑1 (MQ=255)
                                                                                                                                                                                                          gCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACGCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGACTCGTCGTCTGGTTGACGTGGCGCAGGACCTGGTGGTTACCGAAGACGATTGTGGTACCCATGAAGGTATCATGATGACTCCGGTTATCGAGGGTGGTGACGTTAAAGAGCCGCTGCGCGATCGCGTa                                                                                           >  2:166530/1‑301 (MQ=255)
                                                                                                                                                                                                           cTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACGCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGACTCGTCGTCTGGTTGACGTGGCGCAGGACCTGGTGGTTACCGAAGACGATTGTGGTACCCATGCAGGTATCATGATGACTCCGGTTATCGANGGTGGGGACGTTAAAGAGCCGCTGCGCGATCGCGTa                                                                                           >  1:49714/1‑300 (MQ=255)
                                                                                                                                                                                                                   tATGCGTTGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACGCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGACTCGTCGTCTGGTTGACGTGGCGCAGGACCTGGTGGTTACCGAAGACGATTGTGGTACCCATGAAGGTATCATGATGACTCCGGTTATCGAGGGTGGTGACGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCg                                                                                   <  2:66201/300‑1 (MQ=255)
                                                                                                                                                                                                                      gCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACGCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGACTCGTCGTCTGGTTGACGTGGCGCAGGACCTGGTGGTTACCGAAGACGATTGTGGTACCCATGAAGGTATCATGATGAATCCGGTTATCGAGGGTGGTGACGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTa                                                                               >  2:256032/1‑301 (MQ=255)
                                                                                                                                                                                                                       cGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACGCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGACTCGTCGTCTGGTTGACGTGGCGCAGGACCTGGTGGTTACCGAAGACGATTGTGGTACCCATGAAGGTATCATGATGACTCCGGTTATCGAGGGTGGTGACGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTaa                                                                              >  2:120256/1‑301 (MQ=255)
                                                                                                                                                                                                                          ggTCTAAAGGCGAAGCCGGATGGCTCCATCATCGAAACGACAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACGCCAGGACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGACGGATACCGCACTGAAGACTGCGACCTCCGGTTACCTGACACGTCGTCTGGTTGACGTGGCGCAGGACCTGGTGGTTACCGAAGACGATTGTGGTACCCATGAAGGTATCATGATGACTCCGGTTATCGAGGGTGGTGACGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAAct                                                                            <  2:70642/300‑1 (MQ=255)
                                                                                                                                                                                                                            ctgtaTGGCGAAGCCGGATGGCTCCTTCAGCGAAACGCCAATCACCTCGACCTTCCGTGAAGGTCTGAACGTACGCCAGTACTTCATCTCCGCCCACTGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGACTCGTCGTCTGGTTGACGTGGCGCAGGACCTGGTGGTTACCGAAGACGATTGTGGTACCCATGAAGGTATCATGATGACTCCGGTTATCGAGGGTGGTGACGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAActgc                                                                          <  1:95714/296‑1 (MQ=255)
                                                                                                                                                                                                                             cTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACGCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGACTCGTCGTCTGGTTGACGTGGCGCAGGACCTGGTGGTTACCGAAGACGATTGTGGTACCCATGAAGGTATCATGATGACTCCGGTTATCGAGGGTGGTGACGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAActg                                                                           >  1:39019/1‑298 (MQ=255)
                                                                                                                                                                                                                                 tGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACGCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGACTCGTCGTCTGGTTGACGTGGCGCAGGACCTGGTGGTTACCGAAGACGATTGTGGTACCCATGAAGGTATCATGATGACTCCGGTTATCGAGGGTGGTGACGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGa                                                                    <  2:39020/301‑1 (MQ=255)
                                                                                                                                                                                                                                                  acaagcattggaacgacgagcACCGCGAAATCACGGGAAGGGCTGAACGTAAGGCGTCACTGCAGTTCCAAGCGCGGTGCTCGTAAAGGTCTGGCGGATACCGCACGGAAAACTGCGAACTCCGGTTACCTGGATCGTCGTCTGGTTGACGGGGCGCAGGACCTGGTGGTTACCGAAGACGATTGTGGTACCCATGAAGGTATCATGATGACTCCGGTGAGCGAGGGTGGTGACGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGGCGTGTAACTGCTGAAGACGTTCTGAAGCCGGGTa                                                   <  1:202300/281‑1 (MQ=255)
                                                                                                                                                                                                                                                                 ccAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACGCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGACTCGTCGTCTGGTTGACGTGGCGCAGGACCTGGTGGTTACCGAAGACGATTGTGGTACCCATGAAGGTATCATGATGACTCCGGTTATCGAGGGTGGTGACGTTAAAGAGCCGCTGCGCGATCGCGTACTGTGTCGTGTACATGCTGAATACGTTCTGAAGCCGGGTACTGCTGATCTCCt                                      >  1:93870/1‑299 (MQ=255)
                                                                                                                                                                                                                                                                           cgAACTTCCGTGAAGGTCTGAACGTACGCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGACTCGTCGTCTGGTTGACGTGGCGCAGGACCTGGTGGTTACCGAAGACGATTGTGGTACCCATGAAGGTATCATGATGACTCCGGTTATCGAGGGTGGTGACGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGAAGCCGGGTACTGCTGATATCCTCGTTCCGCGCa                           <  2:49716/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                 tCCGTGAAGGTCTGAACGTACGCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGACTCGTCGTCTGGTTGACGTGGCGCAGGACCTGGTGGTTACCGAAGACGATTGTGGTACCCATGAAGGTATCATGATGACTCCGGTTATCGAGGGTGGTGACGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGAAGCCGGGTACTGCTGATATCCTCGTTCCGCGCAACACg                      >  2:217819/1‑299 (MQ=255)
                                                                                                                                                                                                                                                                                              gAACGTACGCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGACTCGTCGTCTGGTTGACGTGGCGCAGGACCTGATGGTTACCGAAGACGACTGTGGTACCCATGAAGGTATCATGATGACTCCGGTTATCGAGGGTGGTGACGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGACGACGTTCTGAAGCCGGGTACTGCTGATATTCTCGTTCCGCGCCACACGCTGCTGCACGACCAg       >  2:220546/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                                   tACGCCAGTACTTCAGCTCCACCCACGGTGCTCGTAACGGGCTGGCGGATACCGCACAGAAAACTGCGAGCGCCGGTTAGCTGACTCGACGGCTGGTAAACGTGGCGCAGGACCTGGTGGTTACCGAAGACGATTGTGGTACCCATGAAGGTATCATGATGACTCCGGTTATCGATGGTGGTGACGTTAATGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGAAGCCGGGTACTGCTGATATCCTCGTTCCGCGCAACACGCTGCTGCACGAACAgtggtg  <  2:239962/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                                                                         
TGGTAACTGCGGGCGAACGCTACAACAAAGTTATCGATATCTGGGCTGCGGCGAACGATCGTGTATCCAAAGCGATGATGGATAACCTGCAAACTGAAACCGTGATTAACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGACTCGTCGTCTGGTTGACGTGGCGCAGGACCTGGTGGTTACCGAAGACGATTGTGGTACCCATGAAGGTATCATGATGACTCCGGTTATCGAGGGTGGTGACGTTAAAGAGCCGCTGCGCGATCGCGTACTGGGTCGTGTAACTGCTGAAGACGTTCTGAAGCCGGGTACTGCTGATATCCTCGTTCCGCGCAACACGCTGCTGCACGAACAGTGGTG  >  NC_000913/4187364‑4187955

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