Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,583,584 T→A 7.7% I877I (ATT→ATA acrD → aminoglycoside/multidrug efflux system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,583,5840TA 7.7% 995.2 / 55.8 304I877I (ATT→ATAacrDaminoglycoside/multidrug efflux system
Reads supporting (aligned to +/- strand):  ref base T (136/144);  new base A (12/12);  total (148/156)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 315 total aligned reads displayed.

AAGCCGCACCGGGGGTCAGTACCGGTACGGCGATGGATATTATGGAATCGTTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCTGTGTCTGGCTGCGTTGTATGAAAGCTGGTCGGTGCCGTTCTCGGTAATGCTGGTCGTGCCGCTGGGGGTAATCGGCGCGCTGCTGGCAACCTGGATGCGCGGGCTGGAAAACGACGTTTACTTCCAGGT  >  CP009273/2583436‑2583733
                                                                                                                                                    |                                                                                                                                                     
aaGCCGCACCGGGGGTCAGTACCGGTACGGCGATGGATATTATGGAATCGTTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCAttt                                                                                                                                                      <  1:1017149/150‑1 (MQ=255)
 aGCCGCACCGGGGGTCAGTACCGGTACGGCGATGGATATTATGGAATCGTTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTc                                                                                                                                                     <  2:698293/150‑1 (MQ=255)
 aGCCGCACCGGGGGTCAGTACCGGTACGGCGATGGATATTATGGAATCGTTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTc                                                                                                                                                     <  2:538505/150‑1 (MQ=255)
  gCCGCACCGGGGGTCAGTACCGGTACGGCGATGGATATTATGGAATCGTTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTcc                                                                                                                                                    <  2:4032581/150‑1 (MQ=255)
  gCCGCACCGGGGGTCAGTACCGGTACGGCGATGGATATTATGGAATCGTTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTcc                                                                                                                                                    <  2:3846452/150‑1 (MQ=255)
   ccGCACCGGGGGTCAGTACCGGTACGGCGATGGATATTATGGAATCGTTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCt                                                                                                                                                   <  2:5123825/150‑1 (MQ=255)
     gCACCGGGGGTCTGTACCGGTACGGTGTTGTATATTATGGAATCTTTAGTGAAGCAGCTGCCATACGGCTTTGGTCTGGAGTGGACGTCGATGTCTTATCAGGAGCGTCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTg                                                                                                                                                 <  2:2985888/150‑1 (MQ=255)
     gCACCGGGGGTCAGTACCGGTACGGCGATGGATATTATGGAATCGTTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTg                                                                                                                                                 <  2:3749305/150‑1 (MQ=255)
          gggggTCAGTACCGGTACGGCGATGGATATTATGGAATCGTTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCtggt                                                                                                                                            <  1:706793/150‑1 (MQ=255)
          gggggTCAGTACCGGTACGGCGATGGATATTATGGAATCGTTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCtggt                                                                                                                                            >  2:2314557/1‑150 (MQ=255)
          gggggTCAGTACCGGTACGGCGATGGATATTATGGAATCGTTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCtggt                                                                                                                                            >  2:4877183/1‑150 (MQ=255)
           ggggTCAGTACCGGTACGGCGATGGATATTATGGAATCGTTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCtggtg                                                                                                                                           <  1:2353562/150‑1 (MQ=255)
              gTCAGTACCGGTACGGCGATGGATATTATGGAATCGTTAGTGAAGCAGCTTCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCATGAGCGGCTTTCCTGCGCGCAGGCTCCGTCGCTGTACGCCAATTCCTTGCTGGTGGTa                                                                                                                                        >  1:1817213/1‑150 (MQ=255)
              gTCAGTACCGGTACGGCGATGGATATTATGGAATCGTTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTa                                                                                                                                        >  1:2052289/1‑150 (MQ=255)
              gTCAGTACCGGTACGGCGATGGATATTATGGAATCGTTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTa                                                                                                                                        <  1:1807721/150‑1 (MQ=255)
              gTCAGTACCGGTACGGCGATGGATATGATGGACTCGTTAGTGAAGCAGCTGCCACACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTa                                                                                                                                        <  1:1957587/150‑1 (MQ=255)
                cAGTACCGGTACGGCGATGGATATTATGGAATCGTTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTAtt                                                                                                                                      <  1:2887916/150‑1 (MQ=255)
                cAGTACCGGTACGGCGATGGATATTATGGAATCGTTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTAtt                                                                                                                                      <  1:4813856/150‑1 (MQ=255)
                 aGTACCGGTACGGCGATGGATATTATGGAATCGTTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTc                                                                                                                                     <  2:2152113/150‑1 (MQ=255)
                 aGTACCGGTACGGCGATGGATATTATGGAATCGTTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTc                                                                                                                                     >  1:3609120/1‑150 (MQ=255)
                   tACCGGTACGGCGATGGATATTATGGAATCGTTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCt                                                                                                                                   >  2:3954080/1‑150 (MQ=255)
                   tACCGGTACGGCGATGGATATTATGGAATCGTTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCt                                                                                                                                   <  1:3350259/150‑1 (MQ=255)
                    aCCGGTACGGCGATGGATATTATGGAATCGTTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCtg                                                                                                                                  >  2:4514925/1‑150 (MQ=255)
                     ccGGTACGGCGATGGATATTATGGAATCGTTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCtgt                                                                                                                                 <  1:3662040/150‑1 (MQ=255)
                     ccGGTACGGCGATGGATATTATGGAATCGTTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATATCCTTGCTGGTGGTATTCCtgt                                                                                                                                 <  2:4755118/150‑1 (MQ=255)
                      cGGTACGGCGATGGATATTATGGAATCGTTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCtgtg                                                                                                                                >  1:2304219/1‑150 (MQ=255)
                         tACGGCGATGGATATTATGGAATCGTTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCTGTGTCt                                                                                                                             >  2:2686412/1‑150 (MQ=255)
                          aCGGCGATGGATATTATGGAATCGTTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCTGTGTCTg                                                                                                                            >  2:746146/1‑150 (MQ=255)
                             gCGATGGATATTATGGAATCGTTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCTGTGTCTGGCt                                                                                                                         <  2:4907853/150‑1 (MQ=255)
                              cGATGGATATTATGGAATCGTTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCTGTGTCTGGCTg                                                                                                                        <  1:2837020/150‑1 (MQ=255)
                                aTGGATATTATGGAATCGTTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCTGTGTCTGGCTGCg                                                                                                                      >  2:3233045/1‑150 (MQ=255)
                                aTGGATATTATGGAATCGTTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCTGTGTCTGGCTGCg                                                                                                                      >  2:4377586/1‑150 (MQ=255)
                                      attatGGAATCGTTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCTGTGTCTGGCTGCGTTGTAt                                                                                                                >  1:625206/1‑150 (MQ=255)
                                       ttatGGAATCGTTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCTGTGTCTGGCTGCGTTGTATg                                                                                                               <  2:2035880/150‑1 (MQ=255)
                                        tatGGAATCGTTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCTGTGTCTGGCTGCGTTGTATGa                                                                                                              >  1:979256/1‑150 (MQ=255)
                                        tatGGAATCGTTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCTGTGTCTGGCTGCGTTGTATGa                                                                                                              <  1:5080813/150‑1 (MQ=255)
                                          tGGAATCGTTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCTGTGTCTGGCTGCGTTGTATGaaa                                                                                                            >  2:3582590/1‑150 (MQ=255)
                                           ggAATCGTTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCTGTGTCTGGCTGCGTTGTATGAAAg                                                                                                           <  1:1496577/150‑1 (MQ=255)
                                                cGTTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCTGTGTCTGGCTGCGTTGTATGAAAGCTGGt                                                                                                      <  1:3842690/150‑1 (MQ=255)
                                                 gTTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCTGTGTCTGGCTGCGTTGTATGAAAGCTGGTc                                                                                                     >  2:1274438/1‑150 (MQ=255)
                                                 gTTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCTGTGTCTGGCTGCGTTGTATGAAAGCTGGTc                                                                                                     <  1:882824/150‑1 (MQ=255)
                                                  ttAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCTGTGTCTGGCTGCGTTGTATGAAAGCTGGTCg                                                                                                    <  2:809996/150‑1 (MQ=255)
                                                  ttAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCGGTGTCTGGCTGCGTTGTATGAAAGCTGGTCg                                                                                                    >  2:3882684/1‑150 (MQ=255)
                                                   tAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCTGTGTCTGGCTGCGTTGTATGAAAGCTGGTCgg                                                                                                   <  1:4829788/150‑1 (MQ=255)
                                                   tAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCTGTGTCTGGCTGCGTTGTATGAAAGCTGGTCgg                                                                                                   >  2:1705573/1‑150 (MQ=255)
                                                    aGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCTGTGTCTGGCTGCGTTGTATGAAAGCTGGTCGGt                                                                                                  <  2:1622077/150‑1 (MQ=255)
                                                           cagcTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCTGTGTCTGGCTGCGTTGTATGAAAGCTGGTCGGTGCCGTtc                                                                                           >  1:1192855/1‑150 (MQ=255)
                                                            agcTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCTGTGTCTGGCTGCGTTGTATGAAAGCTGGTCGGTGCCGTtct                                                                                          <  1:1363883/150‑1 (MQ=255)
                                                            agcTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCTGTGTCTGGCTGCGTTGTATGAAAGCTGGTCGGTGCCGTtct                                                                                          <  2:3343389/150‑1 (MQ=255)
                                                             gcTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCTGTGTCTGGCTGCGTTGTATGAAAGCTGGTCGGTGCCGTtctc                                                                                         <  1:342481/150‑1 (MQ=255)
                                                             gcTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCTGTGTCTGGCTGCGTTGTATGAAAGCTGGTCGGTGCCGTtctc                                                                                         <  2:3776870/150‑1 (MQ=255)
                                                              cTGCCACACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCTGTGTCTGGCTGCGTTGTATGAAAGCTGGTCGGTGCCGTTCTCg                                                                                        <  2:2383154/150‑1 (MQ=255)
                                                               tGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCTGTGTCTGGCTGCGTTGTATGAAAGCTGGTCGGTGCCGTTCTCgg                                                                                       <  2:3088741/150‑1 (MQ=255)
                                                                   aaaCGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCTGTGTCTGGCTGCGTTGTATGAAAGCTGGTCGGTGCCGTTCTCGGTAAt                                                                                   <  1:2555284/150‑1 (MQ=255)
                                                                    aaCGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCTGTGTCTGGCTGCGTTGTATGAAAGCTGGTCGGTGCCGTTCTCGGTAATg                                                                                  >  1:4267748/1‑150 (MQ=255)
                                                                       ggCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCTGTGTCTGGCTGCGTTGTATGAAAGCTGGTCGGTGCCGTTCTCGGTAATGCTg                                                                               <  1:1524349/150‑1 (MQ=255)
                                                                        gCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCTGTGTCTGGCTGCGTTGTATGAAAGCTGGTCGGTGCCGTTCTCGGTAATGCTgg                                                                              <  2:1754785/150‑1 (MQ=255)
                                                                        gCTGTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCTGTGTCTGGCTGCGTTGTATGAAAGCTGGTCGGTGCCGTTCTCGGTAATGCTgg                                                                              >  1:2988537/1‑150 (MQ=255)
                                                                         cTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCTGTGTCTGGCTGCGTTGTATGAAAGCTGGTCGGTGCCGTTCTCGGTAATGCTGGt                                                                             >  1:2311504/1‑150 (MQ=255)
                                                                          tttGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCTGTGTCTGGCTGCGTTGTATGAAAGCTGGTCGGTGCCGTTCTCGGTAATGCTGGTc                                                                            <  2:3507090/150‑1 (MQ=255)
                                                                          tttGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCTGTGTCTGGCTGCGTTGTATGAAAGCTGGTCGGTGCCGTTCTCGGTAATGCTGGTc                                                                            <  1:3271615/150‑1 (MQ=255)
                                                                           ttGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCTGTGTCTGGCTGCGTTGTATGAAAGCTGGTCGGTGCCGTTCTCGGTAATGCTGGTCg                                                                           >  1:3906415/1‑150 (MQ=255)
                                                                                       ggACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCTGTGTCTGGCTGCGTTGTATGAAAGCTGGTCGGTGCCGTTCTCGGTCATGCTGGTCGTGCCGCTggggg                                                               >  2:4411786/1‑150 (MQ=255)
                                                                                         aCGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCTGTGTCTGGCTGCGTTGTATGAAAGCTGGTCGGTGCCGTTCTCGGTAATGCTGGTCGTGCCGCTGGGGGTa                                                             <  2:4978223/150‑1 (MQ=255)
                                                                                          cGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCTGTGTCTGGCTGCGTTGTATGAAAGCTGGTCGGTGCCGTTCTCGGTAATGCTGGTCGTGCCGCTGGGGGTaa                                                            >  2:729440/1‑150 (MQ=255)
                                                                                           ggCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCTGTGTCTGGCTGCGTTGTATGAAAGCTGGTCGGTGCCGTTCTCGGTAATGCTGGTCGTGCCGCTGGGGGTAAt                                                           >  2:188408/1‑150 (MQ=255)
                                                                                            gCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCTGTGTCTGGCTGCGTTGTATGAAAGCTGGTCGGTGCCGTTCTCGGTAATGCTGGTCGTGCCGCTGGGGGTAATc                                                          <  1:3042176/150‑1 (MQ=255)
                                                                                            gCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCTGTGTCTGGCTGCGTTGTATGAAAGCTGGTCGGTGCCGTTCTCGGTAATGCTGGTCGTGCCGCTGGGGGTAATc                                                          <  1:593667/150‑1 (MQ=255)
                                                                                             cGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCTGTGTCTGGCTGCGTTGTATGAAAGCTGGTCGGTGCCGTTCTCGGTAATGCTGGTCGTGCCGCTGGGGGTAATCg                                                         >  2:3276932/1‑150 (MQ=255)
                                                                                               aTGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCTGTGTCTGGCTGCGTTGTATGAAAGCTGGTCGGTGCCGTTCTCGGTAATGCTGGTCGTGCCGCTGGGGGTAATCGgc                                                       <  1:5005523/150‑1 (MQ=255)
                                                                                               aTGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGAGCTGTACGCCATTTCCTTGCTGGTGGTATTCCTGTGTCTGGCTGCTTTGTATGAAAGCTGGTCTGTGCCGTTCTCGGTAATGCTGGTCGTGCCGCTGGGTGTAATCGgc                                                       <  1:1446149/150‑1 (MQ=255)
                                                                                                 gTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATATCCTTGCTGGTGGTATTCCTGTGTCTGGCTGCGTTGTATGAAAGCTGGTCGGTGCCGTTCTCGGTAATGCTGGTCGTGCCGCTGGGGGTAATCGgcgc                                                     >  2:2031761/1‑150 (MQ=255)
                                                                                                   cGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCTGTGTCTGGCTGCGTTGTATGAAAGCTGGTCGGTGCCGTTCTCGGTAATGCTGGTCGTGCCGCTGGGGGTCATCGgcgcgc                                                   >  2:4573725/1‑150 (MQ=255)
                                                                                                     tATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCTGTGTCTGGCTGCGTTGTATGAAAGCTGGTCGGTGCCGTTCTCGGTAATGCTGGTCGTGCCGCTGGGGGTAATCGGCGCgctg                                                 <  1:1396574/150‑1 (MQ=255)
                                                                                                        cAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCTGTGTCTGGCTGCGTTGTATGAAAGCTGGTCGGTGCCGTTCTCGGTAATGCTGGTCGTGCCGCTGGGGGTAATCGGCGCGctgctg                                              >  1:2632629/1‑150 (MQ=255)
                                                                                                             gCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCTGTGTCTGGCTGCGTTGTATGAAAGCTGGTCGGTGCCGTTCTCGGTAATGCTGGTCGTGCCGCTGGGGGTAATCGGCGCGCTGCTGGCAAc                                         <  1:2933228/150‑1 (MQ=255)
                                                                                                             gCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCTGTGTCTGGCTGCGTTGTATGAAAGCTGGTCGGTGCCGTTCTCGGTAATGCTGGTCGTGCCGCTGGGGGTAATCGGCGCGCTGCTGGCAAc                                         <  1:4053243/150‑1 (MQ=255)
                                                                                                               ggCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTCCCTTGCTGGTGGTATTCCTGTGTCTGGCTGCGTTGTATGAAAGCTGGTCGGTGCCGTTCTAGGTAATGCTGGTCGTGCCGCTGGGGGTAATCGGCGCGCTGCTGGCAACCt                                       <  1:4002307/150‑1 (MQ=255)
                                                                                                                   ttCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCTGTGTCTGGCTGCGTTGTATGAAAGCTGGTCGGTGCCGTTCTCGGTAATGCTGGTCGTGCCGCTGGGGGTAATCGGCGCGCTGCTGGCAACCTGGAt                                   <  2:425736/150‑1 (MQ=255)
                                                                                                                    tCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCTGTGTCTGGCTGCGTTGTATGAAAGCTGGTCGGTGCCGTTCTCGGTAATGCTGGTCGTGCCGCTGGGGGTAATCGGCGCGCTGCTGGCAACCTGGATg                                  >  2:2482861/1‑150 (MQ=255)
                                                                                                                      cGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCTGTGTCTGGCTGCGTTGTATGAAAGCTGGTCGGTGCCGTTCTCGGTAATGCTGGTCGTGCCGCTGGGGGTAATCGGCGCGCTGCTGGCAACCTGGATgcg                                <  1:2169168/150‑1 (MQ=255)
                                                                                                                       ggCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCTGTGTCTGGCTGCGTTGTATGAAAGCTGGTCGGTGCCGTTCTCGGTAATGCTGGTCGTGCCGCTGGGGGTAATCGGCGCGCTGCTGGCAACCTGGATgcgc                               >  1:2260625/1‑150 (MQ=255)
                                                                                                                         cgcgcAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCTGTGTCTGGCTGCGTTGTATGAAAGCTGGTCGGTGCCGTTCTCGGTAATGCTGGTCGTGCCGCTGGGGGTAATCGGCGCGCTGCTGGCAACCTGGATGCGCgg                             <  1:3545800/150‑1 (MQ=255)
                                                                                                                         cgcgcAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCTGTGTCTGGCTGCGTTGTATGAAAGCTGGTCGGTGCCGTTCTCGGTAATGCTGGTCGTGCCGCTGGGGGTAATCGGCGCGCTGCTGGCAACCTGGATGCGCgg                             >  1:5132203/1‑150 (MQ=255)
                                                                                                                          gcgcAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCTGTGTCTGGCTGCGTTGTATGAAAGCTGGTCGGTGCCGTTCTCGGTAATGCTGGTCGTGCCGCTGGGGGTAATCGGCGCGCTGCTGGCAACCTGGATGCGCggg                            <  2:168791/150‑1 (MQ=255)
                                                                                                                           cgcAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCTGTGTCTGGCTGCGTTGTATGAAAGCTGGTCGGTGCCGTTCTCGGTAATGCTGGTCGTGCCGCTGGGGGTAATCGGCGCGCTGCTGGCAACCTGGATGCGCGGGc                           <  1:4434487/150‑1 (MQ=255)
                                                                                                                            gcAGGCTCCGGCGCTGTACGCCATATCCTTGCTGGTGGTATTCCTGTGTCTGGCTGCGTTGTATGAAAGCTGGTCGGTGCCGTTCTCGGTAATGCTGGTCGTGCCGCTGGGGGTAATCGGCGCGCTGCTGGCAACCTGGATGCGCGGGCt                          <  2:371434/150‑1 (MQ=255)
                                                                                                                               ggCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCTGTGTCTGGCTGCGTTGTATGAAAGCTGGTCGGTGCCGTTCTCGGTAATGCTGGTCGTGCCGCTGGGGGTAATCGGCGCGCTGCTGGCAACCTGGATGCGCGGGCTGGa                       <  1:1753797/150‑1 (MQ=255)
                                                                                                                                gCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCTGTGTCTGGCTGCGTTGTATGAAAGCTGGTCGGTGCCGTTCTCGGTAATGCTGGTCGTGCCGCTGGGGGTAATCGGCGCGCTGCTGGCAACCTGGATGCGCGGGCTGGaa                      >  1:1083253/1‑150 (MQ=255)
                                                                                                                                   ccGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCTGTGTCTGGCTGCGTTGTATGAAAGCTGGTCGGTGCCGTTCTCGGTAATGCTGGTCGTGCCGCTGGGGGTAATCGGCGCGCTGCTGGCAACCTGGATGCGCGGGCTGGAAAac                   >  2:354729/1‑150 (MQ=255)
                                                                                                                                     ggCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCTGTGTCTGGCTGCGTTGTATGAAAGCTGGTCGGTGCCGTTCTCGGTAATGCTGGTCGTGCCGCTGGGGGTAATCGGCGCGCTGCTGGCAACCTGGATGCGCGGGCTGGAAAacga                 >  2:249964/1‑150 (MQ=255)
                                                                                                                                       cgcTGTACGCCATATCCTTGCTGGTGGTATTCCTGTGTCTGGCTGCGTTGTATGAAAGCTGGTCGGTGCCGTTCTCGGTAATGCTGGTCGTGCCGCTGGGGGTAATCGGCGCGCTGCTGGCAACCTGGATGCGCGGGCTGGAAAacgacg               <  1:3783180/150‑1 (MQ=255)
                                                                                                                                           gTACGCCATATCCTTGCTGGTGGTATTCCTGTGTCTGGCTGCGTTGTATGAAAGCTGGTCGGTGCCGTTCTCGGTAATGCTGGTCGTGCCGCTGGGGGTAATCGGCGCGCTGCTGGCAACCTGGATGCGCGGGCTGGAAAACGACGTTTa           >  1:4601299/1‑150 (MQ=255)
                                                                                                                                              cGCCATTTCCTTGCTGGTGGTATTCCTGTGTCTGGCTGCGTTGTATGAAAGCTGGTCGGTGCCGTTCTCGGTAATGCTGGTCGTGCCGCTGGGGGTAATCGGCGCGCTGCTGGCAACCTGGATGCGCGGGCTGGAAAACGACGTTTACtt        <  1:4687113/150‑1 (MQ=255)
                                                                                                                                               gCCATTTCCTTGCTGGTGGTATTCCTGTGTCTGGCTGCGTTGTATGAAAGCTGGTCGGTGCCGTTCTCGGTAATGCTGGTCGTGCCGCTGGGGGTAATCGGCGCGCTGCTGGCAACCTGGATGCGCGGGCTGGAAAACGACGTTTACTTc       <  2:1503592/150‑1 (MQ=255)
                                                                                                                                                ccAATTCCTTGCTGGTGGTATTCCTGTGTCTGGCTGCGTTGTATGAAAGCTGGTCGGTGCCGTTCTCGGTAATGCTGGTCGTGCCGCTGGGGGTAATCGGCGCGCTGCTGGCAACCTGGATGCGCGGGCTGGAAAACGACGTTTACTTcc      >  2:4643681/1‑150 (MQ=255)
                                                                                                                                                 cATTTCCTTGCTGGTGGTATTCCTGTGTCTGGCTGCGTTGTATGAAAGCTGGTCGGTGCCGTTCTCGGTAATGCTGGTCGTGCCGCTGGGGGTAATCGGCGCGCTGCTGGCAACCTGGATGCGCGGGCTGGAAAACGACGTTTACTTCCa     <  2:1133671/150‑1 (MQ=255)
                                                                                                                                                 cATTTCCTTGCTGGTGGTATTCCTGTGTCTGGCTGCGTTGTATGAAAGCTGGTCGGTGCCGTTCTCGGTAATGCTGGTCGTGCCGCTGGGGGTAATCGGCGCGCTGCTGGCAACCTGGATGCGCGGGCTGGAAAACGACGTTTACTTCCa     <  2:2894941/150‑1 (MQ=255)
                                                                                                                                                   tttCCTTGCTGGTGGTATTCCTGTGTCTGGCTGCGTTGTATGAAAGCTGGTCGGTGCCGTTCTCGGTAATGCTGGTCGTGCCGCTGGGGGTAATCGGCGCGCTGCTGGCAACCTGGATGCGCGGGCTGGAAAACGACGTTTACTTCCAgg   >  1:4186493/1‑150 (MQ=255)
                                                                                                                                                    ttCCTTGCTGGTGGTATTCCTGTGTCTGGCTGCGTTGTATGAAAGCTGGTCGGTGCCGTTCTCGGTAATGCTGGTCGTGCCGCTGGGGGTAATCGGCGCGCTGCTGGCAACCTGGATGCGCGGGCTGGAAAACGACGTTTACTTCCAGGt  <  2:2540452/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
AAGCCGCACCGGGGGTCAGTACCGGTACGGCGATGGATATTATGGAATCGTTAGTGAAGCAGCTGCCAAACGGCTTTGGTCTGGAGTGGACGGCGATGTCGTATCAGGAGCGGCTTTCCGGCGCGCAGGCTCCGGCGCTGTACGCCATTTCCTTGCTGGTGGTATTCCTGTGTCTGGCTGCGTTGTATGAAAGCTGGTCGGTGCCGTTCTCGGTAATGCTGGTCGTGCCGCTGGGGGTAATCGGCGCGCTGCTGGCAACCTGGATGCGCGGGCTGGAAAACGACGTTTACTTCCAGGT  >  CP009273/2583436‑2583733

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