Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,388,902 G→T 93.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,9020GT93.8% 898.8 / 42.0 274A558A (GCG→GCTmutLmethyl‑directed mismatch repair protein
Reads supporting (aligned to +/- strand):  ref base G (9/8);  new base T (145/112);  total (154/120)
Fisher's exact test for biased strand distribution p-value = 8.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

Only 100 of 274 total aligned reads displayed.

ATTGGGTATTGATTTCCAGTCAGATGCACAGCATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACATCC  >  CP009273/4388755‑4389048
                                                                                                                                                   |                                                                                                                                                  
aTTGGGTATTGATTTCCAGTCAGATGCACAGCATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCa                                                                                                                                                  >  2:5111316/1‑150 (MQ=255)
aTTGGGTATTGATTTCCAGTCAGATGCACAGCATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCa                                                                                                                                                  >  1:1574355/1‑150 (MQ=255)
 ttGGGTATTGATTTCCAGTCAGATGCACAGCATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAg                                                                                                                                                 <  1:1693248/150‑1 (MQ=255)
  tGGGTATTGATTTCCAGTCAGATGCACAGCATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGt                                                                                                                                                >  1:2428207/1‑150 (MQ=255)
   gggTATTGATTTCCAGTCAGATGCACAGCATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTg                                                                                                                                               <  2:553653/150‑1 (MQ=255)
   gggTATTGATTTCCAGTCAGATGCACAGCATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTg                                                                                                                                               >  1:3385647/1‑150 (MQ=255)
      tATTGATTTCCAGTCAGATGCACAGCATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGAt                                                                                                                                            >  2:472578/1‑150 (MQ=255)
      tATTGATTTCCAGTCAGATGCACAGCATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGAt                                                                                                                                            <  1:5056550/150‑1 (MQ=255)
      tATTGATTTCCAGTCAGATGCACAGCATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGAt                                                                                                                                            <  1:1449388/150‑1 (MQ=255)
      tATTGATTTCCAGTCAGATGCACAGCATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGAt                                                                                                                                            >  2:2572550/1‑150 (MQ=255)
      tATTGATTTCCAGTCAGATGCACAGCATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGAt                                                                                                                                            <  1:1941183/150‑1 (MQ=255)
       aTTGATTTCCAGTCAGATGCACAGCATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGAtt                                                                                                                                           <  1:2788143/150‑1 (MQ=255)
                cAGTCAGATGCACAGCATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAAt                                                                                                                                  >  2:2657556/1‑150 (MQ=255)
                cAGTCAGATGCACAGCATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTCGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAAt                                                                                                                                  >  2:4792544/1‑150 (MQ=255)
                 aGTCAGATGCACAGCATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATc                                                                                                                                 <  2:3604512/150‑1 (MQ=255)
                 aGTCAGATGCACAGCATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATc                                                                                                                                 >  1:3973690/1‑150 (MQ=255)
                 aGTCAGATGCACAGCATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATc                                                                                                                                 >  1:5164393/1‑150 (MQ=255)
                  gTCAGATGCACAGCATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCt                                                                                                                                >  1:4833418/1‑150 (MQ=255)
                  gTCAGATGCACAGCATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCt                                                                                                                                <  1:3202307/150‑1 (MQ=255)
                    cAGATGCACAGCATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCtga                                                                                                                              >  2:3887923/1‑150 (MQ=255)
                    cAGATGCACAGCATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCtga                                                                                                                              <  1:4655713/150‑1 (MQ=255)
                      gATGCACAGCATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCATGAAATCtgatg                                                                                                                            >  1:2065632/1‑150 (MQ=255)
                      gATGCACAGCATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCtgatg                                                                                                                            <  2:340617/150‑1 (MQ=255)
                       aTGCACAGCATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCggatga                                                                                                                           >  2:4165048/1‑150 (MQ=255)
                         gCACAGCATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGc                                                                                                                         >  2:4111085/1‑150 (MQ=255)
                         gCACAGCATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGc                                                                                                                         >  1:1079339/1‑150 (MQ=255)
                            caGCATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGaa                                                                                                                      <  1:1874613/150‑1 (MQ=255)
                              gCATGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACa                                                                                                                    <  1:1498777/150‑1 (MQ=255)
                                aTGTGACCATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATg                                                                                                                  >  2:2373348/1‑150 (MQ=255)
                                       cATCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTg                                                                                                           >  1:1037966/1‑150 (MQ=255)
                                         tCAGGGCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGt                                                                                                         >  2:343842/1‑150 (MQ=255)
                                            gggCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCaa                                                                                                      >  1:4865692/1‑150 (MQ=255)
                                            gggCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCaa                                                                                                      >  1:3742674/1‑150 (MQ=255)
                                             ggCAGTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAAt                                                                                                     <  2:4454246/150‑1 (MQ=255)
                                                 gTGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGca                                                                                                 >  2:545044/1‑150 (MQ=255)
                                                  tGCCTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGcac                                                                                                >  1:273486/1‑150 (MQ=255)
                                                    ccTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAg                                                                                              <  1:3129325/150‑1 (MQ=255)
                                                    ccTTTACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAg                                                                                              <  2:2428207/150‑1 (MQ=255)
                                                        tACCCTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCa                                                                                          <  2:3540231/150‑1 (MQ=255)
                                                            cTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAAc                                                                                      <  1:1248188/150‑1 (MQ=255)
                                                            cTTACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAAc                                                                                      >  1:5128020/1‑150 (MQ=255)
                                                             ttACGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAAcc                                                                                     >  2:4658773/1‑150 (MQ=255)
                                                               aCGCCAACAAACTTTACAAATCTTGATTCCTGAACTGAGAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACAGGCGCAGTGGTCAATGGCACAGGCCATAACCct                                                                                   <  1:4146093/150‑1 (MQ=255)
                                                                cGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCctg                                                                                  >  2:4409635/1‑150 (MQ=255)
                                                                cGCCAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCctg                                                                                  <  2:1483735/150‑1 (MQ=255)
                                                                  ccAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCctgct                                                                                <  1:4210027/150‑1 (MQ=255)
                                                                  ccAACAAAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCctgct                                                                                >  2:2671622/1‑150 (MQ=255)
                                                                   caacaaAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCctgctg                                                                               >  2:1928779/1‑150 (MQ=255)
                                                                    aacaaAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTgg                                                                              >  1:1391399/1‑150 (MQ=255)
                                                                      caaAATTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCg                                                                            >  2:2742436/1‑150 (MQ=255)
                                                                       aaaaTTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCgg                                                                           >  1:888651/1‑150 (MQ=255)
                                                                        aaaTTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGa                                                                          <  1:3212753/150‑1 (MQ=255)
                                                                          aTTTACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACg                                                                        >  2:2238862/1‑150 (MQ=255)
                                                                           tttACAAATCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGt                                                                       <  1:3061884/150‑1 (MQ=255)
                                                                                  aTCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACgg                                                                <  2:202851/150‑1 (MQ=255)
                                                                                   tCTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGt                                                               <  1:3746306/150‑1 (MQ=255)
                                                                                    cTTGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGtt                                                              <  2:4833418/150‑1 (MQ=255)
                                                                                      tGATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTAt                                                            <  1:3270083/150‑1 (MQ=255)
                                                                                       gATTCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATg                                                           <  2:11849/150‑1 (MQ=255)
                                                                                         ttCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTc                                                         >  2:3830029/1‑150 (MQ=255)
                                                                                          tCCTGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTcc                                                        >  1:341506/1‑150 (MQ=255)
                                                                                             tGAACTGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCa                                                     >  1:3006577/1‑150 (MQ=255)
                                                                                                  tGATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTtg                                                >  1:3502568/1‑150 (MQ=255)
                                                                                                   gATAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTtgt                                               >  2:4565357/1‑150 (MQ=255)
                                                                                                     tAGGCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGa                                             <  2:3020032/150‑1 (MQ=255)
                                                                                                       ggCTACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGaaa                                           <  1:1904326/150‑1 (MQ=255)
                                                                                                          tACCTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAAcg                                        <  2:3993836/150‑1 (MQ=255)
                                                                                                            ccTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAAcgcc                                      <  1:2687336/150‑1 (MQ=255)
                                                                                                            ccTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAAcgcc                                      >  2:5012940/1‑150 (MQ=255)
                                                                                                            ccTGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAAcgcc                                      <  1:2035840/150‑1 (MQ=255)
                                                                                                              tGGCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAAcgccgc                                    >  1:2881843/1‑150 (MQ=255)
                                                                                                               ggCGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACgccgcc                                   >  1:2316404/1‑150 (MQ=255)
                                                                                                                 cGAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCgg                                 <  2:3586695/150‑1 (MQ=255)
                                                                                                                  gAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGTTTATGTCCGCAACTTGTGAAAACGCCGCCggg                                >  1:4732058/1‑150 (MQ=255)
                                                                                                                  gAAGCAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCggg                                >  2:4296227/1‑150 (MQ=255)
                                                                                                                    aGCTGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGAGGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGggtg                              <  2:486970/150‑1 (MQ=255)
                                                                                                                      cAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGggtggt                            >  2:2062593/1‑150 (MQ=255)
                                                                                                                      cAGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGggtggt                            >  1:335678/1‑150 (MQ=255)
                                                                                                                       aGTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTc                           >  1:4529805/1‑150 (MQ=255)
                                                                                                                        gTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATg                                                           >  2:3298772/1‑117 (MQ=255)
                                                                                                                        gTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATg                                                           <  1:3298772/117‑1 (MQ=255)
                                                                                                                        gTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCt                          <  2:366482/150‑1 (MQ=255)
                                                                                                                        gTCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACCAAATCTGATGAGCGAACATGCGCAGTGGACAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCt                          >  2:876556/1‑150 (MQ=255)
                                                                                                                         tCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTg                         <  1:5111316/150‑1 (MQ=255)
                                                                                                                         tCCGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTg                         <  2:2792003/150‑1 (MQ=255)
                                                                                                                          ccGTATTCGAACCTGGCAATATTGCGCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGt                        >  2:1700481/1‑150 (MQ=255)
                                                                                                                           cGTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGtt                       <  2:4150854/150‑1 (MQ=255)
                                                                                                                            gTATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTa                      >  1:1356946/1‑150 (MQ=255)
                                                                                                                             tATTCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTAc                     >  2:1023768/1‑150 (MQ=255)
                                                                                                                             tATTAGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTAc                     >  2:493914/1‑150 (MQ=255)
                                                                                                                               ttCGAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACaa                   >  1:2516262/1‑150 (MQ=255)
                                                                                                                                  gAACCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCt                <  2:1794551/150‑1 (MQ=255)
                                                                                                                                    aCCTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGt              >  2:935128/1‑150 (MQ=255)
                                                                                                                                      cTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGTTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTg            >  1:948304/1‑150 (MQ=255)
                                                                                                                                      cTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTg            >  1:5019391/1‑150 (MQ=255)
                                                                                                                                      cTGGCAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTg            <  1:1713311/150‑1 (MQ=255)
                                                                                                                                          cAATATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGAttt        <  1:4588168/150‑1 (MQ=255)
                                                                                                                                           aaTATTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTa       >  2:4595903/1‑150 (MQ=255)
                                                                                                                                            atatTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTAc      >  2:1505891/1‑150 (MQ=255)
                                                                                                                                              atTGCTCAGTGGATTGCACGAAATCTGATGAGCGAACATGCGCAGTGGTCAATGGCACAGGCCATAACCCTGCTGGCGGACGTGGAACGGTTATGTCCGCAACTTGTGAAAACGCCGCCGGGTGGTCTGTTACAATCTGTTGATTTACAt    <  2:4401988/150‑1 (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: