Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,565,805 A→T 7.1% 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 7.1% 206.4 / 7.2 71C578S (TGC→AGC) mutLDNA mismatch repair protein MutL
Reads supporting (aligned to +/- strand):  ref base A (36/29);  new base T (2/3);  total (39/32)
Fisher's exact test for biased strand distribution p-value = 6.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.53e-01

AGTTTGTCCAGATCGTCCAGCCCCATCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAGGGTCTCGG  >  NC_002947/5565657‑5565954
                                                                                                                                                    |                                                                                                                                                      
aGTTTGTCCAGATCGTCCAGCCCCATCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCaa                                                                                                                                                       <  2:606116/150‑1 (MQ=255)
 gTTTGTCCAGATCGTCCAGCCCCATCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAg                                                                                                                                                      <  1:694926/150‑1 (MQ=255)
  tttGTCCAGATCGTCCAGCCCCATCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGc                                                                                                                                                     >  1:541300/1‑150 (MQ=255)
  tttGTCCAGATCGTCCAGCCCCATCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGc                                                                                                                                                     <  2:1843868/150‑1 (MQ=255)
  attGTCCAGATCGTCCAGCCCCATCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGc                                                                                                                                                     >  2:688346/2‑150 (MQ=255)
            tCGTCCAGCCCCATCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCg                                                                                                                                           <  2:844267/150‑1 (MQ=255)
              gTCCATCCCCATCTGGGTCCAGGTAGGACGACCATGTTTTCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGag                                                                                                                                         >  2:407628/1‑150 (MQ=255)
               tCCAGCCCCATCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGagc                                                                                                                                        >  1:1230543/1‑150 (MQ=255)
                 cAGCCCCATCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCTCGGCACCATGGCTTGCCACGGTGCCgtgcgc                                                                                                                                      <  1:1483956/150‑6 (MQ=255)
                  aGCCCCATCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCGGGTCGGTCATGTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGagcagc                                                                                                                                     >  1:723628/1‑150 (MQ=255)
                  aGCCCCATCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGagcagc                                                                                                                                     >  1:405758/1‑150 (MQ=255)
                   gCCCCATCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGACACCGTGGCAAGCCATGGTGCCGag                                                                                                                                         <  2:204875/145‑1 (MQ=255)
                   gCCCCATCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGACACCGTGGCAAGCCATGGTGCCGag                                                                                                                                         >  1:204875/1‑145 (MQ=255)
                      ccATCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCATGGCTTGCCACGGTGCCGtgcgcgccaa                                                                                                                                 >  2:1808211/1‑143 (MQ=38)
                                  aGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATg                                                                                                                                                 <  1:1651684/122‑1 (MQ=255)
                                  aGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATg                                                                                                                                                 >  2:1651684/1‑122 (MQ=255)
                                    gTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCGGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGa                                                                                                                   >  1:1146490/1‑150 (MQ=255)
                                          cgacCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCGGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGa                                                                                                             <  2:1146490/150‑1 (MQ=255)
                                             ccATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGagca                                                                                                                                       <  2:368322/121‑1 (MQ=255)
                                             ccATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGagca                                                                                                                                       >  1:368322/1‑121 (MQ=255)
                                                       aCTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGt                                                                                                >  1:659841/1‑150 (MQ=255)
                                                       aCTGGCCGCTGCGCTCGGTGTTTTCCATCTCGCGCAGCCCGGCGCTCCCTCCGTGAATTGCCAGACGCCGTGGTGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGt                                                                                                >  1:1007821/1‑150 (MQ=255)
                                                        cTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTa                                                                                               <  2:1230543/150‑1 (MQ=255)
                                                        cTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTa                                                                                               >  1:446949/1‑150 (MQ=255)
                                                          ggCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGGCACCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGCTGCCTTACt                                                                                             <  2:1878863/150‑1 (MQ=255)
                                                                tGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATaa                                                                                       >  1:554435/1‑150 (MQ=255)
                                                                tGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATaa                                                                                       <  2:541300/150‑1 (MQ=255)
                                                                tGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTAGCGCGCACGGCACCGTGGCAAGCCATGCTGCCGAGCAGCTCATTGCGGTACGCCGGAATGCGATCACTGGTGCCGTACTCCATaa                                                                                       >  2:305881/1‑150 (MQ=255)
                                                                  cgcTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAgg                                                                                     >  2:1464097/1‑150 (MQ=255)
                                                                    cTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTc                                                                                   <  2:554435/150‑1 (MQ=255)
                                                                    cTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTc                                                                                   <  1:627225/150‑1 (MQ=255)
                                                                             ttCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAg                                                                                      <  1:536909/138‑1 (MQ=255)
                                                                             ttCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAg                                                                                      >  2:536909/1‑138 (MQ=255)
                                                                                cATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTc                                                                       >  2:757708/1‑150 (MQ=255)
                                                                                  tGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGt                                                                                                     >  2:575873/1‑118 (MQ=255)
                                                                                  tGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGt                                                                                                     <  1:575873/118‑1 (MQ=255)
                                                                                              gggCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTACGCCTGAATGCGATCACTAGTGACCTCCTCAATAAGGTCGCCAAGCACGTGCTGCACCCGGCAGt                                                        >  2:608167/1‑150 (MQ=255)
                                                                                                gCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTg                                                      <  1:407628/150‑1 (MQ=255)
                                                                                                 cGTTCATCTCGGTAATTGCCAGACGCCGTTTGGCTCCCACTGCACCCTAGCAAGCCATGGTTCCGAGCAGCTCATTGATTTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTgg                                                     <  1:546849/150‑1 (MQ=255)
                                                                                                  gTTCATCTCGGGAATTGCCAGACGCCGGTTGGCTCGCACGGCACCGTGGCCAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTgtc                                                    >  1:750473/1‑148 (MQ=255)
                                                                                                  gTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTcc                                                                      >  2:1335856/1‑133 (MQ=255)
                                                                                                  gTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTcc                                                                      <  1:1335856/133‑1 (MQ=255)
                                                                                                   ttCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCACGGTGc                                                                                                                                             <  1:1460784/61‑1 (MQ=255)
                                                                                                   ttCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCACGGTGc                                                                                                                                             >  2:1460784/1‑61 (MQ=255)
                                                                                                    tCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATaa                                                                                       >  1:726259/1‑114 (MQ=255)
                                                                                                    tCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATaa                                                                                       <  2:726259/114‑1 (MQ=255)
                                                                                                      aTCTCGGGAATTGCCAGACGCCGGTTGGTGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGTAGCTCATTGAGGTGCGTCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCa                                                <  2:750473/150‑1 (MQ=255)
                                                                                                        ctcGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCAc                                                                                                         >  1:1756352/1‑92 (MQ=255)
                                                                                                        ctcGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCAc                                                                                                         <  2:1756352/92‑1 (MQ=255)
                                                                                                         tcGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAATCACGTCCTGCACCAGGCGGTTGGCTTCAGcc                                             >  2:1385921/1‑150 (MQ=255)
                                                                                                         tcGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGcc                                             >  2:1692431/1‑150 (MQ=255)
                                                                                                          cGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCt                                            >  1:625612/1‑150 (MQ=255)
                                                                                                           gggAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTg                                           >  2:921585/1‑150 (MQ=255)
                                                                                                            ggAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAg                                                                                                                             >  1:1071492/1‑68 (MQ=255)
                                                                                                            ggAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAg                                                                                                                             <  2:1071492/68‑1 (MQ=255)
                                                                                                            ggAATTGCCAGACGCCGGTTGGCGCGCACGGAACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGc                                          >  1:1826193/1‑150 (MQ=255)
                                                                                                              aaTTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCtt                                        <  1:765551/150‑1 (MQ=255)
                                                                                                                                 gcgcgcACGGCACCGTGGTAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTTGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTgg                     <  1:815292/150‑1 (MQ=255)
                                                                                                                                 gcgcgcACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATaa                                                                                       <  2:247310/85‑1 (MQ=255)
                                                                                                                                 gcgcgcACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATaa                                                                                       >  1:247310/1‑85 (MQ=255)
                                                                                                                                  gctgctcGGCACCATGGCTTGCCACGGTGCCGTGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGc                    <  1:1158033/144‑1 (MQ=255)
                                                                                                                                  cgcgcACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGtt                                                       >  2:592416/1‑115 (MQ=255)
                                                                                                                                  cgcgcACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGtt                                                       <  1:592416/115‑1 (MQ=255)
                                                                                                                                    tgctcGGCACCATGGCTTGCCACGGTGCCGTGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGa                  >  1:1226082/5‑150 (MQ=255)
                                                                                                                                         ggCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTc                                                                       <  1:1077347/93‑1 (MQ=255)
                                                                                                                                         ggCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTc                                                                       >  2:1077347/1‑93 (MQ=255)
                                                                                                                                          gCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGATGGGATCTGGCGAATCGcc            >  1:591746/1‑150 (MQ=255)
                                                                                                                                           cACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTCCTCCCTAACAACGGCAAGCACGCCTGATACTTGGCGGTTTGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCa           <  1:1692431/150‑1 (MQ=255)
                                                                                                                                           cACCATGGCTTGCCACGGTGCCGTGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGc                                              <  1:679847/115‑1 (MQ=255)
                                                                                                                                            aCCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAg          >  2:1029657/1‑150 (MQ=255)
                                                                                                                                            aCCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAg          <  2:405758/150‑1 (MQ=255)
                                                                                                                                                  gCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAGGGtctc    >  2:1910133/1‑150 (MQ=255)
                                                                                                                                                    aaGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAGGGTCTCgg  >  2:422734/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
AGTTTGTCCAGATCGTCCAGCCCCATCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAGGGTCTCGG  >  NC_002947/5565657‑5565954

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