Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,512,306 G→C 6.8% I13M (ATC→ATG 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,3060GC 6.8% 466.1 / 20.5 133I13M (ATC→ATGyhfZputative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (56/67);  new base C (4/5);  total (60/73)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.56e-01

Only 100 of 137 total aligned reads displayed.

TATAGGGCAAGGGCATTGCACATACCACGTTGTTGATATCCACATGAGTCAGCAATGCTTTGTTATCCATCTCGACCAGATAACTGCCATTGCGCCCACGGCGTTCAATCCGTATCGCTCCGCTTGATTCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATTACTTTGTAGGCCGGAT  >  NC_000913/3512056‑3512577
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                               
tataGGGCAAGGGCATTGCACATACCACGTTGTTGATATCCACATGAGTCAGCAATGCTTTGTTATCCATCTCGACCAGATAACTGCCATTGCGCCCACGGCGTTCAATCCGTATCGCTCCGCTTGATTCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCCCATCCATGGTTTTCAATCGATTACCGAACTTTCCACCCAACAAATTACGGGCCAGGCTCGAGATGACGCCGCCATCTTTTTTGaaa                                                                                                                                                                                                                                                         >  4:129747/1‑273 (MQ=255)
          gggcaTTGCACATACCACGTTGTTGATATCCACATGAGTCAGCAATGCTTTGTTATCCATCTCGACCAGATAACTGCCATTGCGCCCACGGCGTTCAATCCGTATCGCTCCGCTTGATTCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCCTCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACAACGCCTTCTTTTTTGATAAActgtcac                                                                                                                                                                                                                                                >  4:65811/1‑268 (MQ=255)
               ttGCACATACCACGTTGTTGATATCCACATGAGTCAGCAATTCTTTGTTATCCATCTCGACCAGATAACTGCCATTGCGCCCACGGCGTTCAATCCGTATCGCTCCGCTTGATTCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATCACGGGCCCGCGTCGTGATCACGACGCCTTCTTTTTTGATAAACGTTAGACGCata                                                                                                                                                                                                                                          >  2:55950/1‑275 (MQ=255)
               ttGCACATACCACGGTGTTGATATCCACATGAGTCAGCAATGCTTTGTTATCCATCTCGACCAGATAACTGCCATTGCGCCCACGGCGTTCAATCCGTATCGCTCCGCTTGATTCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCat                                                                                                                                                                                                                                           <  2:65811/274‑1 (MQ=255)
                       aCCACGTTGTTGATATCCACATGAGTCAGCAATGCTTTGTTATCCATCTCGACCAGATAACTGCCATTGCGCCCACGGCGTTCAATCCGTATCGCTCCGCTTGATTCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAAttt                                                                                                                                                                                                                                  >  4:111636/1‑275 (MQ=255)
                                      tCCACATGAGTCAGCAATGCTTTGTTATCCATCTCGACCAGATAACTGCCATTGCGCCCACGGCGTTCAATCCGTATCGCTCCGCTTGATTCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGa                                                                                                                                                                                                                      <  2:129747/272‑1 (MQ=255)
                                                                      ctcGACCAGATAACTGCCATTGCGCCCACGGCGTTCAATCCGTATCGCTCCGCTTGATTCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGgcgc                                                                                                                                                                                   <  2:30290/275‑1 (MQ=255)
                                                                                                        tCAATCCGTATCGCTCCGCTTGATTCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGtgat                                                                                                                                                                                                                                                                               <  5:1435106/149‑1 (MQ=255)
                                                                                                           aTCCGTATCGCTCCGCTTGATTCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATgac                                                                                                                                                                                                                                                                            <  5:65435/149‑1 (MQ=255)
                                                                                                           aTCCGTATCGCTCCGCTTGATTCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATgac                                                                                                                                                                                                                                                                            >  7:748734/1‑149 (MQ=255)
                                                                                                              cGTATCGCTCCGCTTGATTCCAGCGTTTTCACCGCGGCCTGTGTCATTCCAACGGATGAACTGCATTCATCTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATgacggc                                                                                                                                                                                                                                                                         <  8:1005271/149‑3 (MQ=255)
                                                                                                              cGTATCGCTCCGCTTGATTCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATgacgac                                                                                                                                                                                                                                                                         <  7:642203/149‑1 (MQ=255)
                                                                                                                 aTCGCTCCGCTTGATTCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGcc                                                                                                                                                                                                                                                                      >  7:41334/1‑149 (MQ=255)
                                                                                                                 aTCGCTCCGCTTGATTCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGcc                                                                                                                                                                                                                                                                      <  5:1107363/149‑1 (MQ=255)
                                                                                                                    gCTCCGCTTGATTCCAGCGTTTTCGACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCcttc                                                                                                                                                                                                                                                                   >  5:1251487/1‑149 (MQ=255)
                                                                                                                    gCTCCGCTTGATTCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTCATCACGACGctggc                                                                                                                                                                                                                                                                   <  8:1038427/149‑5 (MQ=255)
                                                                                                                      tCCGCTTGATTCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCcttctt                                                                                                                                                                                                                                                                 >  6:502508/1‑149 (MQ=255)
                                                                                                                       ccGCTTGATTCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCttt                                                                                                                                                                                                                                                                <  8:758519/149‑1 (MQ=255)
                                                                                                                        cGCTTGATTCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCtttt                                                                                                                                                                                                                                                               >  5:1311879/1‑149 (MQ=255)
                                                                                                                        cGCTTGATTCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCtttt                                                                                                                                                                                                                                                               <  8:1155163/149‑1 (MQ=255)
                                                                                                                        cGCTTGATTCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCtttt                                                                                                                                                                                                                                                               <  8:68474/149‑1 (MQ=255)
                                                                                                                         gCTTGATTCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCttttt                                                                                                                                                                                                                                                              <  7:359455/149‑1 (MQ=255)
                                                                                                                               ttCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATaa                                                                                                                                                                                                                                                        >  8:722616/1‑149 (MQ=255)
                                                                                                                                tCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATaaa                                                                                                                                                                                                                                                       <  7:786497/149‑1 (MQ=255)
                                                                                                                                    gCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGtt                                                                                                                                                                                                                                                   >  6:514595/1‑149 (MQ=255)
                                                                                                                                    gCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTCATCACGACGct                                                                                                                                                                                                                                                                      <  8:96907/130‑2 (MQ=255)
                                                                                                                                        tttCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGAc                                                                                                                                                                                                                                               <  5:241505/149‑1 (MQ=255)
                                                                                                                                           cAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCa                                                                                                                                                                                                                                            >  8:1022743/1‑149 (MQ=255)
                                                                                                                                           cAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCa                                                                                                                                                                                                                                            <  6:1210701/149‑1 (MQ=255)
                                                                                                                                            aaCGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCat                                                                                                                                                                                                                                           >  8:562076/1‑149 (MQ=255)
                                                                                                                                                   ccTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAAt                                                                                                                                                                                                                                    <  7:923234/149‑1 (MQ=255)
                                                                                                                                                   ccTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTCTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAAt                                                                                                                                                                                                                                    <  5:756770/149‑1 (MQ=255)
                                                                                                                                                      gCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTccc                                                                                                   <  2:111636/275‑1 (MQ=255)
                                                                                                                                                           aGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTa                                                                                                                                                                                                                            >  5:302775/1‑149 (MQ=255)
                                                                                                                                                           aGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTa                                                                                                                                                                                                                            >  7:74721/1‑149 (MQ=255)
                                                                                                                                                           aGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTa                                                                                                                                                                                                                            >  7:923031/1‑149 (MQ=255)
                                                                                                                                                             gCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCttttt                                                                                                                                                                                                                                                              <  5:465870/113‑1 (MQ=255)
                                                                                                                                                             gCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCttttt                                                                                                                                                                                                                                                              >  6:465870/1‑113 (MQ=255)
                                                                                                                                                              ccAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAAt                                                                                                                                                                                                                         <  5:630318/149‑1 (MQ=255)
                                                                                                                                                               cAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAAtt                                                                                                                                                                                                                        <  7:581914/149‑1 (MQ=255)
                                                                                                                                                                aaatgATGAACGGGTTTCATTTTCCAGCGCATTTATGGGTTTCAATCGATTACCGCACTTTTCACCCCACAAATAACGGGCCAGCGTCGTTATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCaga                                                                                         <  4:336050/271‑1 (MQ=255)
                                                                                                                                                                  cGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGaa                                                                                                                                                                                                                     <  5:853722/149‑1 (MQ=255)
                                                                                                                                                                  cGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGaa                                                                                                                                                                                                                     <  7:722616/149‑1 (MQ=255)
                                                                                                                                                                  cGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGaa                                                                                                                                                                                                                     <  6:302775/149‑1 (MQ=255)
                                                                                                                                                                       gAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATAttt                                                                                                                                                                                                                <  8:415507/149‑1 (MQ=255)
                                                                                                                                                                       gAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATAttt                                                                                                                                                                                                                <  7:568091/149‑1 (MQ=255)
                                                                                                                                                                         aCGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTtat                                                                                                                                                                                                              <  5:985385/149‑1 (MQ=255)
                                                                                                                                                                              attcattTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAAt                                                                                                                                                                                                                                    <  5:585988/122‑1 (MQ=255)
                                                                                                                                                                              attcattTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCtt                                                                                                                                                                                                         >  7:721737/1‑149 (MQ=255)
                                                                                                                                                                                tcattTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCa                                                                                                                                                                                                       <  8:721737/149‑1 (MQ=255)
                                                                                                                                                                                 cattTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAg                                                                                                                                                                                                      <  7:1230418/149‑1 (MQ=255)
                                                                                                                                                                                  attTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAGCGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAgg                                                                                                                                                                                                     >  7:802165/1‑149 (MQ=255)
                                                                                                                                                                                   tttGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGa                                                                                                                                                                                                    >  5:187540/1‑149 (MQ=255)
                                                                                                                                                                                    ttGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGaa                                                                                                                                                                                                   >  6:1407906/1‑149 (MQ=255)
                                                                                                                                                                                    ttGCCAGCTCATCTATGGTTTACAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGaa                                                                                                                                                                                                   >  5:1237056/1‑149 (MQ=255)
                                                                                                                                                                                        cAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTCATCACGACGctg                                                                                                                                                                                                                                                                     <  5:70464/79‑3 (MQ=255)
                                                                                                                                                                                          gCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAAc                                                                                                                                                                                                                                                      <  8:575145/92‑1 (MQ=255)
                                                                                                                                                                                          gCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAAc                                                                                                                                                                                                                                                      >  7:575145/1‑92 (MQ=255)
                                                                                                                                                                                          gCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCat                                                                                                                                                                                                                                           >  5:403634/1‑103 (MQ=255)
                                                                                                                                                                                          gCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGaa                                                                                                                                                                                             <  8:1075238/149‑1 (MQ=255)
                                                                                                                                                                                          gCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGaa                                                                                                                                                                                             <  6:36030/149‑1 (MQ=255)
                                                                                                                                                                                            tCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGa                                                                                                                                                                                           >  7:819630/1‑149 (MQ=255)
                                                                                                                                                                                               tCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCtttttt                                                                                                                                                                                                                                                             >  6:311508/1‑80 (MQ=255)
                                                                                                                                                                                               tCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCtttttt                                                                                                                                                                                                                                                             <  5:311508/80‑1 (MQ=255)
                                                                                                                                                                                                 tATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAgg                                                                                                                                                                                      >  7:249264/1‑149 (MQ=255)
                                                                                                                                                                                                  aTGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGgc                                                                                                                                                                                     >  6:143511/1‑149 (MQ=255)
                                                                                                                                                                                                      ttttCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCaa                                                                                                                                                                                 >  6:1314677/1‑149 (MQ=255)
                                                                                                                                                                                                      ttttCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCaa                                                                                                                                                                                 >  7:661970/1‑149 (MQ=255)
                                                                                                                                                                                                         tCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGAGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAATGGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGAAAAATTGCCTTATGCGCTCCGCTTATc                                                >  4:143990/1‑275 (MQ=255)
                                                                                                                                                                                                           aaTCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCg                                                                                                                                                                            >  7:917762/1‑149 (MQ=255)
                                                                                                                                                                                                            aTCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCga                                                                                                                                                                           >  5:892436/1‑149 (MQ=255)
                                                                                                                                                                                                             tCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCgag                                                                                                                                                                          >  6:166209/1‑149 (MQ=255)
                                                                                                                                                                                                               gATTACCGCACTTTTCACCCAACAAATAGCGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGAt                                                                                                                                                                        <  8:802165/149‑1 (MQ=255)
                                                                                                                                                                                                                aTTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGAtg                                                                                                                                                                       <  5:460188/149‑1 (MQ=255)
                                                                                                                                                                                                                aTTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGAtg                                                                                                                                                                       <  5:460189/149‑1 (MQ=255)
                                                                                                                                                                                                                  tACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGAtgtg                                                                                                                                                                     <  5:1309658/149‑1 (MQ=255)
                                                                                                                                                                                                                   aCCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTgg                                                                                                                                                                    <  8:251795/149‑1 (MQ=255)
                                                                                                                                                                                                                      gCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGaaa                                                                                                                                                                 <  7:478981/149‑1 (MQ=255)
                                                                                                                                                                                                                       cACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAAc                                                                                                                                                                <  5:1407906/149‑1 (MQ=255)
                                                                                                                                                                                                                            ttCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGTGAGATGTGGAAACGGCGa                                                                                                                                                           <  5:1314677/149‑1 (MQ=255)
                                                                                                                                                                                                                            ttCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGa                                                                                                                                                           <  6:1341986/149‑1 (MQ=255)
                                                                                                                                                                                                                            ttCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGa                                                                                                                                                           <  6:260448/149‑1 (MQ=255)
                                                                                                                                                                                                                               aCCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGc                                                                                                                                                        >  6:61851/1‑149 (MQ=255)
                                                                                                                                                                                                                                 ccAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCAc                                                                                                                                                      >  8:1412973/1‑149 (MQ=255)
                                                                                                                                                                                                                                 ccAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCAc                                                                                                                                                      >  8:302601/1‑149 (MQ=255)
                                                                                                                                                                                                                                  caacaaATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCa                                                                                                                                                                                                                                            <  5:173974/62‑1 (MQ=255)
                                                                                                                                                                                                                                  caacaaATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACt                                                                                                                                                     >  6:787192/1‑149 (MQ=255)
                                                                                                                                                                                                                                   aacaaATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACtt                                                                                                                                                    >  6:985516/1‑149 (MQ=255)
                                                                                                                                                                                                                                   aacaaATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACtt                                                                                                                                                    >  8:144987/1‑149 (MQ=255)
                                                                                                                                                                                                                                   aacaaATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACtt                                                                                                                                                    >  6:683697/1‑149 (MQ=255)
                                                                                                                                                                                                                                    acaaATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTc                                                                                                                                                   <  7:717703/149‑1 (MQ=255)
                                                                                                                                                                                                                                      aaaTAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTcac                                                                                                                                                 >  6:323275/1‑149 (MQ=255)
                                                                                                                                                                                                                                         tAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTCACAtt                                                                                                                                              <  5:514595/149‑1 (MQ=255)
                                                                                                                                                                                                                                                 cAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGa                                                                                                                                      <  8:654299/149‑1 (MQ=255)
                                                                                                                                                                                                                                                  ggCGTCGTCATCACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGat                                                                                                                                     >  5:53122/2‑149 (MQ=255)
                                                                                                                                                                                                                                                   gCGTCGTCATCACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGAtgt                                                                                                                                                                      <  6:858093/115‑1 (MQ=255)
                                                                                                                                                                                                                                                     gtcgtGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGattat                                                                                                                                  >  8:752234/1‑149 (MQ=255)
                                                                                                                                                                                                                                                     gtcgtGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGattat                                                                                                                                  <  8:917762/149‑1 (MQ=255)
                                                                                                                                                                                                                                                      tcgtGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGattatt                                                                                                                                 <  6:693644/149‑1 (MQ=255)
                                                                                                                                                                                                                                                       tttgaagaCGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAAGTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGGCCAACAGCGAGATGTAGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATTACTTTGTAGGCCGGAt  <  1:6561/273‑1 (MQ=255)
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                               
TATAGGGCAAGGGCATTGCACATACCACGTTGTTGATATCCACATGAGTCAGCAATGCTTTGTTATCCATCTCGACCAGATAACTGCCATTGCGCCCACGGCGTTCAATCCGTATCGCTCCGCTTGATTCCAGCGTTTTCAACGCGGCCTGCGTCAGGCCAACGGATGAACGGCATTCATTTGCCAGCTCATCTATGGTTTTCAATCGATTACCGCACTTTTCACCCAACAAATAACGGGCCAGCGTCGTGATGACGACGCCTTCTTTTTTGATAAACGTTCGACGCATAATAAATTTTCAGTAAATTGAATATTTATATCTTCAGGAATTTGAAGATAAGGCGCAACAGCGAGATGTGGAAACGGCGAGGCACTTCACATTTTTTCGGATTATTGCAAATGAGGATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATTACTTTGTAGGCCGGAT  >  NC_000913/3512056‑3512577

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