Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,388,439 T→A 52.8% V404E (GTG→GAG)  mutL → methyl‑directed mismatch repair protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,388,4390TA52.8% 23.3 / 260.5 178V404E (GTG→GAG) mutLmethyl‑directed mismatch repair protein
Reads supporting (aligned to +/- strand):  ref base T (45/39);  new base A (49/44);  total (94/84)
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.24e-01

Only 100 of 181 total aligned reads displayed.

CGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGG  >  CP009273/4388293‑4388586
                                                                                                                                                  |                                                                                                                                                   
  tGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTa                                                                                                                                                  >  1:1012108/1‑148 (MQ=255)
  tGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTa                                                                                                                                                  >  1:1011446/1‑148 (MQ=255)
  tGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTa                                                                                                                                                  <  2:3459511/148‑1 (MQ=255)
  tGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTa                                                                                                                                                  <  1:89799/148‑1 (MQ=255)
   ggcggcGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTAt                                                                                                                                                 <  2:3767998/148‑1 (MQ=255)
    gcggcgGGGCACAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATc                                                                                                                                                >  1:259150/1‑148 (MQ=255)
       gcggGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGcc                                                                                                                                             <  2:3380940/148‑1 (MQ=255)
          gggCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGc                                                                                                                                          <  1:1744134/148‑1 (MQ=255)
            gcgcAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCtt                                                                                                                                        <  1:1526855/148‑1 (MQ=255)
            gcgcAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCtt                                                                                                                                        <  1:2770308/148‑1 (MQ=255)
            gcgcAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAAAAGCAAGGTGAAGTGTATCGCCAGCtt                                                                                                                                        <  1:1526010/148‑1 (MQ=255)
             cgcAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCttt                                                                                                                                       <  2:1829704/148‑1 (MQ=255)
               cAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTg                                                                                                                                     >  1:3151975/1‑148 (MQ=255)
               cAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTg                                                                                                                                     >  1:2179112/1‑148 (MQ=255)
               cAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTg                                                                                                                                     <  2:751705/148‑1 (MQ=255)
                 aTCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCa                                                                                                                                   >  1:494020/1‑148 (MQ=255)
                  tCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGTCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCaa                                                                                                                                  <  2:1871436/148‑1 (MQ=255)
                  tCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCaa                                                                                                                                  >  2:738132/1‑148 (MQ=255)
                  tCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCaa                                                                                                                                  <  1:1384568/148‑1 (MQ=255)
                               ccGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGa                                                                                                                     <  2:1104169/148‑1 (MQ=255)
                                     gCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCaaa                                                                                                               >  2:324274/1‑148 (MQ=255)
                                     gCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCaaa                                                                                                               >  2:1659540/1‑148 (MQ=255)
                                      cTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCaaaa                                                                                                              >  1:226339/1‑148 (MQ=255)
                                      cTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCaaaa                                                                                                              <  1:1245849/148‑1 (MQ=255)
                                        cGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAAt                                                                                                            <  1:2504837/148‑1 (MQ=255)
                                          tGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTa                                                                                                          >  2:2944757/1‑148 (MQ=255)
                                            aGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTaaa                                                                                                        <  1:1421888/148‑1 (MQ=255)
                                            aGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTaaa                                                                                                        <  2:1059101/148‑1 (MQ=255)
                                                 gTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAgcgc                                                                                                   >  1:2160635/1‑148 (MQ=255)
                                                  tAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGcc                                                                                                  >  2:2669665/1‑148 (MQ=255)
                                                     cTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGa                                                                                               >  2:1885818/1‑148 (MQ=255)
                                                     cTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGa                                                                                               >  2:3292907/1‑148 (MQ=255)
                                                      tCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGaa                                                                                              <  2:2407467/148‑1 (MQ=255)
                                                      tCCGCGCTACACTCCTGCGCCAGCATCAGGAAGTCGTCCGGCGGCCCCCTGGCCGAATGCGCAGCCAGGATACCACAAACCGCAAGTTGAAGTGTATCGCAAGCTTTTGAAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGaa                                                                                              <  2:677654/148‑1 (MQ=255)
                                                       ccGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAAc                                                                                             >  2:1952583/1‑148 (MQ=255)
                                                          cgcTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGc                                                                                          >  1:861460/1‑148 (MQ=255)
                                                          cgcTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGc                                                                                          >  1:2776112/1‑148 (MQ=255)
                                                          cgcTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGc                                                                                          >  1:1837893/1‑148 (MQ=255)
                                                             tACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAgg                                                                                       >  2:1357091/1‑148 (MQ=255)
                                                              acacTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGa                                                                                      >  1:3422887/1‑148 (MQ=255)
                                                               cacTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGaa                                                                                     >  1:2988578/1‑148 (MQ=255)
                                                                acTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAAc                                                                                    <  2:866361/148‑1 (MQ=255)
                                                                acTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAAc                                                                                    >  1:2656148/1‑148 (MQ=255)
                                                                  tCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCt                                                                                  >  1:2795499/1‑148 (MQ=255)
                                                                   ccTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTg                                                                                 >  1:3071218/1‑148 (MQ=255)
                                                                   ccTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTg                                                                                 >  1:3070819/1‑148 (MQ=255)
                                                                   ccTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTg                                                                                 >  2:3470201/1‑148 (MQ=255)
                                                                                 aGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACCGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAAc                                                                   <  1:3074979/148‑1 (MQ=255)
                                                                                  ggCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACa                                                                  <  2:1717322/148‑1 (MQ=255)
                                                                                   gCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAg                                                                 >  1:1779288/1‑148 (MQ=255)
                                                                                    cAGTGGTCCGGCGGCCCCCTGGCGGAATGCGCAGCCAGGCTACCACAAACAGCACGGTGAAGAGGAGCGCCAGATTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGt                                                                <  1:3391524/148‑1 (MQ=255)
                                                                                     aGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCACGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTc                                                               <  2:2248766/148‑1 (MQ=255)
                                                                                       tcgtcCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCag                                                             >  2:1427712/1‑148 (MQ=255)
                                                                                        cgtcCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCaga                                                            >  2:3591943/1‑148 (MQ=255)
                                                                                           ccGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAACCGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGtt                                                         <  2:223443/148‑1 (MQ=255)
                                                                                             ggCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGtttt                                                       <  2:1288874/148‑1 (MQ=255)
                                                                                               cTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTgg                                                     <  2:3071218/148‑1 (MQ=255)
                                                                                               cTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTgg                                                     <  2:3070819/148‑1 (MQ=255)
                                                                                               cTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTgg                                                     <  2:1902273/148‑1 (MQ=255)
                                                                                                tGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGt                                                    >  1:1367301/1‑148 (MQ=255)
                                                                                                    cccTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCggg                                                >  1:98355/1‑148 (MQ=255)
                                                                                                     ccTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGt                                               >  2:3498475/1‑148 (MQ=255)
                                                                                                     ccTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGt                                               >  2:956170/1‑148 (MQ=255)
                                                                                                      cTGGCCGACTGCGCAGCCAGGCTACCAGACACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTa                                              <  2:483353/148‑1 (MQ=255)
                                                                                                      cTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCATGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTa                                              <  2:2647203/148‑1 (MQ=255)
                                                                                                      cTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTa                                              <  2:3602320/148‑1 (MQ=255)
                                                                                                      cTGGCCGAATGCGCAGCCAGGCTACCACAACCAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTa                                              <  1:2366312/148‑1 (MQ=255)
                                                                                                         gCCGAATGCGCAGCCAGGCTACCAGAAACAGCACGGTGAAGCGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTg                                           <  2:2276128/148‑1 (MQ=255)
                                                                                                            gAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACt                                        <  1:2379684/148‑1 (MQ=255)
                                                                                                             aaTGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTa                                       <  2:579134/148‑1 (MQ=255)
                                                                                                             aaTGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTa                                       <  2:3522927/148‑1 (MQ=255)
                                                                                                             aaTGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTa                                       >  1:1460282/1‑148 (MQ=255)
                                                                                                               tGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATc                                     >  1:1225358/1‑148 (MQ=255)
                                                                                                               tGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATc                                     >  1:2538473/1‑148 (MQ=255)
                                                                                                                 cgcAGTCAGGCAACCACAACCAGCCTGGTGATGTTTATCGCCAGTTTTGCCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGt                                   <  1:1291728/148‑1 (MQ=255)
                                                                                                                  gcAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTc                                  <  2:494020/148‑1 (MQ=255)
                                                                                                                     gCCTGGCAACCAGAAACCCAAAGTTTAAGAGTTTCGCCAGTTTTGGCAAACGCCCGCGCCGATGAAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCAt                               <  1:681064/148‑1 (MQ=255)
                                                                                                                     gCCCGGATACCAGAAACAGAAAGGTGAAGTGTATCGCCAGATTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCTCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCAt                               <  2:1638429/148‑1 (MQ=255)
                                                                                                                         ggCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCg                           <  1:1565038/148‑1 (MQ=255)
                                                                                                                           ccaCCATAAACCGCACGGTGATGTGTACCGCAAGCTTTTCCAACCGCCCGCGCCGAGGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGAc                         <  1:17091/146‑1 (MQ=255)
                                                                                                                             aCCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGACCCTGCACTTGCGGCGAACAGTCAGAGGTTTGGTCGGGTACTGACTATCGTCCATTCCGACtg                       >  1:506214/1‑148 (MQ=255)
                                                                                                                                aGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGc                    >  1:586271/1‑148 (MQ=255)
                                                                                                                                aGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCAATCCGACTGTGc                    >  1:1989581/1‑148 (MQ=255)
                                                                                                                                 gAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCg                   >  1:3167699/1‑148 (MQ=255)
                                                                                                                                  aaaCAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGt                  >  1:3217639/1‑148 (MQ=255)
                                                                                                                                    ccAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTg                <  1:378574/147‑1 (MQ=255)
                                                                                                                                    accgCACGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTg                <  1:378564/145‑1 (MQ=255)
                                                                                                                                    acacCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTg                <  1:3642978/148‑1 (MQ=255)
                                                                                                                                    aCAGCACGGTGAAGTGTATCGCCTGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTg                <  1:3078367/148‑1 (MQ=255)
                                                                                                                                    aCAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTg                >  1:2207625/1‑148 (MQ=255)
                                                                                                                                    aCAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTg                <  1:3182044/148‑1 (MQ=255)
                                                                                                                                     cAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGc               <  1:2087595/148‑1 (MQ=255)
                                                                                                                                        cAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTgg            >  1:3622110/1‑148 (MQ=255)
                                                                                                                                        cAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTgg            >  2:2877045/1‑148 (MQ=255)
                                                                                                                                        cAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTgg            >  2:2039600/1‑148 (MQ=255)
                                                                                                                                        cAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTgg            >  1:268087/1‑148 (MQ=255)
                                                                                                                                             tGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAgcgc       >  1:1839706/1‑148 (MQ=255)
                                                                                                                                              gAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGcgcg      <  1:956170/148‑1 (MQ=255)
                                                                                                                                                  tgtATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACgg  >  1:697069/1‑148 (MQ=255)
                                                                                                                                                  tgtATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACgg  >  1:1070970/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGG  >  CP009273/4388293‑4388586

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: