Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,553,463 A→C 6.1% E80A (GAA→GCA)  hemF → coproporphyrinogen III oxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,553,4630AC 6.1% 235.7 / 3.3 85E80A (GAA→GCA) hemFcoproporphyrinogen III oxidase
Reads supporting (aligned to +/- strand):  ref base A (36/43);  new base C (3/3);  total (39/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGTCGAAGATAGTTGGCAGCGCGAAGCTGGCGGCGGCGGGCGTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCG  >  NC_000913/2553316‑2553612
                                                                                                                                                   |                                                                                                                                                     
ttGTCGAAGATAGTTGGCAGCGCGAAGCTGGCGGCGGCGGGCGTAGTCGGGTGTTGCGTAATGGTGGTGGTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAAc                                                                                                                                                     <  1:335569/150‑1 (MQ=255)
   tCGAAGATAGTTGGCAGCGCGAAGCTGGCGGCGGCGGGCGTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTg                                                                                                                                                  <  1:682190/150‑1 (MQ=255)
      aaGATAGTTGGCAGCGCGAAGCTGGCGGCGGCGGGCGTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCg                                                                                                                                               >  1:480977/1‑150 (MQ=255)
        gATAGTTGGCAGCGCGAAGCTGGCGGCGGCGGGCGTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCggg                                                                                                                                             <  1:506228/150‑1 (MQ=255)
          tAGTTGGCAGCGCGAAGCTGGCGGCGGCGGGCGTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGgcg                                                                                                                                           <  2:763673/150‑1 (MQ=255)
          tAGTTGGCAGCGCGAAGCTGGCGGCGGCGGGCGTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGgcg                                                                                                                                           <  1:986513/150‑1 (MQ=255)
            cttGGCAGCGCGAAGCTGGCGGCGGCGGGCGTAGTCGGGTGCTGCGTAATGGTGGTGATTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGCAAGTTCCGGgcgat                                                                                                                                         <  1:160253/149‑3 (MQ=255)
             ttGGCAGCGCGAAGCTGGCGGCGGCGGGCGTAGTCGGGTGTTGCGTAATGGTGGTGGTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGCAAGTTCCGGgcgatg                                                                                                                                        <  2:242795/150‑4 (MQ=255)
             ttGGCAGCGCGAAGCTGGCGGCGGCGGGCGTAGTCGGGTGTTGAGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTATTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAg                                                                                                                                        <  2:230409/150‑1 (MQ=255)
              tGGCAGCGCGAAGCTGGCGGCGGCGGGCGTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGt                                                                                                                                       >  2:1067615/1‑150 (MQ=255)
                  aGCGCGAAGCTGGCGGCGGCGGGCGTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCg                                                                                                                                   <  2:985468/150‑1 (MQ=255)
                    cgcgAAGCTGGCGGCGGCGGGCGTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAg                                                                                                                                 <  2:558492/150‑1 (MQ=255)
                      cgAAGCTGGCGGCGGCGGGCGTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGc                                                                                                                               >  2:772238/1‑150 (MQ=255)
                      cgAAGCTGGCGGCGGCGGGCGTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCCGAACTTGCCGGGCGCAGTTTggcggc                                                                                                                               >  2:523663/1‑150 (MQ=255)
                      cgAAGCTGGCGGCGGCGGGCGTAGGCGGGTGTTGCGTAATGGTGGGGTTTTCGAACAGGCAGGCGTCAACTTTACGCATGGCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGc                                                                                                                               <  1:843337/150‑1 (MQ=255)
                           cTGGCGGCGGCGGGCGTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGCGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATgg                                                                                                                          >  1:1206846/1‑150 (MQ=255)
                            tGGCGGCGGCGGGCGTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCAAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATggg                                                                                                                         <  2:313789/150‑1 (MQ=255)
                             ggcggcggcGGGCGTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGc                                                                                                                        >  1:422164/1‑150 (MQ=255)
                              gcggcggcgGGCGTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCg                                                                                                                       <  1:746694/150‑1 (MQ=255)
                                  cggcggGCGTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCTATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTc                                                                                                                   <  1:578454/150‑1 (MQ=255)
                                    gcgggcgTAGTCGGGTGTTGCGTAATGGTGGGGTGTTCGAACAGGCAGGCGTCAACTTTACGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCAc                                                                                                                 <  2:109576/150‑1 (MQ=255)
                                     cgggcgTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACt                                                                                                                >  2:755009/1‑150 (MQ=255)
                                     cgggcgTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCAACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACt                                                                                                                >  1:352285/1‑150 (MQ=255)
                                     cgggcgGAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCGTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACt                                                                                                                <  2:1013658/150‑1 (MQ=255)
                                       ggcgTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCTCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTgg                                                                                                              >  2:662121/1‑150 (MQ=255)
                                       ggcgTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTgg                                                                                                              <  2:859119/150‑1 (MQ=255)
                                        gcgTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGt                                                                                                             <  1:599171/150‑1 (MQ=255)
                                        gcgTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGt                                                                                                             >  1:989930/1‑150 (MQ=255)
                                        gcgTAGTCGGGTGTGGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGGCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGt                                                                                                             <  1:1045814/150‑1 (MQ=255)
                                             gTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGt                                                                                                                                       >  2:1301774/1‑119 (MQ=255)
                                             gTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGt                                                                                                                                       <  1:1301774/119‑1 (MQ=255)
                                              tCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCa                                                                                                       <  2:1134276/150‑1 (MQ=255)
                                              tCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCa                                                                                                       <  2:104841/150‑1 (MQ=255)
                                                    gtTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCa                                                                                                 <  1:1076331/150‑1 (MQ=255)
                                                     ttGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCAt                                                                                                <  2:670200/150‑1 (MQ=255)
                                                      tGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATa                                                                                               <  2:1307932/150‑1 (MQ=255)
                                                       gCGTAAAGGTGGTGGTGTCGAACAGGCAGGCGTCAACTTCTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCTTTTCACTGGTAGTGCATCCGCATaa                                                                                              <  1:199842/150‑1 (MQ=255)
                                                          tAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGCACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAAccc                                                                                           >  1:68910/1‑150 (MQ=255)
                                                          aaaTGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCAGCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCTATGGGCGTTTCACTCGTAGTGCATCCGCATAAccc                                                                                           >  2:794692/2‑150 (MQ=255)
                                                                 ggtgTGTTCGAACAGGCAGGCTTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGtt                                                                                    <  2:1055526/150‑1 (MQ=255)
                                                                    gtTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCCCCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCTTTACCCGTATGTTccc                                                                                 >  2:499855/1‑150 (MQ=255)
                                                                        tCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCccacca                                                                             <  2:712115/150‑1 (MQ=255)
                                                                           aaCAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGcc                                                                          >  1:1273226/1‑150 (MQ=255)
                                                                           aaCAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGcc                                                                          >  2:665349/1‑150 (MQ=255)
                                                                               ggcaggcGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCAcgc                                                                      <  1:1083281/150‑1 (MQ=255)
                                                                                   ggcGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAAt                                                                  >  1:1309653/1‑150 (MQ=255)
                                                                                   ggcGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCACATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCCCTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCCCGCGAAt                                                                  >  2:676927/1‑150 (MQ=255)
                                                                                     cGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAAtgt                                                                >  1:13249/1‑150 (MQ=255)
                                                                                      gTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAAtgtg                                                               >  1:753421/1‑150 (MQ=255)
                                                                                      gTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAAtgtg                                                               >  1:1142224/1‑150 (MQ=255)
                                                                                       tCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTTCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTTGTAGTGCATCCGCCTAAACCGGTTGTTCCCCCCAGCCACGCCAATGAGc                                                              >  2:177367/1‑150 (MQ=255)
                                                                                           cTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGtt                                                          >  2:862624/1‑150 (MQ=255)
                                                                                             tctCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGtgtt                                                        >  1:1266321/3‑150 (MQ=255)
                                                                                              ttCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGCACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGttttt                                                       <  2:870196/150‑1 (MQ=255)
                                                                                                cGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTa                                                     <  2:906643/150‑1 (MQ=255)
                                                                                                 gCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGATTCGAGGCGATGGGCGATTCACTGGTAGTGCATCCGCATAACCCGTATGTTCACACCAGCCACGCGAAGGTGCGGTTTTTGAt                                                    >  1:246645/1‑150 (MQ=255)
                                                                                                    tGTCCTCGGTGCGGCGATGCGTGCTTCCGCCACCGCTCCTCGCCCGGAACATGCCGTGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGc                                                 <  1:49240/150‑1 (MQ=255)
                                                                                                     gTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGcc                                                >  1:28296/1‑150 (MQ=255)
                                                                                                     gTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAAATTGCCGCGCGGAGTTTCGAGGCGATGGGCGTTTCACTGGGAGTGCATCCGCAGCACCCGTAGGTTCCCCCCAGCCCCGCGACTGTGCGGTTTTTGATTGcc                                                >  1:227403/1‑150 (MQ=255)
                                                                                                      tCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCg                                               <  2:917598/150‑1 (MQ=255)
                                                                                                        cACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGACCTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTAGTGCCGaa                                             >  1:1240442/1‑150 (MQ=255)
                                                                                                         aCGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGaaa                                            <  1:665349/150‑1 (MQ=255)
                                                                                                          cGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGaaaa                                           >  1:347177/1‑150 (MQ=255)
                                                                                                              gAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCg                                       <  1:951980/150‑1 (MQ=255)
                                                                                                              gAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTAGCCGGGGGCAGGTTCGAGGCGATGGGCGTGTCACTGGTAGTGCATCCGCAAAACCCGTATGTTCCCACGCGCCACGCGAATGTGCGCTTTATTATAGCCGAAAAACCg                                       >  1:194015/1‑150 (MQ=255)
                                                                                                                ggCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCGTCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCggg                                     <  1:772238/150‑1 (MQ=255)
                                                                                                                  cGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTg                                   <  1:92853/150‑1 (MQ=255)
                                                                                                                  cGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTg                                   <  1:1233989/150‑1 (MQ=255)
                                                                                                                    aTGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGcc                                 <  2:20700/150‑1 (MQ=255)
                                                                                                                    aTGCCTGCTTCCGCCACCGCTCATAGCCCGGAACTTTCCGTGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCGACTCGAATGTGCGGTTTTTTATTGCCGATAAACCGGTTGcc                                 <  2:733992/150‑1 (MQ=255)
                                                                                                                       ccTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCCCCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGAt                              >  2:1305958/1‑150 (MQ=255)
                                                                                                                       ccTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCAAGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGAt                              >  2:669767/1‑150 (MQ=255)
                                                                                                                          gCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATccc                           <  2:723140/150‑1 (MQ=255)
                                                                                                                           cTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGATTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACCCGCAGGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATCCCg                          >  1:702058/1‑150 (MQ=255)
                                                                                                                              ccGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCt                       <  1:328311/150‑1 (MQ=255)
                                                                                                                                  cACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCtggtg                   <  1:371634/150‑1 (MQ=255)
                                                                                                                                   aCCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCtggtgg                  <  1:119296/150‑1 (MQ=255)
                                                                                                                                   aCCGCGCATCGCCCGGACAGTGCCGGGCGCCGTTTCGAGGCGATGGGCGTTTCACTGGGAGGGCATCCGCATATCCCGTATGTTCGCCCCAGCCACGCGAATGGGCGGTTTTTTATTGCCGAAAAACCGGGGGCCGATCCCGTCtggtgg                  <  1:707051/150‑1 (MQ=255)
                                                                                                                                       cTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTg              >  2:557879/1‑150 (MQ=255)
                                                                                                                                         ctgcGCCCGGCAAGTTCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGc            >  2:384617/4‑150 (MQ=255)
                                                                                                                                         ctgcGCCCGGCAAGTTCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGc            >  2:339158/4‑150 (MQ=255)
                                                                                                                                         cATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCCTCCGCATAACCCGTATGTTCCCTCCAGCCACGCGAATGTGCGGTTGTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGc            >  1:190363/1‑150 (MQ=255)
                                                                                                                                         cATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGc            >  2:1035482/1‑150 (MQ=255)
                                                                                                                                              cccGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTgg       >  2:8243/1‑150 (MQ=255)
                                                                                                                                              cccGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTgg       <  1:1079844/150‑1 (MQ=255)
                                                                                                                                                   aaCTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGAGCGGTGTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCg  <  2:655484/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TTGTCGAAGATAGTTGGCAGCGCGAAGCTGGCGGCGGCGGGCGTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCG  >  NC_000913/2553316‑2553612

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