Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,875,096 G→T 100% intergenic (‑80/+47) leuE ← / ← dmlR leucine efflux protein/transcriptional activator of dmlA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,875,0960GT98.2% 714.2 / 3.0 221intergenic (‑80/+47)leuE/dmlRleucine efflux protein/transcriptional activator of dmlA
Reads supporting (aligned to +/- strand):  ref base G (4/0);  new base T (98/119);  total (102/119)
Fisher's exact test for biased strand distribution p-value = 4.39e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 221 total aligned reads displayed.

CCAACACAATAAAAATGGCCCCAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGCTGACGTCGCCAGCCTTGA  >  CP009273/1874950‑1875241
                                                                                                                                                  |                                                                                                                                                 
ccAACACAATAAAAATGGCCCCAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGATaa                                                                                                                                                 >  2:1976123/1‑149 (MQ=255)
 cAACACAATAAAAATGGCCCCAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGATaaa                                                                                                                                                <  1:2975416/149‑1 (MQ=255)
  aaCACAATAAAAATGGCCCCAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGATaaaa                                                                                                                                               >  1:1440468/1‑149 (MQ=255)
  aaCACAATAAAAATGGCCCCAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGATaaaa                                                                                                                                               <  2:553311/149‑1 (MQ=255)
   acacAATAAAAATGGCCCCAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGATAAAAt                                                                                                                                              >  2:183666/1‑149 (MQ=255)
    cacaATAAAAATGGCCCCAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGATAAAATc                                                                                                                                             <  1:3405271/149‑1 (MQ=255)
      caATAAAAATGGCCCCAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGATAAAATCAg                                                                                                                                           >  1:3044042/1‑149 (MQ=255)
      caATAAAAATGGCCCCAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGATAAAATCAg                                                                                                                                           <  2:561350/149‑1 (MQ=255)
      caATAAAAATGGCCCCAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGATAAAATCAg                                                                                                                                           <  2:2019757/149‑1 (MQ=255)
      caATAAAAATGGCCCCAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGATAAAATCAg                                                                                                                                           <  2:2776950/149‑1 (MQ=255)
          acaaaTGGCCACAACCAGTTAGTTCCAGTAATTCATAACCCCGTATTCATCTAACACATTAACCTCTTTAATTATCTTTCGATCTTGAGCGATTAATTGTGATTATGCTTACCTATCTGTATAGGCTTAGAAAGGATAATATCTTGttt                                                                                                                                       <  2:3709043/147‑1 (MQ=255)
          aaaaaTGGCCCCACCCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGATAAAATCAGGttt                                                                                                                                       <  2:839842/149‑1 (MQ=255)
           aaaaTGGCCCCAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGATAAAATCAGGtttt                                                                                                                                      <  2:693467/149‑1 (MQ=255)
            aaaTGGCCCCAACCAGATAGGTCCAGTAATTCAGAACCCAGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGATAAAATCAGGTTTTa                                                                                                                                     >  1:4011442/1‑149 (MQ=255)
             aaTGGCCCCAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGATAAAATCAGGTTTTaa                                                                                                                                    <  2:2658908/149‑1 (MQ=255)
             aaTGGCCCCAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGATAAAATCAGGTTTTaa                                                                                                                                    >  1:614916/1‑149 (MQ=255)
              aTGGCCCCAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGATAAAATCAGGTTTTAAc                                                                                                                                   <  1:733455/149‑1 (MQ=255)
               tGGCCCCAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGATAAAATCAGGTTTTAAcc                                                                                                                                  <  2:449869/149‑1 (MQ=255)
               tGGCCCCAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGATAAAATCAGGTTTTAAcc                                                                                                                                  <  1:2583448/149‑1 (MQ=255)
                  ccccAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGATAAAATCAGGTTTTAACCTGa                                                                                                                               <  2:2688674/149‑1 (MQ=255)
                  ccccAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGATAAAATCAGGTTTTAACCTGa                                                                                                                               >  1:2188175/1‑149 (MQ=255)
                  ccccAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGATAAAATCAGGTTTTAACCTGa                                                                                                                               <  1:1564913/149‑1 (MQ=255)
                    ccAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGATAAAATCAGGTTTTAACCTGata                                                                                                                             <  2:335669/149‑1 (MQ=255)
                    ccAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGATAAAATCAGGTTTTAACCTGata                                                                                                                             <  1:57524/149‑1 (MQ=255)
                    ccAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGATAAAATCAGGTTTTAACCTGata                                                                                                                             <  2:1310240/149‑1 (MQ=255)
                     cAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGATAAAATCAGGTTTTAACCTGatat                                                                                                                            <  2:1032014/149‑1 (MQ=255)
                      aaCCAGATAGGTCCATTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGATAAAATCAGGTTTTAACCTGATATc                                                                                                                           <  2:1649086/149‑1 (MQ=255)
                          aGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGATAAAATCAGGTTTTAACCTGATATCAAcc                                                                                                                       <  1:2666908/149‑1 (MQ=255)
                               ggTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGATAAAATCAGGTTTTAACCTGATATCAACCCGATa                                                                                                                  >  2:4035901/1‑149 (MQ=255)
                                 tCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGATAAAATCAGGTTTTAACCTGATATCAACCCGATAAt                                                                                                                <  1:2326271/149‑1 (MQ=255)
                                 tCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGATAAAATCAGGTTTTAACCTGATATCAACCCGATAAt                                                                                                                <  2:2862789/149‑1 (MQ=255)
                                   cAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGATAAAATCAGGTTTTAACCTGATATCAACCCGATAATTg                                                                                                              >  2:61705/1‑149 (MQ=255)
                                     gTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGATAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGaa                                                                                                            <  1:1804783/149‑1 (MQ=255)
                                        aTTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGATAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCa                                                                                                         <  1:869976/149‑1 (MQ=255)
                                        aTTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGATAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCa                                                                                                         <  1:31911/149‑1 (MQ=255)
                                        aTTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGATAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCa                                                                                                         >  1:211353/1‑149 (MQ=255)
                                         ttCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGATAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCAt                                                                                                        <  2:947469/149‑1 (MQ=255)
                                             gAACCCCGTATTCAGCGTACACATTAAACTCTTTAATTTTATTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGATAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAAt                                                                                                    <  1:2401099/149‑1 (MQ=255)
                                                 cccGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCTGCTTAGAAAGGATAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAgg                                                                                                >  2:2229649/1‑149 (MQ=255)
                                                 cccGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGATAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAgg                                                                                                <  1:1173440/149‑1 (MQ=255)
                                                 cccGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGATAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAgg                                                                                                <  1:2752426/149‑1 (MQ=255)
                                                  ccGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGATAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGc                                                                                               >  1:1967439/1‑149 (MQ=255)
                                                          aGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGATAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCa                                                                                       <  1:2695053/149‑1 (MQ=255)
                                                           gCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGATAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAg                                                                                      >  1:2572181/1‑149 (MQ=255)
                                                           gCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGATAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAg                                                                                      <  1:1517337/149‑1 (MQ=255)
                                                            cGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGATAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGt                                                                                     <  1:433097/149‑1 (MQ=255)
                                                                  caTTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGATAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGaaaaa                                                                               >  1:3567290/1‑149 (MQ=255)
                                                                    ttAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGATAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAAtt                                                                             <  1:3770019/149‑1 (MQ=255)
                                                                     tAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGATAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTc                                                                            <  1:2075840/149‑1 (MQ=255)
                                                                      aaCCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGATAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCg                                                                           >  2:2323812/1‑149 (MQ=255)
                                                                       tccTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGATAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCgg                                                                          <  1:1122996/148‑1 (MQ=255)
                                                                            tttAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGATAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGt                                                                     >  1:354852/1‑149 (MQ=255)
                                                                            tttAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGATAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGt                                                                     <  2:2048670/149‑1 (MQ=255)
                                                                               aaTTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGATAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCt                                                                  <  2:3180526/149‑1 (MQ=255)
                                                                                aTTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGATAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTc                                                                 >  1:1499971/1‑149 (MQ=255)
                                                                                aTTATCTTTCGATCAAGCGCGATTAAAGGTTAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGATAAAAACAGGTTTTAACCTGATATCAACACGATAATAGAATCATTAATCAGGCATGATCCAGAGAAAAATTCGGTTAGTGCTc                                                                 >  1:62786/1‑149 (MQ=255)
                                                                                 ttATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGATAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCg                                                                >  2:2814717/1‑149 (MQ=255)
                                                                                 ttATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGATAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCg                                                                >  2:1251255/1‑149 (MQ=255)
                                                                                 ttATCTTTCGATCATGCGCGATTAAAGGTGAAGATGCTAACCAATCTGTCGCGGCTTAGAAAGGATAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGCAACATTCGGGTAGTGCTCg                                                                >  2:2575156/1‑149 (MQ=255)
                                                                                    tCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGATAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCa                                                             >  2:121443/1‑149 (MQ=255)
                                                                                      tttCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGATAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCTGGTAGTGCTCGGCaaa                                                           >  2:3001421/1‑149 (MQ=255)
                                                                                          gATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGATAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATAc                                                       <  1:3213609/149‑1 (MQ=255)
                                                                                          gATCATGCGCGATTAAAGGTGAAGATGCTAACCAATCTGTAGCGGCTTAGAAAGGATAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATAc                                                       <  2:2880870/149‑1 (MQ=255)
                                                                                          gATAATGCGCGAATAAAGGTGTATATGCTAACCACTCCGTAGCGGCTTAGAAATTATAAAATCAGGTTTTCACCTGATATCAACCCGATAATTGAATCATTAATCAGGGATGCTCCAGTTTAAACTTCGGGTAGTGCTCGGCAAAATAc                                                       <  2:3309484/149‑1 (MQ=255)
                                                                                           aTCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGATAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACt                                                      >  2:1888528/1‑149 (MQ=255)
                                                                                           aTCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGATAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGAAGTGCTCGGAAAAATACt                                                      >  2:3291527/1‑149 (MQ=255)
                                                                                            tCATGCGCGATTAAATGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGATAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTg                                                     <  2:4011442/149‑1 (MQ=255)
                                                                                               tGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTATAAAGGATAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGAAAAATACTGGCg                                                  <  2:3576075/149‑1 (MQ=255)
                                                                                                  gcgATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGATAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTaa                                               >  1:1748329/1‑149 (MQ=255)
                                                                                                      ttAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGATAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAAt                                           >  1:2819305/1‑149 (MQ=255)
                                                                                                      ttAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGATAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAAt                                           <  2:463380/149‑1 (MQ=255)
                                                                                                        aaaGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGATAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTc                                         <  2:2621287/149‑1 (MQ=255)
                                                                                                          aGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGATAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTa                                       >  1:1784718/1‑149 (MQ=255)
                                                                                                            gTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGATAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTAcc                                     >  1:636341/1‑149 (MQ=255)
                                                                                                            gTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGATAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTAcc                                     <  2:3409390/149‑1 (MQ=255)
                                                                                                              gAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGATAATATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTACAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGt                                   <  1:730815/149‑1 (MQ=255)
                                                                                                               aaTATGCTAACCAATCTGTAGCGGCTTAGAAAGGATAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGtt                                  <  1:2147853/149‑1 (MQ=255)
                                                                                                               aaTATGCTAACCAATCTGTAGCGGCTTAGAAAGGATAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGtt                                  <  1:4043681/149‑1 (MQ=255)
                                                                                                                   tGCTAACCAATCTGTAGCGGCTTAGAAAGGATAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATcc                              >  1:3093668/1‑149 (MQ=255)
                                                                                                                   tGCTAACCAATCTGTAGCGGCTTAGAAAGGATAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATcc                              >  1:669997/1‑149 (MQ=255)
                                                                                                                    gCTAACCAATCTGTAGCGGCTTAGAAAGGATAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCg                             >  2:1864917/1‑149 (MQ=255)
                                                                                                                     cTAACCAATCTGTAGCGGCTTAGAAAGGATAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGc                            >  1:163100/1‑149 (MQ=255)
                                                                                                                     cTAAACAATCTGTAGCGGCTTAGAAAGCATAAAATCAGGATTTAACCTGATATCAACCCGATAAATGAATCATTAATCAGTCATGCTCCAGAGAAAAATACGGGTAGTGATCTGCAAAATACTGGCTTAAAAATTATACCGTTATCCGc                            >  2:137689/1‑149 (MQ=255)
                                                                                                                      tAACCAATCTGTAGCGGCTTAGAAAGGATAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCa                           <  2:614916/149‑1 (MQ=255)
                                                                                                                      tAACCAATCTGTAGCGGCTTAGAAAGGATAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCa                           >  1:892203/1‑149 (MQ=255)
                                                                                                                       aaCCAATCTGTAGCGGCTTAGAAAGGATAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCAc                          >  1:3093986/1‑149 (MQ=255)
                                                                                                                        aCCAATCTGTAGCGGCTTAGAAAGGATAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACt                         >  1:1292587/1‑149 (MQ=255)
                                                                                                                        aCCAATCTGTAGCGGCTTAGAAAGGATAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACt                         >  1:3582284/1‑149 (MQ=255)
                                                                                                                               tGTAGCGGCTTAGAAAGGATAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGCTg                  <  1:1877279/149‑1 (MQ=255)
                                                                                                                               tGTAGCGGCTTAGAAAGGAGAAAATTAGGTTTTAACCCGATATCAACCCGATAATTGAATAATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGCTg                  >  2:2717573/1‑149 (MQ=255)
                                                                                                                                gTAGCGGCTTAGAAAGGATAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGCTGa                 <  2:2401743/149‑1 (MQ=255)
                                                                                                                                  aGCGGCTTAGAAAGGATAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGCTGACg               >  1:1404938/1‑149 (MQ=255)
                                                                                                                                  aGCGGCTTAGAAAGGATAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGCTGACg               >  1:1033586/1‑149 (MQ=255)
                                                                                                                                      gCTTAGAAAGGATAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGCTGACGTCGc           >  2:4052248/1‑149 (MQ=255)
                                                                                                                                      gCTTAGAAAGGATAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGCTGACGTCGc           >  2:85073/1‑149 (MQ=255)
                                                                                                                                       cTTAGAAAGGATAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGCTGACGTCGcc          >  2:2208099/1‑149 (MQ=255)
                                                                                                                                       cTTAGAAAGGATAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGCTGACGTCGcc          <  2:930209/149‑1 (MQ=255)
                                                                                                                                        ttAGAAAGGATAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGCTGACGTCGCCa         >  2:3973999/1‑149 (MQ=255)
                                                                                                                                        ttAGAAAGGATAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGCTGACGTCGCCa         <  1:1977278/149‑1 (MQ=255)
                                                                                                                                           gAAAGGATAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGCTGACGTCGCCAGcc      >  1:3454129/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
CCAACACAATAAAAATGGCCCCAACCAGATAGGTCCAGTAATTCAGAACCCCGTATTCAGCGAACACATTAACCTCTTTAATTATCTTTCGATCATGCGCGATTAAAGGTGAATATGCTAACCAATCTGTAGCGGCTTAGAAAGGAGAAAATCAGGTTTTAACCTGATATCAACCCGATAATTGAATCATTAATCAGGCATGCTCCAGTGAAAAATTCGGGTAGTGCTCGGCAAAATACTGGCGTAAAAATTCTACCGTTATCCGCACTTTCGCTGACGTCGCCAGCCTTGA  >  CP009273/1874950‑1875241

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