Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,565,805 A→T 5.5% C578S (TGC→AGC)  mutL ← DNA mismatch repair protein MutL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,565,8050AT 5.5% 288.0 / 5.3 90C578S (TGC→AGC) mutLDNA mismatch repair protein MutL
Reads supporting (aligned to +/- strand):  ref base A (42/43);  new base T (3/2);  total (45/45)
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.20e-01

GTCCAGATCGTCCAGCCCCATCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAGGGTCTCGG  >  NC_002947/5565662‑5565954
                                                                                                                                               |                                                                                                                                                       
gTCCAGATCGTCCAGCCCCATCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCAt                                                                                                                                                   <  1:282583/150‑1 (MQ=255)
 tctaGATCGTCCAGCCCCATCTGGGTCCAGGTAGGACGACCATCGTTGCACTGGCCGCTGCCCCCGGTGTTTTCCACGTCGCGCAGCAGAGCGTTCATCTCGGGAATTGCCAGACTACCGTTCGCGCGCACGGCACCGTGGCAAGCCATg                                                                                                                                                  <  2:581673/147‑1 (MQ=255)
 tCCAGATCGTCCAGCCCCATCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGGCGCGCAGCCGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATg                                                                                                                                                  <  2:393103/150‑1 (MQ=255)
     gATCGTCCAGCCCCATCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGc                                                                                                                                              >  1:1225369/1‑150 (MQ=255)
        cGTCCAGCCCCATCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCATGGCTTGCCACGGTGCCgt                                                                                                                                           <  1:1896652/150‑2 (MQ=255)
             aGCCCCATCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGagcagc                                                                                                                                      >  2:1568707/1‑150 (MQ=255)
              gCCCCATCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGagcgc                                                                                                                                       <  2:2075199/148‑3 (MQ=255)
              gCCCCATCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGagcgc                                                                                                                                       >  1:2075199/1‑146 (MQ=255)
               ccccATCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTc                                                                                                                                    >  1:1553575/1‑150 (MQ=255)
                 ccATCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCAt                                                                                                                                  >  2:1060271/1‑150 (MQ=255)
                 ccATCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCAt                                                                                                                                  >  2:253180/1‑150 (MQ=255)
                 ccATCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCAt                                                                                                                                  >  2:1510049/1‑150 (MQ=255)
                 ccATCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCAt                                                                                                                                  >  2:307665/1‑150 (MQ=255)
                  cATCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCt                                                                                                                                     >  2:648504/1‑146 (MQ=255)
                  cATCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCt                                                                                                                                     <  1:648504/146‑1 (MQ=255)
                  cATCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCAtt                                                                                                                                 <  1:1287782/150‑1 (MQ=255)
                   aTCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTg                                                                                                                                >  2:2066594/1‑150 (MQ=255)
                    tCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGa                                                                                                                               >  1:259552/1‑150 (MQ=255)
                     cTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCNCGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAg                                                                                                                              <  1:2051513/150‑1 (MQ=255)
                     cTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAg                                                                                                                              <  1:2066594/150‑1 (MQ=255)
                         gTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTAGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTgc                                                                                                                          <  2:489562/150‑1 (MQ=255)
                          tCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGagca                                                                                                                                        <  2:524493/135‑1 (MQ=255)
                          tCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGagca                                                                                                                                        >  1:524493/1‑135 (MQ=255)
                           ccAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTgcgc                                                                                                                        <  2:259552/150‑1 (MQ=255)
                           ccAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCACCTCAATGAGCTgctc                                                                                                                        <  1:894156/150‑3 (MQ=255)
                            cAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGcc                                                                                                                       <  2:691987/150‑1 (MQ=255)
                             aGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGCTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCt                                                                                                                      <  1:709307/150‑1 (MQ=255)
                                   gacgacCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGagc                                                                                                                                         <  2:297504/125‑1 (MQ=255)
                                   gacgacCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGagc                                                                                                                                         >  1:297504/1‑125 (MQ=255)
                                   gacgacCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGCTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGc                                                                                                               >  1:1769739/1‑150 (MQ=255)
                                      gacCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGAt                                                                                                            >  1:8932/1‑150 (MQ=255)
                                       acCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCt                                                                                                                      >  2:1893377/1‑140 (MQ=255)
                                       acCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCt                                                                                                                      <  1:1893377/140‑1 (MQ=255)
                                          aTGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACt                                                                                                        >  1:108290/1‑150 (MQ=255)
                                              ttGCACTGGCCGCTTCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGCAATTGCCCGACGCCGGTTGGCACGCACGCCACCGTGGCAAGCCTTGGGGCAGAGCAGCTAATTGAAGTGCGACTTAATGCGATCCCGGGtc                                                                                                    >  2:10768/1‑149 (MQ=255)
                                                gCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGcc                                                                                                  >  1:1002303/1‑150 (MQ=255)
                                                 cACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCg                                                                                                 <  1:1060271/150‑1 (MQ=255)
                                                  aCTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGt                                                                                                >  2:2011658/1‑150 (MQ=255)
                                                    tGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTAc                                                                                              <  1:307665/150‑1 (MQ=255)
                                                           tGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATac                                                                                       <  2:1553575/150‑2 (MQ=255)
                                                            gcgcTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAg                                                                                      <  1:1510049/150‑1 (MQ=255)
                                                              gcTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGt                                                                                    <  1:175184/150‑1 (MQ=255)
                                                               ttCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTc                                                                                   <  2:510298/149‑1 (MQ=255)
                                                                tCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCg                                                                                  >  1:1431529/1‑150 (MQ=255)
                                                                            aTGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTcc                                                                      >  2:2343221/1‑150 (MQ=255)
                                                                             tGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCt                                                                     >  1:2092419/1‑150 (MQ=255)
                                                                                ctcgcAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGCAGTGGCCCCTCTGTCCATGACTTGGGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCa                                                                  <  1:1253528/148‑1 (MQ=255)
                                                                                 gcgcAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCAc                                                                 >  2:1979646/1‑150 (MQ=255)
                                                                                 gcgcAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCAc                                                                 >  2:250646/1‑150 (MQ=255)
                                                                                  cgcAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGt                                                                                                     >  2:1657706/1‑113 (MQ=255)
                                                                                  cgcAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGt                                                                                                     <  1:1657706/113‑1 (MQ=255)
                                                                                   gcagcaGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCa                                                               >  2:222430/1‑150 (MQ=255)
                                                                                   gcagcaGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGCTCGGCAAGCACGTCCTGCAcag                                                               >  2:1253528/1‑148 (MQ=255)
                                                                                          ggCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATaa                                                                                       <  2:2353444/119‑1 (MQ=255)
                                                                                          ggCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATaa                                                                                       >  1:2353444/1‑119 (MQ=255)
                                                                                          ggCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGtt                                                        <  2:1225369/150‑1 (MQ=255)
                                                                                          ggCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGtt                                                        >  2:1544758/1‑150 (MQ=255)
                                                                                             gTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGc                                                     >  2:1029622/1‑150 (MQ=255)
                                                                                             gTTCATCTCGGGAATTGCCAGACGCCGGCTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGc                                                     <  2:1769739/150‑1 (MQ=255)
                                                                                                  tctcGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAg                                                >  1:1996698/1‑150 (MQ=255)
                                                                                                  tctcGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAg                                                >  1:1419078/1‑150 (MQ=255)
                                                                                                    tcGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGcc                                              >  1:402739/1‑150 (MQ=255)
                                                                                                    tcGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGcc                                              >  2:1827674/1‑150 (MQ=255)
                                                                                                         aaTTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCtt                                         <  1:1121772/150‑1 (MQ=255)
                                                                                                         aaTTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCtt                                         <  1:1827674/150‑1 (MQ=255)
                                                                                                           ttGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGa                                       <  1:253180/150‑1 (MQ=255)
                                                                                                            tGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCAACAGGCGGTTGGCTTCAGCCTGCTTGAg                                      <  1:2011658/150‑1 (MQ=255)
                                                                                                              ccAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGcc                                                                                                  >  2:784058/1‑88 (MQ=255)
                                                                                                              ccAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGcc                                                                                                  <  1:784058/88‑1 (MQ=255)
                                                                                                              ccAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCa                                    <  1:1183696/150‑1 (MQ=255)
                                                                                                                  aCGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTc                                                                                            <  2:567606/90‑1 (MQ=255)
                                                                                                                  aCGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTc                                                                                            >  1:567606/1‑90 (MQ=255)
                                                                                                                    gCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCg                              >  2:680626/1‑150 (MQ=255)
                                                                                                                          tGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCt                       <  1:1544758/150‑1 (MQ=255)
                                                                                                                             cgcgcACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGTGATCTGGc                    <  2:1002303/150‑1 (MQ=255)
                                                                                                                             cgcgcACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGc                    <  1:680626/150‑1 (MQ=255)
                                                                                                                             cgactaCGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGc                    >  2:2401720/6‑150 (MQ=255)
                                                                                                                               tgctcGGCACCATGGCTTGCCACGGTGCCGTGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCtt                                         >  1:679009/5‑128 (MQ=38)
                                                                                                                               tgctcGGCACCATGGCTTGCCACGGTGCCGTGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCtt                                         <  2:679009/124‑1 (MQ=38)
                                                                                                                               tgctcGGCACCATGGCTTGCCACGGTGCCGTGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGa                  >  1:1559014/5‑150 (MQ=255)
                                                                                                                               tgctcGGCACCATGGCTTGCCACGGTGCCGTGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGa                  >  1:1391532/5‑150 (MQ=255)
                                                                                                                                  ctggCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCATACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATc               <  1:2330649/148‑1 (MQ=255)
                                                                                                                                  aCGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATc               <  2:1505474/150‑1 (MQ=255)
                                                                                                                                    ggCACCGTGGCAAGCCATGGTGCCGAGCAGCTGATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGc             >  1:1538855/1‑150 (MQ=255)
                                                                                                                                      cACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCa           <  2:8932/150‑1 (MQ=255)
                                                                                                                                        ccGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAgg         >  2:247432/1‑150 (MQ=255)
                                                                                                                                              cAAGCCATGGTGCCGCGCAGCTACTAGATGGTCGCCTGAATGCGATCACTGGTGCCGTACTCCATAGGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAGGGTCTCg   <  1:1209667/150‑1 (MQ=255)
                                                                                                                                              cAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGAGACCCTGGCGATTCGCCAgatcccg    <  1:1383249/150‑7 (MQ=255)
                                                                                                                                              cAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCCGCAcc                                                                <  2:1412947/90‑1 (MQ=255)
                                                                                                                                              cAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCCGCAcc                                                                >  1:1412947/1‑90 (MQ=255)
                                                                                                                                               aaGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAGGGTCTCgg  >  1:1342471/1‑150 (MQ=255)
                                                                                                                                               aaGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAGGGTCTCgg  <  1:1029622/150‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                       
GTCCAGATCGTCCAGCCCCATCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAGGGTCTCGG  >  NC_002947/5565662‑5565954

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