Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,565,806 A→T 7.0% A577A (GCT→GCA 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,8060AT 7.0% 208.7 / 6.5 71A577A (GCT→GCAmutLDNA mismatch repair protein MutL
Reads supporting (aligned to +/- strand):  ref base A (37/29);  new base T (2/3);  total (39/32)
Fisher's exact test for biased strand distribution p-value = 6.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.72e-01

AGTTTGTCCAGATCGTCCAGCCCCATCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAGGGTCTCGGG  >  NC_002947/5565657‑5565955
                                                                                                                                                     |                                                                                                                                                      
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                                                                                                      <  1:575873/118‑1 (MQ=255)
                                                                                  tGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGt                                                                                                      >  2:575873/1‑118 (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                                                                       <  1:1335856/133‑1 (MQ=255)
                                                                                                  gTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTcc                                                                       >  2:1335856/1‑133 (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                                                        <  1:592416/115‑1 (MQ=255)
                                                                                                                                  cgcgcACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGtt                                                        >  2:592416/1‑115 (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)
                                                                                                                                                     aGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAGGGTCTCggg  <  2:1007821/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
AGTTTGTCCAGATCGTCCAGCCCCATCTGGGTCCAGGTAGGACGACCATGGTTGCACTGGCCGCTGCGCTCGGTGTTTTCCATGTCGCGCAGCAGGGCGTTCATCTCGGGAATTGCCAGACGCCGGTTGGCGCGCACGGCACCGTGGCAAGCCATGGTGCCGAGCAGCTCATTGAGGTGCGCCTGAATGCGATCACTGGTGCCGTACTCCATAAGGTCGGCAAGCACGTCCTGCACCAGGCGGTTGGCTTCAGCCTGCTTGAGCAAGGCCGGGATCTGGCGAATCGCCAGGGTCTCGGG  >  NC_002947/5565657‑5565955

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: