Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,810,947 A→G 5.6% I69T (ATC→ACC)  mutM ← formamidopyrimidine/5‑formyluracil/ 5‑hydroxymethyluracil DNA glycosylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,810,9470AG 5.6% 223.0 / 2.3 71I69T (ATC→ACC) mutMformamidopyrimidine/5‑formyluracil/ 5‑hydroxymethyluracil DNA glycosylase
Reads supporting (aligned to +/- strand):  ref base A (32/35);  new base G (2/2);  total (34/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGTCCACAGCCAGGCACCAAAGCGGCGCGGATCGGTGTAGCGCAGCACTTTGCCGTTGCTCATCACCAAATCCACATGGTCATGCTTTTCAGGGGGAAGTTCTTCTGGAAGGATGCGCAGGCTGCCAGACATCCCTAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCGCTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGCACCACTGCATGAAGAATGGTTGCAC  >  NC_000913/3810801‑3811094
                                                                                                                                                  |                                                                                                                                                   
tGGTCCACAGCCAGGCACCAAAGCGGCGCGGATCGGTGTAGCGCAGCACTTTGCCGTTGCTCATCACCAAATCCACATGGTCATGCTTTTCAGGGGGAAGTTCTTCTGGAAGGATGCGCAGGCTGCCAGACATCCCTAAATGAATGAt                                                                                                                                                    <  7:187080/148‑1 (MQ=255)
   tCCACAGCCAGGCTCCAAAGCGGCGCGGATCGGTGTAGCGCAGCACTTTGCCGTTGCTCATCCCCAAATGCACATGGGCATGCTTTTCAGGGGGAAGTTCTTCTGGAAGGATGCGCAGGCTGCCAGACATCCCTAAATGAATGATAATc                                                                                                                                                <  4:342387/149‑1 (MQ=255)
   tCCACAGCCAGGCACCAAAGCGGCGCGGATCGGTGTAGCGCAGCACTTTGCCGTTGCTCATCACCAAATCCACATGGTCATGCTTTTCAGGGGGAAGTTCTTCTGGAAGGATGCGCAGGCTGCCAGACATCCCTAAATGAATGATAATc                                                                                                                                                <  1:141193/149‑1 (MQ=255)
             ggCACCAAAGCGGCGCGGATCGGTGTAGCGCAGCACTTTGCCGTTGCTCATCACCAAATCCACATGGTCATGCTTTTCAGGGGGAAGTTCTTCTGGAAGGATGCGCAGGCTGCCAGACATCCCTAAATGAATGATAATCCAGCCCTCAg                                                                                                                                      <  2:293611/149‑1 (MQ=255)
                aCCAAAGCGGCGCGGATCGGTGTAGCGCAGCACTTTGCCGTTGCTCATCACCAAATCCACATGGTCATGCTTTTCAGGGGGAAGTTCTTCTGGAAGGATGCGCAGGCTGCCAGACATCCCTAAATGAATGGTAATCCAGCCCTCAgg                                                                                                                                     >  6:106762/1‑147 (MQ=255)
                aCCAAAGCGGCGCGGATCGGTGTAGCGCAGCACTTTGCCGTTGCTCATAACCAAATCCACATGGTCATGCTTTTCAGGGGGAAGTTCTTCTGGAAGGATGCGCAGGCTGCCAGACATCCCTAAATGAATGGTAATCCAGCCCTCAgg                                                                                                                                     >  6:113027/1‑147 (MQ=255)
                       cggcgCGGATCGGTGTAGCGCAGCACTTTGCCGTTGCTCATCACCAAATCCACATGGTCATGCTTTTCAGGGGGAAGTTCTTCTGGAAGGATGCGCAGGCTGCCAGACATCCCTAAATGAATGATAATCCAGCCCTCAGGCAGCTCcag                                                                                                                            <  2:403395/149‑1 (MQ=255)
                         gcgcGGATCGGTGTAGCGCAGCACTTTGCCGTTGCTCATCACCAAATCCACATGGTCATGCTTTTCAGGGGGAAGTTCTTCTGGAAGGATGCGCAGGCTGCCAGACATCCCTAAATGAATGATAATCCAGCCCTCAGGCAGCTCcagca                                                                                                                          <  6:202578/149‑1 (MQ=255)
                           gcgGATCGGTGTAGCGCAGCACTTTGCCGTTGCTCATCACCAAATCCACATGGTCATGCTTTTCAGGGGGAAGTTCTTCTGGAAGGATGCGCAGGCTGCCAGACATCCCTAAATGAATGATAATCCAGCCCTCAGGCAGCTCcagcag                                                                                                                         >  1:13158/1‑148 (MQ=255)
                                  ggTGTAGCGCAGCACTTTGCCGTTGCTCATCACCAAATCCACATGGTCATGCTTTTCAGGGGGAAGTTCTTCTGGAAGGATGCGCAGGCTGCCAGACATCCCTAAATGAATGATAATCCAGCCCTCAGGCAGTTCCAGCAGCAGatat                                                                                                                  <  4:88753/148‑1 (MQ=255)
                                   gtgtAGCGCAGCACTTTGCCGTTGCTCATCACCAAATCCACATGGTCATGCTTTTCAGGGGGAAGTTCTTCTGGAAGGATGAGCAGGCTGCCAGACATCCCTAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATAtt                                                                                                                 >  5:296058/1‑148 (MQ=255)
                                      tAGCGCAGCACTTTGCCGTTGCTCATCACCAAATCCACATGGTCATGCTTTTCAGGGGGAAGTTCTTCTGGAAGGATGCGCAGGCTGCCAGACATCCCTAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGc                                                                                                             <  2:378578/149‑1 (MQ=255)
                                          gcagcaCTTTGCCGTTGCTCATCACCAAATCCACATGGTCATGCTTTTCAGGGGGAAGTTCTTCTGGAAGGATGCGCAGGCTGCCAGACATCCCTAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCccgc                                                                                                         >  8:241654/1‑149 (MQ=255)
                                                cTTTGCCGTTGCTCATCACCAAATCCACATGGTCATGCTTTTCAGGGGGAAGTTCTTCTGGAAGGATGCGCAGGCTGCCAGACATCCCTAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTg                                                                                                    >  2:347382/1‑148 (MQ=255)
                                                   tGCCGTTGCTCATCACCAAATCCACATGGTCATGCTTTTCAGGGGGAAGTTCTTCTGGAAGGATGCGCAGGCTGCCAGACATCCCTAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCa                                                                                                  <  1:347382/147‑1 (MQ=255)
                                                            tcatcaCCAAATCCACATGGTCATGCTTTTCAGGGGGAAGTTCTTCTGGAAGGATGCGCAGGCTGCCAGACATCCCTAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGa                                                                                                                     >  7:389520/1‑119 (MQ=255)
                                                            tcatcaCCAAATCCACATGGTCATGCTTTTCAGGGGGAAGTTCTTCTGGAAGGATGCGCAGGCTGCCAGACATCCCTAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGa                                                                                                                     <  8:389520/119‑1 (MQ=255)
                                                              atcaCCAAATCCACATGGTCATGCTTTTCAGGGGGAAGTTCTTCTGGAAGGATGCGCAGGCTTCCAGACATCCCTAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACtg                                                                                      <  6:421778/148‑1 (MQ=255)
                                                              atcaCCAAATCCACATGGTCATGCTTTTCAGGGGGAAGTTCTTCTGGAAGGATGCGCAGGCTGCCAGACATCCCTAAATGAATGATAATCCAGCCCTCAg                                                                                                                                      >  2:202040/1‑100 (MQ=255)
                                                              atcaCCAAATCCACATGGTCATGCTTTTCAGGGGGAAGTTCTTCTGGAAGGATGCGCAGGCTGCCAGACATCCCTAAATGAATGATAATCCAGCCCTCAg                                                                                                                                      <  1:202040/100‑1 (MQ=255)
                                                              atcaCCAAATCCACATGGTCATGCTTTTCAGGGGGAAGTTCTTCTGGAAGGATGCGCAGGCTGCCAGACATCCCTAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACtgg                                                                                     >  1:435101/1‑149 (MQ=255)
                                                              atcaCCAAATCCACATGGTCATGCTTTTCAGGGGGAAGTTCTTCTGGAAGGATCCGCAGGCTGCCAGACATCCCTAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACtgg                                                                                     <  7:456000/149‑1 (MQ=255)
                                                               tcaCCAAATCCACATGGTCATGCTTTTCAGGGGGAAGTTCTTCTGGAAGGATGCGCAGGCTGCCAGACATCCCTAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACtggt                                                                                    >  5:366528/1‑149 (MQ=255)
                                                                 aCCAAATCCACATGGTCATGCTTTTCAGGGGGAAGTTCTTCTGGAAGGATGCGCAGGCTGCCAGACATCCCTAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTa                                                                                                               >  7:341653/1‑120 (MQ=255)
                                                                 aCCAAATCCACATGGTCATGCTTTTCAGGGGGAAGTTCTTCTGGAAGGATGCGCAGGCTGCCAGACATCCCTAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTa                                                                                                               <  8:341653/120‑1 (MQ=255)
                                                                 aCCAAATCCACATGGTCATGCTTTTCAGGGGGAAGTTCTTCTGGAAGGATGCGCAGGCTGCCAGACATCCCTAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACtggttg                                                                                  <  6:358283/149‑1 (MQ=255)
                                                                 aCCAAATCCACATGGTCATGCTTTTCAGGGGGAAGTTCTTCTGGAAGGATGCGCAGGCTGCCAGACATCCCTAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACtggttg                                                                                  >  5:358283/1‑149 (MQ=255)
                                                                  ccAGATCCACATTGTCATGCTTTTCAGGGGGAAGTTCTTCTGGAAGGATGCGCAGGCTGCCAGACATCCCTAAATGAATGATAATCCAGc                                                                                                                                            <  4:46241/90‑1 (MQ=255)
                                                                  ccAAATCCACATGGTCATGCTTTTCAGGGGGAAGTTCTTCTGGAAGGATGCGCAGGCTGCCAGACATCCCTAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACtggttgg                                                                                 >  2:149881/1‑149 (MQ=255)
                                                                  ccAAATCCACATGGTCATGCTTTTCAGGGGGAAGTTCTTCTGGAAGGATGCGCAGGCTGCCAGACATCCCTAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACtggttgg                                                                                 >  8:18090/1‑149 (MQ=255)
                                                                  caaaaTCCACATGTTCATGCTTTACAGGGGGAAGTTCTACTGGAAGGATGCGCAGGCTGCCAGACATCCCTAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGCTATTTAGCCCGCCGCTGCACGCTAAGCACtggttgg                                                                                 >  5:172775/3‑149 (MQ=255)
                                                                    aaaTCCACATGGTCATGCTTTTCAGGGGGAAGTTCTTCTGGAAGGATGCGCAGGCTGCCAGACATCCCTAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTc                                                                               <  2:200024/149‑1 (MQ=255)
                                                                      aTCCACATGGTCATGCTTTTCAGGGGGAAGTTCTTCTGGAAGGATGCGCAGGCTGCCAGACATCCCTAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCGc                                                                             <  7:358531/149‑1 (MQ=255)
                                                                        ccACATGGTCATGCTTTTCAGGGGGAAGTTCTTCTGGAAGGATGCGCAGGCTGCCAGACATCCCTAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCGCt                                                                            >  3:377993/1‑148 (MQ=255)
                                                                         cacaTGGTCATGCTTTTCAGGGGGAAGTTCTTCTGGAAGGATGCGCAGGCTGCCAGACATCCCTAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCGCtt                                                                           >  7:400734/1‑148 (MQ=255)
                                                                              ggTCATGCTTTTCAGGGGGAAGTTCTTCTGGAAGGATGCGCAGGCTGCCAGACATCCCTAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCGCTTAAACgg                                                                     >  8:382283/1‑149 (MQ=255)
                                                                                tCATGCTTTTCAGGGGGAAGTTCTTCTGGAAGGATGCGCAGGCTGCCAGACATCCCTAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCGCTTAAACGGTa                                                                   <  7:382283/149‑1 (MQ=255)
                                                                                  aTGCTTTTCAGGGGGAAGTTCTTCTGGAAGGATGCGCAGGCTGCCAGACATCCCTAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCGCTTAAACgg                                                                     >  7:426005/1‑145 (MQ=255)
                                                                                  aTGCTTTTCAGGGGGAAGTTCTTCTGGAAGGATGCGCAGGCTGCCAGACATCCCTAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCGCTTAAACgg                                                                     >  5:361583/1‑145 (MQ=255)
                                                                                      ttttCAGGGGTAAGTTCTTCTGGAAGGATGCGCAGGCTGCCAGACATCACTAAATGAATGATAATCAAGCCCTCAGGCAGCTCCAGCAGCAGAAATTTAGCCAGCCGCTGCAAGCTAAGCACTGGTTGGACGCTTAAACGGTAGAtctc                                                             >  3:300348/1‑149 (MQ=255)
                                                                                      ttttCAGGGGGAAGTTCTTCTGGAAGGATGCGCAGGCTGCCAGACATCCCTAAATGAATGATAATCCAGCCATCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCGCTTAAACGGTAGAtctc                                                             >  3:381821/1‑149 (MQ=255)
                                                                                         tCAGGGGGAATTTCTTCTGGAAGGTTGCGCAGGCTGCCAGACATGCCTCAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCGCTTAAACTGTAGATCtctt                                                           <  6:172775/148‑1 (MQ=255)
                                                                                          cAGGGGGAAGTTCTTCTGGAAGGATGCGCAGGCTGCCAGACATCCCTAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGATAAGCACTGGTTGGTCGCTTAAACGGTAGATCtcttc                                                          <  4:412392/148‑1 (MQ=255)
                                                                                           aGGGGGAAGTTCTTCTGGAAGGATGCGCAGGCTGCCAGACATCCCTAAATGAATGGTAATCCAGCCCTCAGGCAGCTCCAGCATCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCGCTTAAACGGTAGATCtcttct                                                         <  5:113027/148‑1 (MQ=255)
                                                                                           aGGGGGAAGTTCTTCTGGAAGGATGCGCAGGCTGCCAGACATCCCTAAATGAATGGTAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTGGGTCGCTTAAACGGTAGATCtcttct                                                         <  5:106762/148‑1 (MQ=255)
                                                                                           aGGGGGAAGTTCTTCTGGAAGGATGCGCAGGCTGCCAGACATCCCTAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCGCTTAAACGGTAGATCtcttct                                                         <  7:18090/148‑1 (MQ=255)
                                                                                               ggAAGTTCTTCTGGAAGGATGCGCAGGCTGCCAGACATCCCTAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCGCTTAAACGGTAGATCtctt                                                           >  1:321252/1‑142 (MQ=255)
                                                                                               ggAAGTTCTTCTGGAAGGATGCGCAGGCTGCCAGACATCCCTAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCGCTTAAACGGTAGATCtctt                                                           <  2:321252/142‑1 (MQ=255)
                                                                                               ggAAGTTCTTCTGGAAGGATGCGCAGGCTGCCAGACATCCCTAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCGCTTAAACGGTAGATCTCTTCTGAAAc                                                    >  8:417448/1‑149 (MQ=255)
                                                                                                gAAGTTCTTCTGGAAGGATGCGCAGGCTGCCAGACATCCCTAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCGCTTAAACGGTAGATCTCTTCTGAAAcc                                                   >  1:287119/1‑149 (MQ=255)
                                                                                                   gTTCTTCTGGAAGGATGCGCAGGCTGCCAGACATCCCTAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACtgg                                                                                     <  1:318201/112‑1 (MQ=255)
                                                                                                   gTTCTTCTGGAAGGATGCGCAGGCTGCCAGACATCCCTAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACtgg                                                                                     >  2:318201/1‑112 (MQ=255)
                                                                                                      cttctGGAAGGATGCGCATGCTGCCAGACATCCCTAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCGCTTAAACGGTAGATCTCTTCTGAAACCGGCCa                                              <  6:296058/148‑1 (MQ=255)
                                                                                                       ttctGGAAGGATGCGCAGGCTGCCAGACATCCCTAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCGCTTaaa                                                                        >  7:414448/1‑121 (MQ=255)
                                                                                                       ttctGGAAGGATGCGCAGGCTGCCAGACATCCCTAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCGCTTaaa                                                                        <  8:414448/121‑1 (MQ=255)
                                                                                                          tGGAAGGATGCGCAGGCTGCCAGACATCCCTAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCGCTTAAACGGTa                                                                   <  4:64093/123‑1 (MQ=255)
                                                                                                          tGGAAGGATGCGCAGGCTGCCAGACATCCCTAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCGCTTAAACGGTa                                                                   >  3:64093/1‑123 (MQ=255)
                                                                                                          tGGAAGGATGCGCAGGCTGCCAGACATCCCTAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCGCTTAAACGGTAGATCTCTTCTGAAACCGGCCAgcgc                                          >  7:112522/1‑148 (MQ=255)
                                                                                                            gAAGGATGCGCAGGCTGCCAGACATCCCTAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCGCTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCaa                                        <  6:366528/148‑1 (MQ=255)
                                                                                                            gAAGGATGCGCAGGCTGCCAGACATCCCTAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCGCTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCaa                                        <  8:112522/148‑1 (MQ=255)
                                                                                                               ggATGCGCAGGCTGCCAGACATCCCTAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCGCTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGt                                    <  6:361583/149‑1 (MQ=255)
                                                                                                                    cgcAGGCTGCCAGACATCCCTAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCGCTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGtt                               <  8:400734/149‑1 (MQ=255)
                                                                                                                     gcAGGCTGCCAGACATCCCTAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCGCTTAAACGGTAGATCTCTTCTg                                                        >  1:134081/1‑123 (MQ=255)
                                                                                                                     gcAGGCTGCCAGACATCCCTAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCGCTTAAACGGTAGATCTCTTCTg                                                        <  2:134081/123‑1 (MQ=255)
                                                                                                                             ccAGACATCCCTAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCGCTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGcaccac                      >  5:236932/1‑149 (MQ=255)
                                                                                                                                     tccTAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTTGCCCGCCGCTGCACGCTAAGCACTGGTGGGTCGCTTAAACGGTATATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGCACCACTGCATgaa              <  1:149894/148‑1 (MQ=255)
                                                                                                                                      ccTAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCGCTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGCACCACTGCATgaa              <  4:377993/148‑1 (MQ=255)
                                                                                                                                      ccTAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCGCTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGCACCACTGCATgaa              <  2:435101/148‑1 (MQ=255)
                                                                                                                                       cTAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCGCTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGCACCACTGCATgaag             >  7:329620/1‑148 (MQ=255)
                                                                                                                                       cTAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCGCTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGCACCACTGCATgaag             >  8:460493/1‑148 (MQ=255)
                                                                                                                                        tAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCGCTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGCACCACTGCATgaagaa           >  1:387652/1‑149 (MQ=255)
                                                                                                                                            tgaatgaTAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCGCTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGCACCACTGCATGAAGAATGGt       >  8:166557/1‑149 (MQ=255)
                                                                                                                                               atgaTAATCCAGCCCTCAGGCAGCTCCCGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCGCTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGCACCACTGCATGAAGAATGGTTg     >  4:374634/1‑148 (MQ=255)
                                                                                                                                               atgaTAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCGCTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGCACCACTGCATGAAGAATGGTTGc    <  1:149881/149‑1 (MQ=255)
                                                                                                                                               atgaTAATCCAGACCTCAGGCAGCTACAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTTGTTGGTCGCTTAAACGGTAGATCTCTTCTGAAACCGGCCAGTGCAAGCGTCCGTTGCGCACCTCTGCATGAAGAATGGTTGc    <  4:300348/149‑1 (MQ=255)
                                                                                                                                                  aTAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCGCTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGCACCACTGCATGAAGAATGGTTGCAc  >  4:384408/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
TGGTCCACAGCCAGGCACCAAAGCGGCGCGGATCGGTGTAGCGCAGCACTTTGCCGTTGCTCATCACCAAATCCACATGGTCATGCTTTTCAGGGGGAAGTTCTTCTGGAAGGATGCGCAGGCTGCCAGACATCCCTAAATGAATGATAATCCAGCCCTCAGGCAGCTCCAGCAGCAGATATTTAGCCCGCCGCTGCACGCTAAGCACTGGTTGGTCGCTTAAACGGTAGATCTCTTCTGAAACCGGCCAGCGCAAGCGTCCGTTGCGCACCACTGCATGAAGAATGGTTGCAC  >  NC_000913/3810801‑3811094

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: