Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,388,902 G→T 100% 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,9020GT97.1% 1019.7 / 18.5 306A558A (GCG→GCTmutLmethyl‑directed mismatch repair protein
Reads supporting (aligned to +/- strand):  ref base G (5/4);  new base T (171/125);  total (177/129)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 306 total aligned reads displayed.

ATTGGGTATTGATTTCCAGTCAGATGCACAGCATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCC  >  CP009273/4388755‑4389048
                                                                                                                                                   |                                                                                                                                                  
 ttGGGTATTGATTTCCAGTCAGATGCACAGCATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAg                                                                                                                                                 >  2:4406852/1‑150 (MQ=255)
  tGGGTATTGATTTCCAGTCAGATGCACAGCATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGt                                                                                                                                                <  2:790915/150‑1 (MQ=255)
   gggTATTGATTTCCAGTCAGATGCACAGCATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTg                                                                                                                                               >  1:2971265/1‑150 (MQ=255)
       aTTGATTTCCAGTCAGATGCACAGCATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGAtt                                                                                                                                           >  2:2539638/1‑150 (MQ=255)
          gATTTCCAGTCAGATGCACAGCATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCa                                                                                                                                        >  1:4080357/1‑150 (MQ=255)
          gATTTCCAGTCAGATGCACAGCATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCa                                                                                                                                        <  2:1790073/150‑1 (MQ=255)
            tttCCAGTCAGATGCACAGCATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACg                                                                                                                                      >  2:3766996/1‑150 (MQ=255)
            tttCCAGTCAGATGCACAGCATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACg                                                                                                                                      >  1:2542755/1‑150 (MQ=255)
               ccAGTCAGATGCACAGCATGTGACCATCAGNGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGaaa                                                                                                                                   >  1:4279659/1‑150 (MQ=255)
               ccAGTCAGATGCACAGCATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGaaa                                                                                                                                   >  1:1017389/1‑150 (MQ=255)
                cAGTCAGATGCACAGCATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAAt                                                                                                                                  <  1:3553304/150‑1 (MQ=255)
                   tCAGATGCACAGCATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCtg                                                                                                                               >  1:3286465/1‑150 (MQ=255)
                    cAGATGCACAGCATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCtga                                                                                                                              >  2:2714283/1‑150 (MQ=255)
                     aGATGCACAGCATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCtgat                                                                                                                             >  1:805234/1‑150 (MQ=255)
                      gATGCACAGCATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCtaatg                                                                                                                            <  2:3585483/150‑1 (MQ=255)
                       aTGCACAGCATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCtgatga                                                                                                                           <  2:2620499/150‑1 (MQ=255)
                        tGCACAGCATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAg                                                                                                                          <  1:3847368/150‑1 (MQ=255)
                           acaGCATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGa                                                                                                                       >  1:3781370/1‑150 (MQ=255)
                            caGCATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGaa                                                                                                                      <  2:3035316/150‑1 (MQ=255)
                             aGCATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAAc                                                                                                                     <  2:3003495/150‑1 (MQ=255)
                               cATGTGACCATCAGGGCAGTGCCTTTACCCTTTCGCCAACAAAATTTACATATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACTAAATCTGATGAGCGAACAt                                                                                                                   <  2:1720996/150‑1 (MQ=255)
                               cATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACAt                                                                                                                   <  1:2832681/150‑1 (MQ=255)
                               cATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACAt                                                                                                                   >  2:4384236/1‑150 (MQ=255)
                                aTGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATg                                                                                                                  <  2:3104235/150‑1 (MQ=255)
                                aTGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATg                                                                                                                  >  2:143994/1‑150 (MQ=255)
                                aTGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATg                                                                                                                  <  2:546687/150‑1 (MQ=255)
                                 tgtgACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATgc                                                                                                                 >  2:4259396/1‑150 (MQ=255)
                                  gtgACCATCAGGGCAGTGCCTTTACCCTTACGCAAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATgcg                                                                                                                <  2:4127666/150‑1 (MQ=255)
                                      ccATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGt                                                                                                            >  1:2617163/1‑150 (MQ=255)
                                      ccATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGt                                                                                                            >  1:1524605/1‑150 (MQ=255)
                                      ccATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGt                                                                                                            >  2:2155419/1‑150 (MQ=255)
                                       cATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTg                                                                                                           >  1:502306/1‑150 (MQ=255)
                                       cATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTg                                                                                                           <  1:3385469/150‑1 (MQ=255)
                                          cAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCCGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTc                                                                                                        <  1:1225130/150‑1 (MQ=255)
                                             ggCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAAt                                                                                                     <  2:1983983/150‑1 (MQ=255)
                                               cAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATgg                                                                                                   <  1:1373110/150‑1 (MQ=255)
                                                aGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGc                                                                                                  >  2:882256/1‑150 (MQ=255)
                                                   gCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGcaca                                                                                               >  1:1390787/1‑150 (MQ=255)
                                                   gCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGcaca                                                                                               <  2:636810/150‑1 (MQ=255)
                                                    ccTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAg                                                                                              >  1:347203/1‑150 (MQ=255)
                                                       ttACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGcc                                                                                           <  1:2878495/150‑1 (MQ=255)
                                                             ttACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAAcc                                                                                     >  2:4206824/1‑150 (MQ=255)
                                                             ttACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAAcc                                                                                     >  1:3014110/1‑150 (MQ=255)
                                                                 gCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGACATAACCctgc                                                                                 >  2:2321634/1‑150 (MQ=255)
                                                                  ccAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCctgct                                                                                >  1:1233367/1‑150 (MQ=255)
                                                                  ccAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCctgct                                                                                >  2:2786385/1‑150 (MQ=255)
                                                                   caacaaAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCctgctg                                                                               <  2:2101034/150‑1 (MQ=255)
                                                                    aacaaAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTgg                                                                              <  1:2669705/150‑1 (MQ=255)
                                                                       aaaaTTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCgg                                                                           >  1:58972/1‑150 (MQ=255)
                                                                       aaaaTTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCgg                                                                           >  1:3438645/1‑150 (MQ=255)
                                                                        aaaTTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGa                                                                          <  1:2590498/150‑1 (MQ=255)
                                                                        aaaTTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGa                                                                          >  2:4590722/1‑150 (MQ=255)
                                                                        aaaTTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGa                                                                          <  2:3512499/150‑1 (MQ=255)
                                                                        aaaTTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGa                                                                          >  1:3950650/1‑150 (MQ=255)
                                                                          aTTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACg                                                                        >  2:3303975/1‑150 (MQ=255)
                                                                            ttACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTg                                                                      <  2:3822356/150‑1 (MQ=255)
                                                                                aaaTCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAAc                                                                  >  2:3870825/1‑150 (MQ=255)
                                                                                    cTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGtt                                                              >  1:1214979/1‑150 (MQ=255)
                                                                                       gATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATg                                                           >  2:1376352/1‑150 (MQ=255)
                                                                                       gATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATg                                                           >  1:3907709/1‑150 (MQ=255)
                                                                                        aTTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGt                                                          >  2:1107623/1‑150 (MQ=255)
                                                                                           ccTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCg                                                       <  1:2257304/150‑1 (MQ=255)
                                                                                             tGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCa                                                     >  2:445811/1‑150 (MQ=255)
                                                                                             tGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCa                                                     <  1:994682/150‑1 (MQ=255)
                                                                                               aaCTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAAc                                                   <  1:4035828/150‑1 (MQ=255)
                                                                                                   gATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTtgt                                               <  1:558273/150‑1 (MQ=255)
                                                                                                       ggCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGaaa                                           <  2:1154224/150‑1 (MQ=255)
                                                                                                         cTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAAc                                         <  1:3045153/150‑1 (MQ=255)
                                                                                                             cTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAAcgccg                                     >  1:2013402/1‑150 (MQ=255)
                                                                                                             cTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAAcgccg                                     >  1:4525036/1‑150 (MQ=255)
                                                                                                             cTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAAcgccg                                     >  1:3136664/1‑150 (MQ=255)
                                                                                                             cTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAAcgccg                                     <  1:4292636/150‑1 (MQ=255)
                                                                                                               ggCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACgccgcc                                   <  1:1085149/150‑1 (MQ=255)
                                                                                                                gCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACgccgcca                                  <  2:180327/150‑2 (MQ=255)
                                                                                                                  gAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCggg                                <  1:4347295/150‑1 (MQ=255)
                                                                                                                     gCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGggtgg                             >  1:2741596/1‑150 (MQ=255)
                                                                                                                       aGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTc                           >  2:2037263/1‑150 (MQ=255)
                                                                                                                       aGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTc                           >  1:1596343/1‑150 (MQ=255)
                                                                                                                       aGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTc                           >  2:284860/1‑150 (MQ=255)
                                                                                                                          ccGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGt                        >  2:180461/1‑150 (MQ=255)
                                                                                                                          ccGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGt                        >  1:4544791/1‑150 (MQ=255)
                                                                                                                           cGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGtt                       <  1:3321787/150‑1 (MQ=255)
                                                                                                                             tATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTAc                     >  2:4172884/1‑150 (MQ=255)
                                                                                                                              aTTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACa                    <  1:1771611/150‑1 (MQ=255)
                                                                                                                                tCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAAt                  <  2:3964849/150‑1 (MQ=255)
                                                                                                                                   aaCCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTg               >  2:2901249/1‑150 (MQ=255)
                                                                                                                                    aCCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGt              >  2:3251299/1‑150 (MQ=255)
                                                                                                                                      cTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTg            >  2:723503/1‑150 (MQ=255)
                                                                                                                                         gCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGAtt         >  1:3482342/1‑150 (MQ=255)
                                                                                                                                         gCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGAtt         >  1:246946/1‑150 (MQ=255)
                                                                                                                                          cAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGAttt        >  1:4111684/1‑150 (MQ=255)
                                                                                                                                           aaTATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTa       >  2:164706/1‑150 (MQ=255)
                                                                                                                                           aaTATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTa       <  2:155365/150‑1 (MQ=255)
                                                                                                                                            atagTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTAc      <  1:4126788/150‑1 (MQ=255)
                                                                                                                                             tatTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACa     >  1:2900490/1‑150 (MQ=255)
                                                                                                                                             tatTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACa     >  2:4166255/1‑150 (MQ=255)
                                                                                                                                             tatTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACa     >  1:4231058/1‑150 (MQ=255)
                                                                                                                                              atTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACAt    >  2:4005356/1‑150 (MQ=255)
                                                                                                                                              atTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACAt    >  2:3370130/1‑150 (MQ=255)
                                                                                                                                              atTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACAt    >  2:2470057/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
ATTGGGTATTGATTTCCAGTCAGATGCACAGCATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCC  >  CP009273/4388755‑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: