Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,553,468 G→T 6.0% A82S (GCC→TCC)  hemF → coproporphyrinogen III oxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,553,4680GT 6.0% 245.8 / 2.2 86A82S (GCC→TCC) hemFcoproporphyrinogen III oxidase
Reads supporting (aligned to +/- strand):  ref base G (37/43);  new base T (3/3);  total (40/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

TCGAAGATAGTTGGCAGCGCGAAGCTGGCGGCGGCGGGCGTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACT  >  NC_000913/2553319‑2553615
                                                                                                                                                     |                                                                                                                                                   
tCGAAGATAGTTGGCAGCGCGAAGCTGGCGGCGGCGGGCGTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTg                                                                                                                                                     <  1:682190/150‑1 (MQ=255)
   aaGATAGTTGGCAGCGCGAAGCTGGCGGCGGCGGGCGTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCg                                                                                                                                                  >  1:480977/1‑150 (MQ=255)
     gATAGTTGGCAGCGCGAAGCTGGCGGCGGCGGGCGTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCggg                                                                                                                                                <  1:506228/150‑1 (MQ=255)
       tAGTTGGCAGCGCGAAGCTGGCGGCGGCGGGCGTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGgcg                                                                                                                                              <  1:986513/150‑1 (MQ=255)
       tAGTTGGCAGCGCGAAGCTGGCGGCGGCGGGCGTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGgcg                                                                                                                                              <  2:763673/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:989930/1‑150 (MQ=255)
                                     gcgTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGt                                                                                                                <  1:599171/150‑1 (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                                                                             >  2:665349/1‑150 (MQ=255)
                                                                        aaCAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGcc                                                                             >  1:1273226/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:1142224/1‑150 (MQ=255)
                                                                                   gTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAAtgtg                                                                  >  1:753421/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:1233989/150‑1 (MQ=255)
                                                                                                               cGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTg                                      <  1:92853/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:339158/4‑150 (MQ=255)
                                                                                                                                      ctgcGCCCGGCAAGTTCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGc               >  2:384617/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)
                                                                                                                                                aacaggCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCg     <  1:1247449/145‑1 (MQ=255)
                                                                                                                                                aaCTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGAGCGGTGTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCg     <  2:655484/150‑1 (MQ=255)
                                                                                                                                                  cTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGAc   >  1:457347/1‑150 (MQ=255)
                                                                                                                                                   ttGCGGGGCGCAGTTTCGAGGCTATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGATTGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACt  <  1:1068312/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
TCGAAGATAGTTGGCAGCGCGAAGCTGGCGGCGGCGGGCGTAGTCGGGTGTTGCGTAATGGTGGTGTTTTCGAACAGGCAGGCGTCAACTTTTCGCATGTCCACGGTGAGGCGATGCCTGCTTCCGCCACCGCTCATCGCCCGGAACTTGCCGGGCGCAGTTTCGAGGCGATGGGCGTTTCACTGGTAGTGCATCCGCATAACCCGTATGTTCCCACCAGCCACGCGAATGTGCGGTTTTTTATTGCCGAAAAACCGGGTGCCGATCCCGTCTGGTGGTTTGGCGGTGGCTTCGACT  >  NC_000913/2553319‑2553615

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: ‑