Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,512,309 G→C 7.8% V12V (GTC→GTG yhfZ ← putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,512,3090GC 7.8% 405.3 / 20.8 115V12V (GTC→GTGyhfZputative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (48/58);  new base C (4/5);  total (52/63)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.12e-01

Only 100 of 121 total aligned reads displayed.

ATCCGTATCGCTCCGCTTGATTCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAA  >  NC_000913/3512163‑3512454
                                                                                                                                                  |                                                                                                                                                 
atacGTATCGCTCCGCTTGATTCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATgac                                                                                                                                                 <  4:258956/146‑1 (MQ=255)
aTCCGTATCGCTCCGCTTGATTCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATgac                                                                                                                                                 <  1:65435/149‑1 (MQ=255)
aTCCGTATCGCTCCGCTTGATTCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATgac                                                                                                                                                 >  3:748734/1‑149 (MQ=255)
      aTCGCTCCGCTTGATTCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGcc                                                                                                                                           <  1:1107363/149‑1 (MQ=255)
      aTCGCTCCGCTTGATTCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGcc                                                                                                                                           >  3:41334/1‑149 (MQ=255)
         gCTCCGCTTGATTCCAGCGTTTTCGACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCcttc                                                                                                                                        >  1:1251487/1‑149 (MQ=255)
         gCTCCGCTTGATTCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTCATCACGACGctggc                                                                                                                                        <  4:1038427/149‑5 (MQ=255)
           tCCGCTTGATTCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCcttctt                                                                                                                                      >  1:260448/1‑149 (MQ=255)
           tCCGCTTGATTCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCcttctt                                                                                                                                      >  2:502508/1‑149 (MQ=255)
             cGCTTGATTCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCtttt                                                                                                                                    <  4:68474/149‑1 (MQ=255)
              gCTTGATTCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCttttt                                                                                                                                   <  3:755524/149‑1 (MQ=255)
              gCTTGATTCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCttttt                                                                                                                                   <  3:359455/149‑1 (MQ=255)
               cTTGATTCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCtttttt                                                                                                                                  <  1:1435304/149‑1 (MQ=255)
               cTTGATTCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCtttttt                                                                                                                                  <  2:1390170/149‑1 (MQ=255)
                    ttCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATaa                                                                                                                             >  4:722616/1‑149 (MQ=255)
                     tCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATaaa                                                                                                                            <  3:786497/149‑1 (MQ=255)
                         gCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGtt                                                                                                                        >  2:514595/1‑149 (MQ=255)
                         gCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTCATCACGACGct                                                                                                                                           >  3:96907/1‑129 (MQ=255)
                         gCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTCATCACGACGct                                                                                                                                           <  4:96907/130‑2 (MQ=255)
                             tttCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGAc                                                                                                                    <  1:241505/149‑1 (MQ=255)
                                cAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCa                                                                                                                 <  2:1210701/149‑1 (MQ=255)
                                cAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCa                                                                                                                 >  4:1022743/1‑149 (MQ=255)
                                 aaCGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCat                                                                                                                >  4:562076/1‑149 (MQ=255)
                                        ccTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAAt                                                                                                         <  3:923234/149‑1 (MQ=255)
                                        ccTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTCTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAAt                                                                                                         <  1:756770/149‑1 (MQ=255)
                                                aGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTa                                                                                                 >  3:74721/1‑149 (MQ=255)
                                                aGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTa                                                                                                 >  1:302775/1‑149 (MQ=255)
                                                aGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTa                                                                                                 >  3:923031/1‑149 (MQ=255)
                                                  gCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCttttt                                                                                                                                   <  1:465870/113‑1 (MQ=255)
                                                   ccAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAAt                                                                                              <  1:630318/149‑1 (MQ=255)
                                                    cAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAAtt                                                                                             <  3:581914/149‑1 (MQ=255)
                                                       cGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGaa                                                                                          <  1:853722/149‑1 (MQ=255)
                                                       cGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGaa                                                                                          <  2:302775/149‑1 (MQ=255)
                                                       cGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGaa                                                                                          <  3:722616/149‑1 (MQ=255)
                                                            gAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATAttt                                                                                     <  4:415507/149‑1 (MQ=255)
                                                            gAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATAttt                                                                                     <  3:568091/149‑1 (MQ=255)
                                                                   attcattTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAAt                                                                                                         >  2:585988/1‑122 (MQ=255)
                                                                   attcattTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAAt                                                                                                         <  1:585988/122‑1 (MQ=255)
                                                                   attcattTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCtt                                                                              >  1:1331445/1‑149 (MQ=255)
                                                                   attcattTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCtt                                                                              >  3:721737/1‑149 (MQ=255)
                                                                     tcattTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCa                                                                            <  4:721737/149‑1 (MQ=255)
                                                                       attTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAGCGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAgg                                                                          >  3:802165/1‑149 (MQ=255)
                                                                         ttGCCAGCTCATCTATGGTTTACAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGaa                                                                        >  1:1237056/1‑149 (MQ=255)
                                                                             cAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTCATCACGACGctg                                                                                                                                          <  1:70464/79‑3 (MQ=255)
                                                                             cAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTCATCACGACGctg                                                                                                                                          >  2:70464/1‑77 (MQ=255)
                                                                               gCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAAc                                                                                                                           <  4:575145/92‑1 (MQ=255)
                                                                               gCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAAc                                                                                                                           >  3:575145/1‑92 (MQ=255)
                                                                               gCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCat                                                                                                                <  2:403634/103‑1 (MQ=255)
                                                                               gCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCat                                                                                                                >  1:403634/1‑103 (MQ=255)
                                                                               gCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGaa                                                                  <  4:1075238/149‑1 (MQ=255)
                                                                               gCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGaa                                                                  <  4:569461/149‑1 (MQ=255)
                                                                               gCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGaa                                                                  <  1:1044281/149‑1 (MQ=255)
                                                                               gCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGaa                                                                  <  2:36030/149‑1 (MQ=255)
                                                                                   aTCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATa                                                              >  3:686092/1‑149 (MQ=255)
                                                                                    tCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCtttttt                                                                                                                                  >  2:311508/1‑80 (MQ=255)
                                                                                    tCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCtttttt                                                                                                                                  <  1:311508/80‑1 (MQ=255)
                                                                                      tATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAgg                                                           >  3:249264/1‑149 (MQ=255)
                                                                                      tATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAgg                                                           <  4:1054342/149‑1 (MQ=255)
                                                                                       aTGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGgc                                                          <  2:120844/149‑1 (MQ=255)
                                                                                       aTGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGgc                                                          >  2:143511/1‑149 (MQ=255)
                                                                                        tGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGgcg                                                         >  3:251795/1‑149 (MQ=255)
                                                                                           ttttCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCaa                                                      >  3:661970/1‑149 (MQ=255)
                                                                                           ttttCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCaa                                                      >  2:1314677/1‑149 (MQ=255)
                                                                                                aaTCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCg                                                 >  3:917762/1‑149 (MQ=255)
                                                                                                 aTCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCga                                                >  1:892436/1‑149 (MQ=255)
                                                                                                  tCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCgag                                               >  3:578962/1‑149 (MQ=255)
                                                                                                  tCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCgag                                               <  3:900/149‑1 (MQ=255)
                                                                                                  tCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCgag                                               >  2:166209/1‑149 (MQ=255)
                                                                                                    gATTACCGCACTTTTCACCCAACAAATAGCGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGAt                                             <  4:802165/149‑1 (MQ=255)
                                                                                                     aTTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGAtg                                            <  1:460188/149‑1 (MQ=255)
                                                                                                     aTTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGAtg                                            <  1:460189/149‑1 (MQ=255)
                                                                                                       tACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGAtgtg                                          <  1:1309658/149‑1 (MQ=255)
                                                                                                       tACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGAtgtg                                          <  4:748734/149‑1 (MQ=255)
                                                                                                        aCCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTgg                                         <  4:251795/149‑1 (MQ=255)
                                                                                                           gCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGaaa                                      <  3:478981/149‑1 (MQ=255)
                                                                                                            cACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAAc                                     <  1:1407906/149‑1 (MQ=255)
                                                                                                                 ttCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGTGAGATGTGGAAACGGCGa                                <  1:1314677/149‑1 (MQ=255)
                                                                                                                 ttCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGa                                <  2:260448/149‑1 (MQ=255)
                                                                                                                 ttCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGa                                <  2:1341986/149‑1 (MQ=255)
                                                                                                                    aCCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGc                             >  2:61851/1‑149 (MQ=255)
                                                                                                                      ccAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCAc                           >  4:1412973/1‑149 (MQ=255)
                                                                                                                       caacaaATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCa                                                                                                                 >  2:173974/1‑62 (MQ=255)
                                                                                                                       caacaaATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCa                                                                                                                 <  1:173974/62‑1 (MQ=255)
                                                                                                                       caacaaATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACt                          >  2:787192/1‑149 (MQ=255)
                                                                                                                        aacaaATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACtt                         >  2:683697/1‑149 (MQ=255)
                                                                                                                        aacaaATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACtt                         >  2:985516/1‑149 (MQ=255)
                                                                                                                        aacaaATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACtt                         >  4:144987/1‑149 (MQ=255)
                                                                                                                         acaaATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTc                        <  3:717703/149‑1 (MQ=255)
                                                                                                                           aaaTAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTcac                      >  2:323275/1‑149 (MQ=255)
                                                                                                                            aaTAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTcaca                     <  1:502508/149‑1 (MQ=255)
                                                                                                                              tAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTCACAtt                   <  1:514595/149‑1 (MQ=255)
                                                                                                                                   ggCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTc              >  1:539681/1‑149 (MQ=255)
                                                                                                                                      cAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGa           <  4:654299/149‑1 (MQ=255)
                                                                                                                                      cAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGa           <  4:119524/149‑1 (MQ=255)
                                                                                                                                        gCGTCGTCATCACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGAtgt                                           >  1:858093/1‑115 (MQ=255)
                                                                                                                                        gCGTCGTCATCACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGAtgt                                           <  2:858093/115‑1 (MQ=255)
                                                                                                                                          gtcgtGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGattat       >  4:752234/1‑149 (MQ=255)
                                                                                                                                          gtcgtGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGattat       <  4:917762/149‑1 (MQ=255)
                                                                                                                                           tcgtGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGattatt      <  2:963150/149‑1 (MQ=255)
                                                                                                                                               gatgaCGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCaa  >  4:809252/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
ATCCGTATCGCTCCGCTTGATTCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAA  >  NC_000913/3512163‑3512454

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