Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,367,397 T→C 22.6% I110T (ATC→ACC)  arnC → undecaprenyl phosphate‑L‑Ara4FN transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,367,3970TC22.6% 136.4 / 6.2 78I110T (ATC→ACC) arnCundecaprenyl phosphate‑L‑Ara4FN transferase
Reads supporting (aligned to +/- strand):  ref base T (22/30);  new base C (2/13);  total (34/44)
Fisher's exact test for biased strand distribution p-value = 6.46e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGCCTGTGAATCGTTGGGGAAAGAGTATGAGATCCTGCTGATTGATGACGGCAGTAGCGATAATTCCGCGCATATGCTGGTCGAAGCCTCACAAGCGGAGAACAGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCAGCGATTATGGCGGGATTCAGTCACGTTACTGGCGACTTAATTATTACCCTTGATGCCGATCTCCAGAATCCGCCAGAAGAAATCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACGCCAGAACCGCCAGGACAGCTGGTTTCGTAAAACCGCTTCGAAGATGATTAACCGGCTTATTCAGCGCACCACTGGCAAAGCGATGGGTGATTACGGTTGTATGCTGCGCGCCTATCGCCGTCATATTGTCGATGCGATGTTGCACTGCCATGAACGCAGCACCTTTATCCCGATTC  >  NC_000913/2367163‑2367627
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                          
caGCCTGTGAATCGTTGGGGAAAGAGTATGAGATCCTGCTGATTGATGACGGCAGTAGCGATAATTCCGCGCATTGGCTGGTCGAAGCCTCACAAGCGGAGAACAGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCAGCGATTATGGCGGGATTCAGTCACGTTACTGGCGACTTAATTATTACCCTTGATGCCGATCTCCAGAATCCGCCAGAAGAAATcccc                                                                                                                                                                                                                                       <  2:882663/239‑1 (MQ=255)
     tgtgAATCGTTGGGGAAAGAGTATGAGATCCTGCTGATTGATGACGGCAGTAGCGATAATTCCGCGCATATGCTGGTCGAAGCCTCACAAGCGGAGAACAGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCAGCGATTATGGCGGGGTTTAGTTAAGTTACTTGGGGACTTATTATTTCCCTTTATTCCCGTCTCCAGAATCCGCCCCAAGAAAACCCCCGccc                                                                                                                                                                                                                                  >  1:1219189/1‑238 (MQ=255)
     tgtgAATCGTTGGGGAAAGAGTATGAGATCCTGCTGATTGATGACGGCAGTAGCGATAATTCCGCGCATATGCTGGTCGAAGCCTCACAAGCGGAGAACAGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCAGCGATTATGGCGGGGTTCCGTTCCCGTACTTGCGACTTAATTATTACCCTTGATGCCGGTCTCCCGAATCCGCCAGAAAAAAAccccccccc                                                                                                                                                                                                                                  >  1:54682/1‑238 (MQ=255)
     tgtgAATCGTTGGGGAAAGAGTATGAGATCCTGCTGATTGATGACGGCAGTAGCGATAATTCCGCGCATATGCTGGTCGAAGCCTCACAAGCGGAGAACAGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCAGCGATTATGGCGGCATTCATCTACGTCACAGGTGACTTCATTTTTTGCCCTTGTGCCCATCCCCACAAACCACCAGACGAAAAACcccccc                                                                                                                                                                                                                                   >  1:156309/1‑238 (MQ=255)
               tGGGGAAAGAGTATGAGATCCTGCTGATTGAGGACGGCAGTAGCGATAATTCCGCGCATATGCTGGTCGAAGCCTCACAAGCGGAGAACAGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCAGCGATTATGGCGGGATTCAGTCACGTTACTGGCGACTTAATTATTACCCTTGATGCCGATCTCCAGAATCCGCCAGAAGAAATCCCCCGCCTGGTGGCaaaa                                                                                                                                                                                                                       <  1:921410/239‑1 (MQ=255)
                 gggAAAGAGTATGAGATCCTGCTGATTGATGACGGCAGTAGCGATAATTCCGCGCATATGCTGGTCGAAGCCTCACAAGCGGAGAACAGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCAGCGATTATGGCGGGGTTCAGTCAAGTTACTGGCGACTTAATTATTACACTTGTTGCCGATTTCCCGGATACCCCAAGAGAAAACTCCCGCCCGGTGGCagaaa                                                                                                                                                                                                                      >  1:735655/1‑237 (MQ=255)
                     aggagtattgaaTCCTGCTGCTTGAGGATGGCGGTAGCGATCTCTCCGTCCAGATGCGGTTCGAAGCCTCACAAGAGTACAACACCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATGCAGCGATTATGGCGGGATTCAGTCACGTTACTGGCGCCTTAATTATTACCCTTGATGCCGAGCTCCAGAATCCGCCAGAAGAAATCCCCCGCCTGGTGGCAAAAGCCGAt                                                                                                                                                                                                                 <  2:1002060/228‑1 (MQ=255)
                      agagTATGAGATCCTGCTGATTGATGACGGCAGTAGCGATAATTCCGCGCATATGCTGGTCGAAGCCTCACAAGCGGAGAACAGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCAGCGATTATGGCGGGATTCAGTCACGTTACTGGCGACTTAATTATTACCCTTGATGCCGATCTCCAGAATCCGCCAGAAGAAATCCCCCGCCTGGTGGCAAAAGCCGATg                                                                                                                                                                                                                <  1:975068/239‑1 (MQ=255)
                        agTATGAGATCCTGCTGATTGATGACGGCAGTAGCGATAATTCCGCGCATATGCTGGTCGAAGCCTCACAAGCGGAGAACAGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCAGCGATTATGGCGGGGTTCTGTCACATTACTTGGGACTTAATTATTATCCCTTATGCCCATCTCCAGAATCCGCCAGAAGAAATAACCCCCCCGGCGGCaaaaaaagcggaa                                                                                                                                                                                                              >  2:835396/1‑230 (MQ=255)
                               aaTCCGGCTGAGTGATGACGGCAGTAGCGATAATTCCGCGCATATGCTGGTCGAAGCCTCACAAGCGGAGAACAGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCAGCGATTATGGCGGGATTCAGTCACGTTACTGGCGACTTAATTATTACCCTTGATGCCGATCTCCAGAATCCGCCAGAAGAAATCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTacg                                                                                                                                                                                                       <  2:54682/238‑1 (MQ=255)
                                aTCCTGCTGATTGATGACGGCAGGAGCGATAATTCCGCGAATATGCTGGTCGAAGCCTCACAAGCGGAGAACAGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCAGCGATTATGGCGGGATTCAGTCACGTTACTGGCGACTTAATTATTACCCTTGATGCCGATCTCCAGAATCCGCCAGAAGAAATCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTacg                                                                                                                                                                                                       <  1:556668/238‑1 (MQ=255)
                                 tCCTGTTGATGGATGCGGGAAGTGGGGATAATTCCGCGAAAATGCTGGTCGAAGCCTCCAAAGCGGAGAACAGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCAGCGATTATGGCGGGATTCAGTCACGTTACTGGCGACTTAATTATTACCCTGGATGCCGATCTCCAGAATCCGCCAGAAGAAATCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTacgac                                                                                                                                                                                                     <  2:395186/239‑1 (MQ=255)
                                 tCCTGCTGATTGATGACGGCAGTAGCGATAATTCCGCGCATATGCTGGTCGAAGCCTCACAAGCGGAGAACAGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCAGCGATTATGGCGGGATTCAGTCACGTTACTGGCGACTTAATTATTACCCTTGATGCCGATCTCCAGAATCCGCCAGAAGAAATCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTacgac                                                                                                                                                                                                     <  1:787869/239‑1 (MQ=255)
                                 tCCTGCTGATTGATGACGGCAGTAGCGATAATTCCGCGCATATGCTGGTCGAAGCCTCACAAGCGGAGAACAGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCAGCGATTATGGCGGGATTCAGTCACGTTACTGGCGACTTAATTATTACCCTTGATGCCGATCTCCAGAATCCGCCAGAAGAAATCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTacgac                                                                                                                                                                                                     <  1:401396/239‑1 (MQ=255)
                                 tCCTGCTGATTGATGACGGCAGTAGCGATAATTCCGCGCATATGCTGGTCGAAGCCTCACAAGCGGAGAACAGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCAGCGATTATGGCGGGATTCAGTCACGTTACTGGCGACTTAATTATTACCCTTGATGCCGATCTCCAGAATCCGCCAGAAGAAATCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTacg                                                                                                                                                                                                       <  1:835396/237‑1 (MQ=255)
                                    tgctgATTGATGACGGCAGTAGCGATAATTCCGCGCATATGCTGGTCGAAGCCTCACAAGCGGAGAACAGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCAGCGATTATGGCGGGATTTAGTCACGTTAATTGCGGCTTTAATTTTTACCCTTGTGCCGGTCTTCCTAAACCCCCAGAAAAAAACCCCCCCCCGGGGGGCAAAAACCATGAAAGGTacggaggg                                                                                                                                                                                                  >  1:1000130/1‑234 (MQ=255)
                                    tgctgATTGATGACGGCAGTAGCGATAATTCCGCGCATATGCTGGTCGAAGCCTCACAAGCGGAGAACAGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCAGCGATTATGGCGGGATTCAGTCACGTTACTTGCGGCTTTATTATTACCCTTTATACCGATCTCCAAAAACCCCCCGAAGAAAACCCCCCCCTGGGGGCAAAAACCCACGAAGGGTacgaaggg                                                                                                                                                                                                  >  1:1052704/1‑235 (MQ=255)
                                      ctgATTGATGACGGCAGTAGCGATAATTCCGCGCATATGCTGGTCGAAGCCTCACAAGCGGAGAACAGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCAGCGATTATGGCGGGGTTTAGTCACGTTACTTGCGACTTTATTATTTCCCTTTTTGCCGGTCTCCCGAAAACCCCCGAAGAAAACCCCCGCCCTGTGGCAAAAAACCAAGAAGGGTACAACGGGGt                                                                                                                                                                                                >  1:397091/1‑239 (MQ=255)
                                         attgatGACGGCAGTAGCGATAATTCCGCGCATATGCTGGTCGAAGCCTCACAAGCGGAGAACAGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCAGCGATTATGGCGGGGTTTAGTCACCTTACTTGCGGCCTTATTATTATCCCTGTTTCCGATCTCCCACATCCCCCACAAAGAACCCCCCCCCCCGTGGGCAAAACCGCTGAAgggtcacacgcgggagg                                                                                                                                                                                             >  2:1254879/1‑224 (MQ=255)
                                           tgatgaCGGCAGTAGCGATAATTCCGCGCATATGCTGGTCGAAGCCTCACAAGCGGAGAACAGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCAGCGATTATGGCGGGATTCAGTCACGTTACTGGCGACTTAATTATTACCCTTGATGCCGATCTCCAGAATCCGCCAGAAGAAATCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGa                                                                                                                                                                                           <  1:1067930/239‑1 (MQ=255)
                                              tgaGGGGAGTAGCGTAAAATCCGCGCATTTGTTGTCGGAGGCCTCACAGGGGGAGAACAGCCATATTGTGTTTATTTTTCTTACCCGCAATTACGGGCAACATTCAGCGATTATGGCGGGATTCGGTCACGTTACTGGCGACTTAATTATTCCCCTTGATGCCGATCTCCAGAATCCGCCAGAAGAAATCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTg                                                                                                                                                                                        <  2:1064388/239‑1 (MQ=255)
                                              tgaCGGCAGTAGCGATAATTCCGCGCATTTGCTGGTCGAAGCCTCACAAGCGGAGAACAGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCAGCGATTATGGCGGGATTCAGTCACGTTACTGGCGACTTAATTATTACCCTTGATGCCGATCTCCAGAATCCGCCAGAAGAAATCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTg                                                                                                                                                                                        <  1:553539/239‑1 (MQ=255)
                                              tgaCGGCAGTAGCGATAATTCCGCGCATATGCTGGTCGAAGCCTCACAAGCGGAGAACAGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCAGCGATTATGGCGGGATTCAGTCACGTTACTGGCGACTTAATTATTACCCTTGATGCCGATCTCCAGAATCCGCCAGAAGAAATCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTg                                                                                                                                                                                        <  1:357740/239‑1 (MQ=255)
                                              tgaCGGCAGTAGCGATAATTCCGCGCATATGCGGGTCGAAGCCTCAAAAGCGGAGAACAGCATTATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCAGCGATTATGGCGGGATTCAGTCACGTTACTGGCGACTTAATTATTACCCTTGATGCCGATCTCCAGAATCCGCCAGAAGAAATCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTg                                                                                                                                                                                        <  1:503932/239‑1 (MQ=255)
                                                  ggCAGTAGCGATAATTCCGCGCATATGCTGGTCGAAGCCTCACAAGCGGAGAACAGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCAGCGATTATGGCGGGGTTTAGTCAACTTTCCTGGGGCTTAATTATTATCCTTTTTGCCGATCCTCCGAATACGCCCGAAGAAAACCACCCCCTCGTGGCCAAAGCCGCTGAAACTTTCGACCTGGGTGGGCCTggaa                                                                                                                                                                                     >  2:391819/1‑235 (MQ=255)
                                                                    gcgcATATGCTGGTCGAAGCCTCACAAGCGGAGAACAGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCAGCGATTATGGCGGGATTCAGTCACGTTACTGGCGGCTTTATTATTACCCTTGATGCCGATCTCCAGAATCCGCCCGAAGAAATACCCCCCCCGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTGGGGGCTGGACGCCAGAAACGCCCGGACa                                                                                                                                                                   >  1:1313442/1‑238 (MQ=255)
                                                                        atatGCTGGTAAACGCCCAAAAGGCGGAAAACAGCATTTTTTTTTCTATTTTGCTTAACCGCAATTCCGGGCAACATTCAGCGATTATGGCGGGTTTCAGTCCCGTTACTGGCGACTTAATTATTCCCCTTGATGCCGATCCCCAGAACCCGCCAGAAGAAACCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACGCCAGAACCGCCAGGACAGCTgg                                                                                                                                                            <  2:1313442/239‑1 (MQ=255)
                                                                                  cGAAGCCTCACAAGCGGAGAACAGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCAGCGATTATGGCGGGGTTTAAATACAGTTCTTGCGACTTCAATATTATCCCTTTTTTCGATCTCCCGAAACCCCCCGAAAAAAACACCCCCCCCGGTGCCAAAAACCCTGTAAGGTAACACCGGGGTGGGGGCGTTCGCCAGAaccaccccgcaaaccgggcttgtgaaa                                                                                                                                                    >  2:1047707/1‑218 (MQ=255)
                                                                                      gCCTCACAAGCGGAGAACAGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCAGCGATTATGGCGGGGTTTAGGCACGTTACCGGCGACTTAATTATTAACCTTTATGCCGATCTCCAGAATCCCCCCGAAGAAAACCCCCCCCTGGCGGGAAAAACCGATGAAGGGTACGGCGTGGGAGGGACTGTACGCCCGAACCGCCCGGACAGCTGGGTTTGGAAAACCg                                                                                                                                               >  2:372848/1‑238 (MQ=255)
                                                                                      gCCTCACAAGCGGAGAACAGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCAGCGATTATGGCGGGGTTCATTCACGTTACTGGCGAGTTAATTATTTCCCTTGCTGCCGATCTCCACAAACCGCCCGAAGAAATCCCCCCCCCGCTGGCAAAAAGCGACGAAGGCTACGCCGTGGGAGGGAGAGTCCGCCAGAACACCCAAGACAgacggtgtttgaaaaacg                                                                                                                                                 >  2:1289591/1‑221 (MQ=255)
                                                                                       ccccaaaggCGGAGACCAGCAATATTGTGTCTATTTTGTTTAACCCCAATTACGGGCAACATTCAGCGATTATGGCGGGATTAAGCCCCGTTACTGGCGATTTATTTATTCCCCTTGATGCCGACCCCCAGAATCCGCCAGAAGAAACCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACGCCAGAACCGCCAGGACAGCTGGTTTCGTAAAACCGCt                                                                                                                                             <  2:759033/231‑1 (MQ=255)
                                                                                       acacacaAGGGGAGAACAGCCATTTTGTGTCTATTTTGCTAAACCGAAATTACGGGCAACATTCAGCGATTATGGCGGGATTCAGTCACGTTCCTGGCGACTTAATTATTACCCTTGATGCCGATCTCCAGAATCCGCCAGAAGAAATCCCCCGCCTTGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACGCCAGAACCGCCAGGACAGCTGGTTTCGTAAAACCGCt                                                                                                                                             <  2:397091/236‑1 (MQ=255)
                                                                                                     accAGCAATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCAGCGATTATGGCGGGATTCAGTCACGTTACTGGCGACTTAATTATTCCCCTTGATGCCGATCTCCAGAACCCGCCAGAAGAAATCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACGCCAGAACCGCCAGGACAGCTGGTTTCGTAAAACCGCTTCGAAGATGATTaa                                                                                                                               <  1:372848/237‑1 (MQ=255)
                                                                                                     aaCAGCCATATTGTGTCTATTTTGCTTAACCGCAATGACGGGAAACATTCAGCGATTATGGCGGGATTCAGTCACGTTACTGGCGCCTTAATTATTACCCTTGATGCCGATCTCCAGAATCCGCCAGAAGAAACCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACGCCAGAACCGCCAGGACAGCTGGTTTCGTAAAACCGCTTCGAAGATGATTaa                                                                                                                               <  1:396801/239‑1 (MQ=255)
                                                                                                       cAGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCAGCGATTATGGCGGGGTTTAGTCACGTTACTTGGGGCCTAATTATTACCCTTTGTTCCGATCTCCCAAATCCCCCCGAAGAAATCCCCCCCCCTGTGGCAAAAACCGCTGAAGGGTATCACGCGGGGGGGGCTGTTAGCCCGAAACGCCCGGACAGCTGGGTTTTGAAAAAcacctcggaaggggtttaa                                                                                                                               >  2:701846/1‑220 (MQ=255)
                                                                                                               ttGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCAGCGATTATGGCGGGATTCAGTCACGTTACTGGCGACTTAATTATTACCCCTTATGCCGATCTCCAGAATCCGCCAGAAGAAAACCCCCCCCTGGTGGCAAAAGCCGATGAAGGGTACGACGTGGTAGGGACTGTACGCCAGAACCGCCCGGACAGCTGGGTTCGTAAAACCGCTTCGAAGATGATTAACCGGCTTAtt                                                                                                                     >  1:342140/1‑239 (MQ=255)
                                                                                                                    tCTATTTTGCTTAACCGCAATTACGGGCAACATTCAGCGATTATGGCGGGATTCAGTCACGTTTCTGGCGACTTAATTATTACCCTTGCTTCCGATCTCCAGAATCCGCCAGAAGAAATCCCCCGCCTGGTTGCTAAAACCGATTAAGGTTACGACGGGGGAGGGGCCGTACGCCAGAACCGCCCGGACAGCTGGGTTCGTAAAAACGCTTCGAAGATGATTAACCGGcgtttttcgc                                                                                                                 >  1:520752/1‑229 (MQ=255)
                                                                                                                            gCTTAACCGCAATTACGGGCAACATTCAGCGATTATGGCGGGGTTTAGTCACGTTACTGGCGCCCTAATTATTACCCTTTCTGCCGATCTCCAGAAACCCCCCGCAGAAATTCCCCCCCCCGTGGCCAAAGACGATGAAGGTTTCGGCGCGGGAGGGGCAGTACCCCCGCACCGCCACGAGACCACGGTTTCGTACAAAGCTTCCCTGATGATTTAACGGCTCATTCCGCGCCGCACTg                                                                                                        >  2:6715/1‑239 (MQ=255)
                                                                                                                                      aaTTGCGGGCAACAGTCGGCGTTTGGGGCGGAAATGAGTCACGTACCGGGCAATTTATTTATCCCCTGTAATGCCGACCACCAGAACCCGCCAGAAAAATTCCCCCGCTTGGGGGCAAAAGCTGAAGAAGGTTACGGTGTGGTGGGAACTGAACGCAAAAACCGCCAGGACAGCTGGTTTCGTAAAACCGCTTCGAAGATGATTAACCGGCTTATTCAGCGCACCACTGGCAAAGCGAt                                                                                              <  2:156309/239‑1 (MQ=255)
                                                                                                                                      aaTTACGGGCAACATTCAGCGATTATGGCGGGATTCAGCCACGTTACTGGCGATTTAATTTTTACCCTTGATGCCGATCCCCAGAATCCGCCAGAAGAAACCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACGCCAGAACCGCCAGGACAGCTGGTTTCGTAAAACCGCTTCGAAGATGATTAACCGGCTTATTCAGCGCACCACTGGCAAAGCGAt                                                                                              <  1:150396/239‑1 (MQ=255)
                                                                                                                                        ttACGGGCAACATTCAGCGATTATGGCGGGATTCAGTCACGTTACTGGCGCCTTAATTATTACCCTTGATGCGGATCTCCAGAATCCGCCAGAAGAAATCCCCCGCCTGGTGGCAAAAGCCGATGa                                                                                                                                                                                                               <  1:583896/126‑1 (MQ=255)
                                                                                                                                        ttACGGGCAACATTCAGCGATTAGGGCGGGATTCAGTCACGTTACGGGCGACTTAATTATTCCCCTTGATGCCGATCTCCAGAATCCGCCAGAAGAAACCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACGCCAGAACCGCCAGGACAGCTGGTTTCGTAAAACCGCTTCGAAGATGATTAACCGGCTTATTCAGCGCACCACTGGCAAAGCGATgg                                                                                            <  1:1047707/239‑1 (MQ=255)
                                                                                                                                                aaCATTCAGCGATTATGGCGGAATTCAGTCACGTTACTGGCGACTTAATTATTCCCCTTGATGCCGATCTCCAGAATCCGCCAGAAGAAATCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACGCCAGAACCGCCAGGACAGCTGGTTTCGTAAAACCGCTTCGAAGATGATTAACCGGCTTATTCAGCGCACCACTGGCAAAGCGATGGGTGATTa                                                                                     <  1:701846/238‑1 (MQ=255)
                                                                                                                                                   aTTCAGCGATTATGGCGGGATTCAGTCACGTTACTGGCGACTTAATTATTACCCTTGATGCCGATCTCCAGAATCCGCCAGAAGAAATCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTa                                                                                                                                                                                               <  1:389578/131‑1 (MQ=255)
                                                                                                                                                                   cgggatctagcactgTTACGGGCGCTTAAATTATTGCCCGTGATGCGGGTTCCCAAAACCAGCCAAAAAAAACCCCCCGCCTGGGGGCAAAGGCAGAGAAGGGTTAGGAGGTGGTAGGGACGGTACGCCAGAACCGCCAGGACAGCTGGTTTCGTAAAACCGCTTCGAAGATGATTAACCGGCTTATTCAGCGCACCACTGGCAAAGCGATGGGTGATTACGGTTGTATGCTgcgcgc                                                                   <  2:827248/224‑1 (MQ=255)
                                                                                                                                                                        tccactccgtttATTGGCCATTAAATTATACCCCTTGTCGCCGATCTCCAGAACCCCCCAGAAAAAACCCCCCGCCGGGTGGAAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACGCCAGAACCGCCAGGACAGCTGGTTTCGTAAAACCGCTTCGAAGATGATTAACCGGCTTATTCAGCGCACCACTGGCAAAGCGATGGGTGATTACGGTTGTATGCTGCGCGCCTATCg                                                             <  2:342140/229‑1 (MQ=255)
                                                                                                                                                                                cGTTACTGGCGACTTAATTATTACCCTGGATGCCGATCACCAGAACCCGCAAGAAGAAATCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGGGGTAGGGACTGTACGCCAGAACCGCCAGGACAGCTGGTTTCGTAAAACCGCTTCGAAGATGATTAACCGGCTTATTCAGCGCACCACTGGCAAAGCGATGGGTGATTACGGTTGTATGCTGCGCGCCTATCGCCGTCata                                                     <  1:665118/239‑1 (MQ=255)
                                                                                                                                                                                 gTTACTGGCGACTTAATTATTACCCTTGATGCCGATCTCCAGAATCCGCCAGAAGAAATCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACGCCAGAACCGCAAGGACAGCTGGTTTCGTAAAACCGCTTCGAAGATGATTAACCGGCTTATTCAGCGCACCACTGGCAAAGCGATGGGTGATTACGGTTGTATGCTGCGCGCCTATCGCCGTCatat                                                    >  1:166865/1‑239 (MQ=255)
                                                                                                                                                                                    aCTGGCGACTTAATTATTACCCTTGATGCCGATCTCCAGAATCCGCCAGAAGAAATCCCCCGCCTGGTGGCAAAAGCCGATg                                                                                                                                                                                                                <  1:1446224/82‑1 (MQ=255)
                                                                                                                                                                                    aCTGGCGACTTAATTATTACCCTTGATGCCGATCTCCAGAATCCGCCAGAAGAAATCCCCCGCCTGGTGGCAAAAGCCGATg                                                                                                                                                                                                                >  2:1446224/1‑82 (MQ=255)
                                                                                                                                                                                        ggCGACTTAATTATTACCCTTGATGCCGATCTCCAGAATCCGCCAGAAGAAATCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACGCCAGAACCGCCAGGACAGCTGGTTTCGTAAAACCGCTTCGAAGATGATTAACCGGCTTATTCAGCGCACCACTGGCAAAGCGATGGGTGATTACGGTTGTATGCTGCGCGCCTATCGCCGTCATAtt                                                   >  1:1397058/1‑234 (MQ=255)
                                                                                                                                                                                        ggCGACTTAATTATTACCCTGGATGCCGATCTCCAGAACCCCCAAGAAGAAACCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACGCCAGAACCGCCAGGACAGCTGGTTTCGTAAAACCGCTTCGAAGATGATTAACCGGCTTATTCAGCGCACCACTGGCAAAGCGATGGGTGATTACGGTTGTATGCTGCGCGCCTATCGCCGTCATATTGTcga                                              <  1:1254879/239‑1 (MQ=255)
                                                                                                                                                                                           gACTAAATTATTCCCCTGGATGCCGATCTCCAGAATCCGCCAGAAGAAATCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACGCCAGAACCGCCAGGACAGCTGGTTTCGTAAAACCGCTTCGAAGATGATTAACCGGCTTATTCAGCGCACCACTGGCAAAGCGATGGGTGATTACGGTTGTATGCTGCGCGCCTATCGCCGTCATATTGTcgatg                                            <  1:1424038/238‑1 (MQ=255)
                                                                                                                                                                                            gatttAATTATTCCCCGTGATGCCTATCGCCAGAATCCGCCAAAAGAATTCCCCCGCCTGGTGGCAAAAGCCGAGGAAGGTTACGACGTGGAAGGGACTTAACGCCAGAACCGCCAGGACAGCTGGTTTCGTAAAACCGCTTCGAAGATGATTAACCGGCTTATTCAGCGCACCACTGGCAAAGCGATGGGTGATTACGGTTGTATGCTGCGCGCCTATCGCCGTCATATTGTcgatg                                            <  1:1353169/235‑1 (MQ=255)
                                                                                                                                                                                            aCTTAATTATTACCCTTGATGCCGATCTCCAGAATCCGCCAGAAGAAATCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACGCCAGAACCGCCAGGACAGCTGGTTTCGTAAAACCGCTTCGAAGATGATTAACCGGCTTATTCAGCGCACCACTGGCAAAGCGATGGGTGATTACGGTTGTATGCTGCGCGCCTATCGCCGTCATATTGTcgatgcg                                          >  1:242535/1‑239 (MQ=255)
                                                                                                                                                                                               ttaaaTTATTCCCCTGAATGCGAATCTCCAGAAACGCCAGAAAAAAACCCCCCCCTGGGGGGCAAAGCCCGTTGAAGGTTACGACGGGGTAGGGCCTGTACGCCGGAACCGCCAGGACAGCTGGTTTCGTAAAACCGCTTCGAAGATGATTAACCGGCTTATTCAGCGCACCACTGGCAAAGCGATGGGTGATTACGGTTGTATGCTGCGCGCCTATCGCCGTCATATTGTcgatgcga                                         <  2:1000130/236‑1 (MQ=255)
                                                                                                                                                                                                  cattattACCCTCAATGCAGATCTCCAGAATCCGCCAGAAGAAATCCCCCCCCTGGTGCCAAAAGACGATCAAGGTTACGACGTGGAAGGTACTGTACGCAAGAACGGCCAGGACCGCTGGGTTCGTAAAACCGCTTCGCAGACGTCTAACCGGCTCATTCAGAGCACCACTGACAAAGCGATGGGTGATGCCGGATGAATGCTGCGCTCCGATCGACCACATATGGTCGGTGCGAGGt                                      >  2:1212678/2‑239 (MQ=255)
                                                                                                                                                                                                      tttccccTTGAGGCCAGTCTCAAGAGCCGGCAAGAAGAAATCCCCCGCCTGGTGGCAAAAGCCGAGGAAGGTTGCGACGTGGTAGGGACTGTACGCCAGAACCGCCAGGACAGCTGGTTTCGTAAAACCGCTTCGAAGATGATTAACCGGCTTATTCAGCGCACCACTGGCAAAGCGATGGGTGATTACGGTTGTATGCTGCGCGCCTATCGCCGTCATATTGTCGATGCGATGTTGCa                                  <  2:1303471/235‑1 (MQ=255)
                                                                                                                                                                                                      cttccccTGGAGGCTGATCCCAAGAACCCCCCAGAAGAAACCCCCCGCCTGGTGGAAAAAGCCGATGAGGGTTACGACGTGGGAGGGACTGTACGCCAGAACCGCCAGGACAGCTGGTTTCGTAAAACCGCTTCGAAGATGATTAACCGGCTTATTCAGCGCACCACTGGCAAAGCGATGGGTGATTACGGTTGTATGCTGCGCGCCTATCGCCGTCATATTGTCGATGCGATGTTGCa                                  <  2:271681/235‑1 (MQ=255)
                                                                                                                                                                                                      atgaccctggtgcgcgaccccaaGAACCCGCCGGAAAAAACCCCCCGCTGGGTGGAAAAGGCCGATGAAGGTTAGGACGTGGTAGGGCTTGTACCCCAGAACCGCCAGGACAGCTGGTTTCGTAAAACCGCTTCGAAGATGATTAACCGGCTTATTCAGCGCACCACTGGCAAAGCGATGGGTGATTACGGTTGTATGCTGCGCGCCTATCGCCGTCATATTGTCGATGCGATGTTGCa                                  <  2:534050/217‑1 (MQ=255)
                                                                                                                                                                                                        taCCCTTGAGGCCGATCTCAAGAATCCGCAAGAAGAAATCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACGCCAGAACCGCCAGGACAGCTGGTTTCGTAAAACCGCTTCGAAGATGATTAACCGGCTTATTCAGCGCACCACTGGCAAAGCGATGGGTGATTACGGTTGTATACTGCGCGCCTATCGCCGTCATATTGTCGATGCGATGTTGCACt                                <  1:6715/239‑1 (MQ=255)
                                                                                                                                                                                                         aCCCTTGATGCCGATCTCCAGAATCCGCCAGAAGAAATCCCCCGCCTGGTGGCAAAAGCCGCTGACGGTTACGACGTGGTAGGGACTGTACGCCAGAACCGCCAGGCCAGCTGGTTTCGTAAACCCGCTTCGAAGATGATTAACCGGCTTATTCACCGCACAACTGGCAACGCGAGGGGTGATTACGGTTGTATCATGCGCGCCTATCGCCGTCATATTGGCGGTGCGATGTTGTACt                                >  2:816636/1‑238 (MQ=255)
                                                                                                                                                                                                         aCCCTTGATGCCGATCTCCAGAATCCGCCAGAAGAAATCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACGCCAGAACCGCCAGGACAGCTGGTTTCGTAAAACCGCTTCGAAGATGATTAACCGGCTTATTCAGCGCACCACTGGCAAAGCGATGGGTGATTACGGTTGTATGCTGCGCGCCTATCGCCGTCATATTGTCGATGCGATGTTGCACt                                >  1:271681/1‑238 (MQ=255)
                                                                                                                                                                                                          cccTTGATGCCGATCTCCAGAATCCGCCAGAAGAAATCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACGCCAGAACCGCCAGGACAGCTGGTTTCGTAAAACCGCTTCGAAGATGATTAACCGGCTTATTCAGCGCACCACTGGCAAAGCGATGGGTGATTACGGTTGTATGCTGCGCGCCTATCGCCGTCATATTGTCGATGCGATGTTGCACTg                               >  1:1039859/1‑238 (MQ=255)
                                                                                                                                                                                                             ttGATGCCGATCTCCAGAATCCGCCAGAAGAAATCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACGCCAGAACCGCCAGGACAGCTGGTTTCGTAAAACCGCTTCGAAGATGATTAACCGGCTTATTCAGCGCACCACTGGCAAAGCGATGGGTGATTACGCTTGTATGCTGCGCGCCTATCGCCGTCATATTGACGTTGCGATGATGCACTGCCAt                           >  2:1374266/1‑239 (MQ=255)
                                                                                                                                                                                                              tGATGCCGATCTCCAGAATCCGCCAGAAGAAATCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACGCCAGAACCGCCAGGACAGCTGGTTTCGTAAAACCGCTTCGAAGATGATTAACCGGCTTATTCAGCGCACCACTGGCAAAGCGATGGGTGATTACGGTTGTATGCTGCGCGCCTATCGCCGTCATATTGTCGATGCGATGTTGCACTGCCATg                          >  1:1192647/1‑239 (MQ=255)
                                                                                                                                                                                                              tGATGCCGATCCCCAGAATCCGCCAGAAGAAATCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACGCCAGAACCGCCAGGACAGCTGGTTTCGTAAAACCGCTTCGAAGATGATTAACCGGCTTATTCAGCGCACCACTGGCAAAGCGATGGGTGATTACGGTTGTATGCTGCGCGCCTATCGCCGTCATATTGTCGATGCGATGTTGCACTGCCATg                          <  1:804880/239‑1 (MQ=255)
                                                                                                                                                                                                               gATGCCGATCTCCAGAATCCGCCAGAAGAAATCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACGCCAGAACCGCCAGGACAGCTGGTTTCGTAAAACCGCTTCGAAGATGATTAACCGGCTTATTCAGCGCACCACTGGCAAAGCGATGGGTGATTACGGTTGTATGCTGCGCGCCTATCGCCGTCATATTGTCGATGCGATGTTGCACTGCCATGa                         >  2:502595/1‑239 (MQ=255)
                                                                                                                                                                                                               gATGCCGATCTCCAGAATCCGCCAGAAGAAATCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACGCCAGAACCGCCAGGACAGCTGGTTTCGTAAAACCGCTTCGAAGATGATTAACCGGCTTATTCAGCGCACCACTGGCAAAGCGATGGGTGATTACGGTTGTATGCTGCGCGCCTATCGCCGTCATATTGTCGATGCGATGTTGCACTGCCATGa                         >  1:751975/1‑239 (MQ=255)
                                                                                                                                                                                                               gATGCCGATCTCCAGAATCCGCCAGAAGAAATCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACGCCAGAACCGCCAGGACAGCTGGTTTCGTAAAACCGCTTCGAAGATGATTAACCGGCTTATTCAGCGCACCACTGGCAAAGCGATGGGTGATTACGGTTGTATGCTGCGCGCCTATCGCCGTCATATTGTCGATGCGATGTTGCACTGCCATGa                         >  1:749624/1‑239 (MQ=255)
                                                                                                                                                                                                               gATGCCGATCTCCAGAATCCGCCAGAAGAAATCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACGCCAGAACCGCCAGGACAGCTGGTTTCGTAAAACCGCTTCGAAGATGATTAACCGGCTTATTCAGCGCACCACTGGCAAAGCGATGGGTGATTACGGTTGTATGCTGCGCGCCTATCGCCGTCATATTGTCGATGCGATGTTGCACTGCCATGa                         >  1:565222/1‑239 (MQ=255)
                                                                                                                                                                                                               gATGCCGATCTCCAGAATCCGCCAGAAGAAATCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACGCCAGAACCGCCAGGACAGCTGGTTTCGTAAAACCGCTTCGAAGATGATTAACCGGCTTATTCAGCGCACCACTGGCAAAGCGATGGGTGATTACGGTTGTATGCTGCGCGCCTATCGCCGTCATATTGTCGATGCGATGTTGCACTGCCATGa                         >  1:1128689/1‑239 (MQ=255)
                                                                                                                                                                                                                     gATCCCCAGAACCCGCCAGAAGAAACCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACGCCAGAACCGCCAGGACAGCTGGTTTCGTAAAACCGCTTCGAAGATGATTAACCGGCTTATTCAGCGCACCACTGGCAAAGCGATGGGTGATTACGGTTGTATGCTGCGCGCCTATCGCCGTCATATTGTCGATGCGATGTTGCACTGCCATGAACgcag                   <  1:1444157/239‑1 (MQ=255)
                                                                                                                                                                                                                      accccaaaaacccGCCAGAAAAAACCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACGCCAGAACCGCCAGGACAGCTGGTTTCGTAAAACCGCTTCGAAGATGATTAACCGGCTTATTCAGCGCACCACTGGCAAAGCGATGGGTGATTACGGTTGTATGCTGCGCGCCTATCGCCGTCATATTGTCGATGCGATGTTGCACTGCCATGAACgcag                   <  2:242535/227‑1 (MQ=255)
                                                                                                                                                                                                                                   gCCAGAAGAAATCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACGCCAGAACCGCCAGGACAGCTGGTTTCGTAAAACCGCTTCGAAGATGATTAACCGGCTTATTCAGCGCACCACTGGCAAAGCGATGGGTGATTACGGTTGTATGCTGCGCGCCTATCGCCGTCATATTGTCGATGCGATGTTGCACTGCCATGAACGCAGCACCTTTATCCCGa     >  1:44678/1‑239 (MQ=255)
                                                                                                                                                                                                                                     cAGAAGAAATCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACGCCAGAACCGCCAGGACAGCTGGTTTCGTAAAACCGCTTCGAAGATGATTAACCGGCTTATTCAGCGCACCACTGGCAAAGCGATGGGTGATTACGGTTGTATGCTGCGCGCCTATCGCCGTCATATTGTCGATGCGATGTTGCACTGCCATGAACGCAGCACCTTTATCCCGAtt   <  2:1039859/239‑1 (MQ=255)
                                                                                                                                                                                                                                     aagaagaAACCCCCCGCCGGGTGGCAAAAGCCGATGAAGGTTACGGGGGGGTAGGGACTGTACGCCAGAACCGCCAGGACAGCTGGTTTCGTAAAACCGCTTCGAAGATGATTAACCGGCTTATTCAGCGCACCACTGGCAAAGCGATGGGTGATTACGGTTGTATGCTGCGCGCCTATCGCCGTCATATTGTCGATGCGATGTTGCACTGCCATGAACGCAGCACCTTTATCCCGAtt   <  1:157679/238‑1 (MQ=255)
                                                                                                                                                                                                                                       tgtaaccaTCCCGCGCTGGGAGGAAAAATCCGATGAAGGTGTCGAAGGGTAAGGTACTGAAAGCAAGAACCTCCAGGACAGCTGGTTCCGTAAAACCGCTGCGAAGATGAATAACCGGCTTATTCAGCGCACCACTGGCAAAGCGATGGGTGATTACGGTTGTATGCTGCGCGCCTATCGCCGTCATATTGTCGATGCGATGTTGCACTGCCAGGAACGCAGCACCTTTATCCCGATTc  <  2:735655/232‑1 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                          
CAGCCTGTGAATCGTTGGGGAAAGAGTATGAGATCCTGCTGATTGATGACGGCAGTAGCGATAATTCCGCGCATATGCTGGTCGAAGCCTCACAAGCGGAGAACAGCCATATTGTGTCTATTTTGCTTAACCGCAATTACGGGCAACATTCAGCGATTATGGCGGGATTCAGTCACGTTACTGGCGACTTAATTATTACCCTTGATGCCGATCTCCAGAATCCGCCAGAAGAAATCCCCCGCCTGGTGGCAAAAGCCGATGAAGGTTACGACGTGGTAGGGACTGTACGCCAGAACCGCCAGGACAGCTGGTTTCGTAAAACCGCTTCGAAGATGATTAACCGGCTTATTCAGCGCACCACTGGCAAAGCGATGGGTGATTACGGTTGTATGCTGCGCGCCTATCGCCGTCATATTGTCGATGCGATGTTGCACTGCCATGAACGCAGCACCTTTATCCCGATTC  >  NC_000913/2367163‑2367627

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: