Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,842,195 A→G 5.3% A402A (GCA→GCG ascB → cryptic 6‑phospho‑beta‑glucosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,842,1950AG 5.3% 233.4 / 2.9 76A402A (GCA→GCGascBcryptic 6‑phospho‑beta‑glucosidase
Reads supporting (aligned to +/- strand):  ref base A (38/33);  new base G (2/2);  total (40/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.29e-01

GTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACATATCCGCGCAATGGGCGAAGCGATTGCAGACGGCATTCCGCTGATGGGCTACACCACATGGGGCTGTATTGATTTAGTTTCCGCCTCTACGGGTGAAATGAGCAAACGCTACGGCTTTGTCTTTGTTGACCGTGACGACGCAGGCAACGGTACGCTGA  >  NC_000913/2842050‑2842343
                                                                                                                                                 |                                                                                                                                                    
gTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACATATCCGCGCAATg                                                                                                                                                   <  6:298424/149‑1 (MQ=255)
gTATCACCAGGAATCTGATGTACGAGCGCTATCAGAAGCAGCTGATGCTGGAGGAAAACGGCCTGTGCGCAAAATATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACATATCCGCGCaagg                                                                                                                                                   <  6:298431/149‑3 (MQ=255)
 tATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCGGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACATATCCGCGCAATgg                                                                                                                                                  <  3:26059/149‑1 (MQ=255)
  aTCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACATATCCGCGCAATggg                                                                                                                                                 >  7:200094/1‑149 (MQ=255)
  aTCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAAAGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACATATCCGCGCAATggg                                                                                                                                                 >  1:230214/1‑149 (MQ=255)
     aCCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACATATCCGCGCAATGGGCg                                                                                                                                               >  3:187594/1‑148 (MQ=255)
      ccATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACATATCCGCGCAATGGGCg                                                                                                                                               >  1:159781/1‑147 (MQ=255)
      ccATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATAAGCTACTTACGCGAACATATCCGCGCAATGGGCg                                                                                                                                               >  3:2199/1‑147 (MQ=255)
       cATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACATATCCGCGCAATGGGCGAAg                                                                                                                                            <  4:309986/149‑1 (MQ=255)
       cATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACATATCCGCGCAATGGGCGAAg                                                                                                                                            >  2:329636/1‑149 (MQ=255)
       cATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACATATCCGCGCAATGGGCGAAg                                                                                                                                            <  3:315648/149‑1 (MQ=255)
           aaTATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACATATCCGCGCAATGGGCGAAGCGa                                                                                                                                         <  4:29626/148‑1 (MQ=255)
           aaTATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACATATCCGCGCAATGGGCGAAGCGAt                                                                                                                                        <  1:334335/149‑1 (MQ=255)
               tggtgTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACATATCCGCGCAATGGGCGAAGCGATTGc                                                                                                                                     <  6:341090/145‑1 (MQ=255)
               tgatgTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACATATCCGCGCAATGGGCGAAGCGATTGc                                                                                                                                     <  2:51209/148‑1 (MQ=255)
               tgatgTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACATATCCGCGCAATGGGCGAAGCGATTGc                                                                                                                                     <  6:20372/148‑1 (MQ=255)
                 atgTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACATATCCGCGCAATCGGCGAAGCGATTGCAGa                                                                                                                                  >  3:56468/1‑149 (MQ=255)
                         tCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACATATCCGCGCAATGGGCGAAGCGATTGCAGACGGCATTc                                                                                                                          >  7:271059/1‑149 (MQ=255)
                           gTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACATATCCGCGCAATGGGCGAAGCGATTGCAGACGGCATTCCg                                                                                                                        >  7:112037/1‑149 (MQ=255)
                               tCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACATATCCGCGCAATGGGCGAAGCGATTGCAGACGGCATTCCGCTGa                                                                                                                    <  5:69473/149‑1 (MQ=255)
                                    aGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACATATCCGCGCAATGGGCGAAGCGATTGCAGACGGCATTCCGCTGATGGGc                                                                                                               <  7:268035/149‑1 (MQ=255)
                                     gCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACATATCCGCGCAATGGGCGAAGCGATTGCAGACGGCATTCCGCTGATGGGc                                                                                                               >  1:108955/1‑148 (MQ=255)
                                     gCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACATATCCGCGCAATGGGCGAAGCGATTGCAGACGGCATTCCGCTGATGGGc                                                                                                               >  8:178869/1‑148 (MQ=255)
                                        gCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACATATCCGCGCAATGGGCGAAGCGATTGCAGACGGCATTCCGCTGATGGGCTaca                                                                                                           >  1:28017/1‑149 (MQ=255)
                                            tttCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACATATCCGCGCAATGGGCGAAGCGATTGCAGACGGCATTCCGCTCATGGGCTAcacca                                                                                                        >  1:45601/1‑148 (MQ=255)
                                              tCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACATATCCGCGCAATGGGCGAAGCGATTGCAGACGGCATTCCGCTGATGGGCTACACcaca                                                                                                      <  5:276287/148‑1 (MQ=255)
                                                  gtgtAAAACGGCCTGGTCGCAAATGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACATATCCGGGCAATGGGCGAAGCGATTGCAGACGGCATTCCGCTGATGGGCTACACCACATgggg                                                                                                 <  8:92021/149‑1 (MQ=255)
                                                     gAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACATATCCGCGCAATGGGCGAAGCGATTGCAGACGGCATTCCGCTGATGGGCTACACCACATGGGGCTg                                                                                              >  2:90327/1‑149 (MQ=255)
                                                       aaaCGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACATATCCGCGCAATGGGCGAAGCGATTGCAGACGGCATTCCGCTGATGGGCTACACCACATGGGGCTGTa                                                                                            >  8:315321/1‑149 (MQ=255)
                                                       aaaCGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACATATCCGCGCAATGGGCGAAGCGATTGCAGACGGCATTCCGCTGATGGGCTACACCACATGGGGCTGTa                                                                                            >  8:182765/1‑149 (MQ=255)
                                                         aCGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGCCGACTATCGCATCAGCTACTTACGCGAACATATCCGCGCAATGGGCGAAGCGATTGCAGCCGGCAGTCCGCTGATGGGCTACACCACATGGGGCGGTAtt                                                                                          >  2:235747/1‑149 (MQ=255)
                                                          cGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACATATCCGCGCAATGGGCGAAGCGATTGCAGACGGCATTCCGCTGATGGGCTACACCACATGGGGCTGTATTg                                                                                         <  4:187594/149‑1 (MQ=255)
                                                              cTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACATATCCGCGCAATGGGCGAAGCGATTGCAGACGGCATTCCGCTGATGGGCTACACCCCATGGGGCTGTATTGAtt                                                                                      >  6:40248/1‑148 (MQ=255)
                                                              cTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACATATCCGCGCAATGGGCGAAGCGATTGCAGACGGCATTCCGCTGATGGGCTACACCACATGGGGCTGTATTGAtt                                                                                      >  2:286141/1‑148 (MQ=255)
                                                                 ggCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACATATCCGCGCAATGGGCGAAGCGATTGCAGACGGCATTCCGCTGATGGGCTACACCACATGGGGCTGTATTGATTTa                                                                                    <  7:178869/147‑1 (MQ=255)
                                                                  gcgcAAAAGAGGAATTTGCTGCCTATGGCGAGATTCACGACGACTATCGCATCAGCTACTTACGCGAACATATCCGCGCAATGGGCGAAGCGATTGCAGACGGCATTCCGCTGATGGGCTACACCACATGGGGCTGTATTGATTTAGtt                                                                                 <  2:324902/149‑1 (MQ=255)
                                                                    gcAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACATATCCGCGCAATGGGCGAAGCGATTGCAGACGGCATTCCGCTGATGGGCTACACCACATGGGGCTGTATTGATTTAGTTTc                                                                               <  7:315321/149‑1 (MQ=255)
                                                                    gcAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACATATCCGCGCAATGGGCGAAGCGATTGCAGACGGCATTCCGCTGATGGGCTACACCACATGGGGCTGTATTGATTTAGTTTc                                                                               <  7:223390/149‑1 (MQ=255)
                                                                     cAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACATATCCGCGCAATGGGCGAAGCGATTGCAGACGGCATTCCGCTGATGGGCTACACCACATGGGGCTGTATTGATTTAGTTTcc                                                                              <  8:112037/149‑1 (MQ=255)
                                                                      aaaaGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACATATCCGCGCAATGGGCGAAGCGATTGCAGACGGCATTCCGCTGATGGGCTACACCACATGGGGCTGTATTGATTTAGTTTCCg                                                                             <  2:108955/149‑1 (MQ=255)
                                                                      aaaaGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACATATCCGCGCAATGGGCGAAGCGATTGCAGACGGCATTCCGCTGATGGGCTACACCACATGGGGCTGTATTGATTTAGTTTCCg                                                                             >  1:215452/1‑149 (MQ=255)
                                                                      aaaaGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACATATCCGCGCAATGGGCGAAGCGATTGCAGACGGCATTCCGCTGATGGGCTACACCACATGGGGCTGTATTGATTTAGTTTCCg                                                                             >  7:140502/1‑149 (MQ=255)
                                                                       aacgATGAATTTGCTGCCAATGGCGAGATTACCGACGACTACCGCATCAGCTACTTACGCGAACATATCCGCGCAATGGGCGAAGCGCTTGCAGACGGCATTCCGCTGATGGGCTACACCACATGGGGCTGTATTGATTTAGTTTCCGc                                                                            <  6:108913/146‑1 (MQ=255)
                                                                                   gctgcCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACATATCCGCGCAATGGGCGAAGCGATTGCAGACGGCATTCCGCTGATGGGCTACACCACATGGGGCTGTATTGATTTAGTTTCCGCCTCTACGGGTg                                                                 >  2:157403/1‑148 (MQ=255)
                                                                                       ccAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACATATCCGCGCAATGGGCGAAGCGATTGCAGACGGCATTCCGCTGATGGGCTACCCCACATGGGGCTGTATTGATTTAGTTTCCGCCTCTACGGGTGAAAt                                                             >  5:310735/1‑148 (MQ=255)
                                                                                        cAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACATATCCGCGCAATGGGCGAAGCGATTGCAGACGGCATTCCGCTGATGGGCTACACCACATGGGGCTGTATTGATTTAGTTTCCGCCTCTACGGGTGAAATGa                                                           <  1:90327/149‑1 (MQ=255)
                                                                                             gCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACATATCCGCGCTATGGGCGAAGCGATTGCAGACGGCATTCCGGTGTTGGGCTACACCACATGGGGCTGTATTGATTTAGTTTCCGCCTCTACGGGTGAAATGAGCaaa                                                      <  4:2199/149‑1 (MQ=255)
                                                                                                agaTTAACGACGACTATCGCATCAGCTACTTACGCGAACATATCCGCGCAATGGGCGAAGCGATTGCAGACGGCATTCCGCTGATGGGCTACACCACATGGGGCTGTATTGATTTAGTTTCCGCCTCTACGGGTGAAATGAGCAAACGc                                                   >  7:334772/1‑149 (MQ=255)
                                                                                                 gaTTAACGACGACTATCGCATCAGCTACTTACGCGAACATATCCGCGCAATGGGCGAAGCGATTGCAGACGGCATTCCGCTGATGGGCTACACCACATGGGGCTGTATTGATTTAGTTTCCGCCTCTACGGGTGAAATGAGCAAACGc                                                   >  6:7605/1‑148 (MQ=255)
                                                                                                  aTTAACGACGACTATCGCATCAGCTACTTACGCGAACATATCCGCGCAATGGGCGAAGCGATTGCAGACGGCATTCCGCTGATGGGCTACACCACATGGGGCTGTATTGATTTAGTTTCCGCCTCTACGGGTGAAATGAGCAAACGCTa                                                 <  1:157403/149‑1 (MQ=255)
                                                                                                     aaCGACGACTATCGCATCAGCTACTTACGCGAACATATCCGCGCAATGGGCGAAGCGATTGCAGACGGCATTCCGCTGATGGGCTACACCACATGGGGCTGTATTGATTTAGTTTCCGCCTCTACGGGTGAAATGAGCAAACGCTACg                                               <  2:215452/148‑1 (MQ=255)
                                                                                                      acgacgACTATCGCATCAGCTACTTACGCGAACATATCCGCGCAATGGGCGAAGCGATTGCAGACGGCATTCCGCTGATGGGCTACACCCCATGGGGCTGTATTGATTTAGTTTCCGCCTCTACGGGTGAAATGAGCAAACGCTACGGc                                             >  8:205941/1‑149 (MQ=255)
                                                                                                         acgacTATCGCATCAGCTACTTACGCGAACATATCCGCGCAATGGGCGAAGCGATTGCAGACGGCATTCCGCTGATGGGCTACACCACATGGGGCTGTATTGATTTAGTTTCCGCCTCTACGGGTGAAATGAGCAAACGCTACGGcttt                                          <  5:7605/149‑1 (MQ=255)
                                                                                                         acgacTATCGCATCAGCTACTTACGCGAACATATCCGCGCAATGGGCGAAGCGATTGCAGACGGCATTCCGCTGATGGGCTACACCACATGGGGCTGTATTGATTTAGTTTCCGCCTCTACGGGTGAAATGAGCAAACGCTACGGcttt                                          >  5:2404/1‑149 (MQ=255)
                                                                                                         acgacTATCGCATCAGCTACTTACGCGAACATATCCGCGCAATGGGCGAAGCGATTGCAGACGGCATTCCGCTGATGGGCTACACCACATGGGGCTGTATTGATTTAGTTTCCGCCTCTACGGGTGAAATGAGCAAACGCTACGGcttt                                          <  2:45601/149‑1 (MQ=255)
                                                                                                           gacTATCGCATCAGCTACTTACGCGAACATATCCGCGCAATGGGCGAAGCGATTGCAGACGGCATTCCGCTGATGGGCTACACCACATGGGGCTGTATTGATTTAGTTTCCGCCTCTACGGGTGAAATGAGCAAACGCTACGGctttgt                                        >  6:163915/1‑149 (MQ=255)
                                                                                                            acTATCGCATCAGCTACTTACGCGAACATATCCGCGCGATGGGCGAAGCGATTGCAGACGGCATTCCGCTGATGGGCTACACCACATGGGGCTGTATTGATTTAGTTTCCGCCTCTACGGGTGAAATGAGCAAACGCTACGGctttgtc                                       >  3:305780/1‑149 (MQ=255)
                                                                                                            acTATCGCATCAGCTACTTACGCGAACATATCCGCGCGATGGGCGAAGCGATTGCAGACGGCATTCCGCTGATGGGCTACACCACATGGGGCTGTATTGATTTAGTTTCCGCCTCTACGGGTGAAATGAGCAAACGCTACGGctttgtc                                       >  3:161685/1‑149 (MQ=255)
                                                                                                              tATCGCATCAGCTACTTACGCGAACATATCCGCGCAATGGGCGAAGCGATTGCAGACGGCATTCCGCTGATGGGCTACACCACATGGGGCTGTATTGATTTAGTTTCCGCCTCTACGGGTGAAATGAGCAAACGCTACGGctttgtctt                                     >  5:32723/1‑149 (MQ=255)
                                                                                                              tATCGCATCAGCTACTTACGCGAACATATCCGCGCAATGGGCGAAGCGATTGCAGACGGCATTCCGCTGATGGGCTACACCACATGGGGCTGTATTGATTTAGTTTCCGCAGCTACGGGTGAAATGAGCAAACGCTACGGctttgtctt                                     >  6:4199/1‑149 (MQ=255)
                                                                                                                tCGCATCAGCTACTTACGCGAACATATCCGCGCAATGGGCGAAGCGATTGCAGACGGCATTCCGCTGATGTGCTACACCACAGGGGGCTGTATTGATTTAGTTTCCGCCTCTACGGGTGAAATGAGCAAACGCTACGGctatgtctttc                                   >  2:78447/1‑148 (MQ=255)
                                                                                                                tCGCATCAGCTACTTACGCGAACATATCCGCGCAATGGGCGAAGCGATTGCAGACGGCATTCCGCTGATGGTCTACACCACATGGGGCTGTATTGATTTAGTTTCCGCCTCTACGGGTGAAATGAGCAAACGCTACGGctttgtcttt                                    <  1:235747/148‑1 (MQ=255)
                                                                                                                tCGCATCAGCTACTTACGCGAACATATCCGCGCAATGGGCGAAGCGATTGCAGACGGCATTCCGCTGATGGGCTACACCACATGGGGCTGTATTGATTTAGTTTCCGCCTCTACGGGTGAAATGAGCAAACGCTACGGctttgtcttt                                    <  7:182765/148‑1 (MQ=255)
                                                                                                                tCGCATCAGCTACTTACGCGAACATATCCGCGCAATGGGCGAAGCGATTGCAGACGGCATTCCGCTGATGGGCTACACCACATGGGGCTGTATTGATTTAGTTTCCGCCTCTACGGGTGAAATGAGCAAACGCTACGGctttgtcttt                                    <  1:329636/148‑1 (MQ=255)
                                                                                                                    aTCAGCTACTTACGCGAACATATCCGCGCGATGGGCGAAGCGATTGCAGACGGCATTCCGCTGATGGGCTACACCACATGGGGCTGTATTGATTTAGTTTCCGCCTCTACGGGTGAAATGAGCAAACGCTACGGCTTTGTCTTTGTTGa                               <  4:305780/149‑1 (MQ=255)
                                                                                                                    aTCAGCTACTTACGCGAACATATCCGCGCGATGGGCGAAGCGATTGCAGACGGCATTCCGCTGATGGGCTACACCACATGGGGCTGTATTGATTTAGTTTCCGCCTCTACGGGTGAAATGAGCAAACGCTACGGCTTTGTCTTTGTTGa                               <  4:161685/149‑1 (MQ=255)
                                                                                                                    aTCAGCTACTTACGCGAACATATCCGCGCAATGGGCGAAGCGATTGCAGACGGCATTCCGCTGATGGGCTACCCCACATGGGGCTGTATTGATTTAGTTTCCGCCTCTACGGGTGAAATGAGCAAACGCTACGGCTTTGTCTTTGTTGa                               <  8:271059/149‑1 (MQ=255)
                                                                                                                        gCTACTTACGCGAACATATCCGCGCAATGGGCGAAGCGATTGCAGACGGCATTCCGCTGATGGGCTACACCACATGGGGCTGTATTGATTTAGTTTCCGCCTCTACGGGTGAAATGAGCAAACGCTACGGCTTTGTCTTTGTTGACCGt                           >  6:61899/1‑149 (MQ=255)
                                                                                                                          tACTTACGCGAACACATCCGCGCAATGGGCGAAGCGATTGCAGACGGCATTCCGCTGATGGGCTACACCACATGGGGCTGTATTGATTTAGTTTCCGCCTCTACGGGTGAAATGAGCAAACGCTACGGCTTTGTCTTTGTTGACCGt                           >  2:295965/1‑147 (MQ=255)
                                                                                                                               aCGCGAACATATCCGCGCAATGGGCGAAGCGATTGCAGACGGCATTCCGCTGATGGGCTACACCACATGGGGCTGTATTGATTTAGTTTCCGCCTCTACGGGTGAAATGAGCAAACGCTACGGCTTTGTCTTTGTTGACCGTgacga                      <  1:286141/147‑1 (MQ=255)
                                                                                                                                 gcgAACATATCCGCGCAATGGGCGAAGCGATTGCAGACGGCATTCCGCTGATGGGCTACACCACATGGGGCTGTATTGATTTAGTTTCCGCCTCTACGGGTGAAATGAGCAAACGCTACGGCTTTGTCTTTGTTGACCGTGacgacg                    >  2:23461/1‑147 (MQ=255)
                                                                                                                                   gAACATATCCGCGCAATGGGCGAAGCGATTGCAGACGGCATTCCGCTGATGGGCTACACCACATGGGGCTGTATTGATTTAGTTTCCGCCTCTACGGGTGAAATGAGCAAACGCTACGGCTTTGTCTTTGTTGACCGTGACGACGCAgg                >  8:68908/1‑149 (MQ=255)
                                                                                                                                     aCATATCCGCGCAATGGGCGAAGCGATTGCAGACGGCATTCCGCTGATGGGCTACACCACATGGGGCTGTATTGATTTAGTTTCCGCCTCTACGGGTGAAATGAGCAAACGCTACGGCTTTGTCTTTGTTGACCGTGACGACGCAGGc               >  4:212701/1‑148 (MQ=255)
                                                                                                                                       atatCCGCGCAATGGGCGAAGCGATTGCAGACGGCATTCCGCTGATGGGCTACACCACATGGGGCTGTATTGATTTAGTTTCCGCCTCTACGGGTGAAATGAGCAAACGCTACGGctttgtctttgt                                  <  6:33691/127‑1 (MQ=255)
                                                                                                                                       atatCCGCGCAATGGGCGAAGCGATTGCAGACGGCATTCCGCTGATGGGCTACACCACATGGGGCTGTATTGATTTAGTTTCCGCCTCTACGGGTGAAATGAGCAAACGCTACGGctttgtctttgt                                  >  5:33691/1‑127 (MQ=255)
                                                                                                                                             gcgcAATGGGCGAAGCGATTGCAGACGGCATTCCGCTGATGGGCTACACCACATGGGGCTGTATTGATTTAGTTTCCGCCTCTACGGGTGAAATGAGCAAACGCTACGGCTTTGTCTTTGTTGACCGTGACGACGCAGGCAACGGTACg      >  4:110628/1‑149 (MQ=255)
                                                                                                                                               gcAATGGGCGAAGCGATTGCAGACGGCATTCCGCTGATGGGCTACACCACATGGGGCTGTATTGATTTAGTTTCCGCCTCTACGGGTGAAATGAGCAAACGCTACGGCTTTGTCTTTGTTGACCGTGACGACGCAGGCAACGGTACGCt    <  2:159781/149‑1 (MQ=255)
                                                                                                                                                 aaTGGGCGAAGCGATTGCAGACGGCATTCCGCTGATGGGCTACACCACATGGGGCTGTATTGATTTAGTTTCCGCCTCTACGGGTGAAATGAGCAAACGCTACGGCTTTGTCTTTGTTGACCGTGACGACGCAGGCAACGGTCCGCTGa  >  3:104707/1‑149 (MQ=255)
                                                                                                                                                 aaTGGGCGAAGCGATTGCAGACGGCATTCCGCTGATGGGCTACACCACATGGGGCTGTATTGATTTAGTTTCCGCCTCTACGGGTGAAATGAGCAAACGCTACGGCTTTGTCTTTGTTGACCGTGACGACGCAGGCAACGGTACGCTGa  <  1:149669/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GTATCACCATGAATATGATGTACGATCGTTATCAGAAGCCGCTGTTTCTGGTGGAAAACGGCCTGGGCGCAAAAGATGAATTTGCTGCCAATGGCGAGATTAACGACGACTATCGCATCAGCTACTTACGCGAACATATCCGCGCAATGGGCGAAGCGATTGCAGACGGCATTCCGCTGATGGGCTACACCACATGGGGCTGTATTGATTTAGTTTCCGCCTCTACGGGTGAAATGAGCAAACGCTACGGCTTTGTCTTTGTTGACCGTGACGACGCAGGCAACGGTACGCTGA  >  NC_000913/2842050‑2842343

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: