Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,388,902 G→T 90.8% A558A (GCG→GCT mutL → methyl‑directed mismatch repair protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,388,9020GT90.8% 895.1 / 79.6 316A558A (GCG→GCTmutLmethyl‑directed mismatch repair protein
Reads supporting (aligned to +/- strand):  ref base G (11/18);  new base T (176/111);  total (187/129)
Fisher's exact test for biased strand distribution p-value = 1.75e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 316 total aligned reads displayed.

TGGGTATTGATTTCCAGTCAGATGCACAGCATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCC  >  CP009273/4388757‑4389048
                                                                                                                                                 |                                                                                                                                                  
  ggTATTGATTTCCAGTCAGATGCACAGCATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTgg                                                                                                                                              <  1:3185020/150‑1 (MQ=255)
     aTTGATTTCCAGTCAGATGCACAGCATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGAtt                                                                                                                                           >  1:1596980/1‑150 (MQ=255)
      ttGATTTCCAGTCAGATGCACAGCATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTg                                                                                                                                          >  2:249043/1‑150 (MQ=255)
       tGATTTCCAGTCAGATGCACAGCATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGc                                                                                                                                         >  2:3473004/1‑150 (MQ=255)
         aTTTCCAGTCAGATGCACAGCATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCAc                                                                                                                                       <  1:177915/150‑1 (MQ=255)
               aGTCAGATGCACAGCATGTGACCTTCAGTGCAGTTCCTTTACCCTTTCTACTACAAATTTTACAAATCTTGATTCCTGATCTGAGAGGCTAACTGGCGAAGCAGTCCGTATTAGAACCTGGCTATATTGCTCAGTGGATTGCACGatatg                                                                                                                                 <  2:3928423/150‑2 (MQ=255)
               aGTCAGATGCACAGCATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATc                                                                                                                                 >  1:1990471/1‑150 (MQ=255)
               aGTCAGATGCACAGCATGTGACCATCAGGGCAGTGCCTTTACCCGTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATc                                                                                                                                 <  2:4427928/150‑1 (MQ=255)
                  cAGATGCACAGCATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCtga                                                                                                                              >  1:3393492/1‑150 (MQ=255)
                         acaGCATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGa                                                                                                                       <  2:3030890/150‑1 (MQ=255)
                           aGCATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGacc                                                                                                                     >  1:2509024/1‑148 (MQ=255)
                              aTGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATg                                                                                                                  >  2:4113891/1‑150 (MQ=255)
                                 tgACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATgcgc                                                                                                               <  1:2951723/150‑1 (MQ=255)
                                  gACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCa                                                                                                              >  2:58937/1‑150 (MQ=255)
                                     cATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTg                                                                                                           >  1:1185971/1‑150 (MQ=255)
                                       tCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGt                                                                                                         <  1:3251585/150‑1 (MQ=255)
                                        cAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTc                                                                                                        >  2:2104388/1‑150 (MQ=255)
                                         aGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCa                                                                                                       >  1:2455144/1‑150 (MQ=255)
                                           ggCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAAt                                                                                                     <  1:1797524/150‑1 (MQ=255)
                                            gCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATg                                                                                                    <  1:3575616/150‑1 (MQ=255)
                                            gCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATg                                                                                                    >  2:1847805/1‑150 (MQ=255)
                                              aGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGc                                                                                                  <  1:3799773/150‑1 (MQ=255)
                                               gTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGca                                                                                                 >  1:1906221/1‑150 (MQ=255)
                                                tGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGcac                                                                                                <  1:327106/150‑1 (MQ=255)
                                                 gCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGcaca                                                                                               >  2:1175820/1‑150 (MQ=255)
                                                   cTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAgg                                                                                             >  1:35055/1‑150 (MQ=255)
                                                       aCCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCAt                                                                                         >  2:1268708/1‑150 (MQ=255)
                                                        cccTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATa                                                                                        <  2:426716/150‑1 (MQ=255)
                                                         ccTTACGACAACAAACTGTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTTTTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATaa                                                                                       <  2:3645083/150‑1 (MQ=255)
                                                           ttACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAAcc                                                                                     <  1:2844709/150‑1 (MQ=255)
                                                             aCGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCct                                                                                   >  1:1365588/1‑150 (MQ=255)
                                                              cGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCctg                                                                                  <  1:4169267/150‑1 (MQ=255)
                                                               gCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCAGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCctgc                                                                                 >  2:1939543/1‑150 (MQ=255)
                                                                 caacaaAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTTGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCctgctg                                                                               <  1:3701144/150‑1 (MQ=255)
                                                                 caacaaAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCctgctg                                                                               <  1:3670202/150‑1 (MQ=255)
                                                                     aaaaTTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCgg                                                                           >  2:2144727/1‑150 (MQ=255)
                                                                     aaaaTTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGAAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCgg                                                                           >  2:2831121/1‑150 (MQ=255)
                                                                     aaaaTTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCgg                                                                           <  1:2340036/150‑1 (MQ=255)
                                                                      aaaTTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGa                                                                          >  2:141162/1‑150 (MQ=255)
                                                                       aaTTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGAc                                                                         >  2:3785068/1‑150 (MQ=255)
                                                                        aTTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACg                                                                        >  2:3547655/1‑150 (MQ=255)
                                                                          ttACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTg                                                                      <  2:455837/150‑1 (MQ=255)
                                                                          ttACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTg                                                                      <  1:4113852/150‑1 (MQ=255)
                                                                          ttACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTg                                                                      <  1:777261/150‑1 (MQ=255)
                                                                             cAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGaa                                                                   <  2:3517501/150‑1 (MQ=255)
                                                                              aaaTCTTGATTCCTGAACTGATCGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAAc                                                                  >  1:2297764/1‑150 (MQ=255)
                                                                               aaTCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACg                                                                 <  1:1183840/150‑1 (MQ=255)
                                                                               aaTCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACg                                                                 >  2:941872/1‑150 (MQ=255)
                                                                               aaTCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACg                                                                 >  2:4332398/1‑150 (MQ=255)
                                                                               aaTCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACg                                                                 >  2:3257089/1‑150 (MQ=255)
                                                                               aaTCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACg                                                                 >  1:3892944/1‑150 (MQ=255)
                                                                                 tCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGt                                                               >  1:4221070/1‑150 (MQ=255)
                                                                                 tCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGt                                                               >  1:3857518/1‑150 (MQ=255)
                                                                                  cTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGGTTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGtt                                                              >  1:1501370/1‑150 (MQ=255)
                                                                                  cTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATTCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGtt                                                              <  2:3235195/150‑1 (MQ=255)
                                                                                   ttGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTa                                                             >  1:2332939/1‑150 (MQ=255)
                                                                                    tGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTAt                                                            >  2:3859080/1‑150 (MQ=255)
                                                                                    tGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTAt                                                            >  2:360994/1‑150 (MQ=255)
                                                                                    tGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTAt                                                            >  1:1565016/1‑150 (MQ=255)
                                                                                    agATTCCTGAACTGAAAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATTCGCAGTGTTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTag                                                            >  2:54481/2‑149 (MQ=255)
                                                                                      aTTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGt                                                          >  2:1334663/1‑150 (MQ=255)
                                                                                        tCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGACCATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTcc                                                        >  2:372034/1‑150 (MQ=255)
                                                                                          cTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGc                                                      <  1:2173278/150‑1 (MQ=255)
                                                                                           tGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCa                                                     >  1:2677000/1‑150 (MQ=255)
                                                                                           tGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCa                                                     >  1:2677081/1‑150 (MQ=255)
                                                                                            gAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCaa                                                    >  1:4246910/1‑150 (MQ=255)
                                                                                            gAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCaa                                                    >  2:1686033/1‑150 (MQ=255)
                                                                                               cTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACtt                                                 >  2:51957/1‑150 (MQ=255)
                                                                                                tGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTtg                                                <  2:3846950/150‑1 (MQ=255)
                                                                                                tGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTtg                                                <  2:2222757/150‑1 (MQ=255)
                                                                                                 gATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTtgt                                               <  2:2261477/150‑1 (MQ=255)
                                                                                                     ggCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGaaa                                           <  2:2890620/150‑1 (MQ=255)
                                                                                                       cTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAAc                                         >  2:4492634/1‑150 (MQ=255)
                                                                                                       cTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAAc                                         <  2:3154004/150‑1 (MQ=255)
                                                                                                        tACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAAcg                                        >  2:1905100/1‑150 (MQ=255)
                                                                                                         aCCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAAcgc                                       >  2:1071144/1‑150 (MQ=255)
                                                                                                         aCCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAAcgc                                       <  1:3836545/150‑1 (MQ=255)
                                                                                                           cTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAAcgccg                                     <  1:915592/150‑1 (MQ=255)
                                                                                                             ggCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACgccgcc                                   <  1:1309843/150‑1 (MQ=255)
                                                                                                             ggCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACgccgcc                                   <  2:1422351/150‑1 (MQ=255)
                                                                                                             ggCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACgccgcc                                   <  2:649950/150‑1 (MQ=255)
                                                                                                               cGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCgg                                 >  2:1732109/1‑150 (MQ=255)
                                                                                                               cGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCgg                                 >  1:288684/1‑150 (MQ=255)
                                                                                                                    cAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGggtggt                            >  2:2481467/1‑150 (MQ=255)
                                                                                                                      gTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCt                          >  1:2869664/1‑150 (MQ=255)
                                                                                                                        ccGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGt                        >  1:4337035/1‑150 (MQ=255)
                                                                                                                         cGTATTCGAACCTGGCAATATTGCTCAGTGGAGTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGtt                       <  2:40809/150‑1 (MQ=255)
                                                                                                                               cGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATc                 >  2:3461947/1‑150 (MQ=255)
                                                                                                                               cGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATc                 <  1:4355159/150‑1 (MQ=255)
                                                                                                                                gAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCt                >  2:1518074/1‑150 (MQ=255)
                                                                                                                                gAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAAATTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCt                >  2:1570290/1‑150 (MQ=255)
                                                                                                                                  aCCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGt              <  2:1044071/150‑1 (MQ=255)
                                                                                                                                  aCCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGt              >  1:1002292/1‑150 (MQ=255)
                                                                                                                                    cTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTg            >  1:2099401/1‑150 (MQ=255)
                                                                                                                                     tGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGa           >  2:1723186/1‑150 (MQ=255)
                                                                                                                                      ggCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGAt          >  2:1184074/1‑150 (MQ=255)
                                                                                                                                       gCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGAtt         >  2:869153/1‑150 (MQ=255)
                                                                                                                                         aaTATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTa       >  1:1553170/1‑150 (MQ=255)
                                                                                                                                           tatTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACa     <  1:640752/150‑1 (MQ=255)
                                                                                                                                            atTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACAt    >  2:4233474/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
TGGGTATTGATTTCCAGTCAGATGCACAGCATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCC  >  CP009273/4388757‑4389048

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: