Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 889,044 G→T 28.3% intergenic (+155/‑196) ybjN → / → potF negative regulator of motility; multicopy suppressor of coaA(Ts)/putrescine transporter subunit: periplasmic‑binding component of ABC superfamily

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273889,0440GT28.3% 374.9 / 125.7 178intergenic (+155/‑196)ybjN/potFnegative regulator of motility; multicopy suppressor of coaA(Ts)/putrescine transporter subunit: periplasmic‑binding component of ABC superfamily
Reads supporting (aligned to +/- strand):  ref base G (57/70);  new base T (17/34);  total (74/104)
Fisher's exact test for biased strand distribution p-value = 1.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 182 total aligned reads displayed.

TTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTCTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGCTTTTATATAGCGAGCAGTGCTGGCCGGGAGAAAGTTCTCTTTTCTTACACCGCGCCGATAAAAAATATGCACGTTTATTGCATATCTTTCAGTG  >  CP009273/888897‑889192
                                                                                                                                                   |                                                                                                                                                    
tttCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTCTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGtt                                                                                                                                                    <  2:622438/150‑1 (MQ=255)
 ttCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTCTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCtttt                                                                                                                                                   <  2:1025890/150‑1 (MQ=255)
  tCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTCTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGtttt                                                                                                                                                  <  1:3397984/150‑1 (MQ=255)
   cGAGCCCGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTCTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCtttttt                                                                                                                                                 <  2:3063394/150‑1 (MQ=255)
    gAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTCTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCTtttttt                                                                                                                                                >  1:3374139/1‑150 (MQ=255)
    gAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTCTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGtttttt                                                                                                                                                <  1:1058699/150‑1 (MQ=255)
      gCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTCTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTaa                                                                                                                                              <  1:3738644/150‑1 (MQ=255)
         cGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTCTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATcc                                                                                                                                           <  2:3422138/150‑1 (MQ=255)
         cGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTCTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATcc                                                                                                                                           <  2:3420171/150‑1 (MQ=255)
         cGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTCTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATcc                                                                                                                                           <  1:3810649/150‑1 (MQ=255)
               cGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTCTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCta                                                                                                                                     >  2:1079271/1‑150 (MQ=255)
               cGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGATGTTTTTAACCTTTTCTTAAAGATTATTTAACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCACATATGCATAAAAATATGTTAAATACCGTTTTTTAATCCGAGCtc                                                                                                                                     >  2:2139194/1‑149 (MQ=255)
                 cTGGTGCAGTGGCTGCGCGCTGCAAAATTCTCTGCTGTTTTTAACCTATTCTTAAAGATGATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCtata                                                                                                                                   <  1:2559341/150‑1 (MQ=255)
                    gTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTCTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAATTACCGTTTTTTAATACGAGCTATAGtc                                                                                                                                >  1:1006055/1‑150 (MQ=255)
                    gTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTCTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGtc                                                                                                                                >  1:1006112/1‑150 (MQ=255)
                     tGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTCTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGtct                                                                                                                               <  1:1786923/150‑1 (MQ=255)
                              tGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTCTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTg                                                                                                                      <  2:2801667/150‑1 (MQ=255)
                               gcgcgcTGCAAAATTATCTGCTGTTTTTAACCTTTTCTTAAAGATTATTTCACTTCTCTTGTGTCGAGTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTCAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTgg                                                                                                                     >  1:2116539/1‑150 (MQ=255)
                                cgcgcTGCAAAATTATCTGCTGTTTTTAACCTTTTCTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGc                                                                                                                    <  2:3384639/150‑1 (MQ=255)
                                       cAAAATTATCTGCTGTTTTTAACCTTTTCTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGtt                                                                                                             <  2:2134237/150‑1 (MQ=255)
                                          aaTTATCTGCTGTTTTTAACCTTTTCTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCACATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTAtt                                                                                                          >  2:1024003/1‑150 (MQ=255)
                                            ttATCTGCTGTTTTTAACCTTTTCTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCt                                                                                                        >  2:3544322/1‑150 (MQ=255)
                                             tATCTGCTGTTTTTAACCTTTTCTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCtt                                                                                                       >  2:2119662/1‑150 (MQ=255)
                                             tATCTGCTGTTTTTAACCTTTTCTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCtt                                                                                                       <  2:2623922/150‑1 (MQ=255)
                                                ctgctgTTTTTAACCTTTTCTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCg                                                                                                    <  1:576137/150‑1 (MQ=255)
                                                 tgctgTTTTTAACCTTTTCTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATAACGTTTTTTAATCCGAGCTATAGTCTCAAAACCTGGCTAAAGTTATTCTTGcta                                                                                                   >  2:909514/1‑148 (MQ=255)
                                                    tgTTTTTAACCTTTTCTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGc                                                                                                >  1:3614073/1‑150 (MQ=255)
                                                    tgTTTTTAACCTTTTCTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGc                                                                                                >  2:3886287/1‑150 (MQ=255)
                                                       ttttAACCTTTTCTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGCttt                                                                                             <  1:1133351/150‑1 (MQ=255)
                                                       ttttAACCTTTTCTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGCttt                                                                                             <  2:2731292/150‑1 (MQ=255)
                                                        tttAACCTTTTCTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGCtttt                                                                                            <  2:3488802/150‑1 (MQ=255)
                                                         ttAACCTTTTCTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCTTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGCTTTta                                                                                           <  2:3411840/150‑1 (MQ=255)
                                                         ttAACCTTTTCTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGCTTTta                                                                                           >  2:2172380/1‑150 (MQ=255)
                                                         ttAACCTTTTCTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGCTTTta                                                                                           >  1:2198780/1‑150 (MQ=255)
                                                         ttAACCTTTTCTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGCTTTta                                                                                           >  1:4598956/1‑150 (MQ=255)
                                                            aCCTTTTCTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCTTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGCTTTtatat                                                                                        <  2:3112625/150‑1 (MQ=255)
                                                            aCCTTTTCTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCTTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGCTTTtatat                                                                                        <  2:3437386/150‑1 (MQ=255)
                                                            aCCTTTTCTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGCTTTtatat                                                                                        <  2:2465087/150‑1 (MQ=255)
                                                             ccTTTTCTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCTTTTTTTTATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGCTTTtatata                                                                                       <  2:4625687/150‑1 (MQ=255)
                                                              cTTTTCTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGCTTTTATATAg                                                                                      >  1:3167420/1‑150 (MQ=255)
                                                                 ttCTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCTTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGCTTTTATATAGCGa                                                                                   <  2:3373751/150‑1 (MQ=255)
                                                                 ttCTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGCTTTTATATAGCGa                                                                                   >  2:2373076/1‑150 (MQ=255)
                                                                     tAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCTTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGCTTTTATATAGCGAGCAg                                                                               <  1:4234868/150‑1 (MQ=255)
                                                                     tAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGCTTTTATATAGCGAGCAg                                                                               >  1:1405292/1‑150 (MQ=255)
                                                                     tAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGCTTTTATATAGCGAGCAg                                                                               >  2:2938201/1‑150 (MQ=255)
                                                                      aaaGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCTTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGCTTTTATATAGCGAGCAGt                                                                              >  2:3559952/1‑150 (MQ=255)
                                                                      aaaGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCTTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGCTTTTATATAGCGAGCAGt                                                                              >  2:1148897/1‑150 (MQ=255)
                                                                      aaaGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGCTTTTATATAGCGAGCAGt                                                                              >  1:92680/1‑150 (MQ=255)
                                                                       aaGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGCTTTTATATAGCGAGCAGTg                                                                             <  1:61436/150‑1 (MQ=255)
                                                                        aGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGCTTTTATATAGCGAGCAGTGc                                                                            >  1:2969124/1‑150 (MQ=255)
                                                                         gATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGCTTTTATATAGCGAGCAGTGCt                                                                           <  1:4121436/150‑1 (MQ=255)
                                                                         gATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGCTTTTATATAGCGAGCAGTGCt                                                                           <  1:4122427/150‑1 (MQ=255)
                                                                          attattTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCTTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGCTTTTATATAGCGAGCAGTGCTg                                                                          >  2:1766506/1‑150 (MQ=255)
                                                                           ttattTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGCTTTTATATAGCGAGCAGTGCTgg                                                                         >  2:1665907/1‑150 (MQ=255)
                                                                              tttCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGCTTTTATATAGCGAGCAGTGCTGGCCg                                                                      >  2:2580841/1‑150 (MQ=255)
                                                                                tCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCTTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGCTTTTATATAGCGAGCAGTGCTGGCCggg                                                                    <  1:4017155/150‑1 (MQ=255)
                                                                                 cACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGCTTTTATATAGCGAGCAGTGCTGGCCGGga                                                                   <  2:3559713/150‑1 (MQ=255)
                                                                                 cACTTCTCGTGTGTCGATTTGGCTTTATCACCTAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGCTTTTATATAGCGAGCAGTGCTGGCCGGga                                                                   <  2:4227779/150‑1 (MQ=255)
                                                                                  aCTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCTTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGCTTTTATATAGCGAGCAGTGCTGGCCGGgag                                                                  >  2:4746320/1‑150 (MQ=255)
                                                                                  aCTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCTTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGCTTTTATATAGCGAGCAGTGCTGGCCGGgag                                                                  >  2:4746446/1‑150 (MQ=255)
                                                                                   cTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCTTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGCTTTTATATAGCGAGCAGTGCTGGCCGGgaga                                                                 <  2:1047995/150‑1 (MQ=255)
                                                                                   cTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGCTTTTATATAGCGAGCAGTGCTGGCCGGgaga                                                                 >  2:3863144/1‑150 (MQ=255)
                                                                                     tctcTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGCTTTTATATAGCGAGCAGTGCTGGCCGGGAGaaa                                                               <  2:3112982/150‑1 (MQ=255)
                                                                                       tctTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGCTTTTATATAGCGAGCAGTGCTGGCCGGGAGAAAGt                                                             >  2:908857/1‑150 (MQ=255)
                                                                                        ctTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCTTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGCTTTTATATAGCGAGCAGTGCTGGCCGGGAGAAAGtt                                                            >  2:3992522/1‑150 (MQ=255)
                                                                                         gggtgtCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGCTTTTATATAGCGAGCAGTGCTGGCCGGGAGAAAGTtc                                                           <  1:1816866/148‑1 (MQ=255)
                                                                                           gtgtCGATTTGGCTTTATCACATAGAGCAAATATGCATAAATATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGCTTTTATATAGCGAGCAGTGCTGGCCGGGAGAAAGTtctc                                                         <  1:1574027/150‑1 (MQ=255)
                                                                                             gtCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGCTTTTATATAGCGAGCAGTGCTGGCCGGGAGAAAGTTCTCtt                                                       >  1:361742/1‑150 (MQ=255)
                                                                                              tCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCTTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGCTTTTATATAGCGAGCAGTGCTGGCCGGGAGAAAGTTCTCttt                                                      <  1:4535512/150‑1 (MQ=255)
                                                                                                 aTTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCTTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGCTTTTATATAGCGAGCAGTGCTGGCCGGGAGAAAGTTCTCTTTTCt                                                   <  1:3813890/150‑1 (MQ=255)
                                                                                                  tttGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGCTTTTATATAGCGAGCAGTGCTGGCCGGGAGAAAGTTCTCTTTTCtt                                                  <  2:1848046/150‑1 (MQ=255)
                                                                                                   ttGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGCTTTTATATAGCGAGCAGTGCTGGCCGGGAGAAAGTTCTCTTTTCTTa                                                 <  2:3482615/150‑1 (MQ=255)
                                                                                                    tGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGCTTTTATATAGCGAGCAGTGCTGGCCGGGAGAAAGTTCTCTTTTCTTac                                                <  2:1470737/150‑1 (MQ=255)
                                                                                                             ctcaTAGAGCATATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGCTTTTATATAGCGAGCAGTGCTGGCCGGGAGAAAGTTCTCTTTTCTTACACCGCGCCg                                       <  1:3093195/148‑1 (MQ=255)
                                                                                                             cacaTAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGCTTTTATATAGCGAGCAGTGCTGGCCGGGAGAAAGTTCTCTTTTCTTACACCGCGCCg                                       <  1:4382849/150‑1 (MQ=255)
                                                                                                              acaTAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGCTTTTATATAGCGAGCAGTGCTGGCCGGGAGAAAGTTCTCTTTTCTTGCACCGCGCCGa                                      <  2:1718538/150‑1 (MQ=255)
                                                                                                              acaTAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGCTTTTATATAGCGAGCAGTGCTGGCCGGGAGAAAGTTCTCTTTTCTTACACCGCGCCGa                                      >  2:4536055/1‑150 (MQ=255)
                                                                                                               caTAGAGCAAATATGCATAAAAATTTGTTAAATACCTTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGCTTTTATATAGCGAGCAGTGCTGGCCGGGAGAAAGTTCTCTTTTCTTACACCGCGCCGAt                                     <  1:2001625/150‑1 (MQ=255)
                                                                                                               caTAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGCTTTTATATAGCGAGCAGTGCTGGCCGGGAGAAAGTTCTCTTTTCTTACACCGCGCCGAt                                     <  1:1723658/150‑1 (MQ=255)
                                                                                                                    agCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGCTTTTATATAGCGAGCAGTGCTGGCCGGGAGAAAGTTCTCTTTTCTTACACCGCGCCGATaaaaa                                <  2:1946390/150‑1 (MQ=255)
                                                                                                                      cAAATATGCATAAAAATTTGTTAAATACCTTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGCTTTTATATAGCGAGCAGTGCTGGCCGGGAGAAAGTTCTCTTTTCTTACACCGCGCCGATAAAAAat                              <  1:203282/150‑1 (MQ=255)
                                                                                                                       aaaTATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGCTTTTATATAGCGAGCAGTGCTGGCCGGGAGAAAGTTCTCTTTTCTTACACCGCGCCGATAAAAAata                             >  1:2483590/1‑150 (MQ=255)
                                                                                                                        aaTATGCATAAAAATTTGTTAAATACCTTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGCTTTTATATAGCGAGCAGTGCTGGCCGGGAGAAAGTTCTCTTTTCTTACACCGCGCCGATAAAAAatag                            <  1:4723399/150‑2 (MQ=255)
                                                                                                                         atatGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGCTTTTATATAGCGAGCAGTGCTGGCCGGGAGAAAGTTCTCTTTTCTTACACCGCGCCGATAAAAAATATg                           >  2:834575/1‑150 (MQ=255)
                                                                                                                           atGCATAAAAATTTGTTAAATACCTTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGCTTTTATATAGCGAGCAGTGCTGGCCGGGAGAAAGTTCTCTTTTCTTACACCGCGCCGATAAAAAATATGCa                         <  1:2182718/150‑1 (MQ=255)
                                                                                                                               aTAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGCTTTTATATAGCGAGCAGTGCTGGCCGGGAGAAAGTTCTCTTTTCTTACACCGCGCCGATAAAAAATATGCACGtt                     >  2:4021172/1‑150 (MQ=255)
                                                                                                                                    aaTTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGCTTTTATATAGCGAGCAGTGCTGGCCGGGAGAAAGTTCTCTTTTCTTACACCGCGCCGATAAAAAATATGCACGTTTATTg                >  2:4283752/1‑150 (MQ=255)
                                                                                                                                        tGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGCTTTTATATAGCGAGCAGTGCTGGCCGGGAGAAAGTTCTCTTTTCTTACACCGCGCCGATAAAAAATATGCACGTTTATTGCata            <  1:4286397/150‑1 (MQ=255)
                                                                                                                                            aaaTACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGCTTTTATATAGCGAGCAGTGCTGGCCGGGAGAAAGTTCTCTTTTCTTACACCGCGCCGATAAAAAATATGCACGTTTATTGCATATCtt        >  2:4021104/1‑150 (MQ=255)
                                                                                                                                            aaaTACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGCTTTTATATAGCGAGCAGTGCTGGCCGGGAGAAAGTTCTCTTTTCTTACACCGCGCCGATAAAAAATATGCACGTTTATTGCATATCtt        >  1:3179888/1‑150 (MQ=255)
                                                                                                                                            aaaTACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGCTTTTATATAGCGAGCAGTGCTGGCCGGGAGAAAGTTCTCTTTTCTTACACCGCGCCGATAAAAAATATGCACGTTTATTGCATATCtt        >  1:3180316/1‑150 (MQ=255)
                                                                                                                                            aaaTACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGCTTTTATATAGCGAGCAGTGCTGGCCGGGAGAAAGTTCTCTTTTCTTACACCGCGCCGATAAAAAATATGCACGTTTATTGCATATCtt        <  1:1024713/150‑1 (MQ=255)
                                                                                                                                             aaTACCTTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGCTTTTATATAGCGAGCAGTGCTGGCCGGGAGAAAGTTCTCTTTTCTTACACCGCGCCGATAAAAAATATGCACGTTTATTGCATATCttt       <  2:1153198/150‑1 (MQ=255)
                                                                                                                                             aaTACCTTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGCTTTTATATAGCGAGCAGTGCTGGCCGGGAGAAAGTTCTCTTTTCTTACACCGCGCCGATAAAAAATATGCACGTTTATTGCATATCttt       <  2:1153103/150‑1 (MQ=255)
                                                                                                                                             aaTACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGCTTTTATATAGCGAGCAGTGCTGGCCGGGAGAAAGTTCTCTTTTCTTACACCGCGCCGATAAAAAATATGCACGTTTATTGCATATCttt       >  1:2461564/1‑150 (MQ=255)
                                                                                                                                              aTACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGCTTTTATATAGCGAGCAGTGCTGGCCGGGAGAAAGTTCTCTTTTCTTACACCGCGCCGATAAAAAATATGCACGTTTATTGCATATCTTTc      >  1:4771814/1‑150 (MQ=255)
                                                                                                                                               tACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGCTTTTATATAGCGAGCAGTGCTGGCCGGGAGAAAGTTCTCTTTTCTTACACCGCGCCGATAAAAAATATGCACGTTTATTGCATATCTTTCa     >  1:2855482/1‑150 (MQ=255)
                                                                                                                                                aCCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGCTTTTATATAGCGAGCAGTGCTGGCCGGGAGAAAGTTCTCTTTTCTTACACCGCGCCGATAAAAAATATGCACGTTTATTGCATATCTTTCAg    >  2:1885671/1‑150 (MQ=255)
                                                                                                                                                aCCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGCTTTTATATAGCGAGCAGTGCTGGCCGGGAGAAAGTTCTCTTTTCTTACACCGCGCCGATAAAAAATATGCACGTTTATTGCATATCTTTCAg    >  1:924354/1‑150 (MQ=255)
                                                                                                                                                  cGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGCTTTTATATAGCGAGCAGTGCTGGCCGGGAGAAAGTTCTCTTTTCTTACACCGCGCCGATAAAAAATATGCACGTTTATTGCATATCTTTCAgtg  >  1:312571/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TTTCGAGCACGCAGTCGCTGGTGCAGTGGCTGCGCGCTGCAAAATTATCTGCTGTTTTTAACCTTTTCTTAAAGATTATTTCACTTCTCTTGTGTCGATTTGGCTTTATCACATAGAGCAAATATGCATAAAAATTTGTTAAATACCGTTTTTTAATCCGAGCTATAGTCTCAAACCCTGGCTAAAGTTATTCTTGCGATGCTTTTATATAGCGAGCAGTGCTGGCCGGGAGAAAGTTCTCTTTTCTTACACCGCGCCGATAAAAAATATGCACGTTTATTGCATATCTTTCAGTG  >  CP009273/888897‑889192

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: