Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 647,422 G→T 100% 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,4220GT98.5% 662.5 / ‑0.1 194intergenic (‑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 (3/0);  new base T (80/111);  total (83/111)
Fisher's exact test for biased strand distribution p-value = 7.67e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 194 total aligned reads displayed.

GTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGAGGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTTTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCAGGTT  >  CP009273/647279‑647568
                                                                                                                                               |                                                                                                                                                  
gTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTatat                                                                                                                                              >  2:3136725/1‑150 (MQ=255)
  tttACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCt                                                                                                                                            <  2:3448990/150‑1 (MQ=255)
  tttACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCt                                                                                                                                            >  1:342957/1‑150 (MQ=255)
    tACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCttt                                                                                                                                          >  2:3009232/1‑150 (MQ=255)
    tACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCttt                                                                                                                                          >  1:3904152/1‑150 (MQ=255)
     aCGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCtttt                                                                                                                                         <  2:789622/150‑1 (MQ=255)
       gcgGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTtc                                                                                                                                       <  2:1865145/150‑1 (MQ=255)
         gggTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGTAAGGGGGTGATCGATGTATATCTTTTtctc                                                                                                                                     <  1:3859959/150‑1 (MQ=255)
            tGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCtt                                                                                                                                  <  1:3425463/150‑1 (MQ=255)
             gAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCttt                                                                                                                                 <  1:775953/150‑1 (MQ=255)
             gAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCttt                                                                                                                                 <  2:3274367/150‑1 (MQ=255)
              aaaaTATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTc                                                                                                                                <  2:42412/150‑1 (MQ=255)
              aaaaTATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTc                                                                                                                                <  2:564251/150‑1 (MQ=255)
              aaaaTATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGAGGTATATCTTTTTCTCCTTTc                                                                                                                                >  1:2887227/1‑150 (MQ=255)
              aaaaTATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTCCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTc                                                                                                                                <  1:140198/150‑1 (MQ=255)
               aaaTATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCg                                                                                                                               <  2:2027464/150‑1 (MQ=255)
               aaaTATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCg                                                                                                                               <  1:2038152/150‑1 (MQ=255)
               aaaTATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCg                                                                                                                               <  1:4139461/150‑1 (MQ=255)
                aaTATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGc                                                                                                                              >  1:306927/1‑150 (MQ=255)
                  tatCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCta                                                                                                                            <  1:3168226/150‑1 (MQ=255)
                   atCATTGCCGAACATAATAAATTGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCATTGTTTCAAACTGATAACCTTTTATTTTCTCTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCtat                                                                                                                           <  1:1931518/150‑1 (MQ=255)
                    tCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCtata                                                                                                                          <  1:2313458/150‑1 (MQ=255)
                      aTTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACa                                                                                                                        <  2:1504118/150‑1 (MQ=255)
                        tGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATc                                                                                                                      <  2:3835954/150‑1 (MQ=255)
                          ccGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCt                                                                                                                    <  2:3124098/150‑1 (MQ=255)
                           cGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTa                                                                                                                   >  1:2761214/1‑150 (MQ=255)
                           cGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTa                                                                                                                   <  1:2953480/150‑1 (MQ=255)
                            gAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTaa                                                                                                                  >  1:1560473/1‑150 (MQ=255)
                               cATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGa                                                                                                               >  1:1181721/1‑150 (MQ=255)
                                 taataaATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGt                                                                                                             <  2:3651697/150‑1 (MQ=255)
                                   ataaataGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTAt                                                                                                           >  1:950009/1‑150 (MQ=255)
                                     aaataGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTAttt                                                                                                         >  2:239684/1‑150 (MQ=255)
                                      aataGTATCCTGATGGTGCATGTTGTTATCGATTTGCTACGAATGTTGTTCATTGTTGCAAACTGATAAACTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATAAATCCTAAGGAGTATTTc                                                                                                        <  2:1707572/150‑1 (MQ=255)
                                       ataGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCg                                                                                                       >  1:614597/1‑150 (MQ=255)
                                            aTCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTg                                                                                                  >  1:2053445/1‑150 (MQ=255)
                                             tCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGa                                                                                                 <  2:3148897/150‑1 (MQ=255)
                                                tGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAAt                                                                                              <  2:839764/150‑1 (MQ=255)
                                                   aGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAAtttt                                                                                           <  1:640476/150‑1 (MQ=255)
                                                   aGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAAtttt                                                                                           >  1:3271600/1‑150 (MQ=255)
                                                        cATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGattt                                                                                      >  2:3350165/1‑150 (MQ=255)
                                                        cATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGattt                                                                                      >  2:3754394/1‑150 (MQ=255)
                                                                tATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTcaca                                                                              >  2:4056846/1‑150 (MQ=255)
                                                                tATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTcaca                                                                              >  1:1582941/1‑150 (MQ=255)
                                                                tATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTcaca                                                                              <  1:1706879/150‑1 (MQ=255)
                                                                tATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGAGTTATTTcaca                                                                              >  2:1846146/1‑150 (MQ=255)
                                                                tATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAAGTTTGATTTATTTcaca                                                                              >  2:2216097/1‑150 (MQ=255)
                                                                 aTCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACAt                                                                             >  2:91695/1‑150 (MQ=255)
                                                                   cGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATag                                                                           >  1:2215826/1‑150 (MQ=255)
                                                                    gATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATaga                                                                          >  1:435643/1‑150 (MQ=255)
                                                                     aTTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATagag                                                                         <  2:1305855/150‑1 (MQ=255)
                                                                     aTTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATagag                                                                         <  2:1305862/150‑1 (MQ=255)
                                                                     aTTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATagag                                                                         >  1:1293356/1‑150 (MQ=255)
                                                                          cAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGt                                                                    >  2:2191880/1‑150 (MQ=255)
                                                                           aaCGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTg                                                                   >  2:3878906/1‑150 (MQ=255)
                                                                           aaCGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTg                                                                   >  2:1781321/1‑150 (MQ=255)
                                                                             cGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGt                                                                 <  1:2577580/150‑1 (MQ=255)
                                                                               aaTGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGttt                                                               >  2:2954905/1‑150 (MQ=255)
                                                                                aTGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGtttt                                                              >  2:3282649/1‑150 (MQ=255)
                                                                                aTGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGtttt                                                              <  2:262735/150‑1 (MQ=255)
                                                                                aTGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGtttt                                                              <  2:3033086/150‑1 (MQ=255)
                                                                                  gttgttCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGttttct                                                            <  2:4141850/150‑3 (MQ=255)
                                                                                     gttCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTtctatt                                                         <  1:2708020/150‑1 (MQ=255)
                                                                                        cAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTtctatttat                                                      <  1:3592455/150‑1 (MQ=255)
                                                                                         aaTGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTtctatttatt                                                     >  1:2906959/1‑150 (MQ=255)
                                                                                         aaTGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTtctatttatt                                                     >  1:3930519/1‑150 (MQ=255)
                                                                                         aaTGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTtctatttatt                                                     >  1:3930490/1‑150 (MQ=255)
                                                                                          aTGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTtctatttattt                                                    <  1:3814516/150‑1 (MQ=255)
                                                                                            gTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGATTTTCTATTTATTTaa                                                  >  1:3513022/1‑150 (MQ=255)
                                                                                             ttGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAAt                                                 <  1:1297382/150‑1 (MQ=255)
                                                                                              ttcAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATg                                                <  2:23501/148‑1 (MQ=255)
                                                                                              tGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATg                                                <  2:23495/150‑1 (MQ=255)
                                                                                                cAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGAt                                              <  1:1607737/150‑1 (MQ=255)
                                                                                                cAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGAt                                              <  1:1502871/150‑1 (MQ=255)
                                                                                                  aaCTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGAttt                                            >  1:1873500/1‑150 (MQ=255)
                                                                                                     tGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTaa                                         >  1:810366/1‑150 (MQ=255)
                                                                                                     tGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTaa                                         >  2:1435227/1‑150 (MQ=255)
                                                                                                      gATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAg                                        <  1:3721322/150‑1 (MQ=255)
                                                                                                      gATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAg                                        <  2:3971985/150‑1 (MQ=255)
                                                                                                         aaCCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGttt                                     >  2:903332/1‑150 (MQ=255)
                                                                                                          aCCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGtttt                                    >  2:2577669/1‑150 (MQ=255)
                                                                                                            cTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGtttttt                                  <  1:1223170/150‑1 (MQ=255)
                                                                                                             ttttATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTtta                                 <  2:539396/150‑1 (MQ=255)
                                                                                                             ttttATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGATAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTtta                                 >  2:1159638/1‑150 (MQ=255)
                                                                                                               ttATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTttaat                               <  1:1076212/150‑1 (MQ=255)
                                                                                                               ttATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTttaat                               <  2:161474/150‑1 (MQ=255)
                                                                                                                tATTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTttaatt                              <  2:2302587/150‑1 (MQ=255)
                                                                                                                 aTTTTCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTttaatta                             >  2:1544870/1‑150 (MQ=255)
                                                                                                                  ttttCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTttaattaa                            >  2:2800505/1‑150 (MQ=255)
                                                                                                                    ttCACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATg                          <  1:1048835/150‑1 (MQ=255)
                                                                                                                      cACTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTa                        <  1:11260/150‑1 (MQ=255)
                                                                                                                        cTTGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAAt                      >  2:1577093/1‑150 (MQ=255)
                                                                                                                         ttGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAAtt                     <  2:812005/150‑1 (MQ=255)
                                                                                                                          tGGGAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTTGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTa                    >  1:3340271/1‑150 (MQ=255)
                                                                                                                            ggAGAAAGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACg                  <  2:1125499/150‑1 (MQ=255)
                                                                                                                                 aaGGGGGTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATg             <  2:4155620/150‑1 (MQ=255)
                                                                                                                                    ggggTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACt          <  2:2716230/150‑1 (MQ=255)
                                                                                                                                     gggTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTc         <  2:3612065/150‑1 (MQ=255)
                                                                                                                                     gggTGATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTc         <  1:4031083/150‑1 (MQ=255)
                                                                                                                                         gATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCAg     <  1:1310608/150‑1 (MQ=255)
                                                                                                                                         gATCGATGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTCTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCAg     <  2:1406289/150‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
GTTTTACGCGGGTGAAAATATCATTGCCGAACATAATAAATAGTATCCTGAAGGTGCATGTTGTTATCGATTTGCAACGAATGTTGTTCAATGTTGCAAACTGATAACCTTTTATTTTCACTTGGGAGAAAGGGGGTGATCGAGGTATATCTTTTTCTCCTTTCGCTATACATCCTAAGGAGTATTTCGGCGTGAAATTTTGATTTATTTCACATAGAGTTAGTGGTTTTTTATTTATTTAATGATTTTAAGTTTTTTAATTAATGTAATTACGAAATGACTCGCAGGTT  >  CP009273/647279‑647568

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: