Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,969,920 G→T 5.4% G87C (GGT→TGT)  mutH → methyl‑directed mismatch repair protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,969,9200GT 5.4% 256.0 / 2.9 74G87C (GGT→TGT) mutHmethyl‑directed mismatch repair protein
Reads supporting (aligned to +/- strand):  ref base G (37/33);  new base T (2/2);  total (39/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.40e-01

GCTGGTTACGCCAGAGAATTTAAAACGCGATAAAGGCTGGATTGGCGTGTTACTGGAGATCTGGCTAGGTGCCAGCGCAGGGAGTAAACCTGAGCAAGATTTTGCTGCTCTGGGCGTGGAACTTAAAACTATCCCTGTGGATAGTCTTGGTCGTCCGCTGGAAACAACATTCGTTTGTGTTGCCCCGTTAACGGGCAATAGCGGGGTGACCTGGGAAACCAGCCACGTGCGCCACAAACTCAAACGCGTACTGTGGATACCGGTTGAAGGCGAGCGCAGCATCCCGCTGGCGCAGCG  >  NC_000913/2969772‑2970068
                                                                                                                                                    |                                                                                                                                                    
gCTGGTTACGCCAGAGAATTTAAAACGCGATAAAGGCTGGATTGGCGTGTTACTGGAGATCTGGCTAGGTGCCAGCGCAGGGAGTAAACCTGAGCAAGATTTTGCTGCTCTGGGCGTGGAACTTAAAACTATCCCTGTGGATAGTCTTg                                                                                                                                                      <  8:101605/149‑1 (MQ=255)
 cTGGTTACGCCAGAGAATTTAAAACGCGATAAAGGCTGGATTGGCGTGTTACTGGAGATCTGGCTAGGTGCCAGCGCAGGGAGTAAACCTGAGCAAGATTTTGCTGCTCTGGGCGTGGAACTTAAAACTATCCCTGTGGATAGTCTTgg                                                                                                                                                     <  4:23251/149‑1 (MQ=255)
 cTGGTTACGCCAGAGAATTTAAAACGCGATAAAGGCTGGATTGGAGTGTTACTGGAGATCTGGCTAGGTGCCAGCGCAGGGAGTAAACCTTAGCAAGATTTTGCTGCTCTGGGCGTGGAACTTAAAACTATCCCTGTGGATAGTCTTgg                                                                                                                                                     <  5:108721/149‑1 (MQ=255)
    gTTACGCCAGAGAATTTAAAACGCGATAAAGGCTGGATTGGCGTGTTACTGGAGATCTGGCTAGGTGCCAGCGCAGGGAGTAAACCTGAGCAAGATTTTGCTGCTCTGGGCGTGGAACTTAAAACTATCCCTGTGGATAGTCTTGgtcg                                                                                                                                                  >  3:284653/1‑149 (MQ=255)
         gCCAGAGAATTTAAAACGCGATAAAGGCTGGATTGGCGTGTTACTGGAGATCTGGCTAGGTGCCAGCGCAGGGAGTAAACCTGAGCAAGATTTTGCTGCTCTGGGCGTGGAACTTAAAACTATCCCTGTGGATAGTCTTGGTCGTCCGc                                                                                                                                             >  1:294406/1‑149 (MQ=255)
           cAGAGAATTTAAAACGCGATAAAGGCTGGATTGGCGTGTTACTGGAGATCTGGCTAGGTGCCAGCGCAGGGAGTAAACCTGAGCAAGATTTTGCTGCTCTGGGCGTGGAACTTAAAACTATCCCTGTGGATAGTCTTGGTCGTCCGCTg                                                                                                                                           >  1:132895/1‑149 (MQ=255)
           cAGAGAATTTAAAACGCGATAAAGGCTGGATTGGCGTGTTACTGGAGATCTGGCTAGGTGCCAGCGCAGGGAGTAAACCTGAGCAAGATTTTGCTGCTCTGGGCGTGGAACTTAAAACTATCCCTGTGGATAGTCTTGGTCGTCCGCTg                                                                                                                                           >  4:245672/1‑149 (MQ=255)
           cAGAGAATTTAAAACGCGATAAAGGCTGGATTGGCGTGTTACTGGAGATCTGGCTAGGTGCCAGCGCAGGGAGTAAACCTGAGCAAGATTTTGCTACTCTGGGCGTGGAACTTAAAACTATCCCTGTGGATAGTCTTGGTCGTCCGCTg                                                                                                                                           >  3:299709/1‑149 (MQ=255)
                      aaaCGCGATAAAGGCTGGATTGGCGTGTTACTGGAGATCTGGCTAGGTGCCAGCGCAGGGAGTAAACCTGAGCAAGATTTTGCTGCTCTGGGCGTGGAACTTAAAACTATCCCTGTGGATAGTCTTGGTCGTCCGCTGGAAACAACAtt                                                                                                                                >  7:239000/1‑149 (MQ=255)
                       aaCGCGATAAAGGCTGGATTGGCGTGTTACTGGAGATCTGGCTAGGTGCCAGCGCAGGGAGTAAACCTGAGCAAGATTTTGCTGCTCTGGGCGTGGAACTTAAAACTATCCCTGTGGATAGTCTTGGTCGTCCGCTGGAAACAACATTc                                                                                                                               <  7:386708/149‑1 (MQ=255)
                         cgcgATAAAGGCTGGATTGGCGTGTTACTGGAGATCTGGCTAGGTGCCAGCGCAGGGAGTAAACCTGAGCAAGATTTTGCTGCTCTGGGCGTGGAACTTAAAACTATCCCTGTGGATAGTCTTGGTCGTCCGCTGGAAACAACATTCGt                                                                                                                             >  7:270967/1‑149 (MQ=255)
                           cgATAAAGGCTGGATTGGCGTGTTACTGGAGATCTGGCTAGGTGCCAGCGCAGGGAGTAAACCTGAGCAAGATTTTGCTGCTCTGGGCGTGGAACTTAAAACTATCCCTGTGGATAGTCTTGGTCGTCCGCTGGAAACAACATTCGttt                                                                                                                           <  6:290169/149‑1 (MQ=255)
                           cgATAAAGGCTGGATTGGCGTGTTACTGGAGATCTGGCTAGGTGCCAGCGCAGGGAGTAAACCTGAGCAAGATTTTGCTGCTCTGGGCGTGGAACTTAAAACTATCCCTGTGGATAGTCTTGGTCGTCCGCTGGAAACAACATTCGttt                                                                                                                           <  7:229553/149‑1 (MQ=255)
                             aTAAAGGCTGGATTGGCGTGTTACTGGAGATCTGGCTAGGTGCCAGCGCAGGGAGTAAACCTGAGCAAGATTTTGCTGCTCTGGGCGTGGAACTTAAAACTATCCCTGTGGATAGTCTTGGTCGTCCGCTGGAAACAACATTCGTTtgt                                                                                                                         <  2:247812/149‑1 (MQ=255)
                             aTAAAGGCTGGATTGGCGTGTTACTGGAGATCTGGCTAGGTGCCAGCGCAGGGAGTAAACCTGAGCAAGATTTTGCTGCGCTGGGCGTGGAACTTAAAACTATCCCTGTGGATAGTCTTGGTCGTCCGCTGGAAACAACATTCGTTtgt                                                                                                                         <  8:361198/149‑1 (MQ=255)
                               aaaGGCTGGATTGGCGTGTTACTGGAGATCTGGCTAGGTGCCAGCGCAGGGAGTAAACCTGAGCAAGATTTTGCTGCTCTGGGCGTGGAACTTAAAACTATCCCTGTGGATAGTCTTGGTCGTCCGCTGGAAACAACATTCGTTTgtgt                                                                                                                       >  6:117075/1‑149 (MQ=255)
                                aaGGCTGGATTGGCGTGTTACTGGAGATCTGGCTAGGTACCAGCGCAGGGAGTAAACCTGAGCAAGATTTTGCTGCTCTGGGCGTGGAACTTAAAACTATCCCTGTGGATAGTCTTGGTCGTCCGCTGGAAACAACATTCGTTTgtgt                                                                                                                       <  3:3750/148‑1 (MQ=255)
                                  ggCTGGATTGGCGTGTTACTGGAGATCTGGCTAGGTGCCAGCGCAGGGAGTAAACCTGAGCAAGATTTTGCTGCTCTGGGCGTGGAACTTAAAACTATCCCTGTGGATAGTCTTGGTCGTCCGCTGGAAACAACATTCGTTTGTGTTGc                                                                                                                    <  2:132895/149‑1 (MQ=255)
                                    cTGGATTGGCGTGTTACTGGAGATCTGGCTAGGTGCCAGCGCAGGGAGTAAACCAGAGCAAGATTTTGCTGCTCTGGGCGTGGAACTTAAAACTATCCCTGTGGATAGTCTTGGTCGTCCGCTGGAAACAACATTCGTTTGTGTTGccc                                                                                                                  >  4:321478/1‑149 (MQ=255)
                                        aTTGGCGTGTTACTGGAGATCTGGCTAGGTGCCAGCGCAGGGAGTAAACCTGAGCAAGATTTTGCTGCTCTGGGCGTGGAACTTAAAACTATCCCTGTGGATAGTCTTGGTCGTCCGCTGGAAACAACATTCGTTTGTGTTGCCCCGtt                                                                                                              <  3:239983/149‑1 (MQ=255)
                                         ttGGCGTGTTACTGGAGATCTGGCTAGGTGCCAGCGCAGGGAGTAAACCTGAGCAAGATTTTGCTGCTCTGGGCGTGGAACTTAAAACTATCCCTCTGGATAGTCTTGGTCGTCCGCTGGAAACAACATTCGTTTGTGTTGCCCCGtt                                                                                                              >  7:226080/1‑148 (MQ=255)
                                           ggCGTGTTACTGGAGATCTGGCTAGGTGCCAGCGCAGGGAGTAAACCTGAGCAAGATTTTGCTGCTCTGGGCGTGGAACTTAAAACTATCCCTGTGGATAGTCTTGGTCGTCCGCTGGAAACAACATTCGTTTGTGTTGCCCCGTTaa                                                                                                            <  8:239000/148‑1 (MQ=255)
                                            gCGTGTTACTGGAGATCTGGCTAGGTGCCAGCGCAGGGAGTAAACCTGAGCAAGATTTTGCTGCTCTGGGCGTGGAACTTAAAACTATCCCTGTGGATAGTCTTGGTCGTCCGCTGGAAACAACATTCGTTTGTGTTGCCCCGTTAAc                                                                                                           <  7:75789/148‑1 (MQ=255)
                                              gtgtTACTGGAGATCTGGCTAGGTGCCAGCGCAGGGAGTAAACCAGAGCAAGATTTTGCTGCTCTGGGCGTGGAACTTAAAACTATCCCTGTGGATAGTCTTGGTCGTCCGCTGGAAACAACATTCGTTTGTGTTGCCCCGTTAACggg                                                                                                        <  3:321478/149‑1 (MQ=255)
                                                   aCTGGAGATCTGGCTAGGTGCCAGCGCAGGGAGTAAACCTGAGCAAGATTTTGCTGCTCTGGGCGTGGAACTTAAAACTATCCCTGTGGATAGTCTTGGTCGTCCGCTGGAAACAACATTc                                                                                                                               <  3:12521/121‑1 (MQ=255)
                                                   aCTGGAGATCTGGCTAGGTGCCAGCGCAGGGAGTAAACCTGAGCAAGATTTTGCTGCTCTGGGCGTGGAACTTAAAACTATCCCTGTGGATAGTCTTGGTCGTCCGCTGGAAACAACATTc                                                                                                                               >  4:12521/1‑121 (MQ=255)
                                                   aCTGGAGATCTGGCTAGGTGCCAGCGCAGGGAGTAAACCTGAGCAAGATTTTGCTGCTCTGGGCGTGGAACTTAAAACTATCCCTGTGGATAGTCTTGGTCGTCCGCTGGAAACAACAGTCGTTTGTGTTGCCCAGTTAACGGGCaa                                                                                                     >  6:291147/1‑147 (MQ=255)
                                                     tGGAGATCTGGCTAGGTGCCAGCGCAGGGAGTAAACCTGAGCAAGATTTTGCTGCTCTGGGCGTGGAACTTAAAACTATCCCTGTGGATAGTCTTGGTCGTCCGCTGGAAACAACATTCGTTTGTGTTGCCCCGTTAACGGGCAATAGc                                                                                                 >  6:50196/1‑149 (MQ=255)
                                                        tgcgcTGGCTAGGTGCCAGCGCAGGGAGTAAACCTGAGCAAGATTTTGCTGCTCTGGGCGTGGAACTTAAAACTATCCCTGTGGATAGTCTTTGTCGTCCGCTGGAAACAACATTCGTTTGTGTTGCCCCGTTAACGGGCAATAGCggg                                                                                              >  6:215237/5‑149 (MQ=255)
                                                        tgcgcTGGCTAGGTGCCAGCGCAGGGAGTAAACCTGAGCAAGATTTTGCTGCTCTGGGCGTGGAACTTAAAACTATCCCTGTGGATAGTCTTTGTCGTCCGCTGGAAACAACATTCGTTTGTGTTGCCCCGTTAACGGGCAATAGCggg                                                                                              >  6:210004/5‑149 (MQ=255)
                                                        cgctCTGGCTAGGTGCCAGCGCAGGGAGTAACCCTGAGCAAGATTTTGCTGCTCTGGGCGTGGAACTTAAAACTATCCCTGTGGATAGTCTTGGTCGTCCGCTGGAAACAACATTCGTTTGTGGTGCCCCGTTAACGGGCAATAGCggg                                                                                              <  7:192874/146‑1 (MQ=255)
                                                        agaTCTGGCTAGGTGCCAGCGCAGGGAGTAAACCTGAGCAAGATTTTGCTGCTCTGGGCGTGGAACTTAAAACTATCCCTGTGGATAGTCTTGGTCGTCCGCTGGAACCAACATTCGTTTGTGTTGCCCCGTTAACGGGCAATAGCggg                                                                                              <  5:341154/149‑1 (MQ=255)
                                                        agaTCTGGCTAGGGGCCAGCGCAGGGAGTAAACCTGAGCAAGATTTTGCTGCTCTGGGCGTGGAACTTAAAACTATCCCTGTGGATAGTCTTGGTCGTCCGCTGGAAACAACATTCGTTTGTGTTGCCCCGTTAACGGGCAATAGCggg                                                                                              <  5:53859/149‑1 (MQ=255)
                                                               gCTAGGTGCCAGCGCAGGGAGTAAACCTGAGCATGATTTTGCTGCTCTGGGCGTGGAACTTAAAACTATCGCTGTGGATAGTCTTGTTCGTCCGCTGGAAACAACATTCGTTTGTGTTGCCCCGTTAACGGGCAATAGCGGGGt                                                                                            <  4:280581/144‑1 (MQ=255)
                                                               gCTAGGTGCCAGCGCAGGGAGTAAACCTGAGCAAGATTTTGCTGCTCTGGGCGTGGAACTTAAAACTATCCCTGTGGATAGTCTTGGTCGTCCGCTGGAAACAACATTCGTTTGTGTTGCCCCGTTAACGGGCAATAGCGGGGTGACCt                                                                                       <  4:284653/149‑1 (MQ=255)
                                                               gCTAGGTGCCAGCGCAGGGAGTAAACCTGAGCAAGATTTTGCTGCTCTGGGCGTGGAACTTAAAACTATCCCTGTGGATAGTCTTGGTCGTCCGCTGGAAACAACATTCGTTTGTGTTGCCCCGTTAACGGGCAATAGCGGGGTGACCt                                                                                       >  4:84284/1‑149 (MQ=255)
                                                                cgaGGTGCCAGCGCAGGGAGGAAACCGGCGCATGAATTTGCAGCTCTGGGCGTGGCACTTAAAACTATCCCTGTGGATATTCTTGGTCGTCCGCTGGAAACAACATTCGTTTGTGTTGCCCCGTTAACGGGCAATAGCGGGGGGACCt                                                                                       <  4:177977/146‑1 (MQ=255)
                                                                cTAGGTGCCAGCGCAGGGAGTAAACCTGAGCAAGATTTTGCTGCTCTGGGCGTGGAACTTAAAACTATCCCTGTGGATAGTCTTGGTCGTCCGCTGGAAACAACATTCGTTTGTGTTGCCCCGTTAACGGGCAATAGCGGGGTGACCTg                                                                                      <  6:367712/149‑1 (MQ=255)
                                                                cTAGGTGCCAGCGCAGGGAGTAAACCTGAGCAAGATTTTGCTGCTCTGGGCGTGGAACTTAAAACTATCCCTGTGGATAGTCTTGGTCGTCCGCTGGAAACAACATTCGTTTGTGTTGCCCCGTTAACGGGCAATAGCGGGGTGACCTg                                                                                      <  5:117075/149‑1 (MQ=255)
                                                                   ggTGCCAGCGCAGGGAGTAAACCTGAGCAAGATTTTGCTGCTCTGGGCGTGGAACTTAAAACTATCCCTGTGGATAGTCTTGGTCGTCCGCTGGAAACAACATTCGTTTGTGTTGCCCCGTTAACGGGCAATAGCGGGGTGACCTgg                                                                                     <  8:34480/147‑1 (MQ=255)
                                                                     tGCCAGCTCAGGGAGTAAACCTCAGCAAGATTTTGCTTCTCTGGGCGTGGAACTTAAAACTATCCCTGTGGATAGTCTTGGTCGTCCGCTGGAAACAACATTCGTTTGTGTTGCCCCGTTAACGGGCAATAGCGGGGTGACCTGGGaaa                                                                                 <  2:294406/149‑1 (MQ=255)
                                                                       ccAGCGCAGGGAGTAAACCTGAGCAAGATTTTTCTGCTCTGGGCGTGGAACTTAAAACTATCCCTGTGGATAGTCTTGGTCGTCCGCTGGAAACAACATTCGTTTGTGTTGCCCCGTTAACGGGCAATAGCGGGGTGACCTGGGAAAcc                                                                               <  2:269030/149‑1 (MQ=255)
                                                                        ccgcgcAGGGAGAAAACCTGAGCAAGATTTAGCTGCTCTGGGCTTGGAACTTAAAACTATACCTGTGGATAGTCTTGGTCGTCCGCTGTAAACAACATTCGTTTGTGTTGCCCCGTTAACGGTCTCTAGCGGGGTGACCTGGGAAAc                                                                                >  6:265378/3‑147 (MQ=255)
                                                                         aGCGCAGGGAGTAAACCTGAGCAAGATTTTGCTGCTCTGGGCGTGGAACTTAAAACTATCCCTGTGGATAGTCTTGGTCGTCCGCTGGAAACAACATTCGTTTGTGTTGCCCCGTTAACGGGCAATAGCGGGGTGACCTGGGAAACCAg                                                                             >  2:307059/1‑149 (MQ=255)
                                                                         aGCGCAGGGAGTAAACCTGAGCAAGATTTTGCTGCTCTGGGCGTGGAACTTAAAACTATCCCTGTGGATAGTCTTGGTCGTCCGCTGGAAACAACATTCGTTTGTGTTGCCCCGTTAACGGGCAATAGCGGGGTGACCTGGGAAACCAg                                                                             >  6:233112/1‑149 (MQ=255)
                                                                         aGCGCAGGGAGTAAACCTGAGCAAGAATTTGCTGCTCTGGGCGGGGAACTTAAAACTATCCCTGTGGCTAGTCTTGGTCGTCCGCTGGAAACAACATTCGTTTGTGTTGCCCCGTTAACGGGCAATAGCGGGGTGACCTGGGAAACCAg                                                                             <  1:62006/149‑1 (MQ=255)
                                                                             cAGGGAGTAAACCTGAGCAAGATTTTGCTGCTCTGGGCGTGGAACTTAAAACTATCCCTGTGGATAGTCTTGGTCGTCCGCTGGAAACAACATTCGTTTGTGTTGCCCCGTTAACGGGCAATAGCGGGGTGACCTGGGAAACCAGCCAc                                                                         <  3:84284/149‑1 (MQ=255)
                                                                                  aGTAAACCTGAGCAAGATTTTGCTGCTCTGGGCGTGGAACTTAAAACTATCCCTGTGGATAGTCTTGGTCGTCCGCTGGAAACAACATTCGTTTGTGTTGCCCCGTTAACGGGCAATAGCGGGGTGACCTGGGAAACCAGCCACGTgcg                                                                    >  3:346674/1‑149 (MQ=255)
                                                                                    tAAACCTGAGCAAGATTTTGCTGCTCTGGGCGTGGAACTTAAAACTATCCCTGTGGATAGTCTTGGTCGTCCGCTGGAAACAACATTCGTTTGTGTTGCCCCGTTAACGGGCAATAGCGGGGTGACCTGGGAAACCAGCCACGTgcgc                                                                   >  7:17510/1‑148 (MQ=255)
                                                                                              cAAGATTTTGCTGCTCTGGGCGTGGAACTTAAAACTATCCCTGTGGATAGTCTTGGTCGTCCGCTGGAAACAACATTCGTTTGTGTTGCCCCGTTAACGGGCAATAGCGGGGTGACCTGGGAAACCAGCCACGTGCGCCACAAACTCa                                                         >  4:107006/1‑148 (MQ=255)
                                                                                                  aTTTTGCTGCTCTGGGCGTGGAACTTAAAACTATCCCTGTGGATAGTCTTGGTCGTCCGCTGGAAACAACATTCGTTTGTGTTGCCCCGTTAACGGGCAATAGCGGGGTGACCTGGGAAACCAGCCACGTGCGCCCCAAACTCAAAcgc                                                    >  4:232590/1‑149 (MQ=255)
                                                                                                   ttttGCTGCTCTGGGCGTGGAACTTAAAACTATCCCTGTGGATAGTCTTGGTCGTCCGCTGGAAACAACATTCGTTTGTGTTGCCCCGTTAACGGGCAATAGCGGGGTGACCTGGGAAACCAGCCACGTGCGCCCCAAACTCAAAcgcg                                                   >  2:286118/1‑149 (MQ=255)
                                                                                                    tttGCTGCTCTGGGCGTGGAACTTAAAACTATCCCTGTGGATAGTCTTGGTCGTCCGCTGGAAACAACATTCGTTTGTGTTGCCCCGTTAACGGGCAATAGCGGGGTGACCTGGGAAACCAGCCACGTGCGCCACAAACTCAAAcgcg                                                   <  7:145306/148‑1 (MQ=255)
                                                                                                    tttGCTGCTCTGGGCGTGGAACTTAAAACTATCCCTGTGGATAGTCTTGGTCGTCCGCTGGAAACAACATTCGTTTGTGTTGCCCCGTTAACGGGCAATAGCGGGGTGACCTGGGAAACCAGCCACGTGCGCCACAAACTCAAAcgcg                                                   >  8:145306/1‑148 (MQ=255)
                                                                                                        ctgctCTGGGCGTGGAACTTAAAACTATCCCTGTGGATAGTCTTGGTCGTCCGCTGGAAACAACATTCGTTTGTGTTGCCCCGTTAACGGGCAATAGCGGGGTGACCTGGGAAACCAGCCACGTGCGCCACAAACTCAAACGCGTACtg                                              <  6:303700/149‑1 (MQ=255)
                                                                                                         tgctCTGGGCGTGGAACTTAAAACTATCCCTGTGGATAGTCTTTGTCGTCCGCTGGAAACAACATTCGTTTGTGTTGCCCCGTTAACGGGCAATAGCGGGGTGACCTGGGAAACCAGCCACGTGCGCCACAAACTCAAACGCGTACtgt                                             <  5:210004/149‑1 (MQ=255)
                                                                                                          gctCTGGGCGTGGAACTTAAAACTATCCCTGTGGATAGTCTTGGTCGTCCGCTGGAAACAACATTCGTTTGTGTTGCCCCGTTAACGGGCAATAGCGGGGTGACCTGGGACACCAGCCACGTGCGCCACAAACTCAAACGCGTACtgtg                                            >  2:51609/1‑149 (MQ=255)
                                                                                                           ctctGGGCGTGGAACTTAAAACTATCCCTGTGGATAGTCTTTGTCGTCCGCTGGAAACAACATTCGTTTGTGTTGCCCCGTTAACGGGCAATAGCGGGGTGACCTGGGAAACCAGCCACGTGCGCCACAAACTCAAACGCGTACtgt                                             <  5:215237/147‑1 (MQ=255)
                                                                                                           ctctGGGCGTGGAACTTAAAACTATCCCTGTGGATAGTCTTGGTCGTCCGCTGGAAACAACATTCGTTTGTGTTGCCCCGTTAACGGGCAATAGCGGGGTGACCTGGGAAACCAGCCACGTGCGCCACAAACTCAAACGCGTACTGTgg                                           <  1:8187/149‑1 (MQ=255)
                                                                                                           ctctGGGCGTGGAACTTAAAACTATCCCTGTGGATAGTCTTGGTCGTCCGCTGGAAACAACATTCGTTTGTGTTGCCCCGTTAACGGGCAATAGCGGGCTGACCTGGGAAACCAGCCACTTGCGCCACAAACTCAAACTCGTACtgtg                                            >  2:370199/1‑148 (MQ=255)
                                                                                                              tGGGCGTGGAACTTAAATCTATCCCTGTGGATAGTCTTGGTCGTCCGCTGGAAACAACATTCGTTTGTGTTGCCCCGTTAACGGGCAATAGCGGGGTGACCTGGGAAACCAGCCACGTGCGCCACAAACTCAAACGCGTACTGTGGa                                          <  1:149312/147‑1 (MQ=255)
                                                                                                                  cGTGGAACTTAAAACTATCCCTGTGGATAGTCTTGGTCGTCCGCTGGAAACAACATTCGTTTGTGTTGCCCCGTTAACGGGCAATAGCGGGGTGACCTGGGAAACCAGCCACGTGCGCCACAAACTCAAACGCGTACTGTGGATACCgg                                    <  8:17510/149‑1 (MQ=255)
                                                                                                                    tGGAACTTAAAACTATCCCTGTGGATAGTCTTGGTCGTCCGCTGGAAACAACATTCGTTTGTGTTGCCCCGTTAACGGGCAATAGCGGGGTGACCTGGGAAACCAGCCACGTGCGCCACAAACTCAAACGCGTACTGTGGATACCGGtt                                  >  3:166780/1‑149 (MQ=255)
                                                                                                                         cTTAAAACTATCCCTGTGGATAGTCTTGGTCGTCCGCTGGAAACAACATTCGTTTGTGTTGCCCCGTTAACGGGCAATAGCGGGGTGACCTGGGAAACCAGCCACGTGCGCCACAAACTCAAACGCGTACTGTGGa                                          <  6:19525/136‑1 (MQ=255)
                                                                                                                         cTTAAAACTATCCCTGTGGATAGTCTTGGTCGTCCGCTGGAAACAACATTCGTTTGTGTTGCCCCGTTAACGGGCAATAGCGGGGTGACCTGGGAAACCAGCCACGTGCGCCACAAACTCAAACGCGTACTGTGGa                                          >  5:19525/1‑136 (MQ=255)
                                                                                                                           tAAAACTATCCCTGTGGATAGTCTTGGTCGTCCGCTGGAAACAACATTCGTTTGTGTTGCCCCGTTAACGGGCAATAGCGGGGTGACCTGGGAAACCAGCCACGTGCGCCACAAACTCAAACGCGTACTGTGGATACCGGTTGAAGGCg                           >  3:143895/1‑149 (MQ=255)
                                                                                                                            aaaaCTATCCCTGTGGATAGTCTTGGTCGTCCGCTGGAAACAACATTCGTTTGTGTTGCCCCGTTAACGGGCAATAGCGGGGTGACCTGGGAAACCAGCCACGTGCGCCACAAACTCAACCGCGTACTGTGGATACCGGTTGAAGGc                            >  2:341151/1‑147 (MQ=255)
                                                                                                                             aaaCTATCCCTGTGGATAGTCTTGGTCGTCCGCTGGAAACAACATTCGTTTGTGTTGCCCCGTTAACGGGCAATAGCGGGGTGACCTGGGAAACCCGCCACGTGCGCCACAAACTCAAACGCGTACTGTGGATACCGGTTGAAGGCGAg                         >  1:220326/1‑149 (MQ=255)
                                                                                                                             aaaCTATCCCTGTGGATAGTCTTGGTCGTCCGCTGGAAACAACATTCGTTTGTGTTGCCCCGTTAACGGGCAATAGCGGGGTGACCTGGGAAACCAGCCACGTGCGCCACAAACTCAAACGCGTACTGTGGATACCGGTTGAAGGCGAg                         >  6:31186/1‑149 (MQ=255)
                                                                                                                                cTATCCCTGTGGATAGTCTTGGTCGTCCGCTGGAAACAACATTCGTTTGTGTTGCCCCGTTAACGGGCAATAGCGGGGTGACCTGGGAAACCAGCCACGTGCGCCACAAACTCAAACGCGTACTGTGGATACCGGTTGAAGGCGAgcg                       >  1:12680/1‑148 (MQ=255)
                                                                                                                                      cTGTGGATAGTCTTGGTCGTCCGCTGGAAACAACATTCGTTTGTGTTGCCCCGTTAACGGGCAATAGCGGGGTGACCTGGGAAACCAGCCACGTGCGCCACAAACTCAAACGCGTACTGTGGATACCGGTTGAAGGCGAGCGCAGCATc                >  8:99308/1‑149 (MQ=255)
                                                                                                                                       tgtgGATAGTCTTGGTCGTCCGCTGGAAACAACATTCGTTTGTGTTGCCCCGTTAACGGGCAATAGCGGGGTGACCTGGGAAACCAGCCACGTGCGCCACAAACTCAAACGCGTACTGTGGATACCGGTTGAAGGCGAGCGCAGCATcc               >  2:194972/1‑149 (MQ=255)
                                                                                                                                         tgGATAGTCTTGGTTGTCCGCTGGAAACAACATTCGTTTGTGTTGCCCCGTTAACGGGCAATAGCGGGGTGACCTGGGAAACCAGCCACGTGCGCCACAAACTCAAACGCGTACTGTGGATACCGGTTGAAGGCGAGCGCAGCATCCCg             >  3:240185/1‑149 (MQ=255)
                                                                                                                                           gATAGTCTTGGTCGTCCGCTGGAAACAACATTCGTTTGTGTTGCCCCGTTAACGGGCAATAGCGGGGTGACCTGGGAAACCAGCCACGTGCGCCACAAACTCAAACGCGTACTGTGGATACCGGTTGAAGGCGAGCGCAGCATCCCGCt           >  5:176166/1‑149 (MQ=255)
                                                                                                                                            aTAGTCTTGGTCGTCCGCTGGAAACAACATTCGTTTGTGTTGCCCCGTTAACGGGCAATAGCGGGGTGACCTGGGAAACCAGCCACGTGCGCCACAAACTCAAACGCGTACTGTGGATACCGGTTGAAGGCGAGCGCAGCATCCCGCTg          >  2:270267/1‑149 (MQ=255)
                                                                                                                                             tAGTCTTGGTCGTCCGCTGGAAACAACATTCGTTTGTGTTGCCCCGTTAACGGGCAATAGCGGGGTGACCTGGGAAACCAGCCACGTGCGCCACAAACTCAAACGCGTACTGTGGATACCGGTTGAAGGCGAGCGCAGCATCCCGCTgg         <  6:176166/149‑1 (MQ=255)
                                                                                                                                                tCTTGGTCGTCCGCTGGAAACAACATTCGTTTGTGTTGCCCCGTTAACGGGCAATAGCGGGGTGACCTGGGAAACCAGCCACGTGCGCCACAAACTCAAACGCGt                                                  <  6:47678/105‑1 (MQ=255)
                                                                                                                                                    ggTCGTCCGCTGGAAACAACATTCGTTTGTGTTGCCCCGTTAACGGGCAATAGCGGGGTGACCTGGGAAACCAGCCACGTGCGCAACAAACTCAAACGCGTACTGTGGATACCGGTTGAAGGCGAGCGCAGCCTCCCGCTGGCGCcgcg  >  8:288769/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GCTGGTTACGCCAGAGAATTTAAAACGCGATAAAGGCTGGATTGGCGTGTTACTGGAGATCTGGCTAGGTGCCAGCGCAGGGAGTAAACCTGAGCAAGATTTTGCTGCTCTGGGCGTGGAACTTAAAACTATCCCTGTGGATAGTCTTGGTCGTCCGCTGGAAACAACATTCGTTTGTGTTGCCCCGTTAACGGGCAATAGCGGGGTGACCTGGGAAACCAGCCACGTGCGCCACAAACTCAAACGCGTACTGTGGATACCGGTTGAAGGCGAGCGCAGCATCCCGCTGGCGCAGCG  >  NC_000913/2969772‑2970068

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