Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,019,343 G→T 36.6% L73L (CTG→CTT)  guaD → guanine deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,019,3430GT36.6% 336.0 / 238.1 232L73L (CTG→CTT) guaDguanine deaminase
Reads supporting (aligned to +/- strand):  ref base G (76/71);  new base T (27/58);  total (103/129)
Fisher's exact test for biased strand distribution p-value = 3.93e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.94e-01

Only 100 of 234 total aligned reads displayed.

ACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAACACACCTTCCCTACTGAACGTCGTTATGAGGATTTAGAGTACGCCCGCGAAATGTC  >  CP009273/3019195‑3019492
                                                                                                                                                    |                                                                                                                                                     
 cccGGAATAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTTAt                                                                                                                                                     <  2:2526594/150‑1 (MQ=255)
 cccGGAAGTGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGAt                                                                                                                                                     <  1:1447136/150‑1 (MQ=255)
   cGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAg                                                                                                                                                   >  2:3332147/1‑150 (MQ=255)
    ggAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGt                                                                                                                                                  <  2:1680528/150‑1 (MQ=255)
      aaGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAATGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTAc                                                                                                                                                >  1:930036/1‑150 (MQ=255)
       agagATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTAcc                                                                                                                                               <  2:919764/150‑1 (MQ=255)
             tGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCtt                                                                                                                                         <  2:1443952/150‑1 (MQ=255)
             tGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCtt                                                                                                                                         >  2:2799886/1‑150 (MQ=255)
              gCCTCTGCTCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATTGTACCGGGCttt                                                                                                                                        <  1:1822914/150‑1 (MQ=255)
              gCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCttt                                                                                                                                        >  2:377509/1‑150 (MQ=255)
                ctctGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGt                                                                                                                                      <  1:842404/150‑1 (MQ=255)
                ctctGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGt                                                                                                                                      <  1:1037713/150‑1 (MQ=255)
                ctctGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCAAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGt                                                                                                                                      <  1:2930287/150‑1 (MQ=255)
                 tctGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTc                                                                                                                                     >  2:1851892/1‑150 (MQ=255)
                    gcgctgcgGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGAt                                                                                                                                  <  2:1788783/150‑1 (MQ=255)
                        tgcgGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATacac                                                                                                                              <  2:2374258/150‑1 (MQ=255)
                           ggTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTTATAGTACCGGGCTTTGTCGATACACata                                                                                                                           >  1:891006/1‑150 (MQ=255)
                            gTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACatat                                                                                                                          <  1:36203/150‑1 (MQ=255)
                            gTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACatat                                                                                                                          >  2:3994665/1‑150 (MQ=255)
                             tttATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTTATAGTACCGGGCTTTGTCGATACACATATc                                                                                                                         <  1:339745/150‑1 (MQ=255)
                              ttATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTTATAGTACCGGGCTTTGTCGATACACATATcc                                                                                                                        >  1:4037110/1‑150 (MQ=255)
                              ttATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACATATcc                                                                                                                        >  2:3162431/1‑150 (MQ=255)
                                 ttGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACAAATATCCAtt                                                                                                                     >  2:463837/1‑150 (MQ=255)
                                   gAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTTATAGTACCGGGCTTTGTCGATACACATATCCATTAt                                                                                                                   <  1:1502785/150‑1 (MQ=255)
                                   gAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTAt                                                                                                                   <  1:4162987/150‑1 (MQ=255)
                                    aGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTTATAGTACCGGGCTTTGTCGATACACATATCCATTATc                                                                                                                  >  1:92195/1‑150 (MQ=255)
                                    aGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTATc                                                                                                                  >  1:2170792/1‑150 (MQ=255)
                                     ggATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTTATAGTACCGGGCTTTGTCGATACACATATCCATTATcc                                                                                                                 <  1:3305224/150‑1 (MQ=255)
                                     ggATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTTATAGTACCGGGCTTTGTCGATACACATATCCATTATcc                                                                                                                 <  1:3305281/150‑1 (MQ=255)
                                      gATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTTATAGTACCGGGCTTTGTCGATACACATATCCATTATCCg                                                                                                                <  1:3509762/150‑1 (MQ=255)
                                          gTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTTATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCaaa                                                                                                            <  2:3483754/150‑1 (MQ=255)
                                          gTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTTATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCaaa                                                                                                            <  2:3483655/150‑1 (MQ=255)
                                          gTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCaaa                                                                                                            <  2:4053748/150‑1 (MQ=255)
                                           tttATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAg                                                                                                           >  1:1533900/1‑150 (MQ=255)
                                            ttattaCTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTTATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGt                                                                                                          <  1:698429/150‑1 (MQ=255)
                                            ttattaCTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGt                                                                                                          <  2:1333519/150‑1 (MQ=255)
                                                 aCTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTGAAAt                                                                                                     >  2:945218/1‑150 (MQ=255)
                                                 aCTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTGAAAt                                                                                                     <  2:1798057/150‑1 (MQ=255)
                                                  cTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTTATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTGAAAtg                                                                                                    >  1:2620606/1‑150 (MQ=255)
                                                  cTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTTATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTGAAAtg                                                                                                    <  1:3472582/150‑1 (MQ=255)
                                                  cTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTGAAAtg                                                                                                    >  2:2618950/1‑150 (MQ=255)
                                                     aTTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTGAAAtggtg                                                                                                 <  1:2531656/150‑1 (MQ=255)
                                                          aCAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCATACTTATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTGAAATGGTGGGGGc                                                                                            <  2:2022161/150‑1 (MQ=255)
                                                          aCAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTGAAATGGTGGGGGc                                                                                            <  1:1942047/150‑1 (MQ=255)
                                                              ggAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTTATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTAt                                                                                        <  1:3222502/150‑1 (MQ=255)
                                                              ggAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTAt                                                                                        >  2:3123839/1‑150 (MQ=255)
                                                                aaaaGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATgg                                                                                      <  1:3696361/150‑1 (MQ=255)
                                                                 aaaGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGt                                                                                     <  1:3138708/150‑1 (MQ=255)
                                                                   aGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGa                                                                                   >  2:2108527/1‑150 (MQ=255)
                                                                    gTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAg                                                                                  >  2:611719/1‑150 (MQ=255)
                                                                     tGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGc                                                                                 >  1:2810578/1‑150 (MQ=255)
                                                                      ggAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTTATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCa                                                                                >  1:3203478/1‑150 (MQ=255)
                                                                        aaTGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAAt                                                                              <  2:733557/150‑1 (MQ=255)
                                                                           ggTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGc                                                                           >  2:3010122/1‑150 (MQ=255)
                                                                            gTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACCTATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCt                                                                          >  2:2175737/1‑150 (MQ=255)
                                                                             ttgggCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTg                                                                         <  2:1236972/147‑1 (MQ=255)
                                                                               tGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGa                                                                       >  2:780398/1‑150 (MQ=255)
                                                                                   gAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTgg                                                                   >  1:3078832/1‑150 (MQ=255)
                                                                                    aaTGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGt                                                                  >  1:3334638/1‑150 (MQ=255)
                                                                                          aaaaCGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATa                                                            >  1:1724556/1‑150 (MQ=255)
                                                                                           aaaCGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATaa                                                           <  2:4026345/150‑1 (MQ=255)
                                                                                           aaaCGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATaa                                                           >  2:3221007/1‑150 (MQ=255)
                                                                                            aaCGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTTATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATaaa                                                          <  2:3584100/150‑1 (MQ=255)
                                                                                              cGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTTATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTGAAATGGTGGGGGCATATGGTGAGCAATTGCTGGAGTGGTTGAATAAaca                                                        <  2:2550371/150‑1 (MQ=255)
                                                                                               ggAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAacac                                                       <  1:1722129/150‑1 (MQ=255)
                                                                                               ggAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAacac                                                       >  1:2029338/1‑150 (MQ=255)
                                                                                                gAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAacaca                                                      <  1:950825/150‑1 (MQ=255)
                                                                                                 aaaGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTTATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAacacac                                                     <  2:3021295/150‑1 (MQ=255)
                                                                                                  aaGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAACACAcc                                                    >  2:2318267/1‑150 (MQ=255)
                                                                                                  aaGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTGAAATGGTGGGGGCCGATGGTGAGCAATTGCTGGAGTGGTTGAATAAACACAcc                                                    >  2:1097785/1‑150 (MQ=255)
                                                                                                   aGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTTATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAACACACCt                                                   <  2:1606768/150‑1 (MQ=255)
                                                                                                    gCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTTATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAACACACCtt                                                  <  1:287336/150‑1 (MQ=255)
                                                                                                    gCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAACACACCtt                                                  >  2:1957446/1‑150 (MQ=255)
                                                                                                    gCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAACACACCtt                                                  <  1:4171059/150‑1 (MQ=255)
                                                                                                     cATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAACACACCTTc                                                 <  2:1553209/150‑1 (MQ=255)
                                                                                                      aTCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAACACACCTTcc                                                <  2:1376550/150‑1 (MQ=255)
                                                                                                      aTCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAACACACCTTcc                                                >  1:3631912/1‑150 (MQ=255)
                                                                                                         aaaTTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTTATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAACACACCTTCCCTa                                             >  2:2821741/1‑150 (MQ=255)
                                                                                                           gttCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAATGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAACACACCGTCCCTACt                                           >  1:1256498/2‑150 (MQ=255)
                                                                                                               cTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTTATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAACACACCTTCCCTACTGAAc                                       <  2:3673343/150‑1 (MQ=255)
                                                                                                               cTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAACACACCTTCCCTACTGAAc                                       >  2:2419552/1‑150 (MQ=255)
                                                                                                                 gACACCATTCGCGTGCGCGACTATCGCGGCAAACTTATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAACACACCTTCCCTACTGAAcgt                                     >  1:1764852/1‑150 (MQ=255)
                                                                                                                   cacCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAACACACCTTCCCTACTGAAcgtcg                                   <  2:2546185/150‑1 (MQ=255)
                                                                                                                   cacCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAACACACCTTCCCTACTGAAcgtcg                                   <  2:3309272/150‑1 (MQ=255)
                                                                                                                   cacCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAACACACCTTCCCTACTGAAcgtcg                                   <  2:3309303/150‑1 (MQ=255)
                                                                                                                      cATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCTATACACATATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAACACACCTTCCCTACTGAACGTCGTTa                                <  1:1166747/150‑1 (MQ=255)
                                                                                                                          cgcgTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAACACACCTTCCCTACTGAACGTCGTTATGAg                            >  1:2444021/1‑150 (MQ=255)
                                                                                                                             gTGCGCGACTATCGCGGCAAACTTATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAACACACCTTCCCTACTGAACGTCGTTATGAGGAt                         >  1:204153/1‑150 (MQ=255)
                                                                                                                             gTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAACACACCGTCCCTACTGAACGTCGTTATGAGGAt                         <  2:3928290/150‑1 (MQ=255)
                                                                                                                              tGCGCGACTATCGCGGCAAACTTATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAACACACCTTCCCTACTGAACGTCGTTATGAGGAtt                        <  2:299478/150‑1 (MQ=255)
                                                                                                                               gcgcGACTATCGCGGCAAACTTATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAACACACCTTCCCTACTGAACGTCGTTATGAGGAttt                       <  2:2095379/150‑1 (MQ=255)
                                                                                                                               gcgcGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAACACACCTTCCCTACTGAACGTCGTTATGAGGAttt                       <  1:3795156/150‑1 (MQ=255)
                                                                                                                                   gACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAACACACCTTCCCTACTGAACGTCGTTATGAGGATTTagag                   >  1:2655735/1‑150 (MQ=255)
                                                                                                                                    aCTATCGCGGCAAACTTATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAACACACCTTCCCTACTGAACGTCGTTATGAGGATTTAGAGt                  <  2:1714711/150‑1 (MQ=255)
                                                                                                                                     cTATCGCGGCAAACTTATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAACACACCTACCCTACTGAACGTCGGTATGAGGATTTAGAGTa                 >  1:406560/1‑150 (MQ=255)
                                                                                                                                      tATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAACACACCTTCCCTACTGAACGTCGTTATGAGGATTTAGAGTAc                >  2:3695488/1‑150 (MQ=255)
                                                                                                                                          gcgGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAACACACCTTCCCTACTGAACGTCGTTATGAGGATTTAGAGTACGccc            <  2:3156396/150‑1 (MQ=255)
                                                                                                                                           cgGCAAACTTATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAACACACCTTCCCTACTGAACGTCGTTATGAGGATTTAGAGTACGCCcg           >  1:1304996/1‑150 (MQ=255)
                                                                                                                                           cgGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAACACACCTTCCCTACTGAACGTCGTTATGAGGATTTAGAGTACGCCcg           <  2:1769875/150‑1 (MQ=255)
                                                                                                                                              cAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAACACACCTTCCCTACTGAACGTCGTTATGAGGATTTAGAGTACGCCCGCGa        >  1:2297095/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
ACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCGGGCTTTGTCGATACACATATCCATTATCCGCAAAGTGAAATGGTGGGGGCCTATGGTGAGCAATTGCTGGAGTGGTTGAATAAACACACCTTCCCTACTGAACGTCGTTATGAGGATTTAGAGTACGCCCGCGAAATGTC  >  CP009273/3019195‑3019492

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: ‑