Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,388,439 T→A 70.0% 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,4390TA70.0% 395.6 / 223.1 275V404E (GTG→GAG) mutLmethyl‑directed mismatch repair protein
Reads supporting (aligned to +/- strand):  ref base T (40/43);  new base A (104/88);  total (144/131)
Fisher's exact test for biased strand distribution p-value = 4.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 282 total aligned reads displayed.

CGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCA  >  CP009273/4388291‑4388588
                                                                                                                                                    |                                                                                                                                                     
 gcgTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAgagt                                                                                                                                                     <  1:445717/150‑1 (MQ=255)
  cgTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGACACAGCAAGGTGAAgtgtc                                                                                                                                                    >  2:5316440/1‑149 (MQ=255)
  cgTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTa                                                                                                                                                    <  2:3558075/150‑1 (MQ=255)
   gTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTAt                                                                                                                                                   >  2:4814034/1‑150 (MQ=255)
    tGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATc                                                                                                                                                  <  2:4003688/150‑1 (MQ=255)
    tGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATc                                                                                                                                                  <  2:2342558/150‑1 (MQ=255)
       cggcggGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGcc                                                                                                                                               >  2:721400/1‑150 (MQ=255)
        ggcggGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCa                                                                                                                                              >  2:1342144/1‑150 (MQ=255)
             ggCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCttt                                                                                                                                         >  1:1704280/1‑150 (MQ=255)
                 cAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCa                                                                                                                                     <  1:2224211/150‑1 (MQ=255)
                  aaTCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCaa                                                                                                                                    >  2:4102056/1‑150 (MQ=255)
                  aaTCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCaa                                                                                                                                    >  1:2548592/1‑150 (MQ=255)
                    tCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAAc                                                                                                                                  <  1:3998425/150‑1 (MQ=255)
                    tCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAAc                                                                                                                                  <  1:2177168/150‑1 (MQ=255)
                    tCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAAc                                                                                                                                  <  1:146710/150‑1 (MQ=255)
                       cTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGcc                                                                                                                               >  2:1731576/1‑150 (MQ=255)
                        tttGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGccc                                                                                                                              <  1:35446/150‑1 (MQ=255)
                         ttGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCcg                                                                                                                             >  1:4626479/1‑150 (MQ=255)
                           gCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCcgcg                                                                                                                           >  2:2100368/1‑150 (MQ=255)
                            cAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCgcgc                                                                                                                          >  1:1905648/1‑150 (MQ=255)
                               aaCCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCTGCATCAGGCAGTCGTCGGGCTGCCCCCTGGCCGTTTGCGCACCCAGGCTACCAGAACCAGCAAGTTGATGAGTTTCGCCAGATTTTGCAAACGCCCGCGCCGa                                                                                                                       <  1:1563113/150‑1 (MQ=255)
                                    gCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCaa                                                                                                                  <  2:3630385/150‑1 (MQ=255)
                                     cAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCaaa                                                                                                                 >  2:3858864/1‑150 (MQ=255)
                                     cAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCaaa                                                                                                                 >  1:3107593/1‑150 (MQ=255)
                                        cTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAAt                                                                                                              <  2:3152272/150‑1 (MQ=255)
                                            tGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTaaa                                                                                                          <  2:1828299/150‑1 (MQ=255)
                                              aGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAg                                                                                                        <  2:2802219/150‑1 (MQ=255)
                                              aGCCGGTAGCTCCGCGCGACACTCCTGCGCCAGCCTCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGACACCGCGAGGTGACGTGGAGCGCCAGCTTTTGCAAACGCCCGCGCCGCTGCACACATTacacg                                                                                                        >  2:621980/1‑146 (MQ=255)
                                               gCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCATAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAgc                                                                                                       <  2:2492084/150‑1 (MQ=255)
                                               gCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAgc                                                                                                       <  2:4516755/150‑1 (MQ=255)
                                               gCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAgc                                                                                                       <  1:1429124/150‑1 (MQ=255)
                                                  ggTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGcc                                                                                                    <  2:1576003/150‑1 (MQ=255)
                                                     aGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGa                                                                                                 >  2:2185051/1‑150 (MQ=255)
                                                       cTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAAc                                                                                               >  1:2952949/1‑150 (MQ=255)
                                                        tCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGACACAGCAAGGGGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGAGGCAAAAATTAAAAGCGCCGGAAcc                                                                                              >  2:3637031/1‑150 (MQ=255)
                                                         ccGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCg                                                                                             >  2:4903716/1‑150 (MQ=255)
                                                           gcgcTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCa                                                                                           <  2:2032110/150‑1 (MQ=255)
                                                           gcgcTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCa                                                                                           >  1:3670753/1‑150 (MQ=255)
                                                           gcgcTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCa                                                                                           <  1:387975/150‑1 (MQ=255)
                                                             gcTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAgg                                                                                         >  2:2861054/1‑150 (MQ=255)
                                                             gcTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAAGTAAAAGCGCCGGAACCGCAgg                                                                                         >  2:4042488/1‑150 (MQ=255)
                                                                acacTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAAc                                                                                      >  2:2427340/1‑150 (MQ=255)
                                                                acacTCCGGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAAc                                                                                      <  2:4208775/150‑1 (MQ=255)
                                                                    tCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGc                                                                                  >  2:4505815/1‑150 (MQ=255)
                                                                              gCATCAGGCAGTCGTCCGGCTGCCCCCGGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTgcggcg                                                                        <  2:5072391/150‑1 (MQ=255)
                                                                                aTCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGaa                                                                      >  2:3828674/1‑150 (MQ=255)
                                                                                 tCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAAc                                                                     <  2:842465/150‑1 (MQ=255)
                                                                                 tCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAAc                                                                     >  1:2686046/1‑150 (MQ=255)
                                                                                   aGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCATAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAg                                                                   <  1:2071018/150‑1 (MQ=255)
                                                                                   aGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGAGGCAAAAATTAAAAGCGCCGGAACCGCAGGACCCTGCACTTGCGGCGAACAg                                                                   >  2:32542/1‑150 (MQ=255)
                                                                                    ggCTGTCGTCGGGCTGCCCCCGGGCCGAATGCGCAGCCTGGCTACCATAAACACCATGGTGAAGTTTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGt                                                                  <  1:3648002/150‑1 (MQ=255)
                                                                                    ggCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGACAGCTTTTGCAAACTCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAga                                                                  >  1:5216556/1‑149 (MQ=255)
                                                                                      cAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCa                                                                >  1:3935185/1‑150 (MQ=255)
                                                                                       aGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCag                                                               >  1:4476272/1‑150 (MQ=255)
                                                                                       aGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCag                                                               >  1:3900464/1‑150 (MQ=255)
                                                                                          cgtcCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGt                                                            >  1:559493/1‑150 (MQ=255)
                                                                                           gtcCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGtt                                                           >  2:4350500/1‑150 (MQ=255)
                                                                                             cgggCTGCCCCCTGGCCTAATGCGCAGCCTGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGtttt                                                         <  1:4159771/148‑1 (MQ=255)
                                                                                               ggCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTgg                                                       >  2:2802477/1‑150 (MQ=255)
                                                                                                gCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGt                                                      <  1:4842042/150‑1 (MQ=255)
                                                                                                  tGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCg                                                    >  1:974677/1‑150 (MQ=255)
                                                                                                  tGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCAGCACTTGCGGCGAACAGTCAGAGTTTTGGTCg                                                    >  1:1296716/1‑150 (MQ=255)
                                                                                                  tGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCg                                                    >  1:5020883/1‑150 (MQ=255)
                                                                                                     ccccTGGCCGAATGCGCAGCCAGGCTACCAGAACCAGCACGGTGAAGAGTATCGCCAGATTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGt                                                 <  1:2441619/150‑1 (MQ=255)
                                                                                                     ccccTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGt                                                 >  2:4003961/1‑150 (MQ=255)
                                                                                                         tGGCCGAATGCGCAGCCTGGCTACCAGAAACCGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTg                                             <  2:2544128/150‑1 (MQ=255)
                                                                                                         tGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTg                                             >  2:3396440/1‑150 (MQ=255)
                                                                                                          ggCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGa                                            >  2:1349946/1‑150 (MQ=255)
                                                                                                          ggCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGa                                            >  2:462216/1‑150 (MQ=255)
                                                                                                          ggCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTGTTGGTCGGGTACTGa                                            >  2:246027/1‑150 (MQ=255)
                                                                                                           gCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGAc                                           <  2:140583/150‑1 (MQ=255)
                                                                                                            ccGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACt                                          >  2:1245856/1‑150 (MQ=255)
                                                                                                              gAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTTAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTAt                                        <  1:1676816/150‑1 (MQ=255)
                                                                                                              gAATGCGCAGCCAGGCTAACAGAAACCGCATGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTAt                                        <  1:885900/150‑1 (MQ=255)
                                                                                                               aaTGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATc                                       >  1:3544335/1‑150 (MQ=255)
                                                                                                                    gcAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCa                                  >  2:949607/1‑150 (MQ=255)
                                                                                                                    gcAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCa                                  <  1:103000/150‑1 (MQ=255)
                                                                                                                     cAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCAt                                 >  2:2995573/1‑150 (MQ=255)
                                                                                                                        cccggCTACCATAACCCGCAAGGTGAAGAGTATCGCCCGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTcc                              <  1:1846670/147‑1 (MQ=255)
                                                                                                                        ccAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTcc                              <  2:1108308/150‑1 (MQ=255)
                                                                                                                        ccAGGCTACCAGAAACAGCAAGGTGAAGAGAAACGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTcc                              <  2:2631156/150‑1 (MQ=255)
                                                                                                                             cTACCAGAACCCGCACGGTGACGAGTATCGCCTGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACtg                         <  1:5077659/150‑1 (MQ=255)
                                                                                                                              tACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACtgt                        >  2:2250026/1‑150 (MQ=255)
                                                                                                                              tACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACtgt                        <  1:638767/150‑1 (MQ=255)
                                                                                                                               aCCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACtgtg                       >  1:3458469/1‑150 (MQ=255)
                                                                                                                                  aGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGt                    >  2:5257023/1‑150 (MQ=255)
                                                                                                                                  aGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGt                    >  2:369460/1‑150 (MQ=255)
                                                                                                                                   gAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGtt                   <  1:4355872/150‑1 (MQ=255)
                                                                                                                                      aCAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCt                <  1:721400/150‑1 (MQ=255)
                                                                                                                                        aGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTgg              >  1:1172079/1‑150 (MQ=255)
                                                                                                                                        aGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTgg              >  1:3313606/1‑150 (MQ=255)
                                                                                                                                          cAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAg            <  2:5074702/150‑1 (MQ=255)
                                                                                                                                           acggTGAAGAGTAACGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAgc           <  1:4076061/148‑1 (MQ=255)
                                                                                                                                           aaGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAgc           >  1:4528520/1‑150 (MQ=255)
                                                                                                                                           aaGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAgc           >  1:2607431/1‑150 (MQ=255)
                                                                                                                                           aaGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAgc           >  2:2762660/1‑150 (MQ=255)
                                                                                                                                            aGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAgcg          >  1:1199749/1‑150 (MQ=255)
                                                                                                                                             ggTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAgcgc         >  2:4931184/1‑150 (MQ=255)
                                                                                                                                                 aaGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACg     <  1:5316440/150‑1 (MQ=255)
                                                                                                                                                 aaGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACg     >  2:5173677/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CGCGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACGGCA  >  CP009273/4388291‑4388588

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