Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 647,422 G→T 92.8% intergenic (‑110/‑269) citC ← / → citA [citrate [pro‑3S]‑lyase] ligase/sensory histidine kinase in two‑component regulatory system with CitB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273647,4220GT92.8% 331.8 / 16.0 105intergenic (‑110/‑269)citC/citA[citrate [pro‑3S]‑lyase] ligase/sensory histidine kinase in two‑component regulatory system with CitB
Reads supporting (aligned to +/- strand):  ref base G (4/4);  new base T (42/55);  total (46/59)
Fisher's exact test for biased strand distribution p-value = 7.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 105 total aligned reads displayed.

AACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGAGGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTTTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCAGG  >  CP009273/647276‑647566
                                                                                                                                                  |                                                                                                                                                
aaCGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGta                                                                                                                                               >  1:4706662/1‑150 (MQ=255)
aaCGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGta                                                                                                                                               <  1:16766/150‑1 (MQ=255)
  cGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGtata                                                                                                                                             <  1:70661/150‑1 (MQ=255)
   gTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTatat                                                                                                                                            <  2:1507980/150‑1 (MQ=255)
   gTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTatat                                                                                                                                            <  2:1508033/150‑1 (MQ=255)
   gTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTatat                                                                                                                                            <  2:1507996/150‑1 (MQ=255)
         cgcgGGTGAAAATATCATTGCCGAACATAATAAATAGTATACTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCttttt                                                                                                                                      <  1:4067804/150‑1 (MQ=255)
          gcgGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTtc                                                                                                                                     <  1:4751465/150‑1 (MQ=255)
          gcgGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTtc                                                                                                                                     <  1:3959391/150‑1 (MQ=255)
           cgGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTtct                                                                                                                                    >  1:666641/1‑150 (MQ=255)
           cgGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTtct                                                                                                                                    <  2:2706658/150‑1 (MQ=255)
                gAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCttt                                                                                                                               <  2:1593922/150‑1 (MQ=255)
                   aaTATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGc                                                                                                                            >  2:4488717/1‑150 (MQ=255)
                            gCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATcc                                                                                                                   <  1:1880178/150‑1 (MQ=255)
                              cGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTa                                                                                                                 <  2:4273450/150‑1 (MQ=255)
                               gAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTaa                                                                                                                >  1:3261944/1‑150 (MQ=255)
                                aaCATAATAAATAGTATCCTGAAGGTGAATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTTATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAg                                                                                                               <  2:3936932/150‑1 (MQ=255)
                                 aCATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAgg                                                                                                              <  2:1527082/150‑1 (MQ=255)
                                 aCATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGAGGTATATCTTTTTCTCCTTTCGCTATACATCCTAAgg                                                                                                              <  1:996948/150‑1 (MQ=255)
                                 aCAGAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATGTGCAACGACTGGTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTGTTCTCCTTTCGCTATACATCCTAAgg                                                                                                              <  2:2660977/150‑1 (MQ=255)
                                    taataaATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGt                                                                                                           >  1:2367053/1‑150 (MQ=255)
                                    taataaATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGt                                                                                                           >  1:2367179/1‑150 (MQ=255)
                                        aaataGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTAttt                                                                                                       <  2:2205446/150‑1 (MQ=255)
                                         aataGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTc                                                                                                      >  1:434184/1‑150 (MQ=255)
                                         aataGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTc                                                                                                      <  2:4269883/150‑1 (MQ=255)
                                          ataGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCg                                                                                                     <  2:1039687/150‑1 (MQ=255)
                                          ataGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCg                                                                                                     <  2:1039650/150‑1 (MQ=255)
                                              tATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGt                                                                                                 >  2:1330947/1‑150 (MQ=255)
                                               aTCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTg                                                                                                <  2:2153787/150‑1 (MQ=255)
                                               aTCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGAGGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTg                                                                                                >  1:3138563/1‑150 (MQ=255)
                                               aTCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGAGGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTg                                                                                                >  2:1298475/1‑150 (MQ=255)
                                                 ccTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGaa                                                                                              >  2:2103720/1‑150 (MQ=255)
                                                      aGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAAtttt                                                                                         >  1:2747626/1‑150 (MQ=255)
                                                       ggTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTg                                                                                        <  1:1902648/150‑1 (MQ=255)
                                                        gTGCATGTTGTGATCGATTTGCAACGAATGTTGTTCAATGGTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGa                                                                                       <  2:1187734/150‑1 (MQ=255)
                                                          gCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATGTCGGCGTGAAATTTTGatt                                                                                     >  2:4370661/1‑150 (MQ=255)
                                                             tgttgtTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGAGGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGatttat                                                                                  >  2:3648838/1‑150 (MQ=255)
                                                                     tCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGAGGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATa                                                                          <  1:1348886/150‑1 (MQ=255)
                                                                        aTTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATagag                                                                       <  1:1813578/150‑1 (MQ=255)
                                                                         tttGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGt                                                                      <  1:850303/150‑1 (MQ=255)
                                                                          ttGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGtt                                                                     <  1:1981417/150‑1 (MQ=255)
                                                                              aaCGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTg                                                                 <  2:1081501/150‑1 (MQ=255)
                                                                              aaCGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGAGGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTg                                                                 >  1:3643352/1‑150 (MQ=255)
                                                                               aCGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTgg                                                                <  2:2541419/150‑1 (MQ=255)
                                                                                  aaTGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGttt                                                             <  1:1109624/150‑1 (MQ=255)
                                                                                       tgttCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTtctat                                                        >  1:4561952/1‑150 (MQ=255)
                                                                                       tgttCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTtctat                                                        >  1:2181034/1‑150 (MQ=255)
                                                                                       tgttCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTtctat                                                        >  1:2181005/1‑150 (MQ=255)
                                                                                         ttCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTtctattt                                                      >  1:812155/1‑150 (MQ=255)
                                                                                         ttCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTtctattt                                                      <  1:2983113/150‑1 (MQ=255)
                                                                                         ttCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTtctattt                                                      <  1:1582303/150‑1 (MQ=255)
                                                                                         ttCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTtctattt                                                      <  1:1582058/150‑1 (MQ=255)
                                                                                           cAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTtttatttat                                                    <  1:4625674/150‑1 (MQ=255)
                                                                                           cAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTtctatttat                                                    <  1:506533/150‑1 (MQ=255)
                                                                                           cAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTtctatttat                                                    <  1:506564/150‑1 (MQ=255)
                                                                                           cAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTtctatttat                                                    <  1:2958942/150‑1 (MQ=255)
                                                                                                ttGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAAt                                               <  1:1005770/150‑1 (MQ=255)
                                                                                                ttGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAAt                                               <  2:765103/150‑1 (MQ=255)
                                                                                                   cAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGAt                                            >  1:859568/1‑150 (MQ=255)
                                                                                                   cAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGAt                                            >  1:441999/1‑150 (MQ=255)
                                                                                                    aaaCTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGAtt                                           <  1:1023477/150‑1 (MQ=255)
                                                                                                    aaaCTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGAtt                                           >  1:4403497/1‑150 (MQ=255)
                                                                                                    aaaCTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGAtt                                           <  1:3409602/150‑1 (MQ=255)
                                                                                                     aaCTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGAttt                                          <  1:4488717/150‑1 (MQ=255)
                                                                                                        tGATAACCTTTTATTTTCTCTTGGGAGACAGGGTGTGCTCGAGGTATATCTTTTTATCCTTTCGCTAGACATCCTAAGGTGTATTTCGGCGTGATATTTTGATTTATTTCAAATTGATTTAGTGGTTTTTTATTTATTTAATGATTTTaa                                       <  2:3278834/150‑1 (MQ=255)
                                                                                                        tGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTTTATTTATTTAATGATTTTaa                                       <  2:906213/150‑1 (MQ=255)
                                                                                                        tGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTTTATTTATTTAATGATTTTaa                                       <  2:2561528/150‑1 (MQ=255)
                                                                                                        tGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTaa                                       >  1:604201/1‑150 (MQ=255)
                                                                                                          aTAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTTATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGt                                     >  2:2986067/1‑150 (MQ=255)
                                                                                                          aTAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGt                                     >  2:614165/1‑150 (MQ=255)
                                                                                                           tAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGtt                                    >  2:3706588/1‑150 (MQ=255)
                                                                                                            aaCCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGttt                                   <  2:2957274/150‑1 (MQ=255)
                                                                                                            aaCCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGttt                                   >  1:3407788/1‑150 (MQ=255)
                                                                                                             aCCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGtttt                                  <  2:1788783/150‑1 (MQ=255)
                                                                                                             aCCTTNTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGtttt                                  >  1:199083/1‑150 (MQ=255)
                                                                                                                 tttATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTttaa                              >  1:548051/1‑150 (MQ=255)
                                                                                                                 tttATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTttaa                              >  1:549416/1‑150 (MQ=255)
                                                                                                                  ttATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTttaat                             >  1:2116861/1‑150 (MQ=255)
                                                                                                                  ttATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTttaat                             >  1:2117196/1‑150 (MQ=255)
                                                                                                                   tATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTttaatt                            >  2:4691315/1‑150 (MQ=255)
                                                                                                                   tATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTttaatt                            >  2:353338/1‑150 (MQ=255)
                                                                                                                   tATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTACACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTttaatt                            >  2:352605/1‑150 (MQ=255)
                                                                                                                     ttttCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTttaattaa                          >  1:780348/1‑150 (MQ=255)
                                                                                                                     ttttCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTttaattaa                          <  1:2762599/150‑1 (MQ=255)
                                                                                                                     ttttCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTttaattaa                          <  1:2825397/150‑1 (MQ=255)
                                                                                                                       ttCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATg                        >  2:2369358/1‑150 (MQ=255)
                                                                                                                       ttCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATg                        >  2:4749874/1‑150 (MQ=255)
                                                                                                                       gtCACATGGGCTAAATGGGGTTAACGATGTATATACTTTTCTCCTTTCCCTAAACATCCTAATGATTATTTCGGCGTGAAATTTTGATTTAATTCACATAGATTTAGTTGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATg                        <  2:4240167/149‑1 (MQ=255)
                                                                                                                         cACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTa                      <  2:4463743/150‑1 (MQ=255)
                                                                                                                         cACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTa                      <  2:694864/150‑1 (MQ=255)
                                                                                                                              gggAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTAc                 <  2:2367053/150‑1 (MQ=255)
                                                                                                                              gggAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTAc                 <  2:2367179/150‑1 (MQ=255)
                                                                                                                                 agaAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGaa              >  2:13922/1‑150 (MQ=255)
                                                                                                                                   aaaGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAAt            >  2:4450850/1‑150 (MQ=255)
                                                                                                                                   aaaGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTAAGAAAt            >  2:4450851/1‑150 (MQ=255)
                                                                                                                                    aaGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATg           <  1:3703225/150‑1 (MQ=255)
                                                                                                                                    aaGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATg           <  1:1046574/150‑1 (MQ=255)
                                                                                                                                      gggggTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTATGTTTTTTATTTAATGTAATTACGAAATGAc         >  2:583794/1‑150 (MQ=255)
                                                                                                                                         ggTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCg      >  1:1563013/1‑150 (MQ=255)
                                                                                                                                             aTCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCAgg  >  1:3274214/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
AACGTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGAGGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTTTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCAGG  >  CP009273/647276‑647566

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: