Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,075,114 C→A 100% intergenic (‑120/‑158) tktA ← / → loiP transketolase 1, thiamine triphosphate‑binding/Phe‑Phe periplasmic metalloprotease, OM lipoprotein; low salt‑inducible; heat shock protein that binds Era

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,075,1140CA95.2% 861.8 / 30.1 272intergenic (‑120/‑158)tktA/loiPtransketolase 1, thiamine triphosphate‑binding/Phe‑Phe periplasmic metalloprotease, OM lipoprotein; low salt‑inducible; heat shock protein that binds Era
Reads supporting (aligned to +/- strand):  ref base C (7/6);  new base A (123/136);  total (130/142)
Fisher's exact test for biased strand distribution p-value = 7.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 272 total aligned reads displayed.

CATTGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAACATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTGCAAACGGCGCATTTTTAGAATAATCCTGACCTTGTGCG  >  CP009273/3074967‑3075260
                                                                                                                                                   |                                                                                                                                                  
cATTGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAAt                                                                                                                                                  <  2:1147016/150‑1 (MQ=255)
cATTGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAAt                                                                                                                                                  >  2:53215/1‑150 (MQ=255)
cATTGGCAAGCTCGTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCTCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAAt                                                                                                                                                  >  2:3659260/1‑150 (MQ=255)
 aTTGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATc                                                                                                                                                 >  1:55363/1‑150 (MQ=255)
    ggCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAACATCCtt                                                                                                                                              <  2:3144817/150‑1 (MQ=255)
    ggCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCtt                                                                                                                                              >  2:1886403/1‑150 (MQ=255)
     gCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCttt                                                                                                                                             >  1:2660719/1‑150 (MQ=255)
       aaGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTg                                                                                                                                           <  1:192099/150‑1 (MQ=255)
       aaGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTg                                                                                                                                           >  1:2109879/1‑150 (MQ=255)
       aaGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTg                                                                                                                                           <  2:1678783/150‑1 (MQ=255)
          ctctTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCg                                                                                                                                        <  1:201937/150‑1 (MQ=255)
          ctctTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCg                                                                                                                                        <  1:4317694/150‑1 (MQ=255)
            ctTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCt                                                                                                                                      >  1:4227555/1‑150 (MQ=255)
                aCGTGTGGACATTTGGACTCCAGCTCGGATGATGAAGGGCTCGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAg                                                                                                                                  <  1:2314410/150‑1 (MQ=255)
                aCGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAg                                                                                                                                  >  1:2890357/1‑150 (MQ=255)
                 cGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTTCATGGAAAATCCTTTTGCCGCTTCAGa                                                                                                                                 <  1:3401348/150‑1 (MQ=255)
                  gTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGaa                                                                                                                                >  2:326149/1‑150 (MQ=255)
                   tGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGaca                                                                                                                               >  2:3822814/1‑148 (MQ=255)
                       gACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAAtctc                                                                                                                           >  2:3276193/1‑150 (MQ=255)
                        aCATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATctct                                                                                                                          <  2:4154887/150‑1 (MQ=255)
                          tttttGACTCCAGATCTGATGATGAAGGGCTCGCCCTTATCGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAATATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTgg                                                                                                                        <  1:341954/149‑1 (MQ=255)
                          aTTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTgg                                                                                                                        <  1:2177252/150‑1 (MQ=255)
                          aTTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTgg                                                                                                                        >  1:3928182/1‑150 (MQ=255)
                           ttttGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGa                                                                                                                       <  1:2605194/150‑1 (MQ=255)
                           ttttGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGa                                                                                                                       <  1:1932605/150‑1 (MQ=255)
                               gACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCTGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCAt                                                                                                                   <  2:1112889/150‑1 (MQ=255)
                               gACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCAt                                                                                                                   >  2:3177019/1‑150 (MQ=255)
                                aCTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATg                                                                                                                  <  2:7242/150‑1 (MQ=255)
                                      gATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCa                                                                                                            <  2:1067435/150‑1 (MQ=255)
                                        tCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAACATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATg                                                                                                          <  2:3994617/150‑1 (MQ=255)
                                          ggATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGtt                                                                                                        <  1:1312103/150‑1 (MQ=255)
                                          ggATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGtt                                                                                                        >  1:4189458/1‑150 (MQ=255)
                                           gatgatGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTg                                                                                                       <  2:82678/150‑1 (MQ=255)
                                           gatgatGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTg                                                                                                       <  2:81076/150‑1 (MQ=255)
                                             tgatgaAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTgcg                                                                                                     >  1:3705970/1‑150 (MQ=255)
                                             tgatgaAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTgcg                                                                                                     >  1:626735/1‑150 (MQ=255)
                                                       cACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTc                                                                                           >  1:3654579/1‑150 (MQ=255)
                                                         cGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGc                                                                                         >  1:2021458/1‑150 (MQ=255)
                                                                aaCGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCg                                                                                  >  1:2671828/1‑150 (MQ=255)
                                                                 aCGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGGAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGc                                                                                 >  1:1843334/1‑150 (MQ=255)
                                                                 aCGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGc                                                                                 >  1:4091574/1‑150 (MQ=255)
                                                                   gACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTg                                                                               <  1:3986853/150‑1 (MQ=255)
                                                                    aCTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTgc                                                                              >  1:301422/1‑150 (MQ=255)
                                                                      ttGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAACATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTgcgc                                                                            >  1:1770074/1‑150 (MQ=255)
                                                                      ttGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTgcgc                                                                            >  2:3681400/1‑150 (MQ=255)
                                                                         acgacAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCttt                                                                         <  2:107275/150‑1 (MQ=255)
                                                                                cgcgTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTtat                                                                  <  2:4361734/150‑1 (MQ=255)
                                                                                cgcgTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTtat                                                                  <  2:4361768/150‑1 (MQ=255)
                                                                                    ttttGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTGTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATactg                                                              <  2:2622067/150‑1 (MQ=255)
                                                                                     tttGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATactga                                                             >  2:3534328/1‑150 (MQ=255)
                                                                                         ggCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATactgagact                                                         >  1:736046/1‑150 (MQ=255)
                                                                                         ggCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATactgagact                                                         <  2:2006262/150‑1 (MQ=255)
                                                                                          gCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATactgagactg                                                        <  2:1033065/150‑1 (MQ=255)
                                                                                               gCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGt                                                   <  1:3714006/150‑1 (MQ=255)
                                                                                                 cGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAACATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCg                                                 <  1:3217893/150‑1 (MQ=255)
                                                                                                 cGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAACATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCg                                                 >  2:1407571/1‑150 (MQ=255)
                                                                                                  ggAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGa                                                >  1:1405764/1‑150 (MQ=255)
                                                                                                    aaaaTTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGAtt                                              >  2:4408858/1‑150 (MQ=255)
                                                                                                    aaaaTTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGAtt                                              <  1:3905493/150‑1 (MQ=255)
                                                                                                     aaaTTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTc                                             >  2:2666693/1‑150 (MQ=255)
                                                                                                          tGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTg                                        <  2:2109879/150‑1 (MQ=255)
                                                                                                           gCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTGc                                       >  1:675668/1‑150 (MQ=255)
                                                                                                            ccAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTGCa                                      >  1:3357737/1‑150 (MQ=255)
                                                                                                            ccAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTGCa                                      >  1:1083317/1‑150 (MQ=255)
                                                                                                               acaaTTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTGCAAAc                                   >  2:3320132/1‑150 (MQ=255)
                                                                                                                caaTTTACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTGCAAACg                                  <  1:3805958/150‑1 (MQ=255)
                                                                                                                   tttACCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTGCAAACGgcg                               <  2:4250009/150‑1 (MQ=255)
                                                                                                                      aCCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTGCAAACGGCGCAt                            >  2:2022331/1‑150 (MQ=255)
                                                                                                                      aCCGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTGCAAACGGCGCAt                            >  2:2359199/1‑150 (MQ=255)
                                                                                                                       ccGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTGCAAACGGCGCAtt                           <  1:3534328/150‑1 (MQ=255)
                                                                                                                       ccGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTGCAAACGGCGCAtt                           <  1:759294/150‑1 (MQ=255)
                                                                                                                        cGCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTGCAAACGGCGCAttt                          >  2:1275278/1‑150 (MQ=255)
                                                                                                                         gCAAGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTGCAAACGGCGCAtttt                         <  2:797905/150‑1 (MQ=255)
                                                                                                                           aaGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTGCAAACGGCGCATTTTTa                       <  2:2660719/150‑1 (MQ=255)
                                                                                                                            aGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTGCAAACGGCGCATTTTTAg                      >  1:1739268/1‑150 (MQ=255)
                                                                                                                            aGCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTGCAAACGGCGCATTTTTAg                      <  2:1590054/150‑1 (MQ=255)
                                                                                                                             gCCGCGCGTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTAATCGCCAGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTGCAAACGGCGCATTTTTAGa                     >  2:3168084/1‑150 (MQ=255)
                                                                                                                               cgcgcgTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAATTCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTGCAAACGGCGCATTTTTAGaat                   <  1:4472765/150‑1 (MQ=255)
                                                                                                                               cgcgcgTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTGCAAACGGCGCATTTTTAGaat                   <  1:2285064/150‑1 (MQ=255)
                                                                                                                                gcgcgTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTGCAAACGGCGCATTTTTAGaata                  >  2:2944890/1‑150 (MQ=255)
                                                                                                                                gcgcgTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTGCAAACGGCGCATTTTTAGaata                  <  1:3526918/150‑1 (MQ=255)
                                                                                                                                gcgcgTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTGCAAACGGCGCATTTTTAGaata                  >  2:4497816/1‑150 (MQ=255)
                                                                                                                                gcgcgTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTGCAAACGGCGCATTTTTAGaata                  <  2:4368219/150‑1 (MQ=255)
                                                                                                                                gcgcgTAATGTACATTGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTAGCATGTTGCGCAATATACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTGAAAACGGCGCATTTTTAGaata                  <  1:4449640/150‑1 (MQ=255)
                                                                                                                                  gcgTCATGTACATGGAAAATCCTATTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTGCAAACGGCGCATTTTTAGaataat                <  1:1604387/150‑1 (MQ=255)
                                                                                                                                   cgTCATGTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTGCAAACGGCGCATTTTTAGAATAATc               <  1:246894/150‑1 (MQ=255)
                                                                                                                                         gTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTGCAAACGGCGCATTTTTAGAATAATCCTGAcc         <  1:1125142/150‑1 (MQ=255)
                                                                                                                                         gTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTGCAAACGGCGCATTTTTAGAATAATCCTGAcc         <  1:2525503/150‑1 (MQ=255)
                                                                                                                                         gTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTGCAAACGGCGCATTTTTAGAATAATCCTGAcc         >  2:40873/1‑150 (MQ=255)
                                                                                                                                         gTACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTGCAAACGGCGCATTTTTAGAATAATCCTGAcc         <  1:3400066/150‑1 (MQ=255)
                                                                                                                                          tACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTGCAAACGGCGCATTTTTAGAATAATCCTGACCt        >  2:2014276/1‑150 (MQ=255)
                                                                                                                                          tACATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTGCAAACGGCGCATTTTTAGAATAATCCTGACCt        >  2:2014231/1‑150 (MQ=255)
                                                                                                                                           aCATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTGCAAACGGCGCATTTTTAGAATAATCCTGACCtt       <  2:3705970/150‑1 (MQ=255)
                                                                                                                                           aCATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTGCAAACGGCGCATTTTTAGAATAATCCTGACCtt       <  2:626735/150‑1 (MQ=255)
                                                                                                                                           aCATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTGCAAACGGCGCATTTTTAGAATAATCCTGACCtt       <  1:3948625/150‑1 (MQ=255)
                                                                                                                                            cATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTGCAAACGGCGCATTTTTAGAATAATCCTGACCTtg      <  1:1785678/150‑1 (MQ=255)
                                                                                                                                            cATGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTGCAAACGGCGCATTTTTAGAATAATCCTGACCTtg      <  1:1785621/150‑1 (MQ=255)
                                                                                                                                             aTGGAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTGCAAACGGCGCATTTTTAGAATAATCCTGACCTtgt     >  1:983545/1‑150 (MQ=255)
                                                                                                                                               ggAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTGCAAACGGCGCATTTTTAGAATAATCCTGACCTTGTGc   <  1:2864490/150‑1 (MQ=255)
                                                                                                                                                gAAAATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTGCAAACGGCGCATTTTTAGAATAATCCTGACCTTGTGCg  <  1:2652382/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CATTGGCAAGCTCTTTACGTGAGGACATTTTGACTCCAGATCGGATGATGAAGGGCACGCCCTTAACGACTTGACGACAGCGCGTTTTGGGCTACGCCGGAAAATTTGCCAACAATTTACCGCAAGCCGCGCGTCATGTACATGGAACATCCTTTTGCCGCTTCAGAAATCTCTGGATCATGCTCGCATGTTGCGCAATCTACTCGCCCGTCCGCTGCGCTTTTCCTTATACTGAGACTGAGCGTCGATTCACCTGCAAACGGCGCATTTTTAGAATAATCCTGACCTTGTGCG  >  CP009273/3074967‑3075260

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: