Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 559,245 G→T 20.5% G154C (GGT→TGT)  sfmF → FimA homolog

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273559,2450GT20.5% 801.0 / 130.5 303G154C (GGT→TGT) sfmFFimA homolog
Reads supporting (aligned to +/- strand):  ref base G (107/133);  new base T (28/35);  total (135/168)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 307 total aligned reads displayed.

ATCCCGCAATGGCACAAACCGTCGCCATCGAATTACGTAATAGCGATCGCTCCCGGCTCGCACTGGGGGAGGCGAGCCCGACTGAGGAAGTAGATGCAAATGGCAATGTCACACTAAACTTTTTTGCCAATTATCGAGCGTTAGCCAGCGGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATATTTATGATCAATTATAATTAATATTATATTAATTCGTATAATTTGGCGTAGTCGATAAGCTCTACAATTGAATGCAAACCTAGCTTGCCATAAATATTAGATTTATGCGCACTA  >  CP009273/559096‑559394
                                                                                                                                                     |                                                                                                                                                     
    cGCAATGGCACAAACCGTCGCCATCGAATTACGTAATAGCGATCGCTCCCGGCTCGCACTGGGGGAGGCGAGCCCGACTGAGGAAGTAGATGCAAATGGCAATGTCACACTAAACTTTTTTGCCAATTATCGAGCGTTAGCCAGCGgtgt                                                                                                                                                   >  2:3428696/1‑150 (MQ=255)
    cGCAATGGCACAAACCGTCGCCATCGAATTACGTAATAGCGATCGCTCCCGGCTCGCACTGGGGGAGGCGAGCCCGACTGAGGAAGTAGATGCAAATGGCAATGTCACACTAAACTTTTTTGCCAATTATCGAGCGTTAGCCAGCGgtgt                                                                                                                                                   <  1:2629002/150‑1 (MQ=255)
      cAATGGCACAACCCGTCGCCATCGAATTACGTAATAGCGATCGCTCCCGGCTCGCACTGGGGGAGGCGAGCCCGACTGAGGAAGTAGATGCAAATGGCAATGTCACACTAAACTTTTTTGCCAATTATCGAGCGTTAGCCAGCGGTGTTc                                                                                                                                                 <  2:688686/150‑1 (MQ=255)
      cAATGGCACAAACCGTCGCCATCGAATTACGTAATAGCGATCGCTCCCGGCTCGCACTGGGGGAGGCGAGCCCGACTGAGGAAGTAGATGCAAATGGCAATGTCACACTAAACTTTTTTGCCAATTATCGAGCGTTAGCCAGCTGTGTTc                                                                                                                                                 >  1:3178856/1‑150 (MQ=255)
      cAATGGCACAAACCGTCGCCATCGAATTACGTAATAGCGATCGCTCCCGGCTCGCACTGGGGGAGGCGAGCCCGACTGAGGAAGTAGATGCAAATGGCAATGTCACACTAAACTTTTTTGCCAATTATCGAGCGTTAGCCAGCGGTGTTc                                                                                                                                                 <  1:3125317/150‑1 (MQ=255)
      cAATGGCACAAACCGTCGCCATCGAATTACGTAATAGCGATCGCTCCCGGCTCGCACTGGGGGAGGCGAGCCCGACTGAGGAAGTAGATGCAAATGGCAATGTCACACTAAACTTTTTTGCCAATTATCGAGCGTTAGCCAGCGGTGTTc                                                                                                                                                 >  2:4958899/1‑150 (MQ=255)
       aaTGGCACAAACCGTCGCCATCGAATTACGTAATAGCGATCGCTCCCGGCTCGCACGGGGGGAGGCGAGCCCGACTGAGGAAGTCGATGCCAATGGCAATGTCACACTACACTTTTTTGCCAATTATCGAGCGTTAGCCAGCGGTGTTCg                                                                                                                                                <  2:3528866/150‑1 (MQ=255)
        aTGGCACAAACCGTCGCCATCGAATTACGTAATAGCGATCGCTCCCGGCTCGCACTGGGGGAGGCGAGCCCGACTGAGGAAGTAGATGCAAATGGCAATGTCACACTAAACTTTTTTGCCAATTATCGAGCGTTAGCCAGCGGTGTTCgg                                                                                                                                               >  1:3085024/1‑150 (MQ=255)
         tGGCACAAACCGTCGCCATCGAATTACGTAATAGCGATCGCTCCCGGCTCGCACTGGGGGAGGCGAGCCCGACTGAGGAAGTAGATGCAAATGGCAATGTCACACTAAACTTTTTTGCCAAGTATCGAGCGTTCGCCAGCTGTGTTCGGc                                                                                                                                              >  2:4443364/1‑150 (MQ=255)
          ggCACAAACCGTCGCCATCGAATTACGTAATAGCGATCGCTCCCGGCTCGCACTGGGGGAGGCGAGCCCGACTGAGGAAGTAGATGCAAATGGCAATGTCACACTAAACTTTTTTGCCAATTATCGAGCGTTAGCCAGCGGTGTTCGGcc                                                                                                                                             >  2:2698750/1‑150 (MQ=255)
               aaaCCGTCGCCATCGAATTACGTAATAGCGATCGCTCCCGGCTCGCACTGGGGGAGGCGAGCCCGACTGAGGAAGTAGATGCAAATGGCAATGTCACACTAAACTTTTTTGCCAATTATCGAGCGTTAGCCAGCGGTGTTCGGCCAGgtg                                                                                                                                        >  1:3972954/1‑150 (MQ=255)
                   cGTCGCCATCGAATTACGTAATAGCGATCGCTCCCGGCTCGCACTGGGGGAGGCGAGCCCGACTGAGGAAGTAGATGCAAATGGCAATGTCACACTAAACTTTTTTGCCAATTATCGAGCGTTAGCCAGCGGTGTTCGGCCAGGTGTGGc                                                                                                                                    <  1:1205217/150‑1 (MQ=255)
                     tCGCCATCGAATTACGTAATAGCGATCGCTCCCGGCTCGCACTGGGGGAGGCGAGCCCGACTGAGGAAGTAGATGCAAATGGCAATGTCACACTAAACTTTTTTGCCAATTATCGAGCGTTAGCCAGCGGTGTTCGGCCAGGTGTGGCGa                                                                                                                                  <  1:1473780/150‑1 (MQ=255)
                       gCCATCGAATTACGTAATAGCGATCGCTCCCGGCTCGCACTGGGGGAGGCGAGCCCGACTGAGGAAGTAGATGCAAATGGCAATGTCACACTAAACTTTTTTGCCAATTATCGAGCGTTAGCCAGCGGTGTTCGGCCAGGTGTGGCGaaa                                                                                                                                <  2:2612495/150‑1 (MQ=255)
                        ccATCGAATTACGTAATAGCGATCGCTCCCGGCTCGCACTGGGGGAGGCGAGCCCGACTGAGGAAGTAGATGCAAATGGCAATGTCACACTAAACTTTTTTGCCAATTATCGAGCGTTAGCCAGCGGTGTTCGGCCAGGTGTGGCGAAAg                                                                                                                               >  2:2148547/1‑150 (MQ=255)
                        ccATCGAATTACGTAATAGCGATCGCTCCCGGCTCGCACTGGGGGAGGCGAGCCCGACTGAGGAAGTAGATGCAAATGGCAATGTCACACTAAACTTTTTTGCCAATTATCGAGCGTTAGCCAGCGGTGTTCGGCCAGGTGTGGCGAAAg                                                                                                                               >  1:4771345/1‑150 (MQ=255)
                           tCGAATTACGTAATAGCGATCGCTCCCGGCTCGCACTGGGGGAGGCGAGCCCGACTGAGGAAGTAGATGCAAATGGCAATGTCACACTAAACTTTTTTGCCAATTATCGAGCGTTAGCCAGCGGTGTTCGGCCAGGTGTGGCGAAAGCgg                                                                                                                            <  1:4194760/150‑1 (MQ=255)
                           tCGAATTACGTAATAGCGATCGCTCCCGGCTCGCACTGGGGGAGGCGAGCCCGACTGAGGAAGTAGATGCAAATGGCAATGTCACACTAAACTTTTTTGCCAATTATCGAGCGTTAGCCAGCGGTGTTCGGCCAGGTGTGGCGAAAGCgg                                                                                                                            >  1:1521834/1‑150 (MQ=255)
                            cGAATTACGTAATAGCGATCGCTCCCGGCTCGCACTGGGGGAGGCGAGCCCGACTGAGGAAGTAGATGCAAATGGCAATGTCACACTAAACTTTTTTGCCAATTATCGAGCGTTAGCCAGCGGTGTTCGGCCAGGTGTGGCGAAAGCGGa                                                                                                                           <  1:2087500/150‑1 (MQ=255)
                             gAATTACGTAATAGCGATCGCTCCCGGCTCGCACTGGGGGAGGCGAGCCCGACTGAGGAAGTAGATGCAAATGGCAATGTCACACTAAACTTTTTTGCCAATTATCGAGCGTTAGCCAGCGGTGTTCGGCCAGGTGTGGCGAAAGCGGAt                                                                                                                          >  2:3751090/1‑150 (MQ=255)
                               aTTACGTAATAGCGATCGCTCCCGACTCGCACTGGGGGAGGCGAGCCCGACTGAGGAAGTAGATGCAAATGGCAATGTCACACTAAACTTTTTTGCCAATTATCGAGCGTTAGCCAGCGGTGTTCGGCCAGGTGTGGCGAAAGCGGATGc                                                                                                                        <  2:4078789/150‑1 (MQ=255)
                                   cGTAAGAGCTATCGCTCCCGGCTCGCACTGGGGGAGGCGAGCCCGACTGAGGAAGTAGATGCAAATGGCAATGTCACACTAAACTTTTTTGCCAATTATCGAGCGTTAGCCAGCGGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGata                                                                                                                    <  1:1161360/150‑1 (MQ=255)
                                     tAATAGCGATCGCTCCCGGCTCGCACTGGGGGAGGCGAGCCCGACTGAGGAAGTAGATGCAAATGGCAATGTCACACTAAACTTTTTTGCCAATTATCGAGCGTTAGCCAGCGGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATAtt                                                                                                                  <  2:3987336/150‑1 (MQ=255)
                                      aaTAGCGCTCGCTCCCGGCTCGCACTGGGGGAGGCGAGCCCGACTGAGGAAGTAGATGCTAATGGCAATGTCACACTAAACTTTTTTGCCAATTATCGAGCGTTAGCCAGCTGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATAttt                                                                                                                 <  2:5287676/150‑1 (MQ=255)
                                      aaTAGCGATCGCTCCCGGCTCGCACTGGGGGAGGCGAGCCCGACTGAGGAAGTAGATGCAAATGGCAATGTCACACTAAACTTTTTTGCCAATTATCGAGCGTTAGCCAGCGGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATAttt                                                                                                                 <  1:4151872/150‑1 (MQ=255)
                                      aaTAGCGATCGCTCCCGGCTCGCACTGGGGGAGGCGAGCCCGACTGAGGAAGTAGATGCAAATGGCAATGTCACACTAAACTTTTTTGCCAATTATCGAGCGTTAGCCAGCGGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATAttt                                                                                                                 <  1:1009614/150‑1 (MQ=255)
                                         aGCGTTCGCTCCCGGCTCGCACTGGGGGAGGCGAGCCCGACTGAGGAAGTAGATGCAAATGGCAATGTCACACTATACTTTTTTGCCAATTATCGAGCGTTAGCCAGCGGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATTTTTATg                                                                                                              <  1:3856065/150‑1 (MQ=255)
                                           cGATCGCTCCCGGCTCGCACTGGGGGAGGCGAGCCCGACTGAGGAAGTAGATGCAAATGGCAATGTCACACTAAACTTTTTTGCCAATTATCGAGCGTTAGCCAGCTGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATATTTATGAt                                                                                                            <  2:3242105/150‑1 (MQ=255)
                                               cGCTCCCGGCTCGCACTGGGGGAGGCGAGCCCGACTGAGGAAGTAGATGCAAATGGCAATGTCACACTAAACTTTTTTGCCAATTATCGAGCGTTAGCCAGCGGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATATTTATGATCAAt                                                                                                        <  1:5114923/150‑1 (MQ=255)
                                                gCTCCCGGCTCGCACTGGGGGAGGCGAGCCCGACTGAGGAAGTAGATGCAAATGGCAATGTCACACTAAACTTTTTTGCCAATTATCGAGCGTTAGCCAGCGGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATATTTATGATCAAtt                                                                                                       >  2:159244/1‑150 (MQ=255)
                                                  tCCCGGCTCGCACTGGGGGAGGCGAGCCCGACTGAGGAAGTAGATGCAAATGGCAATGTCACACTAAACTTTTTTGCCAATTATCGAGCGTTAGCCAGCGGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATATTTATGATCAATtat                                                                                                     <  2:1548225/150‑1 (MQ=255)
                                                  tCCCGGCTCGCACTGGGGGAGGCGAGCCCGACTGAGGAAGTAGATGCAAATGGCAATGTCACACTAAACTTTTTTGCCAATTATCGAGCGTTAGCCAGCGGTGTTAGGCCAGGTGTGGCGAAAGCGGATGCGATATTTATGATCAATtat                                                                                                     >  2:105730/1‑150 (MQ=255)
                                                     cGGCTCGCACTGGGGGAGGCGAGCCCGACTGAGGAAGTAGATGCAAATGGCAATGTCACACTAAACTTTTTTGCCAATTATCGAGCGTTAGCCAGCTGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATATTTATGATCAATTAtaat                                                                                                  >  1:606888/1‑150 (MQ=255)
                                                     cGGCTCGCACTGGGGGAGGCGAGCCCGACTGAGGAAGTAGATGCAAATGGCAATGTCACACTAAACTTTTTTGCCAATTATCGAGCGTTAGCCAGCTGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATATTTATGATCAATTAtaat                                                                                                  <  1:381099/150‑1 (MQ=255)
                                                     cGGCTCGCACTGGGGGAGGCGAGCCCGACTGAGGAAGTAGATGCAAATGGCAATGTCACACTAAACTTTTTTGCCAATTATCGAGCGTTAGCCAGCGGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATATTTATGATCAATTAtaat                                                                                                  <  1:4827015/150‑1 (MQ=255)
                                                      ggCTCGCACTGGGGGAGGCGAGCCCGACTGAGGAAGTAGATGCAAATGGCAATGTCACACTAAACTTTTTTGCCAATTATCGAGCGTTAGCCAGCGGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATATTTATGATCAATTAtaatt                                                                                                 >  2:4799838/1‑150 (MQ=255)
                                                      ggCTCGCACTGGGGGAGGCGAGCCCGACTGAGGAAGTAGATGCAAATGGCAATGTCACACTAAACTTTTTTGCCAATTATCGAGCGTTAGCCAGCGGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATATTTATGATCAATTAtaatt                                                                                                 <  2:3381719/150‑1 (MQ=255)
                                                          cGCACTGGGGGAGGCGAGCCCGACTGAGGAAGTAGATGCAAATGGCAATGTCACACTAAACTTTTTTGCCAATTATCGAGCGTTAGCCAGCTGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATATTTATGATCAATTATAATTAata                                                                                             >  2:1276129/1‑150 (MQ=255)
                                                             aCTGGGGGAGGCGAGCCCGACTGAGGAAGTAGATGCAAATGGCAATGTCACACTAAACTTTTTTGCCAATTATCGAGCGTTAGCCAGCGGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATATTTATGATCAATTATAATTAatatta                                                                                          <  1:5291567/150‑1 (MQ=255)
                                                              cTGGGGGAGGCGAGCCCGACTGAGGAAGTAGATGCAAATGGCAATGTCACACTAAACTTTTTTGCCAATTATCGAGCGTTAGCCAGCGGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATATTTATGATCAATTATAATTAatattat                                                                                         <  2:3069242/150‑1 (MQ=255)
                                                              cTGGGGGAGGCGAGCCCGACTGAGGAAGTAGATGCAAATGGCAATGTCACACTAAACTTTTTTGCCAATTATCGAGCGTTAGCCAGCGGTGTTCGGCCAGGTGTGGCGAAAGCGGATGAGATATTTATGTTCAATTATAATTAatattat                                                                                         >  2:4646631/1‑150 (MQ=255)
                                                                 ggggAGGCGAGCCCGACTGAGGAAGTAGATGCAAATGGCAATGTCACACTAAACTTTTTTGCCAATTATCGAGCGTTAGCCAGCGGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATATTTATGATCAATTATAATTAatattatatt                                                                                      <  2:2814111/150‑1 (MQ=255)
                                                                 ggggAGGCGAGCCCGACTGAGGAAGTAGATGCAAATGGCAATGTCACACTAAACTTTTTTGCCAATTATCGAGCGTTAGCCAGCGGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATATTTATGATCAATTATAATTAatattatatt                                                                                      <  2:1376949/150‑1 (MQ=255)
                                                                     aGGCGAGCCCGACTGAGGAAGTAGATGCAAATGGCAATGTCACACTAAACTTTTTTGCCAATTATCGAGCGTTAGCCAGCTGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATATTTATGATCAATTATAATTAATATTATATTAAtt                                                                                  <  1:1722162/150‑1 (MQ=255)
                                                                     aGGCGAGCCCGACTGAGGAAGTAGATGCAAATGGCAATGTCACACTAAACTTTTTTGCCAATTATCGAGCGTTAGCCAGCGGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATATTTATGATCAATTATAATTAATATTATATTAAtt                                                                                  >  1:2509888/1‑150 (MQ=255)
                                                                      ggCGAGCCCGACTGAGGAAGTAGATGCAAATGGCAATGTCACACTAAACTTTTTTGCCAATTATCGAGCGTTAGCCAGCGGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATATTTATGATCAATTATAATTAATATTATATTAATTc                                                                                 <  1:3372528/150‑1 (MQ=255)
                                                                          aGCCCGACTGAGGAAGTAGATGCAAATGGCAATGTCACACTAAACTTTTTTGCCAATTATCGAGCGTTAGCCAGCTGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATATTTATGATCAATTATAATTAATATTATATTAATTCGtat                                                                             <  2:2671328/150‑1 (MQ=255)
                                                                          aGCCCGACTGAGGAAGTAGATGCAAATGGCAATGTCACACTAAACTTTTTTGCCAATTATCGAGCGTTAGCCAGCTGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATATTTATGATCAATTATAATTAATATTATATTAATTCGtat                                                                             >  1:3737762/1‑150 (MQ=255)
                                                                          aGCCCGACTGAGGAAGTAGATGCAAATGGCAATGTCACACTAAACTTTTTTGCCAATTATCGAGCGTTAGCCAGCTGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATATTTATGATCAATTATAATTAATATTATATTAATTCGtat                                                                             <  1:1941525/150‑1 (MQ=255)
                                                                            cccGACTGAGGAAGTAGATGCAAATGGCAATGTCACACTAAACTTTTTTGCCAATTATCGAGCGTTAGCCAGCTGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATATTTATGATCAATTATAATTAATATTATATTAATTCGTATaa                                                                           <  1:2413731/150‑1 (MQ=255)
                                                                            cccGACTGAGGAAGTAGATGCAAATGGCAATGTCACACTAAACTTTTTTGCCAATTATCGAGCGTTAGCCAGCGGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATATTTATGATCAATTATAATTAATATTATATTAATTCGTATaa                                                                           <  1:4723401/150‑1 (MQ=255)
                                                                                  tGAGGAAGTAGATGCAAATGGCAATGTCACACTAAACTTTTTTGCCAATTATCGAGCGTTAGCCAGCGGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATATTTATGATCAATTATAATTAATATTATATTAATTCGTATAATTTGGc                                                                     <  2:2888666/150‑1 (MQ=255)
                                                                                      gAAGTAGATGCAAATGGCAATGTCACACTAAACTTTTTTGCCAATTATCGAGCGTTAGCCAGCGGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATATTTATGATCAATTATAATTAATATTATATTAATTCGTATAATTTGGCGTAg                                                                 >  2:4289772/1‑150 (MQ=255)
                                                                                      gAAGTAGATGCAAATGGCAATGTCACACTAAACTTTTTTGCCAATTATCGAGCGTTAGCCAGCGGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATATTTATGATCAATTATAATTAATATTATATTAATTCGTATAATTTGGCGTAg                                                                 >  1:3506329/1‑150 (MQ=255)
                                                                                           aGATGCAAATGGCAATGTCACACTAAACTTTTTTGCCAATTATCGAGCGTTAGCCAGCGGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATATTTATGATCAATTATAATTAATATTATATTAATTCGTATAATTTGGCGTAGTCGAt                                                            >  2:3321743/1‑150 (MQ=255)
                                                                                            gATGCAAATGGCAATGTCACACTAAACTTTTTTGCCAATTATCGAGCGTTAGCCAGCTGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATATTTATGATCAATTATAATTAATATTATATTAATTCGTATAATTTGGCGTAGTCGATa                                                           >  1:4429953/1‑150 (MQ=255)
                                                                                               gCAAATGGCAATGTCACACTAAACTTTTTTGCCAATTATCGAGCGTTAGCCAGCGGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATATTTATGATCAATTATAATTAATATTATATTAATTCGTATAATTTGGCGTAGTCGATAAGc                                                        >  2:3655013/1‑150 (MQ=255)
                                                                                               gCAAATGGCAATGTCACACTAAACTTTTTTGCCAATTATCGAGCGTTAGCCAGCGGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATATTTATGATCAATTATAATTAATATTATATTAATTCGTATAATTTGGCGTAGTCGATAAGc                                                        <  1:332442/150‑1 (MQ=255)
                                                                                                 aaaTGGCAATGTCACACTAAACTTTTTTGCCAATTATCGAGCGTTAGCCAGCTGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATATTTATGATCAATTATAATTAATATTATATTAATTCGTATAATTTGGCGTAGTCGATAAGctc                                                      <  2:1206934/150‑1 (MQ=255)
                                                                                                  aaTGGCAATGTCACACTAAACTTTTTTGCCAATTATCGAGCGTTAGCCAGCGGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATATTTATGATCAATTATAATTAATATTATATTAATTCGTATAATTTGGCGTAGTCGATAAGctct                                                     >  1:3042583/1‑150 (MQ=255)
                                                                                                   aTGGCAATGTCACACTAAACTTTTTTGCCAATTATCGAGCGTTAGCCAGCGGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATATTTATGATCAATTATAATTAATATTATATTAATTCGTATAATTTGGCGTAGTCGATAAGCTCTa                                                    >  2:3851405/1‑150 (MQ=255)
                                                                                                     ggCAATGTCACACTAAACTTTTTTGCCAATTATCGAGCGTTAGCCAGCGGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATATTTATGATCAATTATAATTAATATTATATTAATTCGTATAATTTGGCGTAGTCGATAAGCTCTACa                                                  >  2:4839580/1‑150 (MQ=255)
                                                                                                        aaTGTCACACTAAACTTTTTTGCCAATTATCTAGCGTTAGCCAGCGGTGTTCGGCCAGGTGTTGCGAAAGCGGATGCGATATTTATGATCAATTATAATTAATATTATATTAATTCTTATAATTTGGCGTAGTCGATAAGCTCTACAAtt                                               >  1:3230308/1‑150 (MQ=255)
                                                                                                        aaTGTCACACTAAACTTTTTTGCCAATTATCGAGCGTTAGCCAGCGGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATATTTATGATCAATTATAATTAATATTATATTAATTCGTATAATTTGGCGTAGTCGATAAGCTCTACAAtt                                               <  1:1311138/150‑1 (MQ=255)
                                                                                                        aaTGTCACACTAAACTTTTTTGCCAATTATCGAGCGTTAGCCAGCGGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATATTTATGATCAATTATAATTAATATTATATTAATTCGTATAATTTGGCGTAGTCGATAAGCTCTACAAtt                                               <  1:2871714/150‑1 (MQ=255)
                                                                                                        aaTGTCACACTAAACTTTTTTGCCAATTATCGAGCGTTAGCCAGCGGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATATTTATGATCAATTATAATTAATATTATATTAATTCGTATAATTTGGCGTAGTCGATAAGCTCTACAAtt                                               >  2:1757571/1‑150 (MQ=255)
                                                                                                         aTGTCACACTAAACTTTTTTGCCAATTATCGAGCGTTAGCCAGCGGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATATTTATGATCAATTATAATTAATATTATATTAATTCGTATAATTTGGCGTAGTCGATAAGCTCTACAATTg                                              <  2:1918306/150‑1 (MQ=255)
                                                                                                         aTGTCACACTAAACTTTTTTGCCAATTATCGAGCGTTAGCCAGCGGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATATTTATGATCAATTATAATTAATATTATATTAATTCGTATAATTTGGCGTAGTCGATAAGCTCTACAATTg                                              <  1:865057/150‑1 (MQ=255)
                                                                                                         aTGTCACACTAAACTTTTTTGCCAATTATCGAGCGTTAGCCAGCGGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATATTTATGATCAATTATAATTAATATTATATTAATTCGTATAATTTGGCGTAGTCGATAAGCTCTACAATTg                                              <  1:1599418/150‑1 (MQ=255)
                                                                                                           gTCACACTAAACTTTTTTGCCAATTATCGAGCGTTAGCCAGCGGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATATTTATGATCAATTATAATTAATATTATATTAATTCGTATAATTTGGCGTAGTCGATAAGCTCTACAATTGaa                                            <  2:4352381/150‑1 (MQ=255)
                                                                                                           gTCACACTAAACTTTTTTGCCAATTATCGAGCGTTAGCCAGCGGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATAATTATGATCAATTATAATTAATATTATATTAATTCGTATAATTTGGCGTAGTCGATAAGCTCTACAATTGaa                                            <  2:2227554/150‑1 (MQ=255)
                                                                                                            tCACACTAAACTTTTTTGCCAATTATCGAGCGTTAGCCAGCGGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATATTTATGATCAATTATAATTAATATTATATTAATTCGTATAATTTGGCGTAGTCGATAAGCTCTACAATTGAAt                                           <  1:2297647/150‑1 (MQ=255)
                                                                                                            tCACACTAAACTTTTTTGCCAATTATCGAGCGTTAGCCAGCGGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATATTTATGATCAATTATAATTAATATTATATTAATTCGTATAATTTGGCGTAGTCGATAAGCTCTACAATTGAAt                                           <  1:4458270/150‑1 (MQ=255)
                                                                                                             cacaCTAAACTTTTTTGCCAATTATCGAGCGTTAGCCAGCGGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATATTTATGATCAATTATAATTAATATTATATTAATTCGTATAATTTGGCGTAGTCGATAAGCTCTACAATTGAAtt                                          >  2:3174615/1‑149 (MQ=255)
                                                                                                              acacTAAACTTTTTTGCCAATTATCGAGCGTTAGCCAGCGGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATATTTATGATCAATTATAATTAATATTATATTAATTCGTATAATTTGGCGTAGTCGATAAGCTCTACAATTGAATGc                                         >  1:4511655/1‑150 (MQ=255)
                                                                                                                  tAAACTTTTTTGCCAATTATCGAGCGTTAGCCAGCTGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATATTTATGATCAATTATAATTAATATTATATTAATTCGTATAATTTGGCGTAGTCGATAAGCTCTACAATTGAATGCAAAc                                     <  1:678654/150‑1 (MQ=255)
                                                                                                                   aaaCTTTTTTGCCAATTATCGAGCGTTAGCCAGCGGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATATTTATGATCAATTATAATTAATATTATATTAATTCGTATAATTTGGCGTAGTCGATAAGCTCTACAATTGAATGCAAAcc                                    >  2:1203373/1‑150 (MQ=255)
                                                                                                                   aaaCTTTTTTGCCAATTATCGAGCGTTAGCCAGCGGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATATTTATGATCAATTATAATTAATATTATATTAATTCGTATAATTTGGCGTAGTCGATAAGCTCTACAATTGAATGCAAAcc                                    <  1:4913328/150‑1 (MQ=255)
                                                                                                                    aaCTTTTTTGCCAATTATCGAGCGTTAGCCAGCGGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATATTTATGATCAATTATAATTAATATTATATTAATTCGTATAATTTGGCGTAGTCGATAAGCTCTACAATTGAATGCAAACCt                                   <  2:3829991/150‑1 (MQ=255)
                                                                                                                    aaCTTTTTTGCCAATTATCGAGCGTTAGCCAGCGGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATATTTATGATCAATTATAATTAATATTATATTAATTCGTATAATTTGGCGTAGTCGATAAGCTCTACAATTGAATGCAAACCt                                   <  2:1024213/150‑1 (MQ=255)
                                                                                                                     aCTTTTTTGCCAATTATCGAGCGTTAGCCAGCGGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATATTTATGATCAATTATAATTAATATTATATTAATTCGTATAATTTGGCGTAGTCGATAAGCTCTACAATTGAATGCAAACCTa                                  >  2:4543420/1‑150 (MQ=255)
                                                                                                                        tttttGCCAATTATCGAGCGTTAGCCAGCGGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATATTTATGATCAATTATAATTAATATTATATTAATTCGTATAATTTGGCGTAGTCGATAAGCTCTACAATTGAATGCAAACCTAGCt                               <  2:3923497/150‑1 (MQ=255)
                                                                                                                         ttttGCCAATTATCGAGCGTTAGCCAGCTGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATATTTATGATCAATTATAATTAATATTATATTAATTCGTATAATTTGGCGTAGTCGATAAGCTCTACAATTGAATGCAAACCTAGCtt                              <  1:4272883/150‑1 (MQ=255)
                                                                                                                         ttttGCCAATTATCGAGCGTTAGCCAGCGGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATATTTATGATCAATTATAATTAATATTATATTAATTCGTATAATTTGGCGTAGTCGATAAGCTCTACAATTGAATGCAAACCTAGCtt                              <  1:4605441/150‑1 (MQ=255)
                                                                                                                          tttGCCAATTATCGAGCGTTAGCCAGCGGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATATTTATGATCAATTATAATTAATATTATATTAATTCGTATAATTTGGCGTAGTCGATAAGCTCTACAATTGAATGCAAACCTAGCTTg                             <  2:2416413/150‑1 (MQ=255)
                                                                                                                           ttGCCAATTATCGAGCGTTAGCCAGCGGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATATTTATGATCAATTATAATTAATATTATATTAATTCGTATAATTTGGCGTAGTCGATAAGCTCTACAATTGAATGCAAACCTAGCTTGc                            <  2:2012069/150‑1 (MQ=255)
                                                                                                                              ccAATTATCGAGCGTTAGCCAGCTGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATATTTATGATCAATTATAATTAATATTATATTAATTCGTATAATTAGGCGTAGTCGATAAGCTCTACAATTGAATGCAAACCTAGCTTGCCAt                         >  1:4646304/1‑150 (MQ=255)
                                                                                                                                aaTTATCGAGCGTTAGCCAGCGGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATATTTATGATCAATTATAATTAATATTATATTAATTCGTATAATTTGGCGTAGTCGATAAGCTCTACAATTGAATGCAAACCTAGCTTGCCATaa                       >  1:3401936/1‑150 (MQ=255)
                                                                                                                                 aTTATCGAGCGTTAGCCAGCTGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATATTTATGATCAATTATAATTAATATTATATTAATTCGTATAATTTGGCGTAGTCGATAAGCTCTACAATTGAATGCAAACCTAGCTTGCCATaaa                      >  1:1749574/1‑150 (MQ=255)
                                                                                                                                  ttATCGAGCGTTAGCCAGCGGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATATTTATGATCAATTATAATTAATATTATATTAATTCGTATAATTTGGCGTAGTCGATAAGCTCTACAATTGAATGCAAACCTAGCTTGCCATAAat                     <  2:3711183/150‑1 (MQ=255)
                                                                                                                                  ttATCGAGCGTTAGCCAGCGGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATATTTATGATCAATTATAATTAATATTATATTAATTCGTATAATTTGGCGTAGTCGATAAGCTCTACAATTGAATGCAAACCTAGCTTGCCATAAat                     <  1:1424686/150‑1 (MQ=255)
                                                                                                                                  ttATCGAGCGTTAGCCAGCGGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATATTTATGATCAATTATAATTAATATTATATTAATTCGTATAATTTGGCGTAGTCGATAAGCTCTACAATTGAATGCAAACCTAGCTTGCCATAAat                     <  1:2148547/150‑1 (MQ=255)
                                                                                                                                    aTCGAGCGTTAGCCAGCTGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATATTTATGATCAATTATAATTAATCTTATATTAATTCGTATAATTTGGCGTAGTCGATAAGCTCTACAATTGAATGCAAACCTAGCTTGCCATAAatat                   >  2:2024278/1‑150 (MQ=255)
                                                                                                                                      cGAGCGTTAGCCAGCTGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATATTTATGATCAATTATAATTAATATTATATTAATTCGTATAATTTGGCGTAGTCGATAAGCTCTACAATTGAATGCAAACCTAGCTTGCCATAAATATTa                 <  2:745954/150‑1 (MQ=255)
                                                                                                                                       gAGCGTTAGCCAGCTGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATATTTATGATCAATTATAATTAATATTATATTAATTCGTATAATTTGGCGTAGTCGATAAGCTCTACAATTGAATGCAAACCTAGCTTGCCATAAATATTAg                >  2:5244837/1‑150 (MQ=255)
                                                                                                                                           gTTAGCCAGCGGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATATTTATGATCAATTATAATTAATATTATATTAATTCGTATAATTTGGCGTAGTCGATAAGCTCTACAATTGAATGCAAACCTAGCTTGCCATAAATATTAGAttt            <  2:3788022/150‑1 (MQ=255)
                                                                                                                                             tAGCCAGCTGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATATTTATGATCAATTATAATTAATATTATATTAATTCGTATAATTTGGCGTAGTCGATAAGCTCTACAATTGAATGCAAACCTAGCTTGCCATAAATATTAGATTTAt          >  1:2749326/1‑150 (MQ=255)
                                                                                                                                             tAGCCAGCGGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATATTTATGATCAATTATAATTAATATTATATTAATTCGTATAATTTGGCGTAGTCGATAAGCTCTACAATTGAATGCAAACCTAGCTTGCCATAAATATTAGATTTAt          >  2:4717826/1‑150 (MQ=255)
                                                                                                                                             tAGCCAGCGGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATATTTATGATCAATTATAATTAATATTATATTAATTCGTATAATTTGGCGTAGTCGATAAGCTCTACAATTGAATGCAAACCTAGCTTGCCATAAATATTAGATTTAt          <  2:558891/150‑1 (MQ=255)
                                                                                                                                                  aGCGGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATATTTATGATCAATTATAATTAATATTATATTAATTCGTATAATTTGGCGTAGTCGATAAGCTCTACAATTGAATGCAAACCTAGCTTGCCATAAATATTAGATTTATGCGCa     >  1:1856504/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATCCCGCAATGGCACAAACCGTCGCCATCGAATTACGTAATAGCGATCGCTCCCGGCTCGCACTGGGGGAGGCGAGCCCGACTGAGGAAGTAGATGCAAATGGCAATGTCACACTAAACTTTTTTGCCAATTATCGAGCGTTAGCCAGCGGTGTTCGGCCAGGTGTGGCGAAAGCGGATGCGATATTTATGATCAATTATAATTAATATTATATTAATTCGTATAATTTGGCGTAGTCGATAAGCTCTACAATTGAATGCAAACCTAGCTTGCCATAAATATTAGATTTATGCGCACTA  >  CP009273/559096‑559394

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: ‑