Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,334,464 C→A 88.3% S277R (AGC→AGA yrbG → putative calcium/sodium:proton antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,334,4640CA88.3% 451.6 / 61.0 188S277R (AGC→AGAyrbGputative calcium/sodium:proton antiporter
Reads supporting (aligned to +/- strand):  ref base C (12/10);  new base A (91/75);  total (103/85)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.16e-01

Only 100 of 188 total aligned reads displayed.

TAGCGGGGGTTCGCAAAGGTGAAAACGACATTGCTGTCGGAAATATCATTGGCGCAAACATTTTTAATATTGTCATCGTGTTGGGTTTACCCGCGCTGATAACGCCAGGAGAGATTGATCCACTGGCGTACAGTCGTGACTACAGCGTGATGTTGCTGGTGAGCATTATTTTTGCGTTGCTGTGCTGGCGGCGCTCCCCGCAACCGGGCCGTGGTGTAGGGGTATTATTAACTGGCGGATTTATCGTATGGCTGGCGATGTTGTACTGGTTATCGCCAATACTCGTTGAATA  >  CP009273/3334319‑3334610
                                                                                                                                                 |                                                                                                                                                  
tAGCGGGGGTTCGCAAAGGTGAAAACGACATTGCTGTCGGAAATATCATTGGCGCAAACATTTTTAATATTGTCATCGTGTTGGGTTTACCCGCGCTGATAACGCCAGGAGAGATTGATCCACTGGCGTACAGTCGTGACTACAGAGt                                                                                                                                                  >  1:2411410/1‑148 (MQ=255)
tAGCGGGGGTTCGCAAAGGTGAAAACGACATTGCTGTCGGAAATATCATTGGCGCAAACATTTTTAATATTGTCATCGTGTTGGGTTTACCCGCGCTGATAACGCCAGGAGAGATTGATCCACTGGCGTACAGTCGTGACTACAGAGt                                                                                                                                                  >  2:1322916/1‑148 (MQ=255)
  gCGGGGGTTCGCAAAGGTGAAAACGACATTGCTGTCGGAAATATCATTGGCGCAAACATTTTTAATATTGTCATCGTGTTGGGTTTACCCGCGCTGATAACGCCAGGAGAGATTGATCCACTGGCGTACAGTCGTGACTACAGAGTGa                                                                                                                                                >  2:2045381/1‑148 (MQ=255)
    gggggTTCGCAAAGGTGAAAACGACATTGCTGTCGGAAATATCATTGGCGCAAACATTTTTAATATTGTCATCGTGTTGGGTTTACCCGCGCTGATAACGCCAGGAGAGATTGATCCACTGGCGTACAGTCGTGACTACAGAGTGATg                                                                                                                                              <  1:614955/148‑1 (MQ=255)
     ggggTTCGCAAAGGTGAAAACGACATTGCTGTCGGAAATATCATTGGCGCAAACATTTTTAATATTGTCATCGTGTTGGGTTTACCCGCGCTGATAACGCCAGGAGAGATTGATCCACTGGCGTACAGTCGTGACTACAGAGTGATGt                                                                                                                                             >  1:1061562/1‑148 (MQ=255)
     ggggTTCGCAAAGGTGAAAACGACATTGCTGTCGGAAATATCATTGGCGCAAACATTTTTAATATTGTCATCGTGTTGGGTTTACCCGCGCTGATAACGCCAGGAGAGATTGATCCACTGGCGTACAGTCGTGACTACAGAGTGATGt                                                                                                                                             <  2:861023/148‑1 (MQ=255)
            gCAAAGGTGAAAACGACATTGCTGTCGGAAATATCATTGGCGCAAACATTTTTAATATTGTCATCGTGTTGGGTTTACCCGCGCTGATAACGCCAGGAGAGATTGATCCACTGGCGTAGAGTCGTGACTACAGAGTGATGTTGCTGGt                                                                                                                                      <  1:2965802/148‑1 (MQ=255)
            gCAAAGGTGAAAACGACATTGCTGTCGGAAATATCATTGGCGCAAACATTTTTAATATTGTCATCGTGTTGGGTTTACCCGCGCTGATAACGCCAGGAGAGATTGATCCACTGGCGTACAGTCGTGACTACAGAGTGATGTTGCTGGt                                                                                                                                      <  1:1083740/148‑1 (MQ=255)
               aaGGTGAAAACGACATTGCTGTCGGAAATATCATTGGCGCAAACATTTTTAATATTGTCATCGTGTTGGGTTTACCCGCGCTGATAACGCCAGGAGAGATTGATCCACTGGCGTACAGTCGTGACTACAGAGTGATGTTGCTGGTGAg                                                                                                                                   <  2:3104382/148‑1 (MQ=255)
               aaGGTGAAAACGACATTGCTGTCGGAAATATCATTGGCGCAAACATTTTTAATATTGTCATCGTGTTGGGTTTACCCGCGCTGATAACGCCAGGAGAGATTGATCCACTGGCGTACAGTCGTGACTACAGAGTGATGTTGCTGGTGAg                                                                                                                                   >  2:1382467/1‑148 (MQ=255)
                 ggTGAAAACGACATTGCTGTCGGAAATATCATTGGCGCAAACATTTTTAATATTGTCATCGTGTTGGGTTTACCCGCGCTGATAACGCCAGGAGAGATTGATCCACTGGCGTACAGTCGTGACTACAGAGTGATGTTGCTGGTGAGCa                                                                                                                                 >  1:460883/1‑148 (MQ=255)
                 ggTGAAAACGACATTGCTGTCGGAAATATCATTGGCGCAAACATTTTTAATATTGTCATCGTGTTGGGTTTACCCGCGCTGATAACGCCAGGAGAGATTGATCCACTGGCGTACAGTCGTGACTACAGAGTGATGTTGCTGGTGAGCa                                                                                                                                 <  1:2520334/148‑1 (MQ=255)
                 ggTGAAAACGACATTGCTGTCGGAAATATCATTGGCGCAAACATTTTTAATATTGTCATCGTGTTGGGTTTACCCGCGCTGATAACGCCAGGAGAGATTGATCCACTGGCGTACAGTCGTGACTACAGAGTGATGTTGCTGGTGAGCa                                                                                                                                 <  1:2517005/148‑1 (MQ=255)
                 ggTGAAAACGACATTGCTGTCGGAAATATCATTAGCGCAAACATTTTTAATATTGTCATCGTGTTGGGTTTACCCGCGCTGATAACGCCAGGAGAGATTGATCCACTGGCGTACAGTCGTGACTACAGAGTGATGTTGCTGGTGAGCa                                                                                                                                 >  2:516279/1‑148 (MQ=255)
                  gTGAAAACGACATTGCTGTCGGAAATATCATTGGCGCAAACATTTTTAATATTGTCATCGTGTTGGGTTTACCCGCGCTGATAACGCCAGGAGAGATTGATCCACTGGCGTACAGTCGTGACTACAGAGTGATGTTGCTGGTGAGCat                                                                                                                                >  2:1731487/1‑148 (MQ=255)
                   tGAAAACGACATTGCTGTCGGAAATATCATTGGCGCAAACATTTTTAATATTGTCATCGTGTTGGGTTTACCCGCGCTGATAACGCCAGGAGAGATTGATCCACTGGCGTACAGTCGTGACTACAGAGTGATGTTGCTGGTGAGCatt                                                                                                                               <  2:1172317/148‑1 (MQ=255)
                   tGAAAACGACATTGCTGTCGGAAATATCATTGGCGCAAACATTTTTAATATTGTCATCGTGTTGGGTTTACCCGCGCTGATAACGCCAGGAGAGATTGATCCACTGGCGTACAGTCGTGACTACAGAGTGATGTTGCTGGTGAGCatt                                                                                                                               <  2:1171995/148‑1 (MQ=255)
                    gAAAACGACATTGCTGTCGGAAATATCATTGGCGCAAACATTTTTAATATTGTCATCGTGTTGGGTTTACCCGCGCTGATAACGCCAGGAGAGATTGATCCACTGGCGTACAGTCGTGACTACAGAGTGATGTTGCTGGTGAGCatta                                                                                                                              <  1:2647185/148‑1 (MQ=255)
                     aacaCGACATTGCTGTCGGAAATATCATTGGCGCAAACATTTTTAATATTGTCATCGTGTTGGGTTTACCCGCGCTGATAACGCCAGGAGAGATTGATCCACTGGCGTACAGTCGTGACTACAGAGTGATGTTGCTGGTGAGCattat                                                                                                                             <  2:654866/145‑1 (MQ=255)
                     aaaaCGACATTGCTGTCGGAAATATCATTGGCGCAAACATTTTTAATATTGTCATCGTGTTGGGTTTACCCGCGCTGATAACGCCAGGAGAGATTGATCCACTGGCGTACAGTCGTGACTACAGAGTGATGTTGCTGGTGAGCattat                                                                                                                             <  2:1496916/148‑1 (MQ=255)
                     aaaaCGACATTGCTGTCGGAAATATCATTGGCGCAAACATTTTTAATATTGTCATCGTGTTGGGTTTACCCGCGCTGATAACGCCAGGAGAGATTGATCCACTGGCGTACAGTCGTGACTACAGAGTGATGTTGCTGGTGAGCattat                                                                                                                             <  2:939914/148‑1 (MQ=255)
                     aaaaCGACATTGCTGTCGGAAATATCATTGGCGCAAACATTTTTAATATTGTCATCGTGTTGGGTTTACCCGCGCTGATAACGCCAGGAGAGATTGATCCACTGGCGTACAGTCGTGACTACAGAGTGATGTTGCTGGTGAGCattat                                                                                                                             <  2:2592848/148‑1 (MQ=255)
                          gACATTGCTGTCGGAAATATCATTGGCGCAAACATTTTTAATATTGTCATCGTGTTGGGTTTACCCGCGCTGATAACGCCAGGAGAGATTGATCCACTGGCGTACAGTCGTGACTACAGCGTGATGTTGCTGGTGAGCATTATTTTTg                                                                                                                        <  1:3051625/148‑1 (MQ=255)
                           aCATTGCTGTCGGAAATATCATTGGCGCAAACATTTTTAATATTGTCATCGTGTTGGGTTTACCCGCGCTGATAACGCCAGGAGAGATTGATCCACTGGCGTACAGTCGTGACTACAGAGTGATGTTGCTGGTGAGCATTATTTTTGc                                                                                                                       >  1:1042776/1‑148 (MQ=255)
                            cATTGCTGTCGGAAATATCATTGGCGCAAACATTTTTAATATTGTCATCGTGTTGGGTTTACCCGCGCTGATAACGCCAGGAGAGATTGATCCACTGGCGTACAGTCGTGACTACAGAGTGATGTTGCTGGTGAGCATTATTTTTGCg                                                                                                                      <  1:2414328/148‑1 (MQ=255)
                             aTTGCTGTCGGAAATATCATTGGCGCAAACATTTTTAATATTGTCATCGTGTTGGGTTTACCCGCGCTGATAACGCCAGGAGAGATTGATCCACTGGCGTACAGTCGTGACTACAGAGTGATGTTGCTGGTGAGCATTATTTTTGCGt                                                                                                                     >  1:1268633/1‑148 (MQ=255)
                                gCTGTCGGAAATATCATTGGCGCAAACATTTTTAATATTGTCATCGTGTTGGGTTTACCCGCGCTGATAACGCCAGGAGAGATTGATCCACTGGCGTACAGTCGTGACTACAGCGTGATGTTGCTGGTGAGCATTATTTTTGCGTtgc                                                                                                                  <  2:570431/148‑1 (MQ=255)
                                  tGTCGGAAATATCATTGGCGCAAACATTTTTAATATTGTCATCGTGTTGGGTTTACCCGCGCTGATAACGCCAGGAGAGATTGATCCACTGGCGTACAGTCGTGACTACAGCGTGATGTTGCTGGTGAGCATTATTTTTGCGTtgctg                                                                                                                >  2:411911/1‑148 (MQ=255)
                                     cGGAAATATCATTGGCGCAAACATTTTTAATATTGTCATCGTGTTGGGTTTACCCGCGCTGATAACGCCAGGAGAGATTGATCCACTGGCGTACAGTCGTGACTACAGAGTGATGTTGCTGGTGAGCATTATTTTTGCGTtgctgtgc                                                                                                             <  1:516110/148‑1 (MQ=255)
                                     cGGAAATATCATTGGCGCAAACATTTTTAATATTGTCATCGTGTTGGGTTTACCCGCGCTGATAACGCCAGGAGAGATTGATCCACTGGCGTACAGTCGTGACTACAGAGTGATGTTGCTGGTGAGCATTATTTTTGCGTtgctgtgc                                                                                                             >  1:913808/1‑148 (MQ=255)
                                     cGGAAATATCATTGGCGCAAACATTTTTAATATTGTCATCGTGTTGGGTTTACCCGCGCTGATAACGCCAGGAGAGATTGATCCACTGGCGTACAGTCGTGACTACAGAGTGATGTTGCTGGTGAGCATTATTTTTGCGTtgctgtgc                                                                                                             <  1:2864054/148‑1 (MQ=255)
                                         aaTATCATTGGCGCAAACATTTTTAATATTGTCATCGTGTTGGGTTTACCCGCGCTGATAACGCCAGGAGAGATTGATCCACTGGCGTACAGTCGTGACTACAGAGTGATGTTGCTGGTGAGCATTATTTTTGCGTTGCTGTGCTggc                                                                                                         >  1:826898/1‑148 (MQ=255)
                                           tatCATTGGCGCAAACATTTTTAATATTGTCATCGTGTTGGGTTTACCCGCGCTGATAACGCCAGGAGAGATTGATCCACTGGCGTACAGTCGTGACTACAGATTGATGTTTCTGGTGAGCATTATTTTTGCGTTGCTGTGCTggcgg                                                                                                       >  1:1410812/1‑148 (MQ=255)
                                            atCATTGGCGCAAACATTTTTAATATTGTCATCGTGTTGGGTTTACCCGCGCTGATAACGCCAGGAGAGATTGATCCACTGGCGTACAGTCGTGACTACAGAGTGATGTTGCTGGTGAGCATTATTTTTGCGTTGCTGTGCTggcggc                                                                                                      >  2:3183841/1‑148 (MQ=255)
                                            atCATTGGCGCAAACATTTTTAATATTGTCATCGTGTTGGGTTTACCCGCGCTGATAACGCCAGGAGAGATTGATCCACTGGCGTACAGTCGTGACTACAGAGTGATGTTGCTGGTGAGCATTATTTTTGCGTTGCTGTGCTggcggc                                                                                                      >  1:2086535/1‑148 (MQ=255)
                                              cATTGGCGCAAACATTTTTAATATTGTCATCGTGTTGGGTTTACCCGCGCTGATAACGCCAGGAGAGATTGATCCACTGGCGTACAGTCGTGACTACAGAGTGATGTTGCTGGTGAGCATTATTTTTGCGTTGCTGTGCTGGCGgcgc                                                                                                    >  2:2953325/1‑148 (MQ=255)
                                                ttGGCGCAAACATTTTTAATATTGTCATCGTGTTGGGTTTACCCGCGCTGATAACGCCAGGAGAGATTGATCCACTGGCGTACAGTCGTGACTACAGAGTGATGTTGCTGGTGAGCATTATTTTTGCGTTGCTGTGCTGGCGGCGCTc                                                                                                  >  2:2988693/1‑148 (MQ=255)
                                                 tGGCGCAAACATTTTTAATATTGTCATCGTGTTGGGTTTACCCGCGCTGATAACGCCAGGAGAGATTGATCCACTGGCGTACAGTCGTGACTACAGCGTGATGTTGCTGGTGAGCATTATTTTTGCGTTGCTGTGCTGGCGGCGCTcc                                                                                                 <  1:1480901/148‑1 (MQ=255)
                                                 tGGCGCAAACATTTTTAATATTGTCATCGTGTTGGGTTTACCCGCGCTGATAACGCCAGGAGAGATTGATCCACTGGCGTACAGTCGTGACTACAGAGTGATGTTGCTGGTGAGCATTATTTTTGCGTTGCTGTGCTGGCGGCGCTcc                                                                                                 <  1:676055/148‑1 (MQ=255)
                                                      cAAACATTTTTAATATTGTCATCGTGTTGGGTTTACCCGCGCTGATAACGCCCGGAGAGATTGATCCACTGGCGTACAGTCGTGACTACAGAGTGATGTTGCTGGTGAGCATTATTTTTGCGTTGCTGTGCTGGCGGCGCTCCCCGCa                                                                                            >  2:1621099/1‑148 (MQ=255)
                                                      cAAACATTTTTAATATTGTCATCGTGTTGGGTTTACCCGCGCTGATAACGCCAGGAGAGATTGATCCACTGGCGTACAGTCGTGACTACAGAGTGATGTTGCTGGTGAGCATTATTTTTGCGTTGCTGTGCTGGCGGCGCTCCCCGCa                                                                                            >  2:2215999/1‑148 (MQ=255)
                                                          cATTTTTAATATTGTCATCGTGTTGGGTTTACCCGCGCTGATAACGCCAGGAGAGATTGATCCACTGGCGTACAGTCGTGACTACAGCGTGATGTTGCTGGTGAGCATTATTTTTGCGTTGCTGTGCTGGCGGCGCTCCCCGCAACCg                                                                                        >  1:2436052/1‑148 (MQ=255)
                                                          cATTTTTAATATTGTCATCGTGTTGGGTTTACCCGCGCTGATAACGCCAGGAGAGATTGATCCACTGGCGTACAGTCGTGACTACAGAGTGATGTTGCTGGTGAGCATTATTTTTGCGTTGCTGTGCTGGCGGCGCTCCCCGCAACCg                                                                                        >  2:2298204/1‑148 (MQ=255)
                                                           aTTTTTAATATTGTCATCGTGTTGGGTTTACCCGCGCTGATAACGCCAGGAGAGATTGATCCACTGGCGTACAGTCGTGACTACAGAGTGATGTTGCTGGTGAGCATTATTTTTGCGTTGCTGTGCTGGCGGCGCTCCCCGCAACCgg                                                                                       >  2:1845841/1‑148 (MQ=255)
                                                            tttttAATATTGTCATCGTGTTGGGTTTACCCGCGCTGATAACGCCAGGAGAGATTGATCCACTGGCGTACAGTCGTGACTACAGAGTGATGTTGCTGGTGAGCATTATTTTTTCGTTGCTGTGCTGGCGGCGCTCCCCGCAACCggg                                                                                      >  1:1936574/1‑148 (MQ=255)
                                                            tttttAATATTGTCATCGTGTTGGGTTTACCCGCGCTGATAACGCCAGGAGAGATTGATCCACTGGCGTACAGTCGTGACTACAGAGTGATGTTGCTGGTGAGCATTATTTTTGCGTTGCTGTGCTGGCGGCGCTCCCCGCAACCggg                                                                                      >  1:1936587/1‑148 (MQ=255)
                                                             ttttAATATTGTCATCGTGTTGGGTTTACCCGCGCTGATAACGCCAGGAGAGATTGATCCACTGGCGTACAGTCGTGACTACAGAGTGATGTTGCTGGTGAGCATTATTTTTGCGTTGCTGTGCTGGCGGCGCTCCCCGCAACCGGGc                                                                                     >  1:231570/1‑148 (MQ=255)
                                                                   tatTGTCATCGTGTTGGGTTTACCCGCGCTGATAACGCCAGGAGAGATTGATCCACTGGCGTACAGTCGTGACTACAGAGTGATGTTGCTGGTGAGCATTATTTTTGCGTTGCTGTGCTGGCGGCGCTCCCCGCAACCGGGCCgtggt                                                                               >  1:1332623/1‑148 (MQ=255)
                                                                   tatTGTCATCGTGTTGGGTTTACCCGCGCTGATAACGCCAGGAGAGATTGATCCACTGGCGTACAGTCGTGACTACAGAGTGATGTTGCTGGTGAGCATTATTTTTGCGTTGCTGTGCTGGCGGCGCTCCCCGCAACCGGGCCgtggt                                                                               >  1:11803/1‑148 (MQ=255)
                                                                   tatTGTCATCGTGTTGGGTTTACCCGCGCTGATAACGCCAGGAGAGATTGATCCACTGGCGTACAGTCGTGACTACAGAGTGATGTTGCTGGTGAGCATTATTTTTGCGTTGCTGTGCTGGCGGCGCTCCCCGCAACCGGGCCgtggt                                                                               >  1:3177045/1‑148 (MQ=255)
                                                                      tGTCATCGTGTTGGGTTTACCCGCGCTGATAACGCCAGGAGAGATTGATCCACTGGCGTACAGTCGTGACTACAGAGTGATGTTGCTGGTGAGCATTATTTTTGCGTTGCTGTGCTGGCGGCGCTCCCCGCAACCGGGCCGTGGTGTa                                                                            >  2:1712184/1‑148 (MQ=255)
                                                                      tGTCATCGTGTTGGGTTTACCCGCGCTGATAACGCCAGGAGAGATTGATCCACTGGCGTACAGTCGTGACTACAGAGTGATGTTGCTGGTGAGCATTATTTTTGCGTTGCTGTGCTGGCGGCGCTCCCCGCAACCGGGCCGTGGTGTa                                                                            >  2:1712114/1‑148 (MQ=255)
                                                                          aTCGTGTTGGGTTTACCCGCGCTGATAACGCCAGGAGAGATTGATCCACTGGCGTACAGTCGTGACTACAGAGTGATGTTGCTGGTGAGCATTATTTTTGCGTTGCTGTGCTGGCGGCGCTCCCCGCAACCGGGCCGTGGTGTAgggg                                                                        <  2:2378248/148‑1 (MQ=255)
                                                                            cGTGTTGGGTTTACCCGCGCTGATAACGCCAGGAGAGATTGATCCACTGGCGTACAGTCGTGACTACAGAGTGATGTTGCTGGTGAGCATTATTTTTGCGTTGCTGTGCTGGCGGCGCTCCCCGCAACCGGGCCGTGGTGTAGGGGta                                                                      >  2:507105/1‑148 (MQ=255)
                                                                              tgtTGGGTTTACCCGCGCTGATAACGCCAGGAGAGATTGATCCACTGGCGTACAGTCGTGACTACAGAGTGATGTTGCTGGTGAGCATTATTTTTGCGTTGCTGTGCTGGCGGCGCTCCCCGCAACCGGGCCGTGGTGTAGGGGtatt                                                                    >  1:2786023/1‑148 (MQ=255)
                                                                               gtTGGGTTTACCCGCGCTGATAACGCCAGGAGAGATTGATCCACTGGCGTACAGTCGTGACTACAGCGTGATGTTGCTGGTGAGCATTATTTTTGCGTTGCTGTGCTGGCGGCGCTCCCCGCAACCGGGCCGTGGTGTAGGGGtatta                                                                   >  2:180937/1‑148 (MQ=255)
                                                                               gtTGGGTTTACCCGCGCTGATAACGCCAGGAGAGATTGATCCACTGGCGTACAGTCGTGACTACAGAGTGATGTTGCTGGTGAGCATTATTTTTGCGTTGCTGTGCTGGCGGCGCTCCCCGCAACCGGGCCGTGGTGTAGGGGtatta                                                                   <  2:3278295/148‑1 (MQ=255)
                                                                                ttGGGTTTACCCGCGCTGATAACGCCAGGAGAGATTGATCCACTGGCGTACAGTCGTGACTACAGAGTGATGTTGCTGGTGAGCATTATTTTTGCGTTGCTGTGCTGGCGGCGCTCCCCGCAACCGGGCCGTGGTGTAGGGGtattat                                                                  >  2:461755/1‑148 (MQ=255)
                                                                                  gggTTTACCCGCGCTGATAACGCCAGGAGAGATTGATCCACTGGCGTACAGTCGTGACTACAGAGTGATGTTGCTGGTGAGCATTATTTTTGCGTTGCTGTGCTGGCGGCGCTCCCCGCAACCGGGCCGTGGTGTAGGGGTAttatta                                                                <  2:1098769/148‑1 (MQ=255)
                                                                                   ggTTTACCCGCGCTGATAACGCCAGGAGAGATTGATCCACTGGCGTACAGTCGTGACTACAGAGTGATGTTGCTGGTGAGCATTATTTTTGCGTTGCTGTGCTGGCGGCGCTCCCCGCAACCGGGCCGTGGTGTAGGGGTATTATTaa                                                               <  1:1558879/148‑1 (MQ=255)
                                                                                    gTTTACCCGCGCTGATAACGCCAGGAGAGATTGATCCACTGGCGTACAGTCGTGACTACAGAGTGATGTTGCTGGTGAGCATTATTTTTGCGTTGCTGTGCTGGCGGCGCTCCCCGCAACCGGGCCGTGGTGTAGGGGTATTATTAAc                                                              <  1:365515/148‑1 (MQ=255)
                                                                                     tttACCCGCGCTGATAACGCCAGGAGAGATTGATCCACTGGCGTACAGTCGTGACTACAGCGTGATGTTGCTGGTGAGCATTATTTTTGCGTTGCTGTGCTGGCGGCGCTCCCCGCAACCGGGCCGTGGTGTAGGGGTATTATTAACt                                                             >  2:329250/1‑148 (MQ=255)
                                                                                      ttACCCGCGCTGATAACGCCAGGAGAGATTGATCCACTGGCGTACAGTCGTGACTACAGAGTGATGTTGCTGGTGAGCATTATTTTTGCGTTGCTGTGCTGGCGGCGCTCCCCGCAACCGGGCCGTGGTGTAGGGGTATTATTAACTg                                                            <  2:1268443/148‑1 (MQ=255)
                                                                                          ccGCGCTGATAACGCCAGGAGAGATTGATCCACTGGCGTACAGTCGTGACTACAGAGTGATGTTGCTGGTGAGCATTATTTTTGCGTTGCTGTGCTGGCGGCGCTCCCCGCAACCGGGCCGTGGTGTAGGGGTATTATTAACTGGCgg                                                        <  2:410733/148‑1 (MQ=255)
                                                                                           cgcgCTGATAACGCCAGGAGAGATTGATCCACTGGCGTACAGTCGTGACTACAGAGTGATGTTGCTGGTGAGCATTATTTTTGCGTTGCTGTGCTGGCGGCGCTCCCCGCAACCGGGCCGTGGTGTAGGGGTATTATTAACTGGCGGa                                                       <  1:2002613/148‑1 (MQ=255)
                                                                                            gcgcTGATAACGCCAGGAGAGATTGATCCACTGGCGTACAGTCGTGACTACAGAGTGATGTTGCTGGTGAGCATTATTTTTGCGTTGCTGTGCTGGCGGCGCTCCCCGCAACCGGGCCGTGGTGTAGGGGTATTATTAACTGGCGGAt                                                      <  1:3075849/148‑1 (MQ=255)
                                                                                            gcgcTGATAACGCCAGGAGAGATTGATCCACTGGCGTACAGTCGTGACTACAGAGTGATGTTGCTGGTGAGCATTATTTTTGCGTTGCTGTGCTGGCGGCGCTCCCCGCAACCGGGCCGTGGTGTAGGGGTATTATTAACTGGCGGAt                                                      <  1:1857051/148‑1 (MQ=255)
                                                                                              gcTGATAACGCCAGGAGAGATTGATCCACTGGCGTACAGTCGTGACTACAGAGTGATGTTGCTGGTGAGCATTATTTTTGCGTTGCTGTGCTGGCGGCGCTCCCCGCAACCGGGCCGTGGTGTAGGGGTATTATTAACTGGCGGAttt                                                    <  2:1806008/148‑1 (MQ=255)
                                                                                               cTGATAACGCCAGGAGAGATTGATCCACTGGCGTACAGTCGTGACTACAGAGTGATGTTGCTGGTGAGCATTATTTTTGCGTTGCTGTGCTGGCGGCGCTCCCCGCAACCGGGCCGTGGTGTAGGGGTATTATTAACTGGCGGATTTa                                                   <  2:3118572/148‑1 (MQ=255)
                                                                                                  aTAACGCCAGGAGAGATTGATCCACTGGCGTACAGTCGTGACTACAGAGTGATGTTGCTGGTGAGCATTATTTTTGCGTTGCTGTGCTGGCGGCGCTCCCCGCAACCGGGCCGTGGTGTAGGGGTATTATTAACTGGCGGATTTATCg                                                >  2:2125875/1‑148 (MQ=255)
                                                                                                   tAACGCCAGGAGAGATTGATCCACTGGCGTACAGTCGTGACTACAGAGTGATGTTGCTGGTGAGCATTATTTTTGCGTTGCTGTGCTGGCGGCGCTCCCCGCAACCGGGCCGTGGTGTAGGGGTATTATTAACTGGCGGATTTATCGt                                               <  2:999489/148‑1 (MQ=255)
                                                                                                   tAACGCCAGGAGAGATTGATCCACTGGCGTACAGTCGTGACTACAGAGTGATGTTGCTGGTGAGCATTATTTTTGCGTTGCTGTGCTGGCGGCGCTCCCCGCAACCGGGCCGTGGTGTAGGGGTATTATTAACTGGCGGATTTATCGt                                               >  2:1856669/1‑148 (MQ=255)
                                                                                                    aaCGCCAGGAGAGATTGATCCACTGGCGTACAGTCGTGACTACAGAGTGATGTTGCTGGTGAGCATTATTTTTGCGTTGCTGTGCTGGCGGCGCTCCCCGCAACCGGGCCGTGGTGTAGGGGTATTATTAACTGGCGGATTTATCGTa                                              >  2:1921285/1‑148 (MQ=255)
                                                                                                     aCGCCAGGAGAGATTGATCCACTGGCGTACAGTCGTGACTACAGAGTGATGTTGCTGGTGAGCATTATTTTTGCGTTGCTGTGCTGGCGGCGCTCCCCGCAACCGGGCCGTGGTGTAGGGGTATTATTAACTGTCGGATTTATCGTAt                                             >  2:2578873/1‑148 (MQ=255)
                                                                                                       gCCAGGAGAGATTGATCCACTGGCGTACAGTCGTGACTACAGAGTGATGTTGCTGGTGAGCATTATTTTTGCGTTGCTGTGCTGGCGGCGCTCCCCGCAACCGGGCCGTGGTGTAGGGGTATTATTAACTGGCGGATTTATCGTAtgg                                           <  2:2496669/148‑1 (MQ=255)
                                                                                                           ggAGAGATTGATCCACTGGCGTACAGTCGTGACTACAGAGTGATGTTGCTGGTGAGCATTATTTTTGCGTTGCTGTGCTGGCGGCGCTCCCCGCAACCGGGCCGTGGTGTAGGGGTATTATTAACTGGCGGATTTATCGTAtggctgg                                       >  2:494894/1‑148 (MQ=255)
                                                                                                           ggAGAGATTGATCCACTGGCGTACAGTCGTGACTACAGAGTGATGTTGCTGGTGAGCATTATTTTTGCGTTGCTGTGCTGGCGGCGCTCCCCGCAACCGGGCCGTGGTGTAGGGGTATTATTAACTGGCGGATTTATCGTAtggctgg                                       >  1:738152/1‑148 (MQ=255)
                                                                                                             agagaTTGATCCACTGGCGTACAGTCGTGACTACAGAGTGATTTTGCTGGTGAGCATTATTTTTGCGTTGCTGTGCTGGCGGCGCTCCCCGCAACCGGGCCGTGGTGTAGGGGTATTATTAACTGGCGGATTTATCGTATGGCTGGCg                                     >  1:721461/1‑148 (MQ=255)
                                                                                                                gaTTGATCCACTGGCGTACAGTCGTGACTACAGAGTGATGTTGCTGGTGAGCATTATTTTTGCGTTGCTGTGCTGGCGGCGCTCCCCGCAACCGGGCCGTGGTGTAGGGGTATTATTAACTGGCGGATTTATCGTATGGCTGGCGAtg                                  <  1:1673540/148‑1 (MQ=255)
                                                                                                                   tGATCCACTGGCGTACAGTCGTGACTACAGAGTGATGTTGCTGGTGAGCATTATTTTTGCGTTGCTGTGCTGGCGGCGCTCCCCGCAACCGGGCCGTGGTGTAGGGGTATTATTAACTGGCGGATTTATCGTATGGCTGGCGAtgttg                               >  2:2620088/1‑148 (MQ=255)
                                                                                                                    gATCCACTGGCGTACAGTCGTGACTACAGAGTGATGTTGCTGGTGAGCATTATTTTTGCGTTGCTGTGCTGGCGGCGCTCCCCGCAACCGGGCCGTGGTGTAGGGGTATTATTAACTGGCGGATTTATCGTATGGCTGGCGAtgttgt                              <  2:1061562/148‑1 (MQ=255)
                                                                                                                     aTCCACTGGCGTACAGTCGTGACTACAGAGTGATGTTGCTGGTGAGCATTATTTTTGCGTTGCTGTGCTGGCGGCGCTCCCCGCAACCGGGCCGTGGTGTAGGGGTATTATTAACTGGCGGATTTATCGTATGGCTGGCGATGTTGTa                             <  1:405342/148‑1 (MQ=255)
                                                                                                                       ccACTGGCGTACAGTCGTGACTACAGAGTGATGTTGCTGGTGAGCATTATTTTTGCGTTGCTGTGCTGGCGGCGCTCCCCGCAACCGGGCCGTGGTGTAGGGGTATTATTAACTGGCGGATTTATCGTATGGCTGGCGATGTTGTACt                           <  2:2143068/148‑1 (MQ=255)
                                                                                                                       ccACTGGCGTACAGTCGTGACTACAGAGTGATGTTGCGGGTGAGCATTATTTTTGCGTTGCTGTGCTGGCGGCGCTCCCCGCAACCGGGCCGTGGTGTAGGGGTATTATTAACTGGCGGATTTATCGTATGGCTGGCGATGTTGTACt                           <  2:930338/148‑1 (MQ=255)
                                                                                                                        cACTGGCGTACAGTCGTGACTACAGAGTGATGTTGCTGGTGAGCATTATTTTTGCGTTGCTGTGCTGGCGGCGCTCCCCGCAACCGGGCCGTGGTGTAGGGGTATTATTAACTGGCGGATTTATCGTATGGCTGGCGATGTTGTACTg                          <  2:1415310/148‑1 (MQ=255)
                                                                                                                        cACTGGCGTACAGTCGTGACTACAGAGTGATGGTGCTGGTGAGCATTATTTTTGCGTTGCTGTGCTGGCGGCGCTCCCCGCAACCGGGCCGTGGTGTAGGGGTATTATTAACTGGCGGATTTATCGTATGGCTGGCGATGTTGTACTg                          <  2:3058576/148‑1 (MQ=255)
                                                                                                                         aCTGGCGTACAGTCGTGACTACAGAGTGATGTTGCTGGTGAGCATTATTTTTGCGTTGCTGTGCTGGCGGCGCTCCCCGCAACCGGGCCGTGGTGTAGGGGTATTATTAACTGGCGGATTTATCGTATGGCTGGCGATGTTGTACTgg                         <  1:663516/148‑1 (MQ=255)
                                                                                                                         aCTGGCGTACAGTCGTGACTACAGAGTGATGTTGCTGGTGAGCATTATTTTTGCGTTGCTGTGCTGGCGGCGCTCCCCGCAACCGGGCCGTGGTGTAGGGGTATTATTAACTGGCGGATTTATCGTATGGCTGGCGATGTTGTACTgg                         <  1:522100/148‑1 (MQ=255)
                                                                                                                                tACAGTCGTGACTACAGAGTGATGTTGCTGGTGAGCATTATTTTTGCGTTGCTGTGCTGGCGGCGCTCCCCGCAACCGGGCCGTGGTGTAGGGGTATTATTAACTGGCGGATTTATCGTATGGCTGGCGATGTTGTACTGGTTATCGc                  <  1:173371/148‑1 (MQ=255)
                                                                                                                                 aCAGTCGTGACTACAGAGTGATGTTGCTGGTGAGCATTATTTTTGCGTTGCTGTGCTGGCGGCGCTCCCCGCAACCGGGCCGTGGTGTAGGGGGATTATTAACTGGCGGATTTATCGTATGGCTGGCGATGTTGTACTGGTTATCGcc                 >  1:1225759/1‑148 (MQ=255)
                                                                                                                                  cAGTCGTGACTACAGAGTGATGTTGCTGGTGAGCATTATTTTTGCGTTGCTGTGCTGGCGGCGCTCCCCGCAACCGGGCCGTGGTGTAGGGGTATTATTAACTGGCGGATTTATCGTATGGCTGGCGATGTTGTACTGGTTATCGCCa                >  1:477525/1‑148 (MQ=255)
                                                                                                                                   aGTCGTGACTACAGAGTGATGTTGCTGGTGAGCATTATTTTTGCGTTGCTGTGCTGGCGGCGCTCCCCGCAACCGGGCCGTGGTGTAGGGGTATTATTAACTGGCGGATTTATCGTATGGCTGGCGATGTTGTACTGGTTATCGCCaa               >  2:3065403/1‑148 (MQ=255)
                                                                                                                                     tCGTGACTACAGAGTGATGTTGCTGGTGAGCATTATTTTTGCGTTGCTGTGCTGGCGGCGCTCCCCGCAACCGGGCCGTGGTGTAGGGGTATTATTAACTGGCGGATTTATCGTATGGCTGGCGATGTTGTACTGGTTATCGCCAATa             >  2:2571110/1‑148 (MQ=255)
                                                                                                                                          aCTACAGCGTGATGTTGCTGGTGAGCATTATTTTTGCGTTGCTGTGCTGGCGGCGCTCCCCGCAACCGGGCCGTGGTGTAGGGGTATTATTAACTGGCGGATTTATCGTATGGCTGGCGATGTTGTACTGGTTATCGCCAATACTCGt        <  2:2436026/148‑1 (MQ=255)
                                                                                                                                             aCAGCGTGATGTTGCTGGTGAGCATTATTTTTGCGTTGCTGTGCTGGCGGCGCTCCCCGCAACCGGGCCGTGGTGTAGGGGTATTATTAACTGGCGGATTTATCGTATGGCTGGCGATGTTGTACTGGTTATCGCCAATACTCGTTGa     >  1:1824326/1‑148 (MQ=255)
                                                                                                                                             aCAGAGTGATGTTGCTGGTGAGCATTATTTTTGCGTTGCTGTGCTGGCGGCGCTCCCCGCAACCGGGCCGTGGTGTAGGGGTATTATTAACTGGCGGATTTATCGTATGGCTGGCGATGTTGTACTGGTTATCGCCAATACTCGTTGa     >  1:3179153/1‑148 (MQ=255)
                                                                                                                                             aCAGAGTGATGTTGCTGGTGAGCATTATTTTTGCGTTGCTGTGCTGGCGGCGCTCCCCGCAACCGGGCCGTGGTGTAGGGGTATTATTAACTGGCGGATTTATCGTATGGCTGGCGATGTTGTACTGGTTATCGCCAATACTCGTTGa     >  1:161775/1‑148 (MQ=255)
                                                                                                                                             aCAGAGTGATGTTGCTGGTGAGCATTATTTTTGCGTTGCTGTGCTGGCGGCGCTCCCCGCAACCGGGCCGTGGTGTAGGGGTATTATTAACTGGCGGATTTATCGTATGGCTGGCGATGTTGTACTGGTTATCGCCAATACTCGTTGa     <  1:1389926/148‑1 (MQ=255)
                                                                                                                                             aCAGAGTGATGTTGCTGGTGAGCATTATTTTTGCGTTGCTGTGCTGGCGGCGCTCCCCGCAACCGGGCCGTGGTGTAGGGGTATTATTAACTGGCGGATTTATCGTATGGCTGGCGATGTTGTACTGGTTATCGCCAATACTCGTTGa     >  1:2167571/1‑148 (MQ=255)
                                                                                                                                              cAGAGTGATGTTGCTGGTGAGCATTATTTTTGCGTTGCTGTGCTGGCGGCGCTCCCCGCAACCGGGCCGTGGTGTAGGGGTATTATTAACTGGCGGATTTATCGTATGGCTGGCGATGTTGTACTGGTTATCGCCAATACTCGTTGaa    >  2:2030012/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
TAGCGGGGGTTCGCAAAGGTGAAAACGACATTGCTGTCGGAAATATCATTGGCGCAAACATTTTTAATATTGTCATCGTGTTGGGTTTACCCGCGCTGATAACGCCAGGAGAGATTGATCCACTGGCGTACAGTCGTGACTACAGCGTGATGTTGCTGGTGAGCATTATTTTTGCGTTGCTGTGCTGGCGGCGCTCCCCGCAACCGGGCCGTGGTGTAGGGGTATTATTAACTGGCGGATTTATCGTATGGCTGGCGATGTTGTACTGGTTATCGCCAATACTCGTTGAATA  >  CP009273/3334319‑3334610

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: