Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,265,089 T→A 100% intergenic (‑120/‑310) yjcB ← / → yjcC putative inner membrane protein/putative membrane‑anchored cyclic‑di‑GMP phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,265,0890TA96.3% 685.9 / 10.7 215intergenic (‑120/‑310)yjcB/yjcCputative inner membrane protein/putative membrane‑anchored cyclic‑di‑GMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base T (1/7);  new base A (107/100);  total (108/107)
Fisher's exact test for biased strand distribution p-value = 3.52e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01
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 221 total aligned reads displayed.

ACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTCTGCATTCCTGTAGGAAATTA  >  CP009273/4264944‑4265232
                                                                                                                                                 |                                                                                                                                               
     gCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTaa                                                                                                                                          <  2:2188445/148‑1 (MQ=255)
     gCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTaa                                                                                                                                          >  2:696765/1‑148 (MQ=255)
     gCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTaa                                                                                                                                          <  1:1960932/148‑1 (MQ=255)
        aGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAAttt                                                                                                                                       <  1:3481327/148‑1 (MQ=255)
         gtggtgAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCCCGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAAtttt                                                                                                                                      <  2:2529891/148‑1 (MQ=255)
          tggtgAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTtg                                                                                                                                     >  2:1423712/1‑148 (MQ=255)
             tgAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGa                                                                                                                                  >  2:3044601/1‑148 (MQ=255)
             tgAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGa                                                                                                                                  >  2:2433703/1‑148 (MQ=255)
              gAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGaa                                                                                                                                 <  2:2564576/148‑1 (MQ=255)
                gggTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACt                                                                                                                               <  2:2914749/148‑1 (MQ=255)
                   tAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACtata                                                                                                                            <  2:3241065/148‑1 (MQ=255)
                   tAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACtata                                                                                                                            <  2:2895583/148‑1 (MQ=255)
                   tAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACtata                                                                                                                            >  1:807548/1‑148 (MQ=255)
                      ccATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATca                                                                                                                         >  1:2551800/1‑148 (MQ=255)
                        aTCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATcaca                                                                                                                       <  1:567960/148‑1 (MQ=255)
                         tCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACaa                                                                                                                      <  2:431442/148‑1 (MQ=255)
                         tCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACaa                                                                                                                      <  1:857075/148‑1 (MQ=255)
                         tCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACaa                                                                                                                      >  1:2113534/1‑148 (MQ=255)
                          cAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAat                                                                                                                     >  1:2878051/1‑148 (MQ=255)
                          cAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAat                                                                                                                     <  1:636535/148‑1 (MQ=255)
                           aaaCAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAatt                                                                                                                    >  2:1113901/1‑148 (MQ=255)
                            aaCAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAattg                                                                                                                   >  2:2068773/1‑148 (MQ=255)
                            aaCAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAattg                                                                                                                   <  2:1557458/148‑1 (MQ=255)
                             aCAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAattaa                                                                                                                  >  1:3531338/1‑146 (MQ=255)
                               aaTTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAattgatt                                                                                                                >  1:3072138/1‑148 (MQ=255)
                                aTTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAattgattg                                                                                                               >  2:2252109/1‑148 (MQ=255)
                                aTTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAattgattg                                                                                                               >  2:47416/1‑148 (MQ=255)
                                 ttCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGAttgt                                                                                                              <  1:2734112/148‑1 (MQ=255)
                                   ccTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGAttgttt                                                                                                            <  1:2520114/148‑1 (MQ=255)
                                      ctGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATtgtttgtt                                                                                                         <  1:1208674/148‑1 (MQ=255)
                                         aaaaGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGc                                                                                                      <  1:2713865/148‑1 (MQ=255)
                                         aaaaGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGc                                                                                                      >  2:553023/1‑148 (MQ=255)
                                           aaGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCa                                                                                                    >  1:8724/1‑148 (MQ=255)
                                           aaGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCa                                                                                                    >  1:1609361/1‑148 (MQ=255)
                                            agagCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCat                                                                                                   <  2:1562328/148‑1 (MQ=255)
                                            agagCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCat                                                                                                   >  1:2862311/1‑148 (MQ=255)
                                             gagCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCata                                                                                                  >  1:1719597/1‑148 (MQ=255)
                                              agCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCatat                                                                                                 >  1:2301653/1‑148 (MQ=255)
                                              agCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCatat                                                                                                 >  1:2121076/1‑148 (MQ=255)
                                              agCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCatat                                                                                                 >  1:908012/1‑148 (MQ=255)
                                              agCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCatat                                                                                                 >  1:2937500/1‑148 (MQ=255)
                                                ccGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAAGTGATTGTTTGTTAGCCATATTa                                                                                               >  2:3366701/1‑148 (MQ=255)
                                                  gATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAgg                                                                                             >  1:2260786/1‑148 (MQ=255)
                                                  gATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAgg                                                                                             <  1:293398/148‑1 (MQ=255)
                                                   aTGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGc                                                                                            >  1:1387477/1‑148 (MQ=255)
                                                     gCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCg                                                                                          >  1:1287688/1‑148 (MQ=255)
                                                        cTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGa                                                                                       <  1:3441204/148‑1 (MQ=255)
                                                         tGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGAc                                                                                      >  1:1767231/1‑148 (MQ=255)
                                                         tGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGAc                                                                                      >  1:1499559/1‑148 (MQ=255)
                                                          ggACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACt                                                                                     <  2:1979895/148‑1 (MQ=255)
                                                          ggACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATGAGGCCGTGACt                                                                                     >  1:2922336/1‑148 (MQ=255)
                                                           gACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAAttt                                                                                                                                       <  2:2549558/97‑1 (MQ=255)
                                                           gACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTATGCCGTGACtt                                                                                    >  1:938306/1‑148 (MQ=255)
                                                            aCAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACttt                                                                                   >  1:1977349/1‑148 (MQ=255)
                                                            aCAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACttt                                                                                   <  2:1207971/148‑1 (MQ=255)
                                                            aCAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACttt                                                                                   <  2:3137842/148‑1 (MQ=255)
                                                             cAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACtttt                                                                                  >  2:2384740/1‑148 (MQ=255)
                                                               gcgcGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTAt                                                                                <  2:1922748/148‑1 (MQ=255)
                                                                  cggcgTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGc                                                                             >  2:1236061/1‑148 (MQ=255)
                                                                        tgttCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACa                                                                       >  2:2712435/1‑148 (MQ=255)
                                                                         gttCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAg                                                                      <  2:3485456/148‑1 (MQ=255)
                                                                           tCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGAt                                                                    >  2:3089777/1‑148 (MQ=255)
                                                                             aCGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTa                                                                  <  1:730462/148‑1 (MQ=255)
                                                                                ggAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTAtgt                                                               >  2:204243/1‑148 (MQ=255)
                                                                                  aGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTAtgttg                                                             <  2:737173/148‑3 (MQ=255)
                                                                                  aGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTgg                                                             <  1:2612770/148‑1 (MQ=255)
                                                                                  aGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTgg                                                             <  1:1348936/148‑1 (MQ=255)
                                                                                     cGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGttt                                                          <  2:18461/148‑1 (MQ=255)
                                                                                        ggTACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTc                                                       <  2:3074764/148‑1 (MQ=255)
                                                                                          tACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAg                                                     <  1:2019307/148‑1 (MQ=255)
                                                                                          tACACTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAg                                                     <  1:1423712/148‑1 (MQ=255)
                                                                                             acTCTGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGATGTACAGATTATGTGGTTTTTCAGTgg                                                  >  2:3325983/1‑148 (MQ=255)
                                                                                                ctGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCAt                                               >  1:3428820/1‑148 (MQ=255)
                                                                                                ctGGCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCAt                                               <  2:234111/148‑1 (MQ=255)
                                                                                                  ggCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTa                                             <  2:3499901/148‑1 (MQ=255)
                                                                                                   gCTATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTaa                                            >  2:704065/1‑148 (MQ=255)
                                                                                                     tATCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAgg                                          <  2:508550/148‑1 (MQ=255)
                                                                                                      aTCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAggg                                         >  2:1636049/1‑148 (MQ=255)
                                                                                                      aTCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAAGTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAggg                                         >  2:3409969/1‑148 (MQ=255)
                                                                                                       tCGCGGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGc                                        >  1:382240/1‑148 (MQ=255)
                                                                                                        cgcgGGGCTTGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCa                                       >  2:1138510/1‑148 (MQ=255)
                                                                                                                ttGCAGAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTcc                               >  1:3014624/1‑148 (MQ=255)
                                                                                                                     gAAAACAAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCccgccg                          >  1:3114056/1‑148 (MQ=255)
                                                                                                                      aaaacaAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCcgccgc                         <  1:661054/145‑1 (MQ=255)
                                                                                                                       aaacaAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCgccgcc                        >  1:2155100/3‑148 (MQ=255)
                                                                                                                          caAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCctct                     >  1:240975/1‑148 (MQ=255)
                                                                                                                          caAAAATGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCctct                     >  1:3364371/1‑148 (MQ=255)
                                                                                                                              aaTGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTCTGCAt                 >  1:3513475/1‑148 (MQ=255)
                                                                                                                               aTGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTCTGCAtt                >  2:742118/1‑148 (MQ=255)
                                                                                                                               aTGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTCTGCAtt                >  2:2733896/1‑148 (MQ=255)
                                                                                                                               aTGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTCTGCAtt                >  1:1920277/1‑148 (MQ=255)
                                                                                                                               aTGAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTCTGCAtt                >  1:1921555/1‑148 (MQ=255)
                                                                                                                                 gAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTCTGCATTcc              <  1:3119152/148‑1 (MQ=255)
                                                                                                                                 gAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTCTGCATTcc              <  1:1188815/148‑1 (MQ=255)
                                                                                                                                 gAAACACTCTGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTCTGCATTcc              <  1:3227833/148‑1 (MQ=255)
                                                                                                                                         ctGTTTGTATCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTCTGCATTCCTGTAGGaa      <  2:1719597/148‑1 (MQ=255)
                                                                                                                                          tgtttgtaTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTCTGCATTCCTGTAGGaaa     >  2:1692153/1‑148 (MQ=255)
                                                                                                                                           gtttgtatCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTCTGCATTCCTGTAGGAaat    >  1:3337978/1‑148 (MQ=255)
                                                                                                                                            tttgtatCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTCTGCATTCCTGTAGGAaatt   >  2:622977/1‑148 (MQ=255)
                                                                                                                                             ttgtatCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTCTGCATTCCTGTAGGAaatta  <  2:1333374/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
ACCACGCCAGTGGTGAGGGTAGCCATCAAACAATTCCTCTGAAAAGAGCCGATGCCCTGGACAGCGCGGCGTTGTTCACGGGAGGCGTGGTACACTCTGGCTATCGCGGGGCTTGCAGAACACAAAAATGAAACACTCTGTTTGTTTCATTAATTTTGTGAACTATATCACAATTGATTGTTTGTTAGCCATATTAGGCCGTGACTTTTATTGCTGTACAGATTATGTGGTTTTTCAGTGGCATTAAGGGCATATCTTCCCGCCGCCTCTGCATTCCTGTAGGAAATTA  >  CP009273/4264944‑4265232

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: