Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,388,439 T→A 94.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,4390TA94.0% 1001.9 / 39.4 328V404E (GTG→GAG) mutLmethyl‑directed mismatch repair protein
Reads supporting (aligned to +/- strand):  ref base T (12/8);  new base A (159/149);  total (171/157)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.45e-01

Only 100 of 328 total aligned reads displayed.

CGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACG  >  CP009273/4388293‑4388585
                                                                                                                                                  |                                                                                                                                                  
cgTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTa                                                                                                                                                 >  1:3829917/1‑150 (MQ=255)
cgTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTa                                                                                                                                                 <  1:4293419/150‑1 (MQ=255)
  tGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACTCTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCTGGCTACCAGAAACAGCAAGGTGAAGAGTATc                                                                                                                                               <  2:875542/150‑1 (MQ=255)
  tGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATc                                                                                                                                               >  1:3140215/1‑150 (MQ=255)
  tGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATc                                                                                                                                               >  1:1147171/1‑150 (MQ=255)
   ggcggcGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCg                                                                                                                                              >  2:2643179/1‑150 (MQ=255)
   ggcggcGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCg                                                                                                                                              >  1:4448766/1‑150 (MQ=255)
         ggggCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCt                                                                                                                                        >  1:254975/1‑150 (MQ=255)
              gcAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGc                                                                                                                                   >  2:1972942/1‑150 (MQ=255)
                aaTCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCaa                                                                                                                                 <  1:4020064/150‑1 (MQ=255)
                  tCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAAc                                                                                                                               >  1:2904870/1‑150 (MQ=255)
                    aCTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGATACAGTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGc                                                                                                                             <  1:2221723/150‑1 (MQ=255)
                        tGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCcgc                                                                                                                         <  2:2390550/150‑1 (MQ=255)
                        tGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCcgc                                                                                                                         >  1:2955236/1‑150 (MQ=255)
                           aGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGTCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGcc                                                                                                                      >  2:3470556/1‑150 (MQ=255)
                              aCCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGAt                                                                                                                   <  1:4841085/150‑1 (MQ=255)
                               ccGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATg                                                                                                                  <  2:3884413/150‑1 (MQ=255)
                                 ggCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCa                                                                                                                >  1:783359/1‑150 (MQ=255)
                                 ggCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCa                                                                                                                >  1:3739271/1‑150 (MQ=255)
                                  gCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCaa                                                                                                               <  2:1335898/150‑1 (MQ=255)
                                  gCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCaa                                                                                                               <  1:1071526/150‑1 (MQ=255)
                                  gCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCaa                                                                                                               <  2:5161190/150‑1 (MQ=255)
                                    aGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCaaaa                                                                                                             <  2:2446676/150‑1 (MQ=255)
                                        cGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTa                                                                                                         >  2:3035672/1‑150 (MQ=255)
                                        cGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTa                                                                                                         <  1:1121216/150‑1 (MQ=255)
                                         gTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTaa                                                                                                        <  1:4475294/150‑1 (MQ=255)
                                         gTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTaa                                                                                                        >  2:2126076/1‑150 (MQ=255)
                                         gTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGAAAGTTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTaa                                                                                                        <  1:1488518/150‑1 (MQ=255)
                                               cGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAgcgc                                                                                                  <  2:5125480/150‑1 (MQ=255)
                                                   aGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCACGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGa                                                                                              <  2:2198621/150‑1 (MQ=255)
                                                    gCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGaa                                                                                             >  1:3357/1‑150 (MQ=255)
                                                         gcgcTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCa                                                                                        >  2:1831101/1‑150 (MQ=255)
                                                          cgcTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCCGAAACAGCAAGGGGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCCAAAATTAAAAGCGCCGGAACCGCAg                                                                                       >  2:2932093/1‑150 (MQ=255)
                                                          cgcTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAg                                                                                       <  1:1970064/150‑1 (MQ=255)
                                                           gcTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAgg                                                                                      >  2:3751634/1‑150 (MQ=255)
                                                            cTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGAAAGTTTAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGa                                                                                     <  1:2703637/150‑1 (MQ=255)
                                                              acacTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAAc                                                                                   <  1:47883/150‑1 (MQ=255)
                                                              acacTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAAGTAAAAGCGCCGGAACCGCAGGAAc                                                                                   >  2:4845107/1‑150 (MQ=255)
                                                               cacTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAAcc                                                                                  <  1:2838175/150‑1 (MQ=255)
                                                                 cTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTg                                                                                <  1:550526/150‑1 (MQ=255)
                                                                  tCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGc                                                                               <  1:635675/150‑1 (MQ=255)
                                                                   ccTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCa                                                                              >  2:2725199/1‑150 (MQ=255)
                                                                    cTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCAc                                                                             >  1:1841727/1‑150 (MQ=255)
                                                                     tGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACt                                                                            >  2:90617/1‑150 (MQ=255)
                                                                     tGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCTGGAACCGCAGGAACCTGCACt                                                                            <  1:4351182/150‑1 (MQ=255)
                                                                      gcgcCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCACAGGAACCTGCACtt                                                                           >  1:3042428/1‑150 (MQ=255)
                                                                       cgcCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTg                                                                          >  2:1173456/1‑150 (MQ=255)
                                                                         ccAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTgcg                                                                        <  2:371483/150‑1 (MQ=255)
                                                                         ccAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTgcg                                                                        >  1:1581369/1‑150 (MQ=255)
                                                                          cAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTgcgg                                                                       <  1:2307217/150‑1 (MQ=255)
                                                                          cAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTgcgg                                                                       <  1:3226592/150‑1 (MQ=255)
                                                                           aGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTgcggc                                                                      >  2:1962052/1‑150 (MQ=255)
                                                                             cATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGa                                                                    >  2:713417/1‑150 (MQ=255)
                                                                             cATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGa                                                                    <  2:4872728/150‑1 (MQ=255)
                                                                              aTCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGaa                                                                   >  1:1767117/1‑150 (MQ=255)
                                                                               tCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAAc                                                                  >  2:1269496/1‑150 (MQ=255)
                                                                                cAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACa                                                                 <  1:2813471/150‑1 (MQ=255)
                                                                                cAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACa                                                                 >  2:1547499/1‑150 (MQ=255)
                                                                                  ggCAGTCGTCCTGCTGCACCCTGGCCTAATGCTCAGCCATGCTTCCAGAAACAGCTAGGTTATGAGTATCTGCAGCTTTTGCATACGCCCTAGCCGATGCAAAAATTATTAGCGCCGTTTCCGCAGGAACCTGCACTTGCGGCGAACatt                                                               <  2:3034106/150‑3 (MQ=255)
                                                                                   gCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTc                                                              <  1:3652780/150‑1 (MQ=255)
                                                                                    cAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCa                                                             >  2:4015632/1‑150 (MQ=255)
                                                                                     aGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCTCTTGCGGCGAACAGTCag                                                            <  1:2121396/150‑1 (MQ=255)
                                                                                      gtcgtcCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCaga                                                           >  2:2019923/1‑150 (MQ=255)
                                                                                       tcgtcCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCagag                                                          >  2:2309265/1‑150 (MQ=255)
                                                                                               cTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTc                                                  >  2:3309258/1‑150 (MQ=255)
                                                                                               cTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTc                                                  <  2:1626950/150‑1 (MQ=255)
                                                                                                   ccccTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGt                                              >  1:4699520/1‑150 (MQ=255)
                                                                                                    cccTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTa                                             >  2:4352186/1‑150 (MQ=255)
                                                                                                     ccTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTAc                                            <  2:349816/150‑1 (MQ=255)
                                                                                                      cTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACt                                           >  2:3555016/1‑150 (MQ=255)
                                                                                                      cTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACt                                           >  2:4952174/1‑150 (MQ=255)
                                                                                                       tGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTg                                          >  1:1100567/1‑150 (MQ=255)
                                                                                                       tGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTg                                          <  1:2019923/150‑1 (MQ=255)
                                                                                                       tGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTg                                          >  1:2089521/1‑150 (MQ=255)
                                                                                                        ggCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGa                                         >  1:1808635/1‑150 (MQ=255)
                                                                                                             aaTGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATc                                    >  2:3980601/1‑150 (MQ=255)
                                                                                                               tGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGt                                  >  1:2993282/1‑150 (MQ=255)
                                                                                                                 cgcAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTcc                                >  2:3657764/1‑150 (MQ=255)
                                                                                                                  gcAGCCAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCa                               >  2:62135/1‑150 (MQ=255)
                                                                                                                       cAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCg                          >  2:4333727/1‑150 (MQ=255)
                                                                                                                       cAGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATGCCg                          >  2:4742284/1‑150 (MQ=255)
                                                                                                                        aGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGa                         >  1:4648593/1‑150 (MQ=255)
                                                                                                                        aGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGa                         <  2:5102275/150‑1 (MQ=255)
                                                                                                                        aGGCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGa                         >  2:455873/1‑150 (MQ=255)
                                                                                                                         ggCTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCATCCATTCCGAc                        >  1:3507671/1‑150 (MQ=255)
                                                                                                                           cTACCAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACtg                      >  1:3168115/1‑150 (MQ=255)
                                                                                                                               cAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTCCTGACTATCGTCCATTCCGACTGTGCg                  >  2:4915259/1‑150 (MQ=255)
                                                                                                                               cAGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCg                  <  2:4248891/150‑1 (MQ=255)
                                                                                                                                aGAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGt                 >  1:1358286/1‑150 (MQ=255)
                                                                                                                                 gAAACAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGtt                <  2:448746/150‑1 (MQ=255)
                                                                                                                                  aaaCAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTg               >  1:3646295/1‑150 (MQ=255)
                                                                                                                                  aaaCAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTg               <  1:953592/150‑1 (MQ=255)
                                                                                                                                    aCAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCt             <  1:3286866/150‑1 (MQ=255)
                                                                                                                                    aCAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCt             <  1:1627932/150‑1 (MQ=255)
                                                                                                                                     cAGCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTg            <  1:3273747/150‑1 (MQ=255)
                                                                                                                                       gCAAGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGa          <  2:1135955/150‑1 (MQ=255)
                                                                                                                                         aaGGTGAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAgc        <  1:2126187/150‑1 (MQ=255)
                                                                                                                                              gAAGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGAc   >  1:330328/1‑150 (MQ=255)
                                                                                                                                               aaGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACg  <  2:226412/150‑1 (MQ=255)
                                                                                                                                               aaGAGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACg  >  2:4170380/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CGTGGCGGCGGGGCGCAATCACTTTGCAGAACCGGCAGCTCGTGAGCCGGTAGCTCCGCGCTACACTCCTGCGCCAGCATCAGGCAGTCGTCCGGCTGCCCCCTGGCCGAATGCGCAGCCAGGCTACCAGAAACAGCAAGGTGAAGTGTATCGCCAGCTTTTGCAAACGCCCGCGCCGATGCAAAAATTAAAAGCGCCGGAACCGCAGGAACCTGCACTTGCGGCGAACAGTCAGAGTTTTGGTCGGGTACTGACTATCGTCCATTCCGACTGTGCGTTGCTGGAGCGCGACG  >  CP009273/4388293‑4388585

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: