Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 724,813 A→T 100% intergenic (+17/‑226) ybfA → / → rhsC DUF2517 family protein/Rhs family putative polymorphic toxin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273724,8130AT98.9% 607.6 / ‑2.5 185intergenic (+17/‑226)ybfA/rhsCDUF2517 family protein/Rhs family putative polymorphic toxin
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base T (92/91);  total (93/92)
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.83e-01
Rejected as polymorphism: E-value score below prediction cutoff.
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 185 total aligned reads displayed.

CGCTGCCTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATAT  >  CP009273/724666‑724959
                                                                                                                                                   |                                                                                                                                                  
  cTGCCTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCAtt                                                                                                                                                >  1:4810025/1‑150 (MQ=255)
   tGCCTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTg                                                                                                                                               >  2:2524799/1‑150 (MQ=255)
   tGCCTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTg                                                                                                                                               <  2:451268/150‑1 (MQ=255)
    gCCTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGa                                                                                                                                              <  2:3145635/150‑1 (MQ=255)
        ttCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAAc                                                                                                                                          >  1:2771429/1‑150 (MQ=255)
        ttCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAAc                                                                                                                                          <  2:4248261/150‑1 (MQ=255)
          cATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATACACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCg                                                                                                                                        >  2:4997258/1‑150 (MQ=255)
          cATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCg                                                                                                                                        >  1:929417/1‑150 (MQ=255)
          cATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCg                                                                                                                                        >  2:3199368/1‑150 (MQ=255)
          cATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCg                                                                                                                                        <  1:274262/150‑1 (MQ=255)
          cATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCg                                                                                                                                        <  1:274258/150‑1 (MQ=255)
            tGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGcc                                                                                                                                      >  2:706257/1‑150 (MQ=255)
             gCTTTTCAGGACAGACCGCGCCTGCTTTTAAAGCTACCTGCATCGCGTCTTGTCGAAAACCAGCGATAAACCTGTGTGGATGGATGAGGACGAAAAGGCAACTGGTTATTTTTATTGATTGCTGAGTATTGAATTCATTGAAAAACTCCa                                                                                                                                     <  2:4389042/150‑1 (MQ=255)
               ttttCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGc                                                                                                                                   <  2:3697924/150‑1 (MQ=255)
               ttttCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGc                                                                                                                                   <  2:1948977/150‑1 (MQ=255)
               ttttCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGc                                                                                                                                   <  2:3398412/150‑1 (MQ=255)
                  tCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGa                                                                                                                                >  1:1310359/1‑150 (MQ=255)
                   cTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGaa                                                                                                                               >  1:4889581/1‑150 (MQ=255)
                   cTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGAGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGACGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGaa                                                                                                                               <  2:4689561/150‑1 (MQ=255)
                    tGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGaaa                                                                                                                              >  2:1122197/1‑150 (MQ=255)
                    tGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGaaa                                                                                                                              <  2:1612158/150‑1 (MQ=255)
                    tGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGaaa                                                                                                                              <  2:1612193/150‑1 (MQ=255)
                     ggAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAAc                                                                                                                             <  1:589540/150‑1 (MQ=255)
                      gAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACg                                                                                                                            <  1:77652/150‑1 (MQ=255)
                        aaGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCt                                                                                                                          <  1:386332/150‑1 (MQ=255)
                         aGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTg                                                                                                                         <  1:2264452/150‑1 (MQ=255)
                           aCCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGc                                                                                                                       <  1:4984113/150‑1 (MQ=255)
                               cgcCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGtt                                                                                                                   <  2:2870907/150‑1 (MQ=255)
                                gcCCGCTTTTACAGCTACCTGCATCGCGTCGGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCCTTGAAAACCGCCAGCGGAAACGCTGGCGGttt                                                                                                                  <  2:3772152/150‑1 (MQ=255)
                                  ccGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGttttt                                                                                                                <  1:84630/150‑1 (MQ=255)
                                  ccGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGttttt                                                                                                                <  1:369129/150‑1 (MQ=255)
                                   cGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGtttttt                                                                                                               >  1:3750554/1‑150 (MQ=255)
                                      ttttACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTa                                                                                                            >  1:3261854/1‑150 (MQ=255)
                                        ttACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATc                                                                                                          >  2:1410297/1‑150 (MQ=255)
                                         tACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCa                                                                                                         >  1:382483/1‑150 (MQ=255)
                                          aCAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATACACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAg                                                                                                        <  2:4337327/150‑1 (MQ=255)
                                           cAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAgg                                                                                                       <  2:1046248/150‑1 (MQ=255)
                                                  cTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGt                                                                                                <  2:2771429/150‑1 (MQ=255)
                                                       tCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCaga                                                                                           <  1:3639008/150‑1 (MQ=255)
                                                        cgcgTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCagag                                                                                          <  2:4319859/150‑1 (MQ=255)
                                                         gcgTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCagaga                                                                                         <  1:627390/150‑2 (MQ=255)
                                                         gcgTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAgg                                                                                         <  1:2356960/150‑1 (MQ=255)
                                                           gTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGaa                                                                                       <  1:3401772/150‑1 (MQ=255)
                                                            tCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTCTCAGGAAGATGTTCAGAGGact                                                                                      >  2:2094140/1‑148 (MQ=255)
                                                              tGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAatat                                                                                    >  1:1245573/1‑150 (MQ=255)
                                                                  cGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAg                                                                                >  2:4112096/1‑150 (MQ=255)
                                                                    aaaaCCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGAc                                                                              <  2:1215442/150‑1 (MQ=255)
                                                                     aaaCCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACa                                                                             >  1:3081549/1‑150 (MQ=255)
                                                                     aaaCCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACa                                                                             >  2:3194546/1‑150 (MQ=255)
                                                                       aCCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATc                                                                           <  2:2144232/150‑1 (MQ=255)
                                                                        ccAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCg                                                                          <  1:3926569/150‑1 (MQ=255)
                                                                        ccAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCg                                                                          <  2:1536059/150‑1 (MQ=255)
                                                                          aGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGtt                                                                        >  2:787919/1‑150 (MQ=255)
                                                                               tAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATAcc                                                                   >  1:766724/1‑150 (MQ=255)
                                                                                aaaCCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCa                                                                  >  1:1038505/1‑150 (MQ=255)
                                                                                 aaCCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGTTGATTTTTATTGATTCCTGAGTTTTGAATTCATTGAAAACCGCCAGATGAAACGCTGGCGTTTTTTTGTATCAGGAAGATGTTCAGAGGAATCTTTAGAAATCTTTAAACCaa                                                                 >  2:4827898/1‑150 (MQ=255)
                                                                                  aCCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAAt                                                                <  1:1011960/150‑1 (MQ=255)
                                                                                   ccGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATc                                                               <  1:4769350/150‑1 (MQ=255)
                                                                                    cGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCg                                                              <  2:3902609/150‑1 (MQ=255)
                                                                                     ggTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGa                                                             >  1:1471926/1‑150 (MQ=255)
                                                                                          ggatggatCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACat                                                        <  2:2345808/150‑1 (MQ=255)
                                                                                           gatggatCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACata                                                       <  2:237038/150‑1 (MQ=255)
                                                                                             tggatCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACatata                                                     >  2:3856872/1‑150 (MQ=255)
                                                                                               gatCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATAtt                                                   <  1:3763381/150‑1 (MQ=255)
                                                                                               gatCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATAtt                                                   >  2:5015041/1‑150 (MQ=255)
                                                                                                atCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAATATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGACACATATATTc                                                  >  2:4767176/1‑150 (MQ=255)
                                                                                                atCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATGCCAATCGAAACATATATTc                                                  <  2:668371/150‑1 (MQ=255)
                                                                                                atCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTc                                                  <  2:3476849/150‑1 (MQ=255)
                                                                                                 tCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCa                                                 <  2:2457669/150‑1 (MQ=255)
                                                                                                  cAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCAt                                                >  1:958502/1‑150 (MQ=255)
                                                                                                  cAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCAt                                                >  2:514993/1‑150 (MQ=255)
                                                                                                  cAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGGTATACCAATCGAAACATATATTCAt                                                >  1:3895870/1‑150 (MQ=255)
                                                                                                   aGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTTTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATg                                               >  1:490988/1‑150 (MQ=255)
                                                                                                    ggCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGa                                              <  1:3440829/150‑1 (MQ=255)
                                                                                                    ggCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATAGTTATACCAATCGAAACATATATTCATGa                                              >  2:2270297/1‑150 (MQ=255)
                                                                                                     gCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGaa                                             <  1:4691839/150‑1 (MQ=255)
                                                                                                       ccaaaaGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAat                                           >  1:702080/3‑150 (MQ=255)
                                                                                                        gAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCAGGCGGTTTTTTTTATAATGAATATTTTAAAATGAATATTTATACAAAGTTATAACAAAAGAAACATATATTAATAAAata                                          >  1:450430/1‑150 (MQ=255)
                                                                                                          aaaGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAatata                                        >  1:418947/1‑150 (MQ=255)
                                                                                                              gCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAAtatatata                                    >  1:5058020/1‑150 (MQ=255)
                                                                                                                 aCTGGTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCCCCAGCGGAAACGCTGGCGGTTTTTTGTATCATGAAGATGTTAAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTAATGAAATATATATAAat                                 >  1:189802/1‑150 (MQ=255)
                                                                                                                    ggTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATAtt                              >  1:4298446/1‑150 (MQ=255)
                                                                                                                    ggTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATAtt                              >  1:3059379/1‑150 (MQ=255)
                                                                                                                    ggTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATAAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATAACAAACGAAACATATATTAATGAAATATATATAAATAtt                              >  1:1206441/1‑150 (MQ=255)
                                                                                                                     gTGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATAttt                             >  2:427562/1‑150 (MQ=255)
                                                                                                                      tGATTTTTATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATAtttt                            >  1:1435606/1‑150 (MQ=255)
                                                                                                                         tttttATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCt                         <  2:760276/150‑1 (MQ=255)
                                                                                                                         tttttATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCt                         <  1:884710/150‑1 (MQ=255)
                                                                                                                          ttttATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGACATATATATAAATATTTTCCTa                        >  2:3415047/1‑150 (MQ=255)
                                                                                                                          ttttATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTa                        <  1:4085573/150‑1 (MQ=255)
                                                                                                                             tATTGATTGCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAAtt                     >  1:4331973/1‑150 (MQ=255)
                                                                                                                              attgattgCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTg                    >  1:3485560/1‑150 (MQ=255)
                                                                                                                                  attgCTGAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCt                >  2:3603187/1‑150 (MQ=255)
                                                                                                                                   ttgCTGAGTTTTGAATTCATTGAAAACCGCCAGTGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACAACGTTATAACAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTTTTAtt               >  1:1105530/1‑150 (MQ=255)
                                                                                                                                        gAGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTg          >  1:1587965/1‑150 (MQ=255)
                                                                                                                                         aGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGa         <  1:665148/150‑1 (MQ=255)
                                                                                                                                         aGTTTTGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGa         <  2:1336476/150‑1 (MQ=255)
                                                                                                                                           ttttGAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACa       <  1:4022714/150‑1 (MQ=255)
                                                                                                                                               gAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGata   <  1:712954/150‑1 (MQ=255)
                                                                                                                                               gAATTCATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATAAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGata   >  1:1548966/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CGCTGCCTTTCATGCTTTTCTGGAAAGACCGCGCCCGCTTTTACAGCTACCTGCATCGCGTCTGGTCGAAAACCAGCGATAAACCGGTGTGGATGGATCAGGCCGAAAAGGCAACTGGTGATTTTTATTGATTGCTGAGTTTTGAATACATTGAAAACCGCCAGCGGAAACGCTGGCGGTTTTTTGTATCAGGAAGATGTTCAGAGGAATATTTAGACATCGTTATACCAATCGAAACATATATTCATGAAATATATATAAATATTTTCCTAATTGTTCTTATCTGACAGATAT  >  CP009273/724666‑724959

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: