Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,987,104 T→A 5.3% intergenic (+28/‑432) yqeG → / → yqeH putative transporter/putative LuxR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,987,1040TA 5.3% 325.9 / 7.2 95intergenic (+28/‑432)yqeG/yqeHputative transporter/putative LuxR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (46/44);  new base A (3/2);  total (49/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.56e-01

TTCATCATGCCTACGCTGTCAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGG‑ATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTATTTAAGATGCATAACTT  >  NC_000913/2986957‑2987252
                                                                                                                                                    |                                                                                                                                                     
ttCATCATGCCTACGCTGTCAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGG‑ATGGTTCAAAATTGCTa                                                                                                                                                      >  2:563020/1‑149 (MQ=255)
   atcatGCCTACGCTGTCAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGG‑ATGGTTCAAAATTGCTA‑Caa                                                                                                                                                  >  3:89203/1‑149 (MQ=255)
   atcatGCCTACGCTGTCAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGG‑ATGGTTCAAAATTGCTA‑Caa                                                                                                                                                  >  4:176682/1‑149 (MQ=255)
   atcatGCCTACGCTGTCAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGG‑ATGGTTCAAAATTGCTA‑Caa                                                                                                                                                  >  3:187722/1‑149 (MQ=255)
    tcatGCCTACGCTGTCAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGG‑ATGGTTCAAAATTGCTA‑CAAt                                                                                                                                                 <  1:399122/149‑1 (MQ=255)
       tGCCTACGCTGTCAACGTATCTCATCCTGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGG‑ATGGTTCAAAATTGCTA‑CAAtttt                                                                                                                                              <  1:681496/149‑1 (MQ=255)
            aCGCTGTCAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGGAATGGTTCAAAATTG‑TAGCAATTTTGAAc                                                                                                                                          <  2:630962/149‑1 (MQ=255)
            aCGCTGTCAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGG‑ATGGTTCAAAATTGCTA‑CAATTTTGAAcc                                                                                                                                         <  3:343445/149‑1 (MQ=255)
               cTGTCAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGG‑ATGGTTCAAAATTGCTA‑CAATTTTGAAccatt                                                                                                                                      <  2:61485/149‑1 (MQ=255)
                  tCAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGG‑ATGGTTCAAAATTGCTA‑CAATTTTGAAccattcca                                                                                                                                   <  4:353975/149‑1 (MQ=255)
                        tATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGG‑ATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTc                                                                                                                             <  3:372351/149‑1 (MQ=255)
                         aTCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGG‑ATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCa                                                                                                                            <  3:248777/149‑1 (MQ=255)
                         aTCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGG‑ATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCa                                                                                                                            <  1:61762/149‑1 (MQ=255)
                          tctcATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGG‑ATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCat                                                                                                                           >  1:136/1‑149 (MQ=255)
                            tcATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGG‑ATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCataa                                                                                                                         >  1:469112/1‑149 (MQ=255)
                                cccGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGG‑ATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATAtt                                                                                                                     >  3:74723/1‑149 (MQ=255)
                                  cGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGG‑ATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCa                                                                                                                   >  1:99895/1‑149 (MQ=255)
                                     cgcTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGG‑ATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCAttt                                                                                                                <  2:243750/149‑1 (MQ=255)
                                     cgcTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGG‑ATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCAttt                                                                                                                <  2:157766/149‑1 (MQ=255)
                                        ttAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGG‑ATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGtt                                                                                                             >  2:141809/1‑149 (MQ=255)
                                         tAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGG‑ATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGtta                                                                                                            <  3:432122/149‑1 (MQ=255)
                                            aCCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGG‑ATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGttatta                                                                                                         <  1:344001/149‑1 (MQ=255)
                                              ccTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGG‑ATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTattatt                                                                                                       >  3:338397/1‑149 (MQ=255)
                                               cTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGG‑ATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTc                                                                                                      <  1:82752/149‑1 (MQ=255)
                                                tGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGG‑ATGGTTCAAAATTGCTA‑CAATTTTGa                                                                                                                                            <  2:536444/110‑1 (MQ=255)
                                                tGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGG‑ATGGTTCAAAATTGCTA‑CAATTTTGa                                                                                                                                            >  1:536444/1‑110 (MQ=255)
                                                tGGCGATCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGG‑ATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTcc                                                                                                     >  2:464635/1‑149 (MQ=255)
                                                        cATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGG‑ATGGTTCAAAATTGTAG‑CAATTTTGaa                                                                                                                                           >  1:494749/1‑103 (MQ=255)
                                                        cATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGG‑ATGGTTCAAAATTGTAG‑CAATTTTGaa                                                                                                                                           <  2:494749/103‑1 (MQ=255)
                                                         aTCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGG‑ATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCa                                                                                            >  4:629199/1‑149 (MQ=255)
                                                            ggAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGG‑ATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACtt                                                                                         >  2:579176/1‑149 (MQ=255)
                                                             gAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGG‑ATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTtc                                                                                        <  1:346118/149‑1 (MQ=255)
                                                             gAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGG‑ATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTtc                                                                                        >  4:580218/1‑149 (MQ=255)
                                                                aTCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGG‑ATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGtta                                                                                                            >  4:373546/1‑126 (MQ=255)
                                                                aTCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGG‑ATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGtta                                                                                                            <  3:373546/126‑1 (MQ=255)
                                                                 tCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGG‑ATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTtctctc                                                                                    >  3:165427/1‑149 (MQ=255)
                                                                    gATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGG‑ATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTc                                                                                 >  2:335911/1‑149 (MQ=255)
                                                                    gATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGG‑ATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTc                                                                                 >  2:379088/1‑149 (MQ=255)
                                                                      ttACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGG‑ATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCat                                                                               <  1:289482/149‑1 (MQ=255)
                                                                            tGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGG‑ATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTa                                                                         >  1:463063/1‑149 (MQ=255)
                                                                            tGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGG‑ATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTa                                                                         <  2:516617/149‑1 (MQ=255)
                                                                                 gTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGG‑ATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGtatt                                                                    >  4:478945/1‑149 (MQ=255)
                                                                                  tGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGG‑ATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGtatta                                                                   <  1:66752/149‑1 (MQ=255)
                                                                                   gggTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGG‑ATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGtattat                                                                  <  1:531651/149‑1 (MQ=255)
                                                                                   gggTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGG‑ATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGtattat                                                                  <  2:136/149‑1 (MQ=255)
                                                                                   gggTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGG‑ATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGtattat                                                                  >  4:425839/1‑149 (MQ=255)
                                                                                    ggTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGG‑ATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATc                                                                 >  2:337295/1‑149 (MQ=255)
                                                                                     gTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGG‑ATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCa                                                                <  1:318432/149‑1 (MQ=255)
                                                                                       aTCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCGGTTGTGG‑AGGTTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTGATAATATTCATTTGTTATTATTCCACCTTATCACTTCTCTCGTCATATCCTAGTATTATCagt                                                              >  2:697580/1‑149 (MQ=255)
                                                                                        tCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGG‑ATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCagtt                                                             <  1:37972/149‑1 (MQ=255)
                                                                                        tCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGG‑ATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCagtt                                                             <  4:134603/149‑1 (MQ=255)
                                                                                          cTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGG‑ATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCGCTTCTCTCGTCATATCCTAGTATTATCagttag                                                           >  3:620006/1‑149 (MQ=255)
                                                                                          cTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGG‑ATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCagttag                                                           >  3:130071/1‑149 (MQ=255)
                                                                                          cTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGG‑ATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCagttag                                                           >  3:187212/1‑149 (MQ=255)
                                                                                           tgtgCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGG‑ATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCagttagt                                                          >  4:9315/1‑149 (MQ=255)
                                                                                           tgtgCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGG‑ATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCagttagc                                                          <  1:335911/149‑2 (MQ=255)
                                                                                           tgtgCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGG‑ATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATACTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCagttagt                                                          >  4:282707/1‑149 (MQ=255)
                                                                                            gtgCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGG‑ATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCagttagtt                                                         >  1:166602/1‑149 (MQ=255)
                                                                                             tgCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGG‑ATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAgttagtta                                                        <  3:580218/149‑1 (MQ=255)
                                                                                               cGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGG‑ATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTaaa                                                      <  2:463063/149‑1 (MQ=255)
                                                                                               cGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGG‑ATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTaaa                                                      <  3:122550/149‑1 (MQ=255)
                                                                                                gTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGG‑GTGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAac                                                     >  2:631970/1‑149 (MQ=255)
                                                                                                gTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGG‑ATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGCCATATCCTAGTATTATCAGTTAGTTAAac                                                     >  4:534495/1‑149 (MQ=255)
                                                                                                gTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGG‑ATGGTACAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAac                                                     >  4:233159/1‑149 (MQ=255)
                                                                                                  aTCGGAAATGTGCCTTAGCTGATATCAGTTGTGG‑ATGGGTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAacac                                                   >  2:41441/1‑149 (MQ=255)
                                                                                                         aTGTTCTTTAGCTAATATCAGTTGTGG‑ATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAaga                                            >  4:252190/1‑149 (MQ=255)
                                                                                                         aTGTTCTTTAGCTAATATCAGTTGTGG‑ATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAaga                                            >  1:233165/1‑149 (MQ=255)
                                                                                                          tGTTCTTTAGCTAATATCAGTTGTGG‑ATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAagag                                           <  3:138921/149‑1 (MQ=255)
                                                                                                           gTTCTTTAGCTAATATCAGTTGTGG‑ATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAgaga                                          <  1:563020/149‑1 (MQ=255)
                                                                                                            ttCTTTAGCTAATATCAGTTGTGG‑ATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGAt                                         <  3:176682/149‑1 (MQ=255)
                                                                                                            ttCTTTAGCTAATATCAGTTGTGG‑ATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGAt                                         <  1:579176/149‑1 (MQ=255)
                                                                                                            ttCTTTAGCTAATATCAGTTGTGG‑ATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGAt                                         <  4:130071/149‑1 (MQ=255)
                                                                                                             tCTTTAGCTAATATCAGTTGTGG‑ATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATc                                        <  4:421808/149‑1 (MQ=255)
                                                                                                             tCTTTAGCTAATATCAGTTGTGG‑ATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATc                                        <  3:478945/149‑1 (MQ=255)
                                                                                                                ttAGCTAATATCAGTTGTGG‑ATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGa                                     >  4:98477/1‑149 (MQ=255)
                                                                                                                ttAGCTAATATCAGTTGTGG‑ATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGa                                     <  1:352171/149‑1 (MQ=255)
                                                                                                                 tAGCTAATATCAGTTGTGG‑ATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGaa                                    <  2:469112/149‑1 (MQ=255)
                                                                                                                     tAATATCAATTGTGG‑ATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTACCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATc                                <  2:574330/149‑1 (MQ=255)
                                                                                                                      aaTATCAGTTGTGG‑ATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTaaa                                                      >  4:113589/1‑126 (MQ=255)
                                                                                                                      aaTATCAGTTGTGG‑ATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTaaa                                                      <  3:113589/126‑1 (MQ=255)
                                                                                                                      aaTATCAGTTGTGG‑ATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCt                               >  2:58683/1‑149 (MQ=255)
                                                                                                                         atCAGTTGTGG‑ATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTaa                            >  4:226604/1‑149 (MQ=255)
                                                                                                                                   g‑ATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTAt                  >  2:24628/1‑149 (MQ=255)
                                                                                                                                      tGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTAttt                >  3:258066/1‑149 (MQ=255)
                                                                                                                                      tGGTTCAAAATTG‑TAGCAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTAttt                >  4:111666/1‑149 (MQ=255)
                                                                                                                                       ggTTCAAAATTGTAG‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTATTTa               >  1:502784/1‑149 (MQ=255)
                                                                                                                                        gTTCAAAATTGTAG‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATaa                        <  4:454997/139‑1 (MQ=255)
                                                                                                                                        gTTCAAAATTGTAG‑CAATTTTGAACCATTCCATTCTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATaa                        >  3:454997/1‑139 (MQ=255)
                                                                                                                                        gTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTATTTaa              >  1:417089/1‑149 (MQ=255)
                                                                                                                                          tCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTATTTAAGa            <  4:605660/149‑1 (MQ=255)
                                                                                                                                           ccaaaTTGCTA‑CAATTTTGAACCATTCCATTTTTGATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTATTTAAGAt           <  3:556193/147‑1 (MQ=255)
                                                                                                                                            aaaaTTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGa                      >  2:255684/1‑137 (MQ=255)
                                                                                                                                            aaaaTTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGa                      <  1:255684/137‑1 (MQ=255)
                                                                                                                                              aaTTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTATTTAAGATGCa        >  3:627/1‑149 (MQ=255)
                                                                                                                                               aTTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTATTTAAGATGCAt       <  3:252190/149‑1 (MQ=255)
                                                                                                                                                    tA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTATTTAAGATGCATAACtt  >  2:66914/1‑149 (MQ=255)
                                                                                                                                                    tA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTATTTAAGATGCATAACtt  >  4:197133/1‑149 (MQ=255)
                                                                                                                                                    tA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTATTTAAGATGCATAACtt  >  4:643023/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TTCATCATGCCTACGCTGTCAACGTATCTCATCCCGGCGCTTAAACCCTGGCGTTCCATCGGAAATCTGATTACGCTGATCGTGGGTATCCTGTGCGTATCGGTAATGTTCTTTAGCTAATATCAGTTGTGG‑ATGGTTCAAAATTGCTA‑CAATTTTGAACCATTCCATTTTTCATAATATTCATTTGTTATTATTCCACCCTATCACTTCTCTCGTCATATCCTAGTATTATCAGTTAGTTAAACACACAAAGAGATCTGAACATCTTAAATAAGACTATTTAAGATGCATAACTT  >  NC_000913/2986957‑2987252

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