Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,849,031 T→C 5.6% G65G (GGT→GGC sppA → protease IV (signal peptide peptidase)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,849,0310TC 5.6% 222.1 / 2.6 73G65G (GGT→GGCsppAprotease IV (signal peptide peptidase)
Reads supporting (aligned to +/- strand):  ref base T (36/32);  new base C (2/2);  total (39/34)
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.23e-01

GCGTCTGCTGAATTTCGTCCGTGAAATGGTACTTAACCTGTTCTTTATTTTCCTCGTACTGGTTGGTGTGGGGATTTGGATGCAGGTCAGTGGTGGTGATTCGAAAGAAACGGCCAGTCGTGGCGCACTGCTGCTGGACATTTCTGGTGTGATCGTCGATAAACCCGACAGTTCTCAGCGGTTTAGTAAATTAAGCCGCCAGCTGCTTGGTGCCAGTTCCGATCGTCTGCAGGAAAACTCACTGTTTGATATCGTCAACACTATTCGCCAGGCGAAGGACGACCGCAATATCACCG  >  NC_000913/1848884‑1849179
                                                                                                                                                   |                                                                                                                                                    
gCGTCTGCTGAATTTCGTCCGTGAAATGGTACTTAACCTGTTCTTTATTTTCCTCGTACTGGTTGGTGTGGGGATTTGGATGCAGGTCAGTGGTGGTGATTCGAAAGAAACGGCCAGTCGTGGCGCACTGCTGCTGGACATTTCTGgt                                                                                                                                                      <  3:85527/148‑1 (MQ=255)
gCGTCTGCTGAATTTCGTCCGTGAAATGGTACTTAACCTGTTCTTTATTTTCCTCGTACTGGTTGGTGTGGGGATTTGGATGCAGGTCAGTGGTGGTGATTCGAAAGAAACGGCCAGTCGTGGCGCACTGCTGCTGGACATTTCTGgt                                                                                                                                                      <  3:290535/148‑1 (MQ=255)
gCGTCTGCTGAATTTCGTCCGTGAAATGGTACTTAACCTGTTCTTTATTTTCCTCGTACTGGTTGGTGTGGGGATTTGGATGCAGGTCAGTGGTGGTGATTCGAAAGAAACGGCCAGTCGTGGCGCACTGCTGCTGGACATTTCTGgt                                                                                                                                                      <  5:164376/148‑1 (MQ=255)
       ctgAATTTCGTCCGTGAAATGGTACTTAACCTGTTCTTTATTTTCCTCGTACTGGTTGGTGTGGGGATTTGGATGCAGGTCAGTGGTGGTGATTCGAAAGAAACGGCCAGTCGTGGCGCACTGCTGCTGGACATTTCTGGTGTGAtc                                                                                                                                                <  6:6865/147‑1 (MQ=255)
        tgAATTTCGTCCGTGAAATGGTACTTAACCTGTTCTTTATTTTCCTCGTACTGGTTGGTGTGGGGATTTGGATGCAGGTCAGTGGTGGTGATTCGAAAGAAACGGCCAGTCGTGGCGCACTGCTGCTGGACATTTCTGGTGTGAtcgtc                                                                                                                                             >  7:190559/1‑149 (MQ=255)
        tgAATTTCGTCCGTGAAATGGTACTTAACCTGTTCTTTATTTTCCTCGTACTGGTTGGTGTGGGGATTTGGATGCAGGTCAGTGGTGGTGATTCGAAAGAAACGGCCAGTCGTGGCGCACTGCTGCTGGACATTTCTGGTGTGAtcgtc                                                                                                                                             >  3:289327/1‑149 (MQ=255)
        tgAATTTCGTCCGTGAAATGGTACTTAACCTGTTCTTTATTTTCCTCGTACTGGTTGGTGTGGGGATTTGGATGCAGGTCAGTGGTGGTGATTCGAAAGAAACGGCCAGTCGTGGCGCACTGCTGCTGGACATTTCTGGTGTGAtcgtc                                                                                                                                             >  2:269520/1‑149 (MQ=255)
           aTTTCGTCCGTGAAATGGTACTTAACCTGTTCTTTATTTTCCTCGTACTGGTTGGTGTGGGGATTTGGATGCAGGTCAGTGGTGGTGATTCGAAAGAAACGGCCAGTCGTGGCGCACTGCTGCTGGACATTTCTGGTGTGATCGTCGAt                                                                                                                                          <  4:421718/149‑1 (MQ=255)
                   agTGAATTGGTACTTTATCTGTTGTTTATTTTCCTCGTACTGGTTGGTGTGGGGATTTGGATGCTGGTGTGTTTTGGTGATTGGAATGAAAGGGCCAGTGGTGGGGCACTGCTGATGGACATTTCTGGTGTGTTCGTCGATAAACTCGa                                                                                                                                  <  4:190533/148‑1 (MQ=255)
                    gTGAAATGGTACTTAACCTGTTCTTTATTTTCCTCGTACTGGTTGGTGTGGGGATTTGGATGCAGGTCAGTGGTGGTGATTCGAAAGAAACGGCCAGTCGTGGCGCACTGCTGCTGGACATTTCTGGTGTGATCGTCGAt                                                                                                                                          >  7:168443/1‑140 (MQ=255)
                    gTGAAATGGTACTTAACCTGTTCTTTATTTTCCTCGTACTGGTTGGTGTGGGGATTTGGATGCAGGTCAGTGGTGGTGATTCGAAAGAAACGGCCAGTCGTGGCGCACTGCTGCTGGACATTTCTGGTGTGATCGTCGAt                                                                                                                                          <  8:168443/140‑1 (MQ=255)
                         aTGGTACTTAACCTGTTCTTTATTTTCCTCGTACTGGTTGGTGTGGGGATTTGGATGCAGGTCAGTGGTGGTGATTCGAAAGAAACGGCCAGTCGTGGCGCACTGCTGCTGGACATTTCTGGTGTGATCGTCGATAAACCCGACAGTtc                                                                                                                            >  8:122930/1‑149 (MQ=255)
                            gTACTTAACCTGTTCTTTATTTTCCTCGTACTGGTTGGTGTGGGGATTTGGATGCAGGTCAGTGGTGGTGATTCGAAAGAAACGGCCAGTCGTGGCGCACTGCTGCTGGACATTTCTGGTGTGATCGTCGATAAACCCGACAGTtctc                                                                                                                          <  3:307550/148‑1 (MQ=255)
                            gTACTTAACCTGTTCTTTATTTTCCTCGTACTGGTTGGTGTGGGGATTTGGATGCAGGTCAGTGGTGGTGATTCGAAAGAAACGGCCAGTCGTGGCGCACTGCTGCTGGACATTTCTGGTGTGATCGTCGATAAACCCGACAGTTCTCa                                                                                                                         <  7:340010/149‑1 (MQ=255)
                                    ccTGTTCTTTATTTTCCTCGTACTGGTTGGTGTGGGGATTTGGATGCAGGTCAGTGGTGGTGATTCGAAAGAAACGGCCAGTCGTGGCGCACTGCTGCTGGACATTTCTGGTGTGATCGTCGATAAACCCGACAGTTCTCAGCGGTTTa                                                                                                                 <  4:289327/149‑1 (MQ=255)
                                     cTGTTCTTTATTTTCCTCGTACTGGTTGGTGTGGGGATTTGGATGCAGGTCAGTGGTGGTGATTCGAAAGAAACGGCCAGTCGTGGCGCACTGCTGCTGGACATTTCTGGTGTGATCGTCGATAAACCCGACAGTTCTCAGCGGTTTAg                                                                                                                <  7:48505/149‑1 (MQ=255)
                                          cTTTATTTTCCTGGTAGTGGTTGGTGTGGGGATTTGGATGCAGGTCAGTGGTGGTGAGTCGAAAGAAACGGCCAGTCGTGGCGCACTGCTGCTTGACATTTCTGGTGTGATCGTCGATAAACCCGACAGTTCTCAGCGGTTTAGTaaa                                                                                                            <  6:402576/148‑1 (MQ=255)
                                          cTTTATTTTCCTCGTACTGGTTGGTGTGGGGATTTGGATGCAGGTCAGTGGTGGTGATTCGAAAGAAACGGCCAGTCGTGGCGCACTGCTGCTGGACATTTCTGGTGTGATCGTCGATAAACCCGACAGTTCTCAGCGGTTTAGTaaa                                                                                                            <  7:151166/148‑1 (MQ=255)
                                           tttATTTTCCTCGTACTGGTTGGTGTGGGGATTTGGATGCAGGTCAGTGGTGGTGATTCGAAAGAAACGGCCAGTCGTGGCGCACTGCTGCTGGGCATTTCTGGTGTGATCGTCGATAAACCCGACAGTTCTCAGCGGTTTAGTAAAt                                                                                                           >  4:123646/1‑148 (MQ=255)
                                           tttATTTTCCTCGTACTGGTTGGTGTGGGGATTTGGATGCAGGTCAGTGGTGGTGATTCGAAAGAAACGGCCAGTCGTGGCGCACTGCTGCTGGACATTTCTGGTGTGATCGTCGATAAACCCGACAGTTCTCAGCGGTTTAGTAAAt                                                                                                           >  5:242564/1‑148 (MQ=255)
                                           tttATTTTCCTCGTACTGGTTGGTGTGGGGATTTGGATGCAGGTCAGTGGTGGTGATTCGAAAGAAACGGCCAGTCGTGGCGCACTGCTGCTGGACATTTCTGGTGTGATCGTCGATAAACCCGACAGTTCTCAGCGGTTTAGTAAAt                                                                                                           >  1:282262/1‑148 (MQ=255)
                                           tttATTTTCCTCGTACTGGTTGGTGTGGGGATTTGGATGCAGGTCAGTGGTGGTGATTCGAAAGAAACGGCCAGTCGTGGCGCACTGCTGCTGGACATTTCTGGTGTGATCGTCGATAAACCCGACAGTTCTCAGCGGTTTAGTAAAt                                                                                                           >  8:231897/1‑148 (MQ=255)
                                           tttATTTTCCTCGTACTGGTTGGTGTGGGGATTTGGATGCAGGTCAGTGGTGGTGATACGAAAGAAACGGCCAGTCGTGGCGCACTGCTGCTGGAGATTTCTGGTGTGATCGTCGATAAACACGACAGTTCTCAGCGGTTTAGTAAAt                                                                                                           >  5:344491/1‑148 (MQ=255)
                                              aTTTTCCTCGTACTGGTTGGGGTGGGGATTTGGATGCAGGTCAGTGGTGGTGATTCGAAAGAAACGGCCAGTCGTGGCGCACTGCTGCTGGACATTTCTGGTGTGATCGTCGATAAACCCGACAGTTCTCCGCGGTTTAGTAAATTaa                                                                                                        <  1:269520/148‑1 (MQ=255)
                                                     tCGTACTGGTTGGTGTGGGGATTTGGATGCAGGTCAGTGGTGGTGATTCGAAAGAAACGGCCAGTCGTGGCGCACTGCTGCTGGACATTTCTGGTGTGATCGTCGATAAACCCGACAGTTCTCAGCGGTTTAGTAAATTAAGCCGCCAg                                                                                                >  2:319331/1‑149 (MQ=255)
                                                     tCGTACTGGTTGGTGTGGGGATTTGGATGCAGGTCAGTGGTGGTGATTCGAAAGAAACGGCCAGTCGTGGCGCACTGCTGCTGGACATTTCTGGTGTGATCGTCGATAAACCCGACAGTTCTCAGCGGTTTAGTAAATTAAGCCGCCAg                                                                                                >  5:335033/1‑149 (MQ=255)
                                                       gTACTGGTTGGTGTGGGGATTTGTATGCAGGTCAGTGGTGGTGATTCGAAAGAAACGGCCAGTCGCGGCGCACTGCTGCTGGACATTTCTGGTGTGATCGTCGATAAACCCGACAGTTCTCAGCGGTTTa                                                                                                                 >  4:316321/1‑130 (MQ=255)
                                                       gTACTGGTTGGTGTGGGGATTTGGATGCAGGTCAGTGGTGGTGATTCGAAAGAAACGGCCAGTCGCGGCGCACTGCTGCTGGACATTTCTGGTGTGATCGTCGATAAACCCGACAGTTCTCAGCGGTTTa                                                                                                                 <  3:316321/130‑1 (MQ=255)
                                                                   gtgGGGATTTGGATGCAGGGCAGGGGTGGGGATTCGAATGAAACGGCCAGTCGGGGCGCACTGCTGCTGGAAATTTCTTGTGTTATAGTCGATAAACACGACAGTTCTCAGCGGTTTAGTATATTAAGCCGCCAGCTGCTTGGTGCCAg                                                                                  <  8:419329/149‑1 (MQ=255)
                                                                    tgGGGATTTGGATGCAGGTCAGTGGTGGTGATTCGAAAGAAACGGCCAGTCGTGGCGCACTGCTGCTGGACATTTCTGGTGTGATCGTCGATAAACCCGACAGTTCTCAGCGGTTTAGTAAATTAAGCCGCCAGCTGCTTGGTGCCAGt                                                                                 >  4:171444/1‑149 (MQ=255)
                                                                      gggATTTGGATGCAGGTCAGTGGTGGTGATTCGAAAGAAACGGCCAGTCGTGGCGCACTGCTGCTGGACATTTCTGGTGTGATCGTCGATAAACCCGACAGTTCTCAGCGGTTTAGTAAATTAAGCCGCCAGCTGCTTGGTGCCAGTTc                                                                               <  2:228093/149‑1 (MQ=255)
                                                                      gggATTTGGATGCAGGTCAGTGGTGGTGATTCGAAAGAAACGGCCAGTCGTGGCGCACTGCTGCTGGACATTTCTGGTGTGATCGTCGATAAACCCGACAGTTCTCAGCGGTTTAGTAAATTAAGCCGCCAGCTGCTTGGTGCCAGTTc                                                                               >  1:228093/1‑149 (MQ=255)
                                                                        gATTTGGATGCAGGTCAGTGGTGGTGATTCGAAAGAAACGGCCAGTCGTGGCGCACTGCTGCTGGACATTTCTGGTGTGATCGTCGATAAACCCGACAGTTCTCAGCGGTTTa                                                                                                                 >  1:151519/1‑113 (MQ=255)
                                                                        gATTTGGATGCAGGTCAGTGGTGGTGATTCGAAAGAAACGGCCAGTCGTGGCGCACTGCTGCTGGACATTTCTGGTGTGATCGTCGATAAACCCGACAGTTCTCAGCGGTTTa                                                                                                                 <  2:151519/113‑1 (MQ=255)
                                                                        gATTTGGATGCAGGTCAGTGGTGGTGATTCGAAAGAAACGGCCAGTCGTGGCGCACTGCTGCTGGACATTTCTGGTGTGATCGTCGATAAACCCGACAGTTCTCAGCGGTTTAGTAAATTAAGCCGCCAGCTGCTTGGTGCCAGTTcc                                                                              >  3:404842/1‑148 (MQ=255)
                                                                              gATGCAGGTCAGTGGTGGTGATTCGAAAGAAACGGCCAGTCGTGGCGCACTGCTGCTGGTCATTTCTGGTGTGATCGTCGATAAACCCGACAGTTCTCAGCGGTTTAGTAAATTAAGCCGCCAGCTGCTTGGTGCCAGTTCCGATCGTc                                                                       <  3:123646/149‑1 (MQ=255)
                                                                              gATGCAGGTCAGTGGTGGTGATTCGAAAGAAACGGCCAGTCGTGGCGCACTGCTGCTGGACATTTCTGGTGTGATCGTCGATAAACCCGACAGTTCTCAGCGGTTTAGTAAATTAAGCCGCCAGCTGCTTGGTGCCAGTTCCGATCGTc                                                                       <  7:231897/149‑1 (MQ=255)
                                                                                 gCAGGTCAGTGGTGGTGATTCGAAAGAAACGGCCAGTCGTGGCGCACTGCTGCTGGACATTTCTGGTGTGATCGTCGATAAACCCGACAGTTCTCAGCGGTTTAGTAAATTAAGCCGCCAGCTGCTTGGTGCCAGTTCCGATCGTCTg                                                                     >  8:269839/1‑148 (MQ=255)
                                                                                  cAGGTCAGTGGTGGTGATTCGAAAGAAACGGCCAGTCGTGGCGCACTGCTGCTGGACATTTCTGGTGTGATCGTCGATAAACCCGACAGTTCTCAGCGGTTTAGTAAATTAAGCCGCCAGCTGCTTGGTGCCAGTTCCGATCGTCTGCa                                                                   <  1:319331/149‑1 (MQ=255)
                                                                                   aGGTCAGTGGTGGTGATTCGAAAGAAACGGCCAGTCGTGGCGCACTGCTGCTGGACATTTCTGGTGTGATCGTCGATAAACCCGACAGTTCTCAGCGGTTTAGTAAATTAAGCCGCCAGCTGCTTGGTGCCAGTTCCGATCGTCTGCAg                                                                  >  6:328343/1‑149 (MQ=255)
                                                                                     gTCAGTGGTGGTGATTCGAAAGAAACGGCCAGTCGTGGCGCACTGCTGCTGGACATTTCTGGTGTGATCGTCGATAAACCCGACAGTTCTCAGCGGTTTAGTAAATTAAGCCGCCAGCTGCTTGGTGCCAGTTCCGATCGTCTGCAgg                                                                 >  3:374907/1‑148 (MQ=255)
                                                                                      tCAGTGGTGGTGATTCGAAAGAAACTGCCAGTCGTGGCGCACTTCTGCTGGACATTTCTGGTGTGATCGTCGATAAACCCGACAGTTCTCAGCGGTTTAGTaaa                                                                                                            <  5:171627/104‑1 (MQ=255)
                                                                                      tCAGTGGTGGTGATTCGAAAGAAACGGCCAGTCGTGGCGCACTGCTGCTGGACATTTATGGAGTGATCGTCGATAAACCCGACAGTTCTCATCGGTTTAGTaaa                                                                                                            >  6:171627/1‑104 (MQ=255)
                                                                                            gtggtgTTTCGAAAGAAACGGCCTCGCGGGGCGCACTGCTGGTGGACATTTCTGGTGGGCTCGTCGATTAACCCGTCATTTCTCATCGGTTTAGTAAATTAAGCCGCCAGCTGCTTTGTGCCAGGTCCGATAGTCTGGAGGAAAACTCa                                                         <  4:404842/149‑1 (MQ=255)
                                                                                            gtggtgATTCGAAAGAAACGGCCAGTCGTGGCGCACTGCTGCTGGACATTTCTGGTGTGATCGTCGATAAACCCGACAGTTCTCAGCGGTTTAGTAAATTAAGCCGCCAGCTGCTTGGTGCCAGTTCCGATCGTCTGCAGGAAAACTCa                                                         <  2:282262/149‑1 (MQ=255)
                                                                                            gtggtgATTCGAAAGAAACGGCCAGTCGTGGCGCACTGCTGCTGGACATTTCTGGTGTGATCGTCGATAAACCCGACAGTTCTCAGCGGTTTAGTAAATTAAGCCGCCAGCTGCTTGGTGCCAGTTCCGATCGTCTGCAGGAAAACTCa                                                         <  7:122930/149‑1 (MQ=255)
                                                                                            gtggtgATTCGAAAGAAACGGCCAGTCGTGGCGCACTGCTGCTGGACATTTCTGGCGTGATCGTCGATAAACCCGACAGTTCTCAGCGGTTTAGTAAATTAAGCCGCCAGCTGCTTGGTGCCAGTTCCGATCGTCTGCAGGAAAACTCa                                                         >  5:137372/1‑149 (MQ=255)
                                                                                                  aTTCGAAAGAAACGGCCAGTCGTGGCGCACTGCTGCTGGACATTTCTGGTGTGATCGTCGATAAACCCGACAGTTCTCAGCGGTTTAGTAAATTAAGCCGCCAGCTGCTTGGTGCCAGTTCCGATCGTCTGCAGGAAAACTCACTGttt                                                   <  4:374907/149‑1 (MQ=255)
                                                                                                  aTTCGAAAGAAACGGCCAGTCGTGGCGCACTGCTGCTGGACATTTCTGGTGTGATCGTCGATAAACCCGACAGTTCTCAGCGGTTTAGTAAATTAAGCCGCCAGCTGCTTGGTGCCAGTTCCGATCGTCTGCAGGAAAACTCACTGttt                                                   <  7:269839/149‑1 (MQ=255)
                                                                                                  aTTCGAAAGAAACGGCCAGTCGTGGCGCACTGCTGCTGGACATTTCTGGCGTGATCGTCGATAAACCCGACAGTTCTCAGCGGTTTAGTAAATTAAGCCGCCAGCTGCTTGGTGCCAGTTCCGATCGTCTGCAGGAAAACTCACTGtt                                                    <  6:137372/148‑1 (MQ=255)
                                                                                                    tCGAAAGAAACGGCCAGTCGTGGCGCACTGCTGCTGGACATTTCTGGTGTGATCGTCGATAAACCCGACAGTTCTCAGCGGTTTAGTAAATTAAGCCGCCAGCTGCTTGGTGCCAGTTCCGATCGTCTGCAGGAAAACTCACTGTTTGa                                                 >  1:193109/1‑149 (MQ=255)
                                                                                                      gaaagaaaCGGCCAGTCGTGGCGCACTGCTGCTGGACATTTCTGGTGTGATCGTCGATAACCCCGACAGTTCTCAGCGGTTTAGTAAATTAAGCCGCCAGCTGCTTGGTGCCAGTTCCGATCGTCTGCAGGAAAACTCACTGTTTGata                                               >  5:219280/1‑149 (MQ=255)
                                                                                                      gaaagaaaCGGCCAGTCGTGGCGCACTGCTGCTGGACATTTCTGGTGTGATCGTCGATAAACCCGACAGTTCTCAGCGGTTTAGTAAATTAAGCCGCCAGCTGCTTGGTGCCAGTTCCGATCGt                                                                        >  4:20033/1‑124 (MQ=255)
                                                                                                      gaaagaaaCGGCCAGTCGTGGCGCACTGCTGCTGGACATTTCTGGTGTGATCGTCGATAAACCCGACAGTTCTCAGCGGTTTAGTAAATTAAGCCGCCAGCTGCTTGGTGCCAGTTCCGATCGt                                                                        <  3:20033/124‑1 (MQ=255)
                                                                                                      gaaagaaaCGGCCAGTCGTGGCGCACTGCTGCTGGACATTTCTGGTGTGATCGTCGATAAACCCGACAGTTCTCAGCGGTTTAGTAAATTAAGCCGCCAGCTGCTTGGTGCCAGTTCCGATCGTCTGCAGGAAAACTCACTGTTTGata                                               >  8:119697/1‑149 (MQ=255)
                                                                                                      gaaagaaaCGGCCAGTCGTGGCGCACTGCTGCTGGACAATTCTGGTGTGATCGTCGATAAACCCGACAGTTCTCCGCGGTTAAGTAAATTAAGCCGCCAGCTGCGTGGTGCCAGATCCGCTCGGCTGCAGGAAAACTCACTCTTTCata                                               >  5:19029/1‑149 (MQ=255)
                                                                                                        aagaaaCGGCCAGTCGTGGCGCACTGCTGCTGGACATTTCTGGTGTGATCGTCGATAAACCCGACAGTTCTCAGCGGTTTAGTAAACTAAGCCGCCAGCTGCTTGGTGCCAGTTCCGATCGTCTGCAGGAAAACTCACTGTTTGATATc                                             <  7:333689/149‑1 (MQ=255)
                                                                                                          gaaaCGGCCAGTCGTGGCGCACTGCTGCTGGACATTTCTGGTGTGATCGTCGATAAACCCGACAGTTCTCAGCGGTTTAGTAAATTAAGCCGCCAGCTGCTTGGTGCCAGTTCCGATCGTCTGCAGGAAAACTCACTGTTTGATATCGt                                           >  8:128765/1‑149 (MQ=255)
                                                                                                           aaaCGGCCAGTCGTGGCGCACTGCTGCTGGACATTTCTGGTGTGATCGTCGATAAACCCGACAGTTCTCAGCGGTTTAGTAAATTAAGCCGCCAGCTGCTTGGTGCCAGTTCCGATCGTCTGCAGGAAAACTCACTGTTTGATATCGt                                           >  1:26435/1‑148 (MQ=255)
                                                                                                           aaaCGGCCAGTCGTGGCGCACTGCTGCTGGACATTTCTGGTGTGATCGTCGATAAACCCGACAGTTCTCAGCGGTTTAGTAAATTAAGCCGCCAGCTGCTTGGTGCCAGTTCCGATCGTCTGCAGGAAAACTCACTGTTTGATATCGTc                                          <  8:190559/149‑1 (MQ=255)
                                                                                                                  cAGTCGTGGCGCACTGCTGCTGGACATTTCTGGTGTGATCGTCGATAAACCCGACAGTTCTCAGCGGTTTAGTAAATTAAGCCGCCAGCTGCTTGGTGCCAGTTCCGATCGTCTGCAGGAAAACTCACTGTTTGATATCGTCAACACTa                                   >  3:163896/1‑149 (MQ=255)
                                                                                                                       gTGGCGCACTGCTGCTGGACATTTCTGGTGTGATCGTCGATAAACCCGACAGTTCTCAGCGGTTTAGTAAATTAAGCCGCCAGCTGCTTGGTGCCAGTTCCGATCGTCTGCAGGAAAACTCACTGTTTGATATCGTCAACACTATTCGc                              >  3:143329/1‑149 (MQ=255)
                                                                                                                       gTGGCGCACTGCTGCTGGACATTTCTGGCGTGATCGTCGATAAACCCGACAGTTCTCAGCGGt                                                                                                                    <  8:69246/63‑1 (MQ=255)
                                                                                                                       gTGGCGCACTGCTGCTGGACATTTCTGGCGTGATCGTCGATAAACCCGACAGTTCTCAGCGGt                                                                                                                    >  7:69246/1‑63 (MQ=255)
                                                                                                                            gcACTGCTGCTGGACATTTCTGGTGTGATCGTCGATAAACCCGACAGTTCTCAGCGGTTTAGTAAATTAAGCCGCCAGCTGCTTGGTGCCAGTTCCGATCGTCTGCAGGAAAACTCACTGTTTGATATCGTCAACACTATTCGCCAGGc                         >  2:356011/1‑149 (MQ=255)
                                                                                                                                  ctgctgGACATTTCTGGTGTGATCGTCGATAAACCCGACAGTTCTCAGCGGTTTAGTAAATTAAGCCGCCAGCTGCTTGGTGCCAGTTCCGATCGTCTGCAGGAAAACTCACTGTTTGATATCGTCAACACTATTCGCCAGGCGAAgg                    <  7:128765/148‑1 (MQ=255)
                                                                                                                                      tgGACATTTCTGGTGTGATCGTCGATAAACCCGACAGTTCTCAGCGGTTTAGTa                                                                                                              <  7:433901/54‑1 (MQ=255)
                                                                                                                                      tgGACATTTCTGGTGTGATCGTCGATAAACCCGACAGTTCTCAGCGGTTTAGTa                                                                                                              >  8:433901/1‑54 (MQ=255)
                                                                                                                                        gACATTTCTGGTGTGATCGTCGATAAACCCGACAGTTCTCAGCGGTTTAGTAAATTAAGCCGCCAGCTGCTTGGTGCCAGTTCCGATCGTCTGCAGGAAAACTCACTGTTTGATATCGTCAACACTATTCGCCAGGCGAAGGACGAcc              <  3:171444/148‑1 (MQ=255)
                                                                                                                                        gACATTTCTGGTGTGATCGTCGATAAACCCGACAGTTCTCAGCGGTTTAGTAAATTAAGCCGCCAGCTGCTTGGTGCCAGTTCCGATCGTCTGCAGGAAAACTCACTGTTTGATATCGTCAACACTATTCGCCAGGCGAAGGACGACCg             <  7:119697/149‑1 (MQ=255)
                                                                                                                                        gACATTTCTGGTGTGATCGTCGATAAACCCGACAGTTCTCAGCGGTTTAGTAAATTAAGCCGCCAGCTGCTTGGTGCCAGTTCCGATCGTCTGCAGGAAAACTCACTGTTTGATATCGTCAACACTATTCGCCAGGCGAAGGACGACCg             >  1:277117/1‑149 (MQ=255)
                                                                                                                                        gACATTTCTGGTGTGATCGTCGATAAACCCGACAGTACTCAGAGGTTTAGTAAATTAAGCCGCCAGCTGCTTGGTGCCAGTTCCGATCGTCTGCAGGAAAACTCACTGTTTGATATCGTCAACACTATTCGCCAGGCGAAGGACGACCg             >  7:314704/1‑149 (MQ=255)
                                                                                                                                           aTTTCTGGTGTGATCGTCGATAAACCCGACAGTTCTCAGCGGTTTAGTAAATTAAGCCGCCAGCTGCTTGGTGCCAGTTCCGATCGTCTGCAGGAAAACTCACTGTTTGATATCGTCAACACTATTCGCCAGGCGAAGGACGACCGCaa          >  7:156371/1‑149 (MQ=255)
                                                                                                                                               cTGGTGTGATCGTCGATAAACCCGACAGTTCTCAGCGGTTTAGTAAATTAAGCCGCCAGCTGCTTGGTGCCAGTTCCGATCGTCTGCAGGAAAACTCACTGTTTGATATCGTCAACACTATTCGCCAGGCGAAGGACGACCGc            <  1:256281/143‑1 (MQ=255)
                                                                                                                                               cTGGTGTGATCGTCGATAAACCCGACAGTTCTCAGCGGTTTAGTAAATTAAGCCGCCAGCTGCTTGGTGCCAGTTCCGATCGTCTGCAGGAAAACTCACTGTTTGATATCGTCAACACTATTCGCCAGGCGAAGGACGACCGc            >  2:256281/1‑143 (MQ=255)
                                                                                                                                                 ggTGTGATCGTCGATAAACCCGACAGTTCTCAGCGGTTTAGTAAATTAAGCCGCCAGCTGCTTGGTGCCAGTTCCGATCGTCTGCAGGAAAACTCACTGTTTGATATCGTCAACACTATTCGCCAGGCGAAGGACGACCGCAATATCAc    >  3:352695/1‑149 (MQ=255)
                                                                                                                                                  gtgtGATCGTCGATAAACCCGACAGTTCTCAGCGGTTTAGTAAATTAAGCCGCCAGCTGCTTGGTGCCAGTTCCGATCGTCTGCAGGAAAACTCACTGTTTGATATCGTCAACACTATTCGCCAGGCGAAGGACGACCGCAATATc      <  1:356011/146‑1 (MQ=255)
                                                                                                                                                  gtgtGATCGTCGATAAACCCGACAGTTCTCAGCGGTTTAGTAAATTAAGCCGCCAGCTGCTTGGTGCCAGTTCCGATCGTCTGCAGGAAAACTCACTGTTTGATATCGTCAACACTATTCGCCAGGCGAAGGACGACCGCAATATCAcc   >  4:55326/1‑149 (MQ=255)
                                                                                                                                                   tgtgATCGTCGATAAACCCGACAGTTCTCAGCGGTTTAGTAAATTAAGCCGCCAGCTGCTTGGTGCCAGTTCCGATCGTCTGCAGGAAAACTCACTGTTTGATATCGTCAACACTATTCGCCAGGCGAAGGACGACCGCAATATCACCg  >  2:235902/1‑149 (MQ=255)
                                                                                                                                                   tgtgATCGTCGATAAACCCGACAGTTCTCAGCGGTTTAGTAAATTAAGCCGCCAGCTGCTTGGTGCCAGTTCCGATCGTCTGCAGGAAAACTCACTGTTTGATATCGTCAACACTATTCGCCAGGCGAAGGACGACCGCAATATCACCg  <  3:55326/149‑1 (MQ=255)
                                                                                                                                                   tgtgATCGTCGATAAACCCGACAGTTCTCAGCGGTTTAGTAAATTAAGCCGCCAGCTGCTTGGTGCCAGTTCCGATCGTCTGCAGGAAAACTCACTGTTTGATATCGTCAACACTATTCGCCAGGCGAAGGACGACCGCAATATCACCg  >  7:165766/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GCGTCTGCTGAATTTCGTCCGTGAAATGGTACTTAACCTGTTCTTTATTTTCCTCGTACTGGTTGGTGTGGGGATTTGGATGCAGGTCAGTGGTGGTGATTCGAAAGAAACGGCCAGTCGTGGCGCACTGCTGCTGGACATTTCTGGTGTGATCGTCGATAAACCCGACAGTTCTCAGCGGTTTAGTAAATTAAGCCGCCAGCTGCTTGGTGCCAGTTCCGATCGTCTGCAGGAAAACTCACTGTTTGATATCGTCAACACTATTCGCCAGGCGAAGGACGACCGCAATATCACCG  >  NC_000913/1848884‑1849179

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: