Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,552,023 T→C 5.1% intergenic (+32/+375) fdnI → / ← yddM formate dehydrogenase N subunit gamma/putative DNA‑binding transcriptional regulator YddM

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,552,0230TC 5.1% 306.5 / 4.2 98intergenic (+32/+375)fdnI/yddMformate dehydrogenase N subunit gamma/putative DNA‑binding transcriptional regulator YddM
Reads supporting (aligned to +/- strand):  ref base T (41/52);  new base C (3/2);  total (44/54)
Fisher's exact test for biased strand distribution p-value = 6.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

GATTAAAGGGATGATCGAAGGGAAGGTAAGTCGTCGCTGGGCGAAGAAACACCATCCGCGCTGGTATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAA  >  NC_000913/1551874‑1552169
                                                                                                                                                     |                                                                                                                                                  
gATTAAAGGGATGATCGAAGGGAAGGTAAGTCGTCGCTGGGCGAAGAAACACCATCCGCGCTGGTATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCt                                                                                                                                                    >  1:655709/1‑150 (MQ=255)
 aTTAAAGGGATGATCGAAGGGAAGGTAAGTCGTCGCTGGGCGAAGAAACACCATCCGCGCTGGTATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTg                                                                                                                                                   <  2:1753746/150‑1 (MQ=255)
    aaaGGGATGATCGAAGGGAAGGTAAGTCGTCGCTGGGCGAAGAAACACCATCCGCGCTGGTATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTaa                                                                                                                                                >  1:759817/1‑150 (MQ=255)
         gatgatCGAAGGGAAGGTAAGTCGTCGCTGGGCGAAGAAACACCATCCGCGCTGGTATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGgc                                                                                                                                           >  2:1294903/1‑150 (MQ=255)
          atgatCGAAGGGAAGGTAAGTCGTCGCTGGGCGAAGAAACACCATCCGCGCTGGTATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGgcg                                                                                                                                          <  2:444975/150‑1 (MQ=255)
             atCGAAGGGAAGGTAAGTCGTCGCTGGGCGAAGAAACACCATCCGCGCTGGTATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCa                                                                                                                                       <  2:1454207/150‑1 (MQ=255)
             atCGAAGGGAAGGTAAGTCGTCGCTGGGCGAAGAAACACCATCCGCGCTGGTATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATAGCGCTGTAAAAGGCGCCa                                                                                                                                       >  1:1223565/1‑150 (MQ=255)
              tCGAAGGGAAGGTAAGTCGTCGCTGGGCGAAGAAACACCATCCGCGCTGGTATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCAt                                                                                                                                      <  1:1087022/150‑1 (MQ=255)
              tCGAAGGGAAGGTAAGTCGTCGCTGGGCGAAGAAACACCATCCGCGCTGGTATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCAt                                                                                                                                      >  1:1760525/1‑150 (MQ=255)
                gaagggaaggTAAGTCGTCGCTGGGCAAAGAAACACCATCCGCGCTGGTATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCAttt                                                                                                                                    >  1:244909/1‑150 (MQ=255)
                   gggaaggTAAGTCGTCGCTGGGCGAAGAAACACCATCCGCGCTGGTATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCa                                                                                                                                 <  1:1619054/150‑1 (MQ=255)
                    ggaaggTAAGTCGTCGCTGGGCGAAGAAACACCATCCGCGCTGGTATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCCTTTTACAGCGCCATTTTCag                                                                                                                                <  2:1584114/150‑2 (MQ=255)
                     gaaggTAAGTCGTCGCTGGGCGAAGAAACACCATCCGCGCTGGTATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTTAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCata                                                                                                                               >  1:1392962/1‑150 (MQ=255)
                        ggTAAGTCGTCGCTGGGCGAAGAAACACCATCCGCGCTGGTATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGTGCTGTAAAAGGCGCCATTTTCATATTg                                                                                                                            >  1:495648/1‑150 (MQ=255)
                        ggTAAGTCGTCGCTGGGCGAAGAAACACCATCCGCGCTGGTATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTg                                                                                                                            >  1:706310/1‑150 (MQ=255)
                        ggTAAGTCGTCGCTGGGCGAAGAAACACCATCCGCGCTGGTATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTg                                                                                                                            <  2:1169122/150‑1 (MQ=255)
                        ggTAAGTCGTCGCTGGGCGAAGAAACACCATCCGCGCTGGTATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTg                                                                                                                            <  2:1060281/150‑1 (MQ=255)
                          tAAGTCGTCGCTGGGCGAAGAAACACCATCCGCGCTGGTATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTc                                                                                                                                  >  2:1382232/1‑142 (MQ=255)
                          tAAGTCGTCGCTGGGCGAAGAAACACCATCCGCGCTGGTATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTc                                                                                                                                  <  1:1382232/142‑1 (MQ=255)
                            aGTCGTCGCTGGGCGAAGAAACACCATCCGCGCTGGTATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGa                                                                                                                        >  2:982146/1‑150 (MQ=255)
                              tcgtcgCTGGGCGAAGAAACACCATCCGCGCTGGTATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGTGCTGTAAAAGGCGCCATTTTCATATTGTAGACa                                                                                                                      <  2:495648/150‑1 (MQ=255)
                              tcgtcgCTGGGCGAAGAAACACCATCCGCGCTGGTATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACa                                                                                                                      <  1:1214901/150‑1 (MQ=255)
                                  cgCTGGGCGAAGAAACACCATCCGCGCTGGTATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGt                                                                                                                  <  1:1993381/150‑1 (MQ=255)
                                   gCTGGGCGAAGAAACACCATCCGCGCTGGTATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTa                                                                                                                 >  1:48727/1‑150 (MQ=255)
                                     tGGGCGAAGAAACACCATCCGCGCTGGTATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAgg                                                                                                               <  1:1816714/150‑1 (MQ=255)
                                      gggCGAAGAAACACCATCCGCGCTGGTATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGc                                                                                                              >  1:1736554/1‑150 (MQ=255)
                                       ggCGAAGAAACACCATCCGCGCTGGTATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCt                                                                                                             >  2:729567/1‑150 (MQ=255)
                                        gCGAAGAAACACCATCCGCGCTGGTATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCtt                                                                                                            >  1:1214673/1‑150 (MQ=255)
                                         cGAAGAAACACCATCCGCGCTGGTATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCttt                                                                                                           >  2:1234099/1‑150 (MQ=255)
                                         cGAAGAAACACCATCCGCGCTGGTATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCCGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCttt                                                                                                           >  2:345616/1‑150 (MQ=255)
                                          gaagaaACACCATCCGCGCTGGTATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTg                                                                                                          <  1:539464/150‑1 (MQ=255)
                                            agaaACACCATCCGCGCTGGTATCGTGAACTCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGtt                                                                                                        <  2:1876265/150‑1 (MQ=255)
                                            agaaACACCATCCGCGCTGGTATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCGATTTTCATATTGTAGACAACGTAGGCTTTGtt                                                                                                        <  2:292622/150‑1 (MQ=255)
                                            agaaACACCATCCGCGCTGGTATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGtt                                                                                                        <  2:1506219/150‑1 (MQ=255)
                                              aaaCACCATCCGCGCTGGTATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCa                                                                                                      <  1:855530/150‑1 (MQ=255)
                                                acacCATCCGCGCTGGTATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATg                                                                                                    <  1:1246656/150‑1 (MQ=255)
                                                     aTCCGCGCTGGTATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGa                                                                                               <  1:120979/150‑1 (MQ=255)
                                                            cTGGTATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATgcggcg                                                                                        <  1:313491/150‑1 (MQ=255)
                                                            cTGGTATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATgcggcg                                                                                        >  2:1830679/1‑150 (MQ=255)
                                                            cTGGTATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATgcggcg                                                                                        <  1:1294903/150‑1 (MQ=255)
                                                              ggTATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTg                                                                                      <  2:1559858/150‑1 (MQ=255)
                                                                  tCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCACGCCGCATCCGGCATGAAca                                                                                  <  2:1355022/150‑2 (MQ=255)
                                                                    gTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATTCCGGATGCGGCGTGAACGcc                                                                                >  2:273692/1‑150 (MQ=255)
                                                                    gTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGcc                                                                                <  2:759817/150‑1 (MQ=255)
                                                                    gTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGca                                                                                >  1:1035443/1‑149 (MQ=255)
                                                                      gAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCtt                                                                              >  2:1438029/1‑150 (MQ=255)
                                                                       aaaTCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTa                                                                             <  1:1438029/150‑1 (MQ=255)
                                                                        aaTCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTAt                                                                            <  1:44263/150‑1 (MQ=255)
                                                                           cGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATgcgg                                                                                          >  1:1001047/1‑133 (MQ=255)
                                                                           cGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATgcgg                                                                                          <  2:1001047/133‑1 (MQ=255)
                                                                            gagaAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCgg                                                                        >  1:1723427/1‑150 (MQ=255)
                                                                            gagaAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATAAGGCGTTCACGCCGCATCCgg                                                                        <  2:820884/150‑1 (MQ=25)
                                                                             agaAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                       >  2:1811412/1‑150 (MQ=255)
                                                                               aaGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCAt                                                                     <  1:1157367/150‑1 (MQ=255)
                                                                               aaGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCAt                                                                     <  2:244909/150‑1 (MQ=255)
                                                                                 ggCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGa                                                                   <  2:706310/150‑1 (MQ=255)
                                                                                   cAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGaaa                                                                 >  1:1323244/1‑150 (MQ=255)
                                                                                     gAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCCGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAAc                                                               <  1:345616/150‑1 (MQ=255)
                                                                                       aGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAAccc                                                             >  2:866132/1‑150 (MQ=255)
                                                                                        gCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCt                                                            <  1:1331849/150‑1 (MQ=255)
                                                                                         cGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCtt                                                           >  1:1358075/1‑150 (MQ=255)
                                                                                          gAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTc                                                          <  1:51332/150‑1 (MQ=255)
                                                                                          gAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTc                                                          <  1:982146/150‑1 (MQ=255)
                                                                                           aaaaaaGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCa                                                         >  2:1211502/1‑150 (MQ=255)
                                                                                                   gagTGAAGAAGGGATATAATCTCTTTTTAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAAGGCCTTATCCGGCATGAAAACCCTTCAAATCCAAt                                                 <  2:1392962/150‑1 (MQ=255)
                                                                                                   gagTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAATTCCAAt                                                 <  1:1064759/150‑1 (MQ=255)
                                                                                                    agTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATa                                                >  1:1349889/1‑150 (MQ=255)
                                                                                                         agaagGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTg                                           >  1:1647608/1‑150 (MQ=255)
                                                                                                          gaagGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGc                                          >  1:1771036/1‑150 (MQ=255)
                                                                                                          gaagGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGc                                          <  2:1349889/150‑1 (MQ=255)
                                                                                                          gaagGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGc                                          <  1:729567/150‑1 (MQ=255)
                                                                                                            agGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAg                                        <  2:1760525/150‑1 (MQ=255)
                                                                                                              ggATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTg                                      <  2:615853/150‑1 (MQ=255)
                                                                                                                   taATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACgt                                 >  1:1421009/1‑150 (MQ=255)
                                                                                                                        tctTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAgg                            <  2:1281126/150‑1 (MQ=255)
                                                                                                                         ctTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCACTGAACGTGTAGGc                           <  2:1799371/150‑1 (MQ=255)
                                                                                                                          ttttGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGcc                          >  1:369145/1‑150 (MQ=255)
                                                                                                                            ttGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTg                        <  2:48727/150‑1 (MQ=255)
                                                                                                                             tGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGa                       <  2:1472467/150‑1 (MQ=255)
                                                                                                                             tGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGa                       >  2:1121095/1‑150 (MQ=255)
                                                                                                                              gAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGAt                      >  1:920843/1‑150 (MQ=255)
                                                                                                                              gAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGAt                      <  2:1305213/150‑1 (MQ=255)
                                                                                                                               aaCTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATa                     <  1:1094265/150‑1 (MQ=255)
                                                                                                                                aCTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATaa                    >  1:1141316/1‑150 (MQ=255)
                                                                                                                                 cTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGAtgcg                   >  1:1207836/1‑147 (MQ=255)
                                                                                                                                 cTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAg                   >  1:310401/1‑150 (MQ=40)
                                                                                                                                 cTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAg                   <  2:762575/150‑1 (MQ=40)
                                                                                                                                      aGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAGATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAg              >  1:215983/1‑150 (MQ=37)
                                                                                                                                       gCTGAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAgc             <  2:920843/150‑1 (MQ=14)
                                                                                                                                       gCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAgc             <  2:1214673/150‑1 (MQ=35)
                                                                                                                                        atGAAAATGGCGCCTTTTACAGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAgcg            >  2:68161/2‑150 (MQ=11)
                                                                                                                                         tGAAAATGGCGCCTTTTACAGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGATCGCCTTTTACGGCTGGAGACCCGTTCAAATCGAATATATTGGAGTGAACGTGTAGGCATGATAAGCGTAgcgc           >  1:550500/1‑150 (MQ=11)
                                                                                                                                           caaaTGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCa                                                                                                      <  1:1139659/56‑1 (MQ=18)
                                                                                                                                           aaaaTGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCa                                                                                                      >  2:1139659/1‑57 (MQ=38)
                                                                                                                                              aTGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAg      <  1:1234099/150‑1 (MQ=32)
                                                                                                                                                ggCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATc                                                     >  2:985251/1‑101 (MQ=32)
                                                                                                                                                  cgcTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATg                                                                    <  1:1134423/84‑1 (MQ=33)
                                                                                                                                                  cgcTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATg                                                                    >  2:1134423/1‑84 (MQ=33)
                                                                                                                                                  cgcTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCaa  >  2:1808029/1‑150 (MQ=32)
                                                                                                                                                     |                                                                                                                                                  
GATTAAAGGGATGATCGAAGGGAAGGTAAGTCGTCGCTGGGCGAAGAAACACCATCCGCGCTGGTATCGTGAAATCGAGAAGGCAGAAGCGAAAAAAGAGAGTGAAGAAGGGATATAATCTCTTTTGAACTTTAAGCTGAAAATGGCGCTGTAAAAGGCGCCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGAAAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAA  >  NC_000913/1551874‑1552169

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