Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,845,053 A→G 5.2% F473F (TTT→TTC hycE ← hydrogenase 3, large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,845,0530AG 5.2% 240.0 / 2.7 77F473F (TTT→TTChycEhydrogenase 3, large subunit
Reads supporting (aligned to +/- strand):  ref base A (32/41);  new base G (2/2);  total (34/43)
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.98e-01

CATCGGAAACGGTGTTGCCGCGCAGCATGTAGCGCAGGGTCGGCCAGTTCGCGTAGGTCGCGGCACGGCAGCGCCAGCGGTACAGCTTCTGGTTGTCGCCGGTCATGCTCCAGTGGATATCATCGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGGCAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGATACGCACTTTCAGACGGGAAATAACGTCGCA  >  NC_000913/2844911‑2845199
                                                                                                                                              |                                                                                                                                                  
cATCGGAAACGGTGTTGCCGCGCAGCATGTAGCGCAGGGTCGGCCAGTTCGCGTAGGTCGCGGCACGGCAGCGCCAGCGGTACAGCTTCTGGTTGTCGCCGGTCATGCTCCAGTGGATATCATCGCCGCGCGGCGCTTCGGCAAAGcc                                                                                                                                               >  4:311535/1‑148 (MQ=255)
  tCGGAAACGGTGTTGCCGCGCAGCATGTAGCGCAGGGTCGGCCAGTTCGCGTAGGTCGCGGCACGGCAGCGCCAGCGGTACAGCTTCTGGTTGTCGCCGGTCATGCTCCAGTGGATATCATCGCCGCGCGGCGCTTCGGCGAAGCCCAg                                                                                                                                            >  7:341056/1‑149 (MQ=255)
  tCGGAAACGGTGTTGCCGCGCAGCATGTAGCGCAGGGTCGGCCAGTTCGCGTAGGTCGCGGCACGGCAGCGCCAGCGGTACAGCTTCTGGTTGTCGCCGGTCATGCTCCAGTGGATATCATCGCCGCGCGGCGCTTCGGCAAAGCCCAg                                                                                                                                            >  8:279429/1‑149 (MQ=255)
        aCGGTGTTGCCGCGCAGCATGTAGCGCAGGGTCGGCCAGTTCGCGTAGGTCGCGGCACGGCAGCGCCAGCGGTACAGCTTCTGGTTGTCGCCGGTCATGCTCCAGTGGATATCATCGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCaaa                                                                                                                                      <  3:55529/149‑1 (MQ=255)
        aCGGTGTTGCCGCGCAGCATGTAGCGCAGGGTCGGCCAGTTCGCGTAGGTCGCGGCACGGCAGCGCCAGCGGTACAGCTTCTGGTTGTCGCCGGTCATGCTCCAGTGGATATCATCGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCaaa                                                                                                                                      <  1:256477/149‑1 (MQ=255)
             gtTGCCGCGCAGCATGTAGCGCAGGGTCGGCCAGTTCGCGTAGGTCGCGGCACGGCAGCGCCAGCGGTACAGCTTCTGGTTGTCGCCGGTCATGCTCCAGTGGATATCATCGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGt                                                                                                                                 >  4:2146/1‑149 (MQ=255)
              ttGCCGCGCAGCATGTAGCGCAGGGGCGGCCAGTTCGCGTAGGTCGCGGCACGGCAGCGCCAGCGGTACAGCTTCTGGTTGGCGCCGGTCATGCTCCAGTGGATATCATCGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTg                                                                                                                                <  4:225606/149‑1 (MQ=255)
                 ccGCGCAGCATGTAGCGCATGGTCGGCCAGTTCGCGTAGGTCGCGGCACGGCAGCGCCAGCGGTACAGCTTCTGGTTGTCGCCGGTCATGCTCCAGTGGATATCATCGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCg                                                                                                                              >  6:114936/1‑148 (MQ=255)
                 ccGCGCAGCATGTAGCGCAGGGTCGGCCAGTTCGCGTAGGTCGCGGCACGGCAGCGCCAGCGGTACAGCTTCTGGTTGTCGCCGGTCATGCTCCAGTGGATATCATCGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCg                                                                                                                              <  5:58354/148‑1 (MQ=255)
                   gcgcAGCATGTAGCGCAGGGTCGGCCAGTTCGCGTAGGTCGCGGCACGGCAGCGCCAGCGGTACAGCTTCTGGTTGTCGCCGGTCATGCTCCAGTGGATATCATCGCCGCGCTGGGCTTTGCCGAAGCGCCgcgcggc                                                                                                                                      <  6:245773/138‑4 (MQ=25)
                   gcgcAGCATGTAGCGCAGGGTCGGCCAGTTCGCGTAGGTCGCGGCACGGCAGCGCCAGCGGTACAGCTTCTGGTTGTCGCCGGTCATGCTCCAGTGGATATCATCGCCGCGCTGGGCTTTGCCGAAGCGCCgcgcggc                                                                                                                                      >  5:245773/1‑135 (MQ=25)
                    cgcAGCATGTAGCGCAGGGTCGGCCAGTTCGCGTAGGTCGCGGCACGGCAGCGCCAGCGGTACAGCTTCTGGTTGTCGCCGGTCATGCTCCAGTGGATATCATCGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAAt                                                                                                                          <  3:72038/149‑1 (MQ=255)
                       agcaTGTAGCGCAGGGTCGGCCAGTTCGCGTAGGTCGCGGCACGGCAGCGCCAGCGGTACAGCTTCTGGTTGTCGCCGGTCATGCTCCAGTGGATATCATCGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTa                                                                                                                       <  3:2146/149‑1 (MQ=255)
                        gcaTGTAGCGCAGGGTCGGCCAGTTCGCGTAGGTCGCGGCACGGCAGCGCCAGCGGTACAGCTTCTGGTTGTCGCCGGTCATGCTCCAGTGGATATCATCGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAg                                                                                                                      <  7:116930/149‑1 (MQ=255)
                          aTGTAGCGCAGGGTCGGCCAGTTCGCGTAGGTCGCGGCACGGCAGCGCCAGCGGTACAGCTTCTGGTTGTCGCCGGTCATGCTCCAGTGGATATCATCGCCGCGCGGCGCTGCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGt                                                                                                                    <  2:313426/149‑1 (MQ=255)
                           tGTAGCGCAGGGTCGGCCAGTTCGCGTAGGTCGCGGCACGGCAGCGCCAGCGGTACAGCTTCTGGTTGTCGCCGGTCATGCTCCAGTGGATATCATCGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGt                                                                                                                    >  1:310446/1‑148 (MQ=255)
                             tAGCGCAGGGTCGGCCAGTTCGCGTAGGTCGCGGCACGGCAGCGCCAGCGGTACAGCTTCTGGTTGTCGCCGGTCATGCTCCAGTGGATATCATCGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTaaa                                                                                                                 <  6:16092/149‑1 (MQ=255)
                               gcgcAGGGTCGGCCAGTTCGCGTAGGTCGCGGCACGGCAGCGCCAGCGGTACAGCTTCTGGTTGTCGCCGGTCATGCTCCAGTGGATATCATCGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAg                                                                                                                <  7:279429/148‑1 (MQ=255)
                                cgcAGGGTCGGCCAGTTCGCGTAGGTCGCGGCACGGCAGCGCCAGCGGTACAGCTTCTGGTTGTCGCCGGTCATGCTCCAGTGGATATCATCGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGcc                                                                                                              <  5:114936/149‑1 (MQ=255)
                                   aGGGTCGGCCAGTTCGCGTAGGTCGCGGCACGGCAGCGCCAGCGGTACAGCTTCTGGTTGTCGCCGGTCATGCTCCAGTGGATATCATCGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTc                                                                                                           <  2:175371/149‑1 (MQ=255)
                                   aGGGTCGGCCAGTTCGCGTAGGTCGCGGCACGGCAGCGCCAGCGGTACAGCTTCTGGTTGTCGCCGGTCATGCTCCAGTGGATATCATCGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTc                                                                                                           <  6:41807/149‑1 (MQ=255)
                                         ggCCAGTTCGCGTAGGTCGCGGCACGGCAGCGCCAGCGGTACAGCTTCTGGTTGTCGCCGGTCATGCTCCAGTGGATATCATCGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCAt                                                                                                     <  8:3579/149‑1 (MQ=255)
                                          gCCAGTTCGCGTAGGTCGCGGCACGGCAGCGCCAGCGGTACAGCTTCTGGTTGTCGCCGGTCATGCTCCAGTGGATATCATCGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCAt                                                                                                     >  3:103371/1‑148 (MQ=255)
                                                 cgcgTAGGTCGCGGCACGGCAGCGCCAGCGGTACAGCTTCTGGTTGTCGCCGGTCATGCTCCAGTGGATATCATCGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAgg                                                                                                                     <  7:253752/125‑1 (MQ=255)
                                                 cgcgTAGGTCGCGGCACGGCAGCGCCAGCGGTACAGCTTCTGGTTGTCGCCGGTCATGCTCCAGTGGATATCATCGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAgg                                                                                                                     >  8:253752/1‑125 (MQ=255)
                                                       ggTCGCGGCACGGCAGCGCCAGCGGTACAGCTTCTGGTTGTCGCCGGTCATGCTCCAGTGGATATCATCGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCg                                                                                       >  4:21466/1‑149 (MQ=255)
                                                         tCGCGGCACGGCAGCGCCAGGGGTACAGCTTCTGGTTGTCGCCGGTCATGCTCCAGTGGATATCATCGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGGc                                                                                     <  3:32477/149‑1 (MQ=255)
                                                           gcgGCACGGCAGCGCCAGCGGTACAGCTGCTGGTTGTCGCCGGTCATGCTCCAGTGGATATCATCGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGGCAg                                                                                   <  4:243471/149‑1 (MQ=255)
                                                            cggcacggcaGCGCCAGCGGTACAGCTTCTGGTTGTCGCCGGTCATGCTCCAGTGGATATCATCGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGt                                                                                                                                 <  3:176230/102‑1 (MQ=255)
                                                            cggcacggcaGCGCCAGCGGTACAGCTTCTGGTTGTCGCCGGTCATGCTCCAGTGGATATCATCGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGt                                                                                                                                 >  4:176230/1‑102 (MQ=255)
                                                             ggcacggcaGCGCCAGCGGTACAGCTTCTGGTTGTCGCCGGTCATGCTCCAGTGGATATCATCGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGGCAGGt                                                                                 >  7:280238/1‑149 (MQ=255)
                                                                 cggcaGCGCCAGCGGTACAGCTTCTGGTTGTCGCCGGTCATGCTCCAGTGGATATCATCGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGGCAGGTTATc                                                                             <  1:344993/149‑1 (MQ=255)
                                                                  ggcaGCGCCAGCGGTACAGCTTCTGGTTGTCGCCGGTCATGCTCCAGTGGATATCATCGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGGCAGGTTATcc                                                                            <  3:285299/149‑1 (MQ=255)
                                                                     aGCGCCAGCGGTACAGCTTCTGGTTGTCGCCGGTCATGCTCCAGTGGATATCATCGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGGCAGGTTATCCAGa                                                                         <  4:103371/149‑1 (MQ=255)
                                                                     aGCGCCAGCGGTACAGCTTCTGGTTGTCGCCGGTCATGCTCCAGTGGATATCATCGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGGCAGGTTATCCAGa                                                                         <  2:310446/149‑1 (MQ=255)
                                                                        gcCAGCGGTACAGCTTCTGGTTGTCGCCGGTCATGCTCCAGTGGATATCATCGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGGCAGGTTATCCAGACCg                                                                      >  7:260789/1‑149 (MQ=255)
                                                                        gcCAGCGGTACAGCTTCTGGTTGTCGCCGGTCATGCTCCAGTGGATATCATCGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGGCAGGTTATCCAGACCg                                                                      >  1:253711/1‑149 (MQ=255)
                                                                         ccAGCGGTACAGCTTCTGGTTGTCGCCGGTCATGCTCCAGTGGATATCATCGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGGCAGGTTATCCAGACCg                                                                      >  3:270566/1‑148 (MQ=255)
                                                                           aGCGGTACAGCTTCTGGGTGGCGCCGGTCATGCTCCAGTGGACATCATCGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGGCAGGTTATCCGGACCGTAg                                                                   <  8:295647/149‑1 (MQ=255)
                                                                              ggTACAGCTTCTGGTTGTCGCCGGTCATGCTCCAGTGGATATCATCGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGGCAGGTTATc                                                                             >  1:79511/1‑136 (MQ=255)
                                                                              ggTACAGCTTCTGGTTGTCGCCGGTCATGCTCCAGTGGATATCATCGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGGCAGGTTATc                                                                             <  2:79511/136‑1 (MQ=255)
                                                                               gTACAGCTTCTGGTTGTCGCCGGTCATGCTCCAGTGGATATCATAGCCGCGCGGCGCTTCGGCAAAGCCAAGCGCAAAGCGGTCCGGAAAGTAGGTAAAGCCATCCACCATCAGAGGGCCACCCGCCAGGTTATCCAGACCGTAGTCg                                                                >  3:227506/1‑148 (MQ=255)
                                                                                  cAGCTTCTGGTTGTCGCCGGTCATGCTCCAGTGGATATCATCGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGGCAGGTTATCCAGACCGTAGTCGATCa                                                            >  6:134791/1‑149 (MQ=255)
                                                                                  cAGCTTCTGGTGGTCGCCTGTCATGCTCCAGTGGATATCATCGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAATGCGGTGCTGAATGTAGGTAAAGCCTTCCAGCATCAGTGGGCCACCCGGCAGTTTATCCAGACCGTAGGCGATc                                                             <  4:50800/148‑1 (MQ=255)
                                                                                   aGCTTCTGGTTGTCGCCGGTCATGCTCCAGTGGATATCATCGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGGCAGGTTATCCAGACCGTAGTCGATCAt                                                           >  6:106755/1‑149 (MQ=255)
                                                                                        cTGGTTGTCGCCGGTCATGCTCCAGTGGATATCATCGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGGCAGGTTATCCAGACCGTAGTCGATCATGtt                                                        >  5:122828/1‑147 (MQ=255)
                                                                                           gTTGTCGCCGGTCATGCTCCAGTGGATATCATCGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGGCAGGTTATCCAGACCGTAGTCGATCATGTTCAgcg                                                   >  2:100377/1‑149 (MQ=255)
                                                                                           gTTGTCGCCGGTCATGCTCCAGTGGATATCATCGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTACCACCATCAGTGGGCCACCCGGCAGGATATCCAGACCGAAGTCGATAATGTTCAgcg                                                   >  4:208079/1‑149 (MQ=255)
                                                                                                 gCCGGTCATGCTCCAGTGGATATCATCGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGGCAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGta                                              >  4:269078/1‑148 (MQ=255)
                                                                                                  ccGGTCATGCTCCAGTGGATATCATCGCCGCGCGGCGCTTCGTCAAAGCCCAGCGCAAAGCGGTGCTGAATGTAGGTAAAGACTTCCACCATCAGTGGGCCACCCGGCAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGtata                                            <  1:92928/149‑1 (MQ=255)
                                                                                                  ccGGTCATGCTCCAGTGGATATCATCGCCGCGCGGCGCTTCGGCGAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGGCAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGtata                                            <  8:341056/149‑1 (MQ=255)
                                                                                                     gTCATGCTCCAGTGGATATCATCGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGGCAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGAc                                         >  4:53290/1‑149 (MQ=255)
                                                                                                           ctgcAGTGGATATCATCGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCTCCCGGCAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGtt                                   <  1:282461/146‑1 (MQ=255)
                                                                                                           cTCCAGTGGATATCATCGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAACGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGGCAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGt                                    <  8:287119/148‑1 (MQ=255)
                                                                                                            tCCAGTGGATATCATCGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCCCCATCAGTGGGCCACCCGGCAGGTTATCCAGACCGTAGTCg                                                                >  1:332174/1‑119 (MQ=255)
                                                                                                            tCCAGTGGATATCATCGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGGCAGGTTATCCAGCCCGTAGTCg                                                                <  2:332174/119‑1 (MQ=255)
                                                                                                            tCCAGTGGATATCATCGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGGCAGGTTATCCAGAc                                                                        >  4:71995/1‑111 (MQ=255)
                                                                                                            tCCAGTGGATATCATCGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGGCAGGTTATCCAGAc                                                                        <  3:71995/111‑1 (MQ=255)
                                                                                                            tCCAGTGGATATCATCGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGGCAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGtt                                   >  3:155981/1‑148 (MQ=255)
                                                                                                              cAGTGGATATCATCGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGGCAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGa                                 <  1:160162/148‑1 (MQ=255)
                                                                                                              cAGTGGATATCATCGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGGCAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGa                                 >  8:127703/1‑148 (MQ=255)
                                                                                                              cAGTGGATATCATCGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGGCAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGAt                                <  4:270566/149‑1 (MQ=255)
                                                                                                                 tGGATATCATCGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGGCAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGATACg                             >  2:195774/1‑149 (MQ=255)
                                                                                                                    atatCATCGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGGCAg                                                                                   <  5:323744/92‑1 (MQ=255)
                                                                                                                    atatCATCGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGGCAg                                                                                   >  6:323744/1‑92 (MQ=255)
                                                                                                                      atcatcGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGGCAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGATACGCACtt                        >  4:322448/1‑149 (MQ=255)
                                                                                                                        catcGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGGCAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGATACGCACttt                       >  1:147739/1‑148 (MQ=255)
                                                                                                                         atcGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGGCAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGATACGCACTTTc                      <  5:106755/148‑1 (MQ=255)
                                                                                                                              cgcgcgGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGGCAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGATACGCACTTTCAGACgg                <  3:21466/149‑1 (MQ=255)
                                                                                                                               gcgcgGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGGCAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGATACGCACTTTCAGACggg               <  7:127703/149‑1 (MQ=255)
                                                                                                                                          cGGAAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGGCAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGATACGCACTTTCAGACGGGAAATAACGTc     <  3:53290/148‑1 (MQ=255)
                                                                                                                                           ggCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGGCAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGATACGCACTTTCAGACGGGAAATAACGTCg    <  1:100377/148‑1 (MQ=255)
                                                                                                                                           ggCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGGCAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGATACGCACTTTCAGACGGGAAATAACGTCGc   <  4:155981/149‑1 (MQ=255)
                                                                                                                                            gCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGGCAGGTTATCCAGACCGTAGTCGATCATg                                                          >  4:143298/1‑93 (MQ=255)
                                                                                                                                            gCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGGCAGGTTATCCAGACCGTAGTCGATCATg                                                          <  3:143298/93‑1 (MQ=255)
                                                                                                                                            gCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGGCAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGATACGCACTTTCAGACGGGAAATAACGTCg    >  4:315383/1‑147 (MQ=255)
                                                                                                                                            gCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGGCAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGATACGCACTTTCAGACGGGAAATAACGTCGCa  <  3:269078/149‑1 (MQ=255)
                                                                                                                                              aaaGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGGCAGGTTATCCAGACCGTAGTCGATc                                                             <  8:303207/88‑1 (MQ=255)
                                                                                                                                              aaaGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTACCACCATCAGTGGGCCACCCGGCAGGTTATCCAGACCGTAGTCGATc                                                             >  7:303207/1‑88 (MQ=255)
                                                                                                                                              |                                                                                                                                                  
CATCGGAAACGGTGTTGCCGCGCAGCATGTAGCGCAGGGTCGGCCAGTTCGCGTAGGTCGCGGCACGGCAGCGCCAGCGGTACAGCTTCTGGTTGTCGCCGGTCATGCTCCAGTGGATATCATCGCCGCGCGGCGCTTCGGCAAAGCCCAGCGCAAAGCGGTGCGGAATGTAGGTAAAGCCTTCCACCATCAGTGGGCCACCCGGCAGGTTATCCAGACCGTAGTCGATCATGTTCAGCGCGGTATAGACTTCGTTGATACGCACTTTCAGACGGGAAATAACGTCGCA  >  NC_000913/2844911‑2845199

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: