Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,266,318 G→T 100% R307L (CGT→CTT)  yjcC → putative membrane‑anchored cyclic‑di‑GMP phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,266,3180GT96.6% 538.1 / 8.2 177R307L (CGT→CTT) yjcCputative membrane‑anchored cyclic‑di‑GMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base G (2/4);  new base T (80/91);  total (82/95)
Fisher's exact test for biased strand distribution p-value = 6.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 178 total aligned reads displayed.

TTACTCTGGCTACGTATTCGACAAAACTATTTATCTCCCAAACGTAAATTGCAACGCGCCCTCGAAAAACATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACGTTGGCCTGGTGAGCAGGGGCAAATAATGAATCCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGATCTGGGCGATTACCTGGCAACACATGCAGA  >  CP009273/4266173‑4266462
                                                                                                                                                 |                                                                                                                                                
    tctGGCTACGTATTCGACAAAACTATTTATCTCCCAAACGTAAATTGCAACGCGCCCTCGAAAAACATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACTTTGGcc                                                                                                                                            <  2:562205/148‑1 (MQ=255)
     ctGGCTACGTATTCGACAAAACTATTTATCTCCCAAACGTAAATTGCAACGCGCCCTCGAAAAACATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACTTTGGCCt                                                                                                                                           <  2:233011/148‑1 (MQ=255)
     ctGGCTACGTATTCGACAAAACTATTTATCTCCCAAACGTAAATTGCAACGCGCCCTCGAAAAACATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACTTTGGCCt                                                                                                                                           >  1:1589571/1‑148 (MQ=255)
      tGGCTACGTATTCGACAAAACTATTTATCTCCCAAACGTAAATTGCAACGCGCCCTCGAAAAACATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACTTTGGCCTg                                                                                                                                          <  2:1968995/148‑1 (MQ=255)
       ggCTACGTATTCGACAAAACTATTTATCTCCCAAACGTAAATTGCAACGCGCCCTCGAAAAACATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACTTTGGCCTgg                                                                                                                                         <  2:3446998/148‑1 (MQ=255)
       ggCTACGTATTCGACAAAACTATTTATCTCCCAAACGTAAATTGCAACGCGCCCTCGAAAAACATCAACTTGGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACTTTGGCCTgg                                                                                                                                         <  2:2855928/148‑1 (MQ=255)
         cTACGTATTCGAAAAAACTATTTATCTCCCAAACGTAAATTGCAACGCGACCTCGAAAAACATCAACTTTGTCTTTATTACAAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACTTTGGCCTGGTg                                                                                                                                       <  1:1782362/148‑1 (MQ=255)
           aCGTATTCGACAAAACTATTTATCTCCCAAACGTAAATTGCAACGCGCCCTCGAAAAACATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACTTTGGCCTGGTGAg                                                                                                                                     >  1:1821362/1‑148 (MQ=255)
             gTATTCGACAAAACTATTTATCACCAAAACGTAAATTGCAACGCGCCCTCTAAAAACATCAACTTTGTCTTTATTACCAGCAAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACTTTGGCCTGGTGAGCa                                                                                                                                   <  1:82435/148‑1 (MQ=255)
              tATTCGACAAAACTATTTATCTCCCAAACGTAAATTGCAACGCGCCCTCGAAAAACATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACTTTGGCCTGGTGAGCAg                                                                                                                                  >  1:2334640/1‑148 (MQ=255)
                ttCGACAAAACTATTTATCACCCAAACGTAAATTGCAACGCGCCCTCGAAAAACATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACTTTGGCCTGGTGAGCAggg                                                                                                                                >  2:2429369/1‑148 (MQ=255)
                 tCGACAAAACTATTTATCTCCCAAACGTAAATTGCAACGCGCCCTCGAAAAACATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACTTTGGCCTGGTGAGCAgggg                                                                                                                               >  1:2340154/1‑148 (MQ=255)
                 tCGACAAAACTATTTATCTCCCAAACGTAAATTGCAACGCGCCCTCGAAAAACATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACTTTGGCCTGGTGAGCAgggg                                                                                                                               >  1:2340073/1‑148 (MQ=255)
                       aaaCTATTTATCTCCCAAACGTAAATTGCAACGCGCCCTCGAAAAACATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACTTTGGCCTGGTGAGCAGGGGCAaata                                                                                                                         <  2:1006507/148‑1 (MQ=255)
                         aCTATTTATCTCCCAAACGTAAATTGCAACGCGCCCTCGAAAAACATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACTTTGGCCTGGTGAGCAGGGGCAaataat                                                                                                                       >  2:1447839/1‑148 (MQ=255)
                          cTATTTATCTCCCAAACGTAAATTGCAACGCGCCCTCGAAAAACATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACTTTGGCCTGGTGAGCAGGGGCAAATAATg                                                                                                                      >  1:3070796/1‑148 (MQ=255)
                             tttATCTCCCAAACGTAAATTGCAACGCGCCCTCGAAAAACATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACTTTGGCCTGGTGAGCAGGGGCAAATAATGAAt                                                                                                                   <  1:2965776/148‑1 (MQ=255)
                              ttATCTCCCAAACGTAAATTGCAACGCGCCCTCGAAAAACATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAATCGTTGTTACTTTTGCATGGTGAGAAGGGGCAAATAATGAAta                                                                                                                  >  1:1530382/1‑147 (MQ=255)
                                  ctcCCAAACGTAAATTGCAACGCGCCCTCGAAAAACATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACTTTGGCCTGGTGAGCAGGGGCAAATAATGAATCCGGc                                                                                                              >  1:2396313/1‑148 (MQ=255)
                                    cccAAACGTAAATTGCAACGCGCCCTCGAAAAACATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACTTTGGCCTGGTGAGCAGGGGCAAATAATGAATCCGGCag                                                                                                            <  2:1868584/148‑1 (MQ=255)
                                    cccAAACGTAAATTGCAACGCGCCCTCGAAAAACATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACTTTGGCCTGGTGAGCAGGGGCAAATAATGAATCCGGCag                                                                                                            <  2:2527838/148‑1 (MQ=255)
                                    cccAAACGTAAATTGCAACGCGCCCTCGAAAAACATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACTTTGGCCTGGTGAGCAGGGGCAAATAATGAATCCGGCag                                                                                                            >  2:1444356/1‑148 (MQ=255)
                                      cAAACGTAAATTGCAACGCGCCCTCGAAAAACATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACTTTGGCCTGGTGAGCAGGGGCAAATAATGAATCCGGCagag                                                                                                          >  1:836319/1‑148 (MQ=255)
                                        aaCGTAAATTGCAACGCGCCCTCGAAAAACATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACTTTGGCCTGGTGAGCAGGGGCAAATAATGAATCCGGCAGAGtt                                                                                                        <  2:2990053/148‑1 (MQ=255)
                                         aCGTAAATTGCAACGCGCCCTCGAAAAACATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACTTTGGCCTGGTGAGCAGGGGCAAATAATGAATCCGGCAGAGttt                                                                                                       >  1:1825465/1‑148 (MQ=255)
                                         aCGTAAATTGCAACGCGCCCTCGAAAAACATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACTTTGGCCTGGTGAGCAGGGGCAAATAATGAATCCGGCAGAGttt                                                                                                       <  2:2891356/148‑1 (MQ=255)
                                            tAAATTGCAACGCGCCCTCGAAAAACATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACTTTGGCCTGGTGAGCAGGGGCAAATAATGAATCCGGCAGAGTTTAtt                                                                                                    <  1:1158438/148‑1 (MQ=255)
                                              aaTTGCAACGCGCCCTCGAAAAACATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACTTTGGCCTGGTGAGCAGGGGCAAATAATGAATCCGGCAGAGTTTATTcc                                                                                                  >  2:2321601/1‑148 (MQ=255)
                                              aaTTGCAACGCGCCCTCGAAAAACATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACTTTGGCCTGGTGAGCAGGGGCAAATAATGAATCCGGCAGAGTTTATTcc                                                                                                  >  2:3258854/1‑148 (MQ=255)
                                               aTTGCAACGCGCCCTCGAAAAACATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACTTTGGCCTGGTGAGCAGGGGCAAATAATGAATCCGGCAGAGTTTATTCCg                                                                                                 <  2:1253608/148‑1 (MQ=255)
                                                ttGCAACGCGCCCTCGAAAAACATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACTTTGGCCTGGTGAGCAGGGGCAAATAATGAATCCGGCAGAGTTTATTCCGc                                                                                                >  2:420344/1‑148 (MQ=255)
                                                 tGCAACGCGCCCTCGAAAAACATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACGTTGGCCTGGTGAGCAGGGGCAAATAATGAATCCGGCAGAGTTTATTCCGCt                                                                                               >  2:347350/1‑148 (MQ=255)
                                                  gCAACGCGCCCTCGAAAAACATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACTTTGGCCTGGTGAGCAGGGGCAAATAATGAATCCGGCAGAGTTTATTCCGCTg                                                                                              <  2:1770894/148‑1 (MQ=255)
                                                  gCAACGCGCCCTCGAAAAACATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACTTTGGCCTGGTGAGCAGGGGCAAATAATGAATCCGGCAGAGTTTATTCCGCTg                                                                                              <  2:3309584/148‑1 (MQ=255)
                                                      cgcgCCCTCGAAAAACATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACTTTGGCCTGGTGAGCAGGGGCAAATAATGAATCCGGCAGAGTTTATTCCGCTGGCAg                                                                                          >  1:3233397/1‑148 (MQ=255)
                                                         gcCCTCGAAAAACATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACTTTGGCCTGGTGAGCAGGGGCAAATAATGAATCCGGCAGAGTTTATTCCGCTGGCAGaaa                                                                                       <  2:3187408/148‑1 (MQ=255)
                                                          cccTCGAAAAACATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACTTTGGCCTGGTGAGCAGGGGCAAATAATGAATCCGGCAGAGTTTATTCCGCTGGCAGaaaa                                                                                      <  1:1009983/148‑1 (MQ=255)
                                                             tCGAAAAACATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACTTTGGCCTGGTGAGCAGGGGCAAATAATGAATCCGGCAGAGTTTATTCCGCTGGCAGAACagga                                                                                   >  2:2261533/1‑148 (MQ=255)
                                                              cGAAAAACATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACTTTGGCCTGGTGAGCAGGGGCAAATAATGAATCCGGCAGAGTTTATTCCGCTGGCAGAAAaggag                                                                                  <  1:1670216/148‑1 (MQ=255)
                                                               gAAAAACATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACTTTGGCCTGGTGAGCAGGGGCAAATAATGAATCCGGCAGAGTTTATTCCGCTGGCAGAAAaggagg                                                                                 <  2:585973/148‑1 (MQ=255)
                                                                aaaaaCATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACTTTGGCCTGGTGAGCAGGGGCAAATAATGAATCCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAggg                                                                                <  2:802673/148‑1 (MQ=255)
                                                                aaaaaCATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACTTTGGCCTGGTGAGCAGGGGCAAATAATGAATCCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAggg                                                                                <  2:802497/148‑1 (MQ=255)
                                                                aaaaaCATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACTTTGGCCTGGTGAGCAGGGGCAAATAATGAATCCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAggg                                                                                >  2:598156/1‑148 (MQ=255)
                                                                 aaaaCATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACTTTGGCCTGGTGAGCAGGGGCAAATAATGAATCCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAgggg                                                                               <  1:2528029/148‑1 (MQ=255)
                                                                 aaaaCATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACTTTGGCCTGGTGAGCAGGGGCAAATAATGAATCCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAgggg                                                                               <  2:565113/148‑1 (MQ=255)
                                                                 aaaaCATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACTTTGGCCTGGTGAGCAGGGGCAAATAATGAATCCGGCAGAGTTTATTCCGCTGACAGAAAAGGAgggg                                                                               >  2:2608867/1‑148 (MQ=255)
                                                                  aaaCATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACTTTGGCCTGGTGAGCAGGGGCAAATAATGAATCCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGga                                                                              >  2:2583554/1‑148 (MQ=255)
                                                                    aCATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACTTTGGCCTGGTGAGCAGGGGCAAATAATGAATCCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGgatg                                                                            >  2:2874057/1‑148 (MQ=255)
                                                                      aTCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACTTTGGCCTGGTGAGCAGGGGCAAATAATGAATCCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGgatgat                                                                          <  1:3267905/148‑1 (MQ=255)
                                                                       tCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACTTTGGCCTGGTGAGCAGGGGCAAATAATGAATCCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATa                                                                         >  1:2663241/1‑148 (MQ=255)
                                                                        cAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACTTTGGCCTGGTGAGCAGGGGCAAATAATGAATCCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAg                                                                        >  1:3535807/1‑148 (MQ=255)
                                                                         aaCTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACTTTGGCCTGGTGAGCAGGGGCAAATAATGAATCCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGa                                                                       <  1:2888732/148‑1 (MQ=255)
                                                                         aaCTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACTTTGGCCTGGTGAGCAGGGGCAAATAATGAATCCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGa                                                                       >  2:3046632/1‑148 (MQ=255)
                                                                         aaCTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACTTTGGCCTGGTGAGCAGGGACAAATAATGAATCCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGa                                                                       >  1:2051230/1‑148 (MQ=255)
                                                                          aCTTTGACTTTATTACCAGCCAATAATCGAAATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACTTTGGCCTGGTGAGCAGGGGCAAATAATGAATCCGGCAGAGTTTATTCCGCTGGCAGACAAGGAGGGGATGATAGaa                                                                      >  2:296125/1‑148 (MQ=255)
                                                                           cTTTTTCATTTTTACCAGACAATAATCTATAACAAAACACAAAAATGTATCTGCGCTGAAGCGTTGTTACTTTGGCCTGGTGAGAAGGGGCAAATAATGAATCCGGCAGATTTTATTCCGCTGGAAGATAAGGAGGGGATGATAGAAc                                                                     <  1:1320616/148‑1 (MQ=255)
                                                                           cTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACTTTGGCCTGGTGAGCAGGGGCAAATAATGAATCCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAAc                                                                     <  2:497406/148‑1 (MQ=255)
                                                                             ttGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACTTTGGCCTGGTGAGCAGGGGCAAATAATGAATCCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAg                                                                   <  2:156236/148‑1 (MQ=255)
                                                                              tGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACTTTGGCCTGGTGAGCAGGGGCAAATAATGAATCCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGa                                                                  <  1:357438/148‑1 (MQ=255)
                                                                                tCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACTTTGGCCTGGTGAGCAGGGGCAAATAATGAATCCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATa                                                                <  2:1077565/148‑1 (MQ=255)
                                                                                 cTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACTTTGGCCTGGTGAGCAGGGGCAAATAATGAATCCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATaa                                                               <  1:1444356/148‑1 (MQ=255)
                                                                                    tattaCCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACTTTGGCCTGGTGAGCAGGGGCAAATAATGAATCCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTg                                                            <  2:2862772/148‑1 (MQ=255)
                                                                                      ttaCCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACTTTGGCCTGGTGAGCAGGGGCAAATAATGAATCCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGAt                                                          <  1:2392387/148‑1 (MQ=255)
                                                                                        aCCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACTTTGGCCTGGTGAGCAGGGGCAAATAATGAATCCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTa                                                        <  2:2185000/148‑1 (MQ=255)
                                                                                        aCCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACTTTGGCCTGGTGAGCAGGGGCAAATAATGAATCCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTa                                                        <  2:2185090/148‑1 (MQ=255)
                                                                                        aCCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACTTTGGCCTGGTGAGCAGGGGCAAATAATGAATCCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTa                                                        >  2:1408317/1‑148 (MQ=255)
                                                                                         ccAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACTTTGGCCTGGTGAGCAGGGGCAAATAATGAATCCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTAt                                                       <  2:2755134/148‑1 (MQ=255)
                                                                                             ccAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACTTTGGCCTGGTGAGCAGGGGCAAATAATGAATCCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTa                                                   >  1:1393086/1‑148 (MQ=255)
                                                                                             ccAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACTTTGGCCTGGTGAGCAGGGGCAAATAATGAATCCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTa                                                   >  2:1251114/1‑148 (MQ=255)
                                                                                              cAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACTTTGGCCTGGTGAGCAGGGGCAAATAATGAATCCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTAt                                                  <  2:3229040/148‑1 (MQ=255)
                                                                                               aataatCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACTTTGGCCTGGTGAGCAGGGGCAAATAATGAATCCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTAtt                                                 <  1:3103594/148‑1 (MQ=255)
                                                                                                ataatCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACTTTGGCCTGGTGAGCAGGGGCAAATAATGAATCCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTg                                                <  2:237436/148‑1 (MQ=255)
                                                                                                 taatCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACTTTGGCCTGGTGAGCAGGGGCAAATAATGAATCCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGa                                               <  2:1393086/148‑1 (MQ=255)
                                                                                                   atCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACTTTGGCCTGGTGAGCAGGGGCAAATAATGAATCCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATa                                             >  2:1791361/1‑148 (MQ=255)
                                                                                                    tCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACTTTGGCCTGGTGAGCAGGGGCAAATAATGAATCCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATaa                                            >  1:2725868/1‑148 (MQ=255)
                                                                                                     cGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACTTTGGCCTGGTGAGCAGGGGCAAATAATGAATCCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAAt                                           <  1:509878/148‑1 (MQ=255)
                                                                                                       atatCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACTTTGGCCTGGTGAGCAGGGGCAAATAATGAATCCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGt                                         <  2:2340129/148‑1 (MQ=255)
                                                                                                        tatCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACTTTGGCCTGGTGAGCAGGGGCAAATAATGAATCCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTc                                        <  1:2746496/148‑1 (MQ=255)
                                                                                                         atCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACTTTGGCCTGGTGAGCAGGGGCAAATAATGAATCCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCt                                       <  1:943355/148‑1 (MQ=255)
                                                                                                          tCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACTTTGGCCTGGTGAGCAGGGGCAAATAATGAATCCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCtt                                      >  2:2285654/1‑148 (MQ=255)
                                                                                                          tCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACTTTGGCCTGGTGAGCAGGGGCAAATAATGAATCCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCtt                                      >  1:2367742/1‑148 (MQ=255)
                                                                                                          tCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACTTTGGCCTGGTGAGCAGGGGCAAATAATGAATCCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCtt                                      <  1:672189/148‑1 (MQ=255)
                                                                                                           cAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACGTTGGCCTGGTGAGCAGGGGCAAATAATGAATCCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTc                                     <  2:1493376/148‑1 (MQ=255)
                                                                                                               aCAGAAAAATGTATCGGCGCTGAAGCGTTGTTACTTTGGCCTGGTGAGCAGGGGCAAATAATGAATCCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCcgcg                                 >  2:1464803/1‑148 (MQ=255)
                                                                                                                cAGAAAAATGTATCGGCGCTGAAGCGTTGTTACTTTGGCCTGGTGAGCAGGGGCAAATAATGAATCCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGa                                <  2:2305534/148‑1 (MQ=255)
                                                                                                                cAGAAAAATGTATCGGCGCTGAAGCGTTGTTACTTTGGCCTGGTGAGCAGGGGCAAATAATGAATCCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGa                                >  1:1213593/1‑148 (MQ=255)
                                                                                                                cAGAAAAATGTATCGGCGCTGAAGCGTTGTTACTTTGGCCTGGTGAGCAGGGGCAAATAATGAATCCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGa                                >  1:969379/1‑148 (MQ=255)
                                                                                                                   aaaaaTGTATCGGCGCTGAAGCGTTGTTACTTTGGCCTGGTGAGCAGGGGCAAATAATGAATCCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGATCt                             >  1:1056354/1‑148 (MQ=255)
                                                                                                                    aaaaTGTATCGGCGCTGAAGCGTTGTTACTTTGGCCTGGTGAGCAGGGGCAAATAATGAATCCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGATCTg                            <  1:1234981/148‑1 (MQ=255)
                                                                                                                    aaaaTGTATCGGCGCTGAAGCGTTGTTACTTTGGCCTGGTGAGCAGGGGCAAATAATGAATCCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGATCTg                            <  1:1234971/148‑1 (MQ=255)
                                                                                                                         gTATCGGCGCTGAAGCGTTGTTACTTTGGCCTGGTGAGCAGGGGCAAATAATGAATCCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTACGCGATCTGGGCGa                       >  1:3346322/1‑148 (MQ=255)
                                                                                                                          tATCGGCGCTGAAGCGTTGTTACTTTGGCCTGGTGAGCAGGGGCAAATAATGAATCCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGATCTGGGCGAt                      >  2:3112872/1‑148 (MQ=255)
                                                                                                                               gcgcTGAAGCGTTGTTACTTTGGCCTGGTGAGCAGGGGCAAATAATGAATCCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGATCTGGGCGATTACCt                 <  1:2182461/148‑1 (MQ=255)
                                                                                                                                 gcTGAAGCGTTGTTACTTTGGCCTGGTGAGCAGGGGCAAATAATGAATCCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGATCTGGGCGATTACCTgg               <  2:2023952/148‑1 (MQ=255)
                                                                                                                                   tGAAGCGTTGTTACTTTGGCCTGGTGAGCAGGGGCAAATAATGAATCCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGATCTGGGCGATTACCTGGCa             <  2:1091168/148‑1 (MQ=255)
                                                                                                                                    gAAGCGTTGTTACTTTGGCCTGGTGAGCAGGGGCAAATAATGAATCCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGATCTGGGCGATTACCTGGCaa            <  2:92451/148‑1 (MQ=255)
                                                                                                                                        cGTTGTTACTTTGGCCTGGTGAGCAGGGGCAAATAATGAATCCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGATCTGGGCGATTACCTGGCAAcaca        <  2:1821362/148‑1 (MQ=255)
                                                                                                                                        cGTTGTTACTTTGGCCTGGTGAGCAGGGGCAAATAATGAATCCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGATCTGGGCGATTACCTGGCAAcaca        >  2:1449829/1‑148 (MQ=255)
                                                                                                                                              tACGTTGGCCTGGTGAGCAGGGGCAAATAATGAATCCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGATCTGGGCGATTACCTGGCAACACATGCAGa  <  1:347350/148‑1 (MQ=255)
                                                                                                                                              tACGTTGGCCTGGTGAGCAGGGGCAAATAATGAATCCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGATCTGGGCGATTACCTGGCAACACATGCAGa  <  1:2202526/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
TTACTCTGGCTACGTATTCGACAAAACTATTTATCTCCCAAACGTAAATTGCAACGCGCCCTCGAAAAACATCAACTTTGTCTTTATTACCAGCCAATAATCGATATCAAAACAGAAAAATGTATCGGCGCTGAAGCGTTGTTACGTTGGCCTGGTGAGCAGGGGCAAATAATGAATCCGGCAGAGTTTATTCCGCTGGCAGAAAAGGAGGGGATGATAGAACAGATAACTGATTATGTTATTGATAATGTCTTCCGCGATCTGGGCGATTACCTGGCAACACATGCAGA  >  CP009273/4266173‑4266462

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: