Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,714,168 Δ1 bp 5.8% coding (2227/2334 nt) bisC ← biotin sulfoxide reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,714,1680A. 5.8% 369.1 / 14.9 104coding (2227/2334 nt)bisCbiotin sulfoxide reductase
Reads supporting (aligned to +/- strand):  ref base A (54/44);  new base . (2/4);  total (56/48)
Fisher's exact test for biased strand distribution p-value = 4.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.36e-01

Only 100 of 105 total aligned reads displayed.

GATTGCGCATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGACCGCTCCGGCAAGAATTTGC  >  NC_000913/3714021‑3714316
                                                                                                                                                    |                                                                                                                                                    
gATTGCGCATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAg                                                                                                                                                     >  3:361662/1‑149 (MQ=255)
   tGCGCATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTg                                                                                                                                                  <  4:366679/149‑1 (MQ=255)
       cATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTgggaa                                                                                                                                              <  4:620156/149‑3 (MQ=255)
         tGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCTCGCG‑GCTGGGgaat                                                                                                                                            <  4:124512/149‑3 (MQ=255)
         tGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCTCGCG‑GCTGGGgaa                                                                                                                                             <  2:199985/148‑2 (MQ=255)
         tGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGgag                                                                                                                                             >  1:248736/1‑148 (MQ=255)
          ggCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGgaga                                                                                                                                            <  2:629372/148‑1 (MQ=255)
                aTGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTCCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGT‑CGCGAGCTGGGGAGATCTTTgg                                                                                                                                     >  3:554156/1‑149 (MQ=255)
                 tGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTgg                                                                                                                                     <  1:173618/148‑1 (MQ=255)
                 tGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTgg                                                                                                                                     >  2:601840/1‑148 (MQ=255)
                     aCCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGc                                                                                                                                >  4:243653/1‑149 (MQ=255)
                      cccAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGc                                                                                                                                >  2:410875/1‑148 (MQ=255)
                       ccAACTACCCACATGGATTCTGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTgg                                                                                                                                                 >  3:749290/1‑130 (MQ=255)
                       ccAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTgg                                                                                                                                                 <  4:749290/130‑1 (MQ=255)
                        cAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCAc                                                                                                                              >  2:676840/1‑148 (MQ=255)
                         aaCTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACg                                                                                                                             >  2:707717/1‑148 (MQ=255)
                          aCTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACg                                                                                                                             <  1:128600/147‑1 (MQ=255)
                          aCTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACg                                                                                                                             >  2:128600/1‑147 (MQ=255)
                            tACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTc                                                                                                                          <  2:694085/148‑1 (MQ=255)
                               ccACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATCCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTg                                                                                                                      >  3:92025/1‑149 (MQ=255)
                                 acaTGGATTATGAGCTGGCCGGTGGTTTAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTgcg                                                                                                                    >  4:125764/1‑149 (MQ=255)
                                 acaTGGATTATGAGCTGGCCGGTGGTTCACACGCTGTAAGTGTCAGTTCGGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCCCGTTCACTgcg                                                                                                                    >  4:804193/1‑149 (MQ=255)
                                 acaTGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTgcg                                                                                                                    >  3:500751/1‑149 (MQ=255)
                                 acaTGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGCATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTgg                                                                                                                                                 >  3:372111/1‑120 (MQ=255)
                                 acaTGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGCATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTgg                                                                                                                                                 <  4:372111/120‑1 (MQ=255)
                                  caTGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCttt                                                                                                                                       <  1:719391/129‑1 (MQ=255)
                                  caTGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCttt                                                                                                                                       >  2:719391/1‑129 (MQ=255)
                                  caTGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTgcg                                                                                                                    >  1:271996/1‑148 (MQ=255)
                                    tGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCTCGCG‑GCTGGGgaa                                                                                                                                             <  2:544356/121‑2 (MQ=255)
                                    tGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGCTCGCG‑GCTGGGgaa                                                                                                                                             >  1:544356/1‑120 (MQ=255)
                                    tGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGcc                                                                                                                  >  1:471707/1‑148 (MQ=255)
                                      gATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGtt                                                                                                               >  4:332366/1‑149 (MQ=255)
                                       aTTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGtt                                                                                                               >  1:196309/1‑148 (MQ=255)
                                          aTGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTa                                                                                                           <  4:68356/149‑1 (MQ=255)
                                            gAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACa                                                                                                         <  4:500751/149‑1 (MQ=255)
                                            gAGCTGGCCGGTGGTTCAAACGCTGTAAGCGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACa                                                                                                         >  3:737059/1‑149 (MQ=255)
                                                tGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAAt                                                                                                      >  1:534268/1‑148 (MQ=255)
                                                           tCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCgg                                                                                           >  2:845648/1‑148 (MQ=255)
                                                           tCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCgg                                                                                           >  1:557208/1‑148 (MQ=255)
                                                            cAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGt                                                                                          <  1:688249/148‑1 (MQ=255)
                                                               aCGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTaa                                                                                       >  2:222206/1‑148 (MQ=255)
                                                               aCGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTaa                                                                                       >  2:674341/1‑148 (MQ=255)
                                                                cGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTAGTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAAt                                                                                     <  4:737059/149‑1 (MQ=255)
                                                                 gCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATc                                                                                    >  4:75169/1‑149 (MQ=255)
                                                                   tGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATcc                                                                                   <  1:346127/148‑1 (MQ=255)
                                                                   tGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCa                                                                                  >  4:394465/1‑149 (MQ=255)
                                                                   tGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCa                                                                                  <  3:243653/149‑1 (MQ=255)
                                                                     tAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAg                                                                                 >  1:327015/1‑148 (MQ=255)
                                                                     tAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGa                                                                                >  4:140846/1‑149 (MQ=255)
                                                                      aaGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGa                                                                                <  2:248736/148‑1 (MQ=255)
                                                                         tgtCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCg                                                                            >  3:649628/1‑149 (MQ=255)
                                                                          gtCAGTCCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCTCAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAACTACCGTCAGCGGTTAAATCCAGATCCgg                                                                           <  4:554156/149‑1 (MQ=255)
                                                                           tCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAAt                                                                                     >  3:333257/1‑138 (MQ=255)
                                                                            cAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGc                                                                          >  1:101622/1‑148 (MQ=255)
                                                                              gTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCAt                                                                       <  3:241482/149‑1 (MQ=255)
                                                                                   ggACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGcccc                                                                  <  3:781983/149‑1 (MQ=255)
                                                                                     aCCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATg                                                                      <  2:31295/143‑1 (MQ=255)
                                                                                     aCCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATg                                                                      >  1:31295/1‑143 (MQ=255)
                                                                                     aCCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCtt                                                                >  4:599728/1‑149 (MQ=255)
                                                                                        gTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCg                                                              >  1:44792/1‑148 (MQ=255)
                                                                                         ttGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGCGATCTTTGGTCAGCACGTTCACTGCGCCGtttt                                                                                                             <  3:571610/100‑1 (MQ=255)
                                                                                         ttGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGtttt                                                                                                             >  4:571610/1‑100 (MQ=255)
                                                                                           gTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGaa                                                          <  3:332366/149‑1 (MQ=255)
                                                                                           gTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGaa                                                          <  4:732683/149‑1 (MQ=255)
                                                                                                tttCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCa                                                      >  2:469919/1‑148 (MQ=255)
                                                                                                      gCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCAc                                                >  1:755698/1‑148 (MQ=255)
                                                                                                      gCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCAc                                                >  1:152513/1‑148 (MQ=255)
                                                                                                       ccATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGc                                              <  3:602619/149‑1 (MQ=255)
                                                                                                       ccATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGc                                              <  4:102394/149‑1 (MQ=255)
                                                                                                          tGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCa                                            >  2:479514/1‑148 (MQ=255)
                                                                                                               gcgcCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAgg                                          >  1:117453/1‑145 (MQ=255)
                                                                                                               gcgcCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAgg                                          <  2:117453/145‑1 (MQ=255)
                                                                                                                cgcCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTaa                                     <  4:361662/149‑1 (MQ=255)
                                                                                                                 gcCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTaa                                     <  1:410875/148‑1 (MQ=255)
                                                                                                                 gcCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAAt                                    <  3:125764/149‑1 (MQ=255)
                                                                                                                   cGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAAtt                                   >  2:153853/1‑148 (MQ=255)
                                                                                                                    gTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCCAATCACGCCAGGTTTAATTc                                  <  2:237482/148‑1 (MQ=255)
                                                                                                                     tATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTcc                                 >  1:448165/1‑148 (MQ=255)
                                                                                                                     tATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTcc                                 <  2:327015/148‑1 (MQ=255)
                                                                                                                        tACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCg                             >  4:533347/1‑149 (MQ=255)
                                                                                                                         aCCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATc                                                                                    <  3:809916/93‑1 (MQ=255)
                                                                                                                         aCCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATc                                                                                    >  4:809916/1‑93 (MQ=255)
                                                                                                                         aCCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGc                            >  3:327256/1‑149 (MQ=255)
                                                                                                                          cccGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCt                           >  4:549696/1‑149 (MQ=255)
                                                                                                                               acaGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATc                                                                                    >  4:298896/1‑87 (MQ=255)
                                                                                                                               acaGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATc                                                                                    <  3:298896/87‑1 (MQ=255)
                                                                                                                                caGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGa                      <  2:316584/148‑1 (MQ=255)
                                                                                                                                caGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGa                      <  2:44792/148‑1 (MQ=255)
                                                                                                                                 aGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGAcc                    >  4:134057/1‑149 (MQ=255)
                                                                                                                                 aGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGAcc                    <  3:394465/149‑1 (MQ=255)
                                                                                                                                      ttCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCa                                            >  2:407203/1‑120 (MQ=255)
                                                                                                                                      ttCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCa                                            <  1:407203/120‑1 (MQ=255)
                                                                                                                                       tCCCCAGCTCGCG‑GCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGaa                                                          >  3:27109/1‑105 (MQ=255)
                                                                                                                                       tCCCCAGCTCGCG‑GCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGaa                                                          <  4:27109/105‑1 (MQ=255)
                                                                                                                                         cccAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGACCGCTCCgg             >  1:429617/1‑148 (MQ=255)
                                                                                                                                           cAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGACCGCTCCGGCa           <  1:676840/148‑1 (MQ=255)
                                                                                                                                            aGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGt                                                             >  4:818690/1‑97 (MQ=255)
                                                                                                                                            aGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGt                                                             <  3:818690/97‑1 (MQ=255)
                                                                                                                                                cgcgAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGACCGCTCCGGCAAGaa       <  2:152513/148‑1 (MQ=255)
                                                                                                                                                    aGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGACCGCTCCGGCAAGAATTTGc  <  4:649628/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GATTGCGCATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGACCGCTCCGGCAAGAATTTGC  >  NC_000913/3714021‑3714316

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑