Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 334,215 G→T 5.1% T73T (ACC→ACA yahB ← putative DNA‑bindng transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913334,2150GT 5.1% 272.0 / 3.0 79T73T (ACC→ACAyahBputative DNA‑bindng transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (37/38);  new base T (2/2);  total (39/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.01e-01

CTGCTTTTTCAGCACCTGTAACAAGCGTGCGGTGTGTTTGGGCGTATAAAGCAGCTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGTTAACTCAATGCTGCGGGTACTGCGAGTGAACAGCACCACATCCAGCCCCGTCTCCATACGCTTAATGGTGTAGCTAATGGCGGAAGTGGTTAAACCCAACTCTTCT  >  NC_000913/334071‑334356
                                                                                                                                                |                                                                                                                                             
cTGCTTTTTCAGCACCTGTAACAAGCGTGCGGTGTGTTTGGGCGTATAAAGCAGCTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAg                                                                                                                                           >  3:330904/1‑149 (MQ=255)
 tGCTTTTTCAGCCCCTGTAACAAGCGTGCGGTGTGTTTGGGCGTATAAAGCAGCTGATTAATACAAATGCGCACCCGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGc                                                                                                                                          >  5:260509/1‑149 (MQ=255)
   cTTTTTCAGCACCTGTAACAAGCGTGCGGTGTGTTTGGGCGTATAAAGCAGCTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCt                                                                                                                                         <  5:12852/148‑1 (MQ=255)
    tttttCAGCACCTGTAACAAGCGTGCGGTGTGTTTGGGCGTATAAAGCAGCTGATTAATACAAATGAGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCttt                                                                                                                                       >  3:319217/1‑149 (MQ=255)
      tttCAGCACCTGTAACAAGCGTGCGGTGTGTTTGGGCGTATAAAGCAGCTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCttt                                                                                                                                       <  8:290234/147‑1 (MQ=255)
       ttCAGCACCTGTAACAAGCGTGCGGTGTGTTTGGGCGTATAAAGCAGCTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGc                                                                                                                                     <  3:148449/148‑1 (MQ=255)
        tCAGCACCTGTAACAAGCGTGCGGTGTGTTTGGGCGTATAAAGCAGCTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCg                                                                                                                                    >  6:347093/1‑148 (MQ=255)
              ccTGTAACAAGCGTGCGGTGTGTTTGGGCGTATAAAGCAGCTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGaaaaaa                                                                                                                             >  3:187803/1‑149 (MQ=255)
               cTGTAACAAGCGTGCGGTGTGTTTGGGCGTATAAAGCAGCTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGaaaaaa                                                                                                                             >  3:137167/1‑148 (MQ=255)
                   aaCAAGCGTGCGGTGTGTTTGGGCGTATAAAGCAGCTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCg                                                                                                                                               <  7:344980/126‑1 (MQ=255)
                   aaCAAGCGTGCGGTGTGTTTGGGCGTATAAAGCAGCTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCg                                                                                                                                               >  8:344980/1‑126 (MQ=255)
                       aGCGTGCGGTGTGTTTGGGCGTATAAAGCAGCTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGCCCg                                                                                                                    <  8:102287/149‑1 (MQ=255)
                       aGCGTGCGGTGTGTTTGGGCGTATAAAGCAGCTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGCCCg                                                                                                                    <  8:241125/149‑1 (MQ=255)
                         cGTGCGGTGTGTTTGGGCGTATAAAGCAGCTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGa                                                                                                                  <  4:187803/149‑1 (MQ=255)
                               gtgtgtTTGGGCGTATAAAGCAGCTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGt                                                                                                            <  3:325039/149‑1 (MQ=255)
                               gtgtgtTGGGGCGGATCAAGCAGCTGCTTCATACAAAGGCGCACACGCGCTTCTATGCCCGGCGAACTAGTCTCGATTCTGCGTGTGATGGCATCACAATCATTCCGCCGGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGt                                                                                                            <  3:249235/149‑1 (MQ=255)
                                tgtgtTTGGGCGTATAAAGCAGCTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGt                                                                                                            >  1:379484/1‑148 (MQ=255)
                                 gtgtTTGGGCGTATAAAGCAGCTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGtt                                                                                                           >  8:207853/1‑148 (MQ=255)
                                    tttGGGCGTATAAAGCAGCTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGTTAAc                                                                                                        <  3:225746/148‑1 (MQ=255)
                                        ggCGTATAAAGCAGCTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGGCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGTTAACTCAAt                                                                                                   <  2:58213/149‑1 (MQ=255)
                                         gCGTATAAAGCAGCTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGTTAACTCAAtg                                                                                                  >  8:376064/1‑149 (MQ=255)
                                          cGTATAAAGCAGCTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGTTAACTCAAtgc                                                                                                 <  4:137167/149‑1 (MQ=255)
                                              taAAGCAGCTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGTTAACTCAATGCTGCg                                                                                             >  5:121080/1‑149 (MQ=255)
                                                aaGCAGCTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGTTAACTCAATGCTGCggg                                                                                           <  4:330904/149‑1 (MQ=255)
                                                  gcagcTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGTTAACTCAATGCTGCGGGTa                                                                                         <  5:49033/149‑1 (MQ=255)
                                                      cTGATTAATCCAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGTTAACTCAATGCTGCGGGTACTGc                                                                                     <  2:122282/149‑1 (MQ=255)
                                                      cTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTa                                                                                                                                            >  7:343097/1‑94 (MQ=255)
                                                      cTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTa                                                                                                                                            <  8:343097/94‑1 (MQ=255)
                                                      cTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGTTAACTCAATGCTGCGGGTACTGc                                                                                     >  3:150806/1‑149 (MQ=255)
                                                      cTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGTTAACTCAATGCTGCGGGTACTGc                                                                                     <  3:73124/149‑1 (MQ=255)
                                                          ttAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGTTAACTCAATGCTGCGGGTACTGCGAGt                                                                                 >  5:339752/1‑149 (MQ=255)
                                                           tAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTg                                                                                                                                      <  8:383089/95‑1 (MQ=255)
                                                           tAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTg                                                                                                                                      >  7:383089/1‑95 (MQ=255)
                                                             aTACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCATGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGTTAACTCAATGATTCGGGTACTGCGAGTTaa                                                                              >  3:72772/1‑149 (MQ=255)
                                                             aTACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGTTAACTCAATGCTGCGGGTACTGCGAGTGaa                                                                              >  6:208995/1‑149 (MQ=255)
                                                             aTACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGTTAACTCAATGCTGCGGGTACTGCGAGTGaa                                                                              >  2:229074/1‑149 (MQ=255)
                                                              tACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGTTAACTCAATGCTGCGGGTACTGCGAGTGaa                                                                              >  5:193478/1‑148 (MQ=255)
                                                              tACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGTTAACTCAATGCTGCGGGTACTGCGAGTGAAc                                                                             <  1:229074/149‑1 (MQ=255)
                                                                 aaaTGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGTTAACTCAATGCTGCgg                                                                                            <  5:247789/131‑1 (MQ=255)
                                                                 aaaTGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGTTAACTCAATGCTGCGGGTACTGCGAGTGAACAg                                                                           >  6:28962/1‑148 (MQ=255)
                                                                 aaaTGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTATGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGTTAACTCAATGCTGCgg                                                                                            >  6:247789/1‑131 (MQ=255)
                                                                  aaTGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGGGTTTGATGGCATAAAAATCATTCTGCAGGTCGGTAGCTTTGCGGAAAAAATTACGCCCGGACTCCGTTATCTCAATGCTGCGGGTACTGCGAGTGTACAGca                                                                         <  4:319217/149‑1 (MQ=255)
                                                                    tGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCCGCAGGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGTTAACTCAATGCTGCGGGTACTGCGAGTGAACCGccc                                                                        >  2:408958/1‑146 (MQ=255)
                                                                          cacGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGTTAACTCAATGCTGCGGGTACTGCGAGt                                                                                 >  6:341236/1‑133 (MQ=255)
                                                                          cacGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGTTAACTCAATGCTGCGGGTACTGCGAGt                                                                                 <  5:341236/133‑1 (MQ=255)
                                                                               gcTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGTTAACTCAATGCTGCGGGTACTGCGAGTGAACAGCACCACATCCAGccc                                                            <  6:121080/149‑1 (MQ=255)
                                                                                cTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGTTAACTCAATGCTGCGGGTACTGCGAGTGAACAGCACCACATCCAGcccc                                                           <  6:260509/149‑1 (MQ=255)
                                                                                         ccTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGTTAACTCAATGCTGCGGGTACTGCGAGTGAACAGCACCACATCCAGCCCCGTCTCCATa                                                  >  3:333415/1‑149 (MQ=255)
                                                                                          cTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCTGTAGCTTTGCGGAATAAATAACGGCCGGAGTGCGTTAACTCAATGCTGCGGGTACTGCGAGTGAACAGCACCACATCCAGCCCCGTCTCCATAc                                                 >  5:92624/1‑149 (MQ=255)
                                                                                          cTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCTGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGTTAACTCAATGCTGCGGGTACTGCGAGTGAACAGCACCACATCCAGCCCCGTCTCCATAc                                                 >  5:98308/1‑149 (MQ=255)
                                                                                             cGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGTTAACTCAATGCTGCGGGTACTGCGAGTGAACAGCACCACATCCAGCCCCGTCTCCATACGc                                               <  5:347093/148‑1 (MQ=255)
                                                                                              gAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGTTAACTCAATGCTGCGGGTACTGCGAGTGAACAGCACCACATCCAGCCCCGTCTCCATACGCt                                              >  2:211357/1‑148 (MQ=255)
                                                                                               aaaTAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGTTAACTCAATGCTGCGGGTACTGCGAGTGAACAGCACCACATCCAGCCCCGTCTCCATACGCTTa                                            <  5:28962/149‑1 (MQ=255)
                                                                                                aaTAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGTTAACTCAATGCTGCGGGTACTGCGAGTGAACAGCACCACATCCAGCCCCGTCTCCATACGCTTaa                                           <  4:150806/149‑1 (MQ=255)
                                                                                                     tATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGTTAACTCAATGCTGCGGGTACTGCGAGTGAACAGcacc                                                                       >  4:148426/1‑116 (MQ=255)
                                                                                                     tATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGTTAACTCAATGCTGCGGGTACTGCGAGTGAACAGcacc                                                                       <  3:148426/116‑1 (MQ=255)
                                                                                                     tATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGTTAACTCAATGCTGCGGGTACTGCGAGTGAACAGCACCACATCCAGCCCCGTCTCCATACGCTTAATGgtg                                      >  8:114835/1‑149 (MQ=255)
                                                                                                       tCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGTTAACTCAATGCTGCGGGTACTGCGAGTGAACAGCACCACATCCAGCCCCGTCTCCATACGCTTAATGgtgt                                     <  6:339752/148‑1 (MQ=255)
                                                                                                        cGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGTTAACTCAATGCTGCGGGTACTGCGAGTGAACAGCACCCCATCCAGCCCCGTCTCCATACGCTTAATGGTGTAg                                   >  7:327020/1‑149 (MQ=255)
                                                                                                           ttCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGTTAACTCAATGCTGCGGGTACTGCGAGTGAACAGCACCACATCCAGCCCCGTCTCCATACGCt                                              <  2:121582/135‑1 (MQ=255)
                                                                                                           ttCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGTTAACTCAATGCTGCGGGTACTGCGAGTGAACAGCACCACATCCAGCCCCGTCTCCATACGCt                                              >  1:121582/1‑135 (MQ=255)
                                                                                                                     gATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGTTAACTCAATGCTGCGGGTACTGCGAGTGAACAGCACCACATCCAGCCCCGTCTCCATACGCTTAATGGTGTAGCTAATGGCGGaa                       <  4:333415/148‑1 (MQ=255)
                                                                                                                     gATGGCATAAAAATCATTCAGCAGGTCGGGAGCTTTGCGGAAAAAATTACGCCCGGACGCCGTTAACTCAATGCTGCGGGTACTGAGAGTGAACATCACCACATCCAGCCCCGTCTCCATACGCTTAATGGTGTAGCTAATGGCggta                       <  4:72772/148‑3 (MQ=255)
                                                                                                                       tGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGTTAACTCAATGCTGCGGGTACTGCGAGTGAACAGCACCACATCCAGCCCCGTCTCCATACGCTTAATGGTGTAGCTAATGGCGGAAGTg                    >  1:169650/1‑149 (MQ=255)
                                                                                                                        ggCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGTTAACTCAATGCTGCGGGTACTGCGAGTGaa                                                                              <  1:86941/90‑1 (MQ=255)
                                                                                                                        ggCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGTTAACTCAATGCTGCGGGTACTGCGAGTGaa                                                                              >  2:86941/1‑90 (MQ=255)
                                                                                                                         gCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGTTAACTCAATGCTGCGGGTACTGCGAGTGAACAGCACCACATCCAGCCCCGTCTCCATACGCTTAATGGTGTAGCTAATGGCGGAAGTgg                   <  7:376064/148‑1 (MQ=255)
                                                                                                                           aTAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGTTAACTCAATGCTGCGGGTACTGCGAGTGAACAGCACCACATCCAGCCCCGTCTCCATACGCTTAATGgtgt                                     >  3:230476/1‑128 (MQ=255)
                                                                                                                           aTAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGTTAACTCAATGCTGCGGGTACTGCGAGTGAACAGCACCACATCCAGCCCCGTCTCCATACGCTTAATGgtgt                                     <  4:230476/128‑1 (MQ=255)
                                                                                                                           aTAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGTTAACTCAATGCTGCGGGTACTGCGAGTGAACAGCACCACATCCAGCCCCGTCTCCATACGCTTAATGGTGTAGCTAATGGCGGAAGTGGtt                 >  3:189533/1‑148 (MQ=255)
                                                                                                                            tAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGTTAACTCAATGCTGCGGGTACTGCGAGTGAACAGCCCCACATCCAGCCCCGTCTCCATACGCTTAATGGTGTAGCTAATGGCGGAAGTGGTTa                >  3:27062/1‑148 (MQ=255)
                                                                                                                            tAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGTTAACTCAATGCTGCGGGTACTGCGAGTGAACAGCACCACATCCAGCCCCGTCTCCATACGCTTAATGGTGTAGCTAATGGCGGAAGTGGTTaa               <  4:189533/149‑1 (MQ=255)
                                                                                                                            tAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGTTAACTCAATGCTGCGGGTACTGCGAGTGAACAGCACCACATCCAGCCCCGTCTCCATACGCTTAATGGTGTAGCTAATGGCGGAAGTGGTTa                >  3:163603/1‑148 (MQ=255)
                                                                                                                            tAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAAAGCCCGGACTCCGTTAACTCAATGCTGCGGGTACTGCGAGTGAACAGCACCACATCCAGCCCCGTCTCCATACGCTTAATGGTGTAGCTAATGGCGGAAGTGGTTa                >  8:148899/1‑148 (MQ=255)
                                                                                                                              aaaaTCATTCAGCAGGTCTGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGTTAACTCAATGGTGCGGGTACTGCGAGTGAACAGCACCACATCCAGCCCCGTCTGCATACGCATAATGGTGTAGCTAATGGCGGAAGTGGTTAAAc             <  6:92624/149‑1 (MQ=255)
                                                                                                                              aaaaTCATTCAGCAGGTCTGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGTTAACTCAATGCTGCGGGTACTGCGAGTGAACAGCACCACATCCAGCCCCGTCTCCATACGCTTAATGGTGTAGCTAATGGCGGAAGTGGTTAAAc             <  6:98308/149‑1 (MQ=255)
                                                                                                                              aaaaTCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGTTAACTCAATGCTGCGGGTACTGCGAGTGAACAGCACCACATCCAGCCCCGTCTCCATACGCTTAATGGTGTAGCTAATGGCGGAAGTGGTTaaa              <  1:211357/148‑1 (MQ=255)
                                                                                                                               aaaTCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGTTAACTCAATGCTGCGGGTACTGCGAGt                                                                                 <  6:302811/80‑1 (MQ=255)
                                                                                                                                        agcagGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGTTAACTCAATGCTGCGGGTACTGCGAGTGAACAGCACCACATCCAGCCCCGTCTCCATACGCTTAATGGTGTAGCTAATGGCGGAAGTGGTTAAACCCAACtcttc   >  7:49939/1‑149 (MQ=255)
                                                                                                                                         gcagGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGTTAACTCAATGCTGCGGGTACTGCGAGTGAACAGCACCACATCCAGCCCCGTCTCCATACGCTTAATGGTGTATCTAATGGCGTAAGTGGTTAAACCCAACtcttct  >  7:335570/1‑149 (MQ=255)
                                                                                                                                         gcagGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGTTAACTCAATGCTGCGGGTACTGCGAGTGAACAGCACCACATCCAGCCCCGTCTCCATACGCTTAATGGTGTAGCTAATGGCGGAAGTGGTTAAACCCAACtcttct  <  2:169650/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                             
CTGCTTTTTCAGCACCTGTAACAAGCGTGCGGTGTGTTTGGGCGTATAAAGCAGCTGATTAATACAAATGCGCACACGCGCTTCTATGCCCTGCGAAATAGTATCGATTCTGCGTTTGATGGCATAAAAATCATTCAGCAGGTCGGTAGCTTTGCGGAAAAAATAACGCCCGGACTCCGTTAACTCAATGCTGCGGGTACTGCGAGTGAACAGCACCACATCCAGCCCCGTCTCCATACGCTTAATGGTGTAGCTAATGGCGGAAGTGGTTAAACCCAACTCTTCT  >  NC_000913/334071‑334356

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: