Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 1,169,937 A→G 15.4% H92R (CAT→CGT)  U712_05805 → Transcriptional activator protein med

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068811,169,9370AG15.4% 877.4 / 163.3 377H92R (CAT→CGT) U712_05805Transcriptional activator protein med
Reads supporting (aligned to +/- strand):  ref base A (169/149);  new base G (28/30);  total (198/179)
Fisher's exact test for biased strand distribution p-value = 5.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 379 total aligned reads displayed.

TGGGGCACAAAAGGGTATAAAGGATTACTGAATATACAATCCAAATACAATGTGGACGTCTATTATAAAGAAGGCGTTAAAACAGAAGAAGATATCATAAATGCGATTGAGGATTTTCATAAGAGGGGCGTCAATCTGCTCTACGGCCATGGGAGTGAATATGCTGAGGTGTTTAACTTGGTCGGCGAAGATTATCCAGATATGGAATTCGTCATTTCCAATGCGAAAGCGAAAGCTGATAATGTAACAAGTGTCCATTTTAGCGGGGAAGCAATGGGCTTTTTTGGAGGAATGACTG  >  CP006881/1169789‑1170086
                                                                                                                                                    |                                                                                                                                                     
  gggCACAAAAGGGTATAAAGGATTACTGAATATACAATCCAAATACAATGTGGACGTCTATTATAAACAAGGCGTTAAAACAGAAGAAGATATCATAAATGCGATTGAGGATTTTCATAAGAGGGGCGTCAATCTGCTCTACGGCCGTgg                                                                                                                                                    <  1:5070739/150‑1 (MQ=255)
    gCACAAAAGGGTATAAAGGATTACTGAATATACAATCCAAATACAATGTGGACGTCTATTATAAAGAAGGCGTTAAAACAGAAGAAGATATCATAAATGCGATTGAGGATTTTCATAAGAGGGGCGTCAATCTGCTCTACGGCCATGGGa                                                                                                                                                  >  1:3858489/1‑150 (MQ=255)
     cacaAAAGGGTATAAAGGATTACTGAATATACAATCCAAATACAATGTGGACGTCTATTATAAAGAAGGCGTTAAAACAGAAGAAGATATCATAAATGCGATTGAGGATTTTCATAAGAGGGGCGTCAATCTGCTCTACGGCCATGGGAg                                                                                                                                                 <  1:4151275/150‑1 (MQ=255)
     cacaAAAGGGTATAAAGGATTACTGAATATACAATCCAAATACAATGTGGACGTCTATTATAAAGAAGGCGTTAAAACAGAAGAAGATATCATAAATGCGATTGAGGATTTTCATAAGAGGGGCGTCAATCTGCTCTACGGCCATGGGAg                                                                                                                                                 <  2:3310951/150‑1 (MQ=255)
      acaAAAGGGTATAAAGGATTACTGAATATACAATCCAAATACAATGTGGACGTCTATTATAAAGAAGGCGTTAAAACAGAAGAAGATATCATAAATGCGATTGAGGATTTTCATAAGAGGGGCGTCAATCTGCTCTACGGCCATGGGAGt                                                                                                                                                >  2:3703061/1‑150 (MQ=255)
          aaGGGTATAAAGGATTACTGAATATACAATCCAAATACAATGTGGACGTCTATTATAAAGAAGGCGTAAAAACAGAAGAAGATATCATAAATGCGATTGAGGATTTTCATAAGAGGGGCGTCAATCTGCTCTACGGCCATGGGAGTGAat                                                                                                                                            >  1:3372295/1‑150 (MQ=255)
           aGGGTATAAAGGATTACTGAATATACAATCCAAATACAATGTGGACGTCTATTATAAAGAAGGCGTTAAAACAGAAGAAGATATCATAAATGCGATTGAGGATTTTCATAAGAGGGGCGTCAATCTGCTCTACGGCCATGGGAGTGAata                                                                                                                                           >  2:5212788/1‑150 (MQ=255)
            gggTATAAAGGATTACTGAATATACAATCCAAATACAATGTGGACGTCTATTATAAAGAAGGCGTTAAAACAGAAGAAGATATCATAAATGCGATTGAGGATTTTCATAAGAGGGGCGTCAATCTGCTCTACGGCCATGGGAGTGAatat                                                                                                                                          >  2:3634541/1‑150 (MQ=255)
                 taAAGGATTACTGAATATACAATCCAAATACAATGTGGACGTCTATTATAAAGAAGGCGTTAAAACAGAAGAAGATATCATAAATGCGATTGAGGATTTTCATAAGAGGGGCGTCAATCTGCTCTACGGCCATGGGAGTGAATATGCTGa                                                                                                                                     >  2:1698526/1‑150 (MQ=255)
                  aaaGGATTACTGAATATACAATCCAAATACAATGTGGACGTCTATTATAAAGAAGGCGTTAAAACAGAAGAAGATATCATAAATGCGATTGAGGATTTTCATAAGAGGGGCGTCAATCTGCTCTACGGCCATGGGAGTGAATATGCTGAg                                                                                                                                    >  2:4260400/1‑150 (MQ=255)
                    aGGATTACTGAATATACAATCCAAATACAATGTGGACGTCTATTATAAAGAAGGCGTTAAAACAGAAGAAGATATCATAAATGCGATTGAGGATTTTCATAAGAGGGGCGTCAATCTGCTCTACGGCCATGGGAGTGAATATGCTGAGgt                                                                                                                                  >  2:5649737/1‑150 (MQ=255)
                     ggATTACTGAATATACAATCCAAATACAATGTGGACGTCTATTATAAAGAAGGCGTTAAAACAGAAGAAGATATCATAAATGCGATTGAGGATTTTCATAAGAGGAGCGTCAATCTGCTCTACGGCCGTGGCAGTGAATATGCTGAGgtg                                                                                                                                 >  1:2204329/1‑150 (MQ=255)
                         tACTGAATATACAATCCAAATACAATGTGGACGTCTATTATAAAGAAGGCGTTAAAACAGAAGAAGATATCATAAATGCGATTGAGGATTTTCATAAGAGGGGCGTCAATCTGCTCTACGGCCGTGGGAGTGAATATGCTGAGGTGTTTa                                                                                                                             >  2:1227814/1‑150 (MQ=255)
                          aCTGAATATACATTTCAAATACAATGTGGACGTCTATTATAAAGAAGGCGTTAAAACAGAAGAAGATATCATAAATGCGATTGAGGATTTTCATAAGAGGGGCGTCAATCTGCTCTACGGCCATGGGAGTGAATATGCTGAGGTGTTTaa                                                                                                                            <  1:5770053/150‑1 (MQ=255)
                          aCTGAATATACAATCCAAATACAATGTGGACGTCTATTATAAAGAAGGCGTTAAAACAGAAGAAGATATCATAAATGCGATTGAGGATTTTCATAAGAGGGGCGTCAATCTGCTCTACGGCCATGGGAGTGAATATGCTGAGGTGTTTaa                                                                                                                            >  2:3442410/1‑150 (MQ=255)
                              aaTATACAATCCAAATACAATGTGGACGTCTATTATAAAGAAGGCGTTAAAACAGAAGAAGATATCATAAATGCGATTGAGGATTTTCATAAGAGGGGCGTCAATCTGCTCTACGGCCATGGGAGTGAATATGCTGAGGTGTTTAACTTg                                                                                                                        >  2:4020148/1‑150 (MQ=255)
                                  taCAATCCAAATACAATGTGGACGTCTATTATAAAGAAGGCGTTAAAACAGAAGAAGATATCATAAATGCGATTGAGGATTTTCATAAGAGGGGCGTCAATCTGCTCTACGGCCATGGGAGTGAATATGCTGAGGTGTTTAACTTGGTCg                                                                                                                    <  2:4152047/150‑1 (MQ=255)
                                   aCAATCCAAATACAATGTGGACGTCTATTATAAAGAAGGCGTTAAAACAGAAGAAGATATCATAAATGCGATTGAGGATTTTCATAAGAGGGGCGTCAATCTGCTCTACGGCCATGGGAGTGAATATGCTGAGGTGTTTAACTTGGTCgg                                                                                                                   <  1:5419424/150‑1 (MQ=255)
                                    cAATCCAAATACAATGTGGACGTCTATTATAAAGAAGGCGTTAAAACAGAAGAAGATATCATAAATGCGATTGAGGATTTTCATAAGAGGGGCGTCAATCTGATCTACGGCCATGGGAGTGAATATGCTGTGGTGTGTAACTTGGTCGGc                                                                                                                  >  1:2623878/1‑150 (MQ=255)
                                     aaTCCAAATACAATGTGGACGTCTATTATAAAGAAGGCGTTAAAACAGAAGAAGATATCATAAATGCGATTGAGGATTTTCATAAGAGGGGCGTCAATCTGCTCTACGGCCATGGGAGTGAATATGCTGAGGTGTTTAACTTGGTCGGCg                                                                                                                 <  1:5714744/150‑1 (MQ=255)
                                      aTCCAAATACAATGTGGACGTCTATTATAAAGAAGGCGTTAAAACAGAAGAAGATATCATAAATGCGATTGAGGATTTTCATAAGAGGGGCGTCAATCTGCTCTACGGCCATGGGAGTGAATATGCTGAGGTGTTTAACTTGGTCGGCGa                                                                                                                <  2:3335897/150‑1 (MQ=255)
                                      aTCCAAATACAATGTGGACGTCTATTATAAAGAAGGCGTTAAAACAGAAGAAGATATCATAAATGCGATTGAGGATTTTCATAAGAGGGGCGTCAATCTGCTCTACGGCCATGGGAGTGAATATGCTGAGGTGTTTAACTTGGTCGGCGa                                                                                                                >  1:3705157/1‑150 (MQ=255)
                                       tCCAAATACAATGTGGACGTCTATTATAAAGAAGGCGTTAAAACAGAAGAAGATATCATAAATGCGATTGAGGATTTTCATAAGAGGGGCGTCAATCTGCTCTACGGCCATGGGAGTGAATATGCTGAGGTGTTTAACTTGGTCGGCGaa                                                                                                               <  2:2447131/150‑1 (MQ=255)
                                        ccAAATACAATGTGGACGTTTATTATAAAGAAGGCGTTAAAACAGAAGAAGATATCATAAATGCGATTGAGGATTTTCATAAGAGGGGCGTCAATCTGCTCTACGGCCATGGGAGTGAATATGCTGAGGTGTTTAACTTGGTCGGCGAAg                                                                                                              <  2:938397/150‑1 (MQ=255)
                                           aaTACAATGTGGACGTCTATTATAAAGAAGGCGTTAAAACAGAAGAAGATATCATAAATGCGATTGAGGATTTTCATAAGAGGGGCGTCAATCTGCTCTACGGCCATGGGAGTGAATATGCTGAGGTGTTTAACTTGGTCGGCGAAGAtt                                                                                                           <  2:84048/150‑1 (MQ=255)
                                              aCAATGTGGACGTCTATTATAAAGAAGGCGTTAAAACAGAAGAAGATATCATAAATGCGATTGAGGATTTTCATAAGAGGGGCGTCAATCTGCTCTACGGCCATGGGAGTGAATATGCTGAGGTGTTTAACTTGGTCGGCGAAGATTATc                                                                                                        >  2:4516711/1‑150 (MQ=255)
                                                aaTGTGGACGTCTATTATAAAGAAGGCGTTAAAACAGAAGAAGATATCATAAATGCGATTGAGGATTTTCATAAGAGGGGCGTCAATCTGCTCTACGGCCATGGGAGTGAATATGCTGAGGTGTTTAACTTGGTCGGAGAAGATTATCCa                                                                                                      >  2:1559251/1‑150 (MQ=255)
                                                 aTGTGGACGTCTATTATAAAGAAGGCGTTAAAACAGAAGAAGATATCATAAATGCGATTGAGGATTTTCATAAGAGGGGCGTCAATCTGCTCTACGGCCATGGGAGTGAATATGCTGAGGTGTTTAACTTGGTCGGCGAAGATTATCCAg                                                                                                     >  2:3465848/1‑150 (MQ=255)
                                                   gtgGACGTCTATTATAAAGAAGGCGTTAAAACAGAAGAAGATATCATAAATGCGATTGAGGATTTTCATAAGAGGGGCGTCAATCTGCTCTACGGCCATGGGAGTGAATATGCTGAGGTGTTTAACTTGGTCGGCGAAGATTATCCAGat                                                                                                   <  1:2579944/150‑1 (MQ=255)
                                                    tgGACGTCTATTATAAAGAAGGCGTTAAAACAGAAGAAGATATCATAAATGCGATTGAGGATTTTCATAAGAGGGGCGTCAATCTGCTCTACGGCCATGGGAGTGAATATGCTGAGGTGTTTAACTTGGTCGGCGAAGATTATCCAGata                                                                                                  <  1:1698526/150‑1 (MQ=255)
                                                       aCGTCTATTATAAAGAAGGCGTTAAAACAGAAGAAGATATCATAAATGCGATTGAGGATTTTCATAAGAGGGGCGTCAATCTGCTCTACGGCCATGGGAGTGAATATGCTGAGGTGTTTAACTTGGTCGGCGAAGATTATCCAGATATgg                                                                                               <  1:3020869/150‑1 (MQ=255)
                                                        cGTCTATTATAAAGAAGGCGTTAAAACAGAAGAAGATATCATAAATGCGATTGAGGATTTTCATAAGAGGGGCGTCAATCTGCTCTACGGCCATGGGAGTGAATATGCTGAGGTGTTTAACTTGGTCGGCGAAGATTATCCAGATATGGa                                                                                              >  2:3901797/1‑150 (MQ=255)
                                                        cGTCTATTATAAAGAAGGCGTTAAAACAGAAGAAGATATCATAAATGCGATTGAGGATTTTCATAAGAGGGGCGTCAATCTGCTCTACGGCCATGGGAGTGAATATGCTGAGGTGTTTAACTTGGTCGGCGAAGATTATCCAGATATGGa                                                                                              >  1:5335962/1‑150 (MQ=255)
                                                         gTCTATTATAAAGAAGGCGTTAAAACAGAAGAAGATATCATAAATGCGATTGAGGATTTTCATAAGAGGGGCGTCAATCTGCTCTACGGCCATGGGAGTGAATATGCTGAGGTGTTTAACTTGGTCGGCGAAGATTATCCAGATATGGaa                                                                                             >  2:541011/1‑150 (MQ=255)
                                                          tCTATTATAAAGAAGGCGTTAAAACAGAAGAAGATATCATAAATGCGATTGAGGATTTTCATAAGAGGGGCGTCAATCTGCTCTACGGCCATGGGAGTGAATATGCTGAGGTGTTTAACTTGGTCGGCGAAGATTATCCAGATATGGAAt                                                                                            <  1:678971/150‑1 (MQ=255)
                                                             attatAAAGAAGGCGTTAAAACAGAAGAAGATATCATAAATGCGATTGAGGATTTTCATAAGAGGGGCGTCAATCTGCTCTACGGCCATGGGAGTGAATATGCTGAGGTGTTTAACTTGGTCGGCGAAGATTATCCAGATATGGAATTCg                                                                                         >  1:901679/1‑150 (MQ=255)
                                                               tataAAGAAGGCGTTAAAACAGAAGAAGATATCATAAATGCTATTGAGGATTTTCATAAGAGGGGCGTCAATCTGCTGTACGGCCGTGGGAGTGAATATGCTGAGGTGTTTAACTTGGTCGGCGAAGATTATCCAGATATGGAATTCGTc                                                                                       <  1:952258/150‑1 (MQ=255)
                                                                  aaaGAAGGCGTTAAAACAGAAGAAGATATCATAAATGCGATTGAGGATTTTCATAAGAGGGGCGTCAATCTGCTCTACGGCCATGGGAGTGAATATGCTGAGGTGTTTAACTTGGTCGGCGAAGATTATCCAGATATGGAATTCGTCAtt                                                                                    >  1:2298061/1‑150 (MQ=255)
                                                                   aagaagGCGATAAAACAGAAGAAGATATCATAAATGCGATTGAGGATTTTCATAAGAGGGGCGTCAATCTGCTCTACGGCCATGGGAGTGAATATGCTGAGGTGTTTAACTTGGTCGGCGAAGATTATCCAGATATGGAATTCGTCAttt                                                                                   <  2:1409480/150‑1 (MQ=255)
                                                                      aagGCGTTAAAACAGAAGAAGATATCATAAATGCGATTGAGGATTTTCATAAGAGGGGCGTCAATCTGCTCTACGGCCATGGGAGTGAATATGCTGAGGTGTTTAACTTGGTCGGCGAAGATTATCCAGATATGGAATTCGTCATTTCCa                                                                                >  2:5328591/1‑150 (MQ=255)
                                                                      aagGCGTTAAAACAGAAGAAGATATCATAAATGCGATTGAGGATTTTCATAAGAGGGGCGTCAATCTGCTCTACGGCCATGGGAGTGAATATGCTGAGGTGTTTAACTTGGTCGGCGAAGATTATCCAGATATGGAATTCGTCATTTCCa                                                                                <  1:1040985/150‑1 (MQ=255)
                                                                        ggCGTTAAAACAGAAGAAGATATCATAAATGCGATTGAGGATTTTCATAAGAGGGGCGTCAATCTGCTCTACGGCCATGGGAGTGAATATGCTGAGGTGTTTAACTTGGTCGGCGAAGATTATCCAGATATGGAATTCGTCATTTCCAAt                                                                              <  1:4789020/150‑1 (MQ=255)
                                                                        ggCGTTAAAACAGAAGAAGATATCATAAATGCGATTGAGGATTTTCATAAGAGGGGCGTCAATCTGCTCTACGGCCATGGGAGTGAATATGCTGAGGTGTTTAACTTGGTCGGCGAAGATTATCCAGATATGGAATTCGTCATTTCCAAt                                                                              >  1:242940/1‑150 (MQ=255)
                                                                        ggCGTTAAAACAGAAGAAGATATCATAAATGCGATTGAGGATTTTCATAAGAGGGGCGTCAATCTGCTCTACGGCCATGGGAGTGAATATGCTGAGGTGTTTAACTTGGTCGGCGAAGATTATCCAGATATGGAATTCGTCATTTCCAAt                                                                              <  2:5047627/150‑1 (MQ=255)
                                                                          cGTTAATACAGAAGAAGATATCATAAATGCGATTGAGGATTTTCATAAGAGGGGCGTCAATCTGCTCTACGGCCATGGGAGTGAATATGCTGAGGTGTTTAACTTGGTCGGCGAAGATTATCCAGATATGGAATTCGTCATTTCCAATgc                                                                            >  1:5630872/1‑150 (MQ=255)
                                                                          cGTTAAAACAGAAGAAGATATCATAAATGCGATTGAGGATTTTCATAAGAGGGGCGTCAATCTGCTCTACGGCCATGGGAGTGAATATGCTGAGGTGTTTAACTTGGTCGGCGAAGATTATCCAGATATGGAATTCGTCATTTCCAATgc                                                                            <  2:458941/150‑1 (MQ=255)
                                                                          cGTTAAAACAGAAGAAGATATCATAAATGCGATTGAGGATTTTCATAAGAGGGGCGTCAATCTGCTCTACGGCCATGGGAGTGAATATGCTGAGGTGTTTAACTTGGTCGGCGAAGATTATCCAGATATGGAATTCGTCATTTCCAATgc                                                                            >  2:5117878/1‑150 (MQ=255)
                                                                          cGTTAAAACAGAAGAAGATATCATAAATGCGATTGAGGATTTTCATAAGAGGGGCGTCAATCTGCTCTACGGCCATGGGAGTGAATATGCTGAGGTGTTTAACTTGGTCGGCGAAGATTATCCAGATATGGAATTCGTCATTTCCAATgc                                                                            <  2:5506589/150‑1 (MQ=255)
                                                                              aaaaCAGAAGAAGATATCATAAATGCGATTGAGGATTTTCATAAGAGGGGCGTCAATCTGCTCTACGGCCATGGGAGTGAATATGCTGAGGTGTTTAACTTGGTCGGCGAAGATTATCCAGATATGGAATTCGTCATTTCCAATgcgaaa                                                                        <  1:756327/150‑1 (MQ=255)
                                                                               aaaCAGAAGAAGATATCATAAATGCGATTGAGGATTTTCATAAGAGGGGCGTCAATCTGCTCTACGGCCGTGGGAGTGAATATGCTGAGGTGTTTAACTTGGTCGGCGAAGATTATCCAGATATGGAATTCGTCATTTCCAATgcgaaag                                                                       <  1:1227814/150‑1 (MQ=255)
                                                                                  cAGAAGAAGATATCATAAATGCGATTGAGGATTTTCATAAGAGGGGCGTCAATCTGCTCTACGGCCATGGGAGTGAATATGCTGAGGTGTTTAACTTGGTCGGCGAAGATTATCCAGATATGGAATTCGTCATTTCCAATgcgaaagcga                                                                    >  2:1861307/1‑150 (MQ=255)
                                                                                    gaagaagaTATCATAAATGCGATTGAGGATTTTCATAAGAGGGGCGTCAATCTGCTCTACGGCCGTGGGAGTGAATATGCTGAGGTGTTTAACTTGGTCGGCGAAGATTATCCAGATATGGAATTCGTCATTTCCAATgcgaaagcgaaa                                                                  <  2:1613475/150‑1 (MQ=255)
                                                                                     aagaagaTATCATAAATGCGATTGAGGATTTTCATAAGAGGGGCGTCAATCTGCTCTACGGCCATGGGAGTGAATATGCTGAGGTGTTTAACTTGGTCGGCGAAGATTATCCAGATATGGAATTCGTCATTTCCAATGcgaaagcgaaag                                                                 <  2:3009308/150‑1 (MQ=255)
                                                                                      agaagaTATCATAAATGCGATTGAGGATTTTCATAAGAGGGGCGTCAATCTGCTCTACGGCCATGGGAGTGAATATGCTGAGGTGTTTAACTTGGTCGGCGAAGATTATCCAGATATGGAATTCGTCATTTCCAATGCgaaagcgaaagc                                                                >  1:3598429/1‑150 (MQ=255)
                                                                                       gaagaTATCATAAATGCGATTGAGGATTTTCATAAGAGGGGCGTCAATCTGCTCTACGGCCATGGGAGTGAATATGCTGAGGTGTTTAACTTGGTCGGCGAAGATTATCCAGATATGGAATTCGTCATTTCCAATGCGAAAGCGAAAGCt                                                               <  2:2321976/150‑1 (MQ=255)
                                                                                       gaagaTATCATAAATGCGATTGAGGATTTTCATAAGAGGGGCGTCAATCTGCTCTACGGCCATGGGAGTGAATATGCTGAGGTGTTTAACTTGGTCGGCGAAGATTATCCAGATATGGAATTCGTCATTTCCAATGCGAAAGCGAAAGCt                                                               <  2:2702627/150‑1 (MQ=255)
                                                                                        aagaTAGCATAAATGCGATTGAGGATTTTCATAAGAGGGGCGTCAATCTGCTCTACGGCCGTGGGAGTGAATATGCTGAGGTGTTTAACTTGGTCGGCGAAGATTATCCAGATATGGAATTCGTCATTTCCAATGCGAAAGCGAAAGCTg                                                              <  2:2168305/150‑1 (MQ=255)
                                                                                         agaTATCATAAATGCGATTGAGGATTTTCATAAGAGGGGCGTCAATCTGCTCTACGGCCATGGGAGTGAATATGCTGAGGTGTTTAACTTGGTCGGCGAAGATTATCCAGATATGGAATTCGTCATTTCCAATGCGAAAGCGAAAGCTGa                                                             <  1:3415810/150‑1 (MQ=255)
                                                                                          gaTATCATAAATGCGATTGAGGATTTTCATAAGAGGGGCGTCAATCTGCTCTACGGCCATGGGAGTGAATATGCTGAGGTGTTTAACTTGGTCGGCGAAGATTATCCAGATATGGAATTCGTCATTTCCAATGCGAAAGCGAAAGCTGAt                                                            >  1:4438554/1‑150 (MQ=255)
                                                                                             attaTAAATGCGATTGAGGATTTTCATAAGAGGGGCGTCAATCTGCTCTACCGCCATGGGAGTGAATATGCTGAGGTGTTTAACTTGGTCGGCGAAGATTATCTAGATATGGAATACGTCATTTCCAATGCGAAAGCGAAAGCTGATAAt                                                         <  1:3495587/147‑1 (MQ=255)
                                                                                             aacATAAATGCGATTGAGGATTTTCATAAGAGGGGCGTCAATCTGCTCTACGGCCATGGGAGTGAATATGCTGAGGTGTTTAACTTGGTCGGCGAAGATTATCCAGATATGGAATTCGTCATTTCCAATGCGAAAGCGAAAGCTGATAAt                                                         <  1:5576273/148‑1 (MQ=255)
                                                                                              tCATAAATGCGATTGAGGATTTTCATAAGAGGGGCGTCAATCTGCTCTACGGCCATGGGAGTGAATATGCTGAGGTGTTTAACTTGGTCGGCGAAGATTATCCAGATATGGAATTCGTCATTTCCAATGCGAAAGCGAAAGCTGATAATg                                                        >  2:1931461/1‑150 (MQ=255)
                                                                                                 caaaTGCGATTGAGGATTTTCATAAGAGGGGCGTCAATCTGCTCTACGACCATGGGAGTGAATATGCTGAGGTGTTTAACTTGTTCGGCGAAGATTATCCAGATATGGAATTCGTCATTTCCAATGCGAAAGCGGAAGATGATAATGTaa                                                     >  2:3740462/2‑150 (MQ=255)
                                                                                                   aaTGCGATTGAGGATTTTCATAACACGGGCGTCAATCTGCTCTACGGCCATGGGAGTGAATATGCTGAGGTGTTTAACTTGGTCGGCGAAGATTATCCAGATATGGAATTCGTCATTTCCAGTGCGAAAGCGAAACTTGATAATGTTACa                                                   >  1:1141012/1‑150 (MQ=255)
                                                                                                    aTGTGATTGAGGATTTTCATAAGAGGGGCGTCAATCTGCTCTACGGCCATGGGAGTGAATATGCTGAGGTGTTTAACTTGGTCGGCGAAGATTATCCAGATATGGAATTCGTCATTTCCAATGCGAAAGCGAAAGCTGATAATGTAACaa                                                  <  1:3996634/150‑1 (MQ=255)
                                                                                                    aTGCGATTGAGGATTTTCATAAGAGGGGCGTCAATCTGCTCTACGGCCATGGGAGTGAATATGCTGAGGTGTTTAACTTGGTCGGCGAAGATTATCCAGATATGGAATTCGTCATTTCCAATGCGAAAGCGAAAGCTGATAATGTAACaa                                                  >  2:3614948/1‑150 (MQ=255)
                                                                                                       cGATTGAGGATTTTCATAAGAGGGGCGTCAATCTGCTCTACGGCCATGGGAGTGAATATGCTGAGGTGTTTAACTTGGTCGGCGAAGATTATCCAGATATGGAATTCGTCATTTCCAATGCGAAAGCGAAAGCTGATAATGTAACAAgtg                                               <  2:5819807/150‑1 (MQ=255)
                                                                                                       cGATTGAGGATTTTCATAAGAGGGGCGTCAATCTGCTCTACGGCCATGGGAGTGAATATGCTGAGGTGTTTAACTTGGTCGGCGAAGATTATCCAGATATGGAATTCGTCATTTCCAATGCGAAAGCGAAAGCTGATAATGTAACAAgtg                                               <  1:3142825/150‑1 (MQ=255)
                                                                                                         aTTGAGGATTTTCATAAGAGGGGCGTCAATCTGCTCTACGGCCGTGGGAGTGAATATGCTGAGGTGTTTAACTTGGTCGGCGAAGATTATCCAGATATGGAATTCGTCATTTCCAATGCGAAAGCGAAAGCTGATAATGTAACAAGTGTc                                             >  2:1140558/1‑150 (MQ=255)
                                                                                                         aTTGAGGATTTTCATAAGAGGGGCGTCAATCTGCTCTACGGCCATGGGAGTGAATATGCTGAGGTGTTTAACTTGGTCGGCGAAGATTATCCAGATATGGAATTCGTCATTTCCAATGCGAAAGCGAAAGCTGATAATGTAACAAGTGTc                                             <  1:1857851/150‑1 (MQ=255)
                                                                                                          ttGAGGTTTTCCATAAGAGGGGCGTCAATCTGCTCTACGGCCATGGGAGTGAATATGCTGAGGTGTTTAACTTGGTCGGCGAAGATTATCNAGATATGGANTTCGTCATTTCCAATGCGAAAGCNAAAGCTGATAATGTAACAAGNGTcn                                            <  1:4549989/150‑2 (MQ=255)
                                                                                                          ttGAGGATTTTCATAAGAAGGGCGTAAATCTTCTCTACGGCCATGCAAGTAAATATGCTGAGGTGTATAACTAGGTCGGAGAAGGTTATCCAGATATGGAATTCGTCATTTCCAATGCGAAAGCGAAAGCTGATAATGTAACAAGTGTct                                            >  2:47179/1‑149 (MQ=255)
                                                                                                             aGGATTTTCATAAGAGGGGCGTCAATCTGCTCTACGGCCGTGGGAGTGAATATGCTGAGGTGTTTAACTTGGTCGGCGAAGATTATCCAGATATGGAATTCGTCATTTCCAATGCGAAAGCGAAAGCTGATAATGTAACAAGTGTCCAtt                                         >  2:743894/1‑150 (MQ=255)
                                                                                                               gATTTTCATAAGAGGGGCGTCAATCTGCTCTACGGCCATGGGAGTGAATATGCTGAGGTGTTTAACTTGGTCGGCGAAGATTATCCAGATATGGAATTCGTCATTTCCAATGCGAAAGCGAAAGCTGATAATGTAACAAGTGTCCAtttt                                       <  1:3620187/150‑1 (MQ=255)
                                                                                                                aTTTTCATAAGAGGGGCGTCAATCTGCTCTACGGCCATGGGAGTGAATATGCTGAGGTGTTTAACTTGGTCGGCGAAGATTATCCAGATATGGAATTCGTCATTTCCAATGCGAAAGCGAAAGCTGATAATGTAACAAGTGTCCATTTTa                                      >  2:5245383/1‑150 (MQ=255)
                                                                                                                  tttCATAAGAGGGGCGTCAATCTGCTCTACGGCCATGGTAGTGAATATGCTGAGGTGTTTAACTTGGTCGGCGAAGATTATCCAGATATGGAATTCGTCATTTCCAATGCGAAAGCGAAAGCTGATAATGTAACAAGTGTCCATTTTAGc                                    >  1:294755/1‑150 (MQ=255)
                                                                                                                     cATAAGAGGGGCGTCAATCTGCTCTACGGCCGTGGGAGTGAATATGCTGAGGTGTTTAACTTGGTCGGCGAAGATTATCCAGATATGGAATTCGTCATTTCCAATGCGAAAGCGAAAGCTGATAATGTAACAAGTGTCCATTTTAGCggg                                 <  1:3357206/150‑1 (MQ=255)
                                                                                                                        aaGAGGGGCGTCAATCTGCTCTACGGCCATGGGAGTGAATATGCTGAGGTGTTTAACTTGGTCGGCGAAGATTATCCAGATATGGAATTCGTCATTTCCAATGCGAAAGCGAAAGCTGATAATGTAACAAGTGTCCATTTTAGCGGGGaa                              >  2:4122917/1‑150 (MQ=255)
                                                                                                                         agagGGGCGTCAATCTGCTCTACGGCCATGGGAGTGAATATGCTGAGGTGTTTAACTTGGTCGGCGAAGATTATCCAGATATGGAATTCGTCATTTCCAATGCGAAAGCGAAAGCTGATAATGTAACAAGTGTCCATTTTAGCGGGGAAg                             >  2:3129711/1‑150 (MQ=255)
                                                                                                                          gagGGGCGTCAATCTGCTCTACGGCCGTGGGAGTGAATATGCTGAGGTGTTTAACTTGGTCGGCGAAGATTATCCAGATATGGAATTCGTCATTTCCAATGCGAAAGCGAAAGCTGATAATGTAACAAGTGTCCATTTTAGCGGGGAAGc                            <  1:1618415/150‑1 (MQ=255)
                                                                                                                           agGGGCGTCAATCTGCTCTACGGCCGTGGGAGTGAATATGCTGAGGTGTTTAACTTGGTCGGCGAAGATTATCCAGATATGGAATTCGTCATTTCCAATGCGAAAGCGAAAGCTGATAATGTAACAAGTGTCGATTTTAGCGGGGAAGCa                           >  1:3072319/1‑150 (MQ=255)
                                                                                                                           agGGGCGTCAATCTGCTCTACGGCCATGGGAGTGAATATGCTGAGGTGTTTAACTTGGTCGGCGAAGATTATCCAGATATGGAATTCGTCATTTCCAATGCGAAAGCGAAAGCTGATAATGTAACAAGTGTCCATTTTAGCGGGGAAGCa                           >  2:4913820/1‑150 (MQ=255)
                                                                                                                            ggggCGTCAATCTGCTCTACGGCCATGGGAGTGAATATGCTGAGGTGTTTAACTTGGTCGGCGAAGATTATCCAGATATGGAATTCGTCATTTCCAATGCGAAAGCGAAAGCTGATAATGTAACAAGTGTCCATTTTAGCGGGGAAGCaa                          >  1:1173166/1‑150 (MQ=255)
                                                                                                                                   cAATCTGCTCTACGGCCATGGGAGTGAATATGCTGAGGTGTTTAACTTGGTCGGCGAAGATTATCCAGATATGGAATTCGTCATTTCCAATGCGAAAGCGAAAGCTGATAATGTAACAAGTGTCCATTTTAGCGGGGAAGCAATGGGCtt                   >  2:1342455/1‑150 (MQ=255)
                                                                                                                                   cAATCTGCTCTACGACCGTGGGAGTGAATATGCTGAGGTGTTTAACTTGGTCGGCGAAGATTATCCAGATATGGAATTCGTCATTTCCAATGCGAAAGCGAAAGCTGATAATGTAACAAGTGTCCATTTTAGCGGGGAAGCAATGGGCtt                   >  1:1934870/1‑150 (MQ=255)
                                                                                                                                    aaTCTGCTCTACGGCCATGGGAGTGAATATGCTGAGGTGTTTAACTTGGTCGGCGAAGATTATCCAGATATGGAATTCGTCATTTCCAATGCGAAAGCGAAAACTGATAATGTAACAAGTGTCCATTTTAGCGGAGAAGCAATGGGCttt                  >  1:1352768/1‑150 (MQ=255)
                                                                                                                                      tgtGCCCTACGGCCATGGGAGTGAATATGCTGAGGTGTTTAACTTGGTCGGCGAAGATTATCCAGATATGGAATTCGTCATTTCCAATGCGAAAGCGAAAGCTGATAATGTAACAAGTGTCCATTTTAGCGGGGAAGCAATGGCCttttt                <  1:2916006/148‑1 (MQ=255)
                                                                                                                                      tCTGCTCTACGGCCATGGGAGTGAATATGCTGAGGTGTTTAACTTGGTCGGCGAAGATTATCCAGATATGGAATTCGTCATTTCCAATGCGAAAGCGAAAGCTGATAATGTAACAAGTGTCCATTTTAGCGGGGAAGCAATGGGCttttt                <  1:5800519/150‑1 (MQ=255)
                                                                                                                                      tCTGCCCTACGGCCATGGGAGTGAATATGCTGAGGTGTTTAACTTGGTCGGCGAAGATTATCCAGATATGGAATTCGTCATTTCCAATGCGAAAGCGAAAGCTGATAATGTAACAAGTGTCCATTTTAGCGGGGAAGCAATGGGCttttt                <  1:2691702/150‑1 (MQ=255)
                                                                                                                                        tGCTCTACGGCCATGGGAGTGAATGTGCTGAGGTGTTTAACTTGGTCGGCGAAGATTATCCAGATATGGAATTCGTCATTTCCAATGCGAAAGCGAAAGCTGATAATGTAACAAGTGTCCATTTTAGCGGGGAAGCAATGGGCTTTTTTg              >  2:1293833/1‑150 (MQ=255)
                                                                                                                                        tGCTCTACGGCCATGGGAGTGAATATGCTGAGGTGTTTAACTTGGTCGGCGAAGATTATCCAGATATGGAATTCGTCATTTCCAATGCGAAAGCGAAAGCTGATAATGTAACAAGTGTCCATTTTAGCGGGGAAGCAATGGGCTTTTTTg              >  1:4231332/1‑150 (MQ=255)
                                                                                                                                        tGCTCTACGGCCATGGGAGTGAATATGCTGAGGTGTTTAACTTGGTCGGCGAAGATTATCCAGATATGGAATTCGTCATTTCCAATGCGAAAGCGAAAGCTGATAATGTAACAAGTGTCCATTTTAGCGGGGAAGCAATGGGCTTTTTTg              >  2:5281128/1‑150 (MQ=255)
                                                                                                                                         gCTCTACGGCCATGGGAGTGAATATGCTGAGGTGTTTAACTTGGTCGGCGAAGATTATCCAGATATGGAATTCGTCATTTCCAATGCGAAAGCGAAAGCTGATAATGTAACAAGTGTCCATTTTAGCGGGGAAGCAATGGGCTTTTTTgg             <  1:157168/150‑1 (MQ=255)
                                                                                                                                             tACGGCCGTGGGAGTGAATATGCTGAGGTGTTTAACTTGGTCGGCGAAGATTATCCAGATATGGAATTCGTCATTTCCAATGCGAAAGCGAAAGCTGATAATGTAACAAGTGTCCATTTTAGCGGGGAAGCAATGGGCTTTTTTggagga         <  1:5182850/150‑1 (MQ=255)
                                                                                                                                                ggCCATGGGAGTGAATATGCTGAGGTGTTTAACTTGGTCGGCGAAGATTATCCAGATATGGAATTCGTCATTTCCAATGCGAAAGCGAAAGCTGATAATGTAACAAGTGTCCATTTTAGCGGGGAAGCAATGGGCTTTTTTGGAGGAATg      <  1:5245880/150‑1 (MQ=255)
                                                                                                                                                 gCCATGGGAGTGAATATGCTGAGGTGTTTAACTTGGTCGGCGAAGATTATCCAGATATGGAATTCGTCATTTCCAATGCGAAAGCGAAAGCTGATAATGTAACAAGTGTCCATTTTAGCGGGGAAGCAATGGGCTTTTTTGGAGGAATGa     <  1:4604984/150‑1 (MQ=255)
                                                                                                                                                  ccATGGGAGTGAATATGCTGAGGTGTTTAACTTGGTCGGCGAAGATTATCCAGATATGGAATTCGTCATTTCCAATGCGAAAGCGAAAGCTGATAATGTAACAAGTGTCCATTTTAGCGGGGAAGCAATGGGCTTTGTTGGAGGAATGAc    >  2:2065032/1‑150 (MQ=255)
                                                                                                                                                   cATGGGAGTGAATATGCTGAGGTGTTTAACTTGGTCGGCGAAGATTATCCAGATATGGAATTCGTCATTTCCAATGCGAAAGCGAAAGCTGATAATGTAACAAGTGTCCATTTTAGCGGGGAAGCAATGGGCTTTTTTGGAGGAATGAct   <  2:3548737/150‑1 (MQ=255)
                                                                                                                                                   cATGGGAGTGAATATGCTGAGGTGTTTAACTTGGTCGGCGAAGATTATCCAGATATGGAATTCGTCATTTCCAATGCGAAAGCGAAAGCTGATAATGTAACAAGTGTCCATTTTAGCGGGGAAGCAATGGGCTTTTTTGGAGGAATGAct   >  1:3241173/1‑150 (MQ=255)
                                                                                                                                                    aTGGGAGTGAATATGCTGAGGTGTTTAACTTGGTCGGCGAAGATTATCCAGATATGGAATTCGTCATTTCCAATGCGAAAGCGAAAGCTGATAATGTAACAAGTGTCCATTTTAGCGGGGAAGCAATGGGCTTTTTTGGAGGAATGActg  <  1:1294355/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TGGGGCACAAAAGGGTATAAAGGATTACTGAATATACAATCCAAATACAATGTGGACGTCTATTATAAAGAAGGCGTTAAAACAGAAGAAGATATCATAAATGCGATTGAGGATTTTCATAAGAGGGGCGTCAATCTGCTCTACGGCCATGGGAGTGAATATGCTGAGGTGTTTAACTTGGTCGGCGAAGATTATCCAGATATGGAATTCGTCATTTCCAATGCGAAAGCGAAAGCTGATAATGTAACAAGTGTCCATTTTAGCGGGGAAGCAATGGGCTTTTTTGGAGGAATGACTG  >  CP006881/1169789‑1170086

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