Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 216,081 C→T 6.9% intergenic (+102/+98) nlpE → / ← yaeF lipoprotein NlpE/peptidase C92 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913216,0810CT 6.9% 221.3 / 4.6 72intergenic (+102/+98)nlpE/yaeFlipoprotein NlpE/peptidase C92 family protein
Reads supporting (aligned to +/- strand):  ref base C (34/33);  new base T (3/2);  total (37/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.45e-01

GGAAATTTTACCCCAACCAGGATTGCAGTAGTTTGGGGCAGTAACCCGTCTTGAGACAGAAACAAACGCAAAACAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGACCCGCCTTGATATAAATCCCAGGC  >  NC_000913/215936‑216227
                                                                                                                                                 |                                                                                                                                                  
ggAAATTTTACCCCAACCAGGATTGCAGTAGTTTGGGGCAGTAACCCGTCTTGAGACAGAAACAAACGCAAAACAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATcc                                                                                                                                                >  1:1233707/1‑150 (MQ=255)
   aaTTTTACCCCATCCAGGATTGCAGTAGGTTGTGGCAGTAACCCGTCTTGAGACAGAAACAAACGGAAAACAGGCGAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTGTGTAGGCCGGATAAGGCGTGTACGCCGCATCAGGc                                                                                                                                             <  2:653159/150‑1 (MQ=38)
   aaTTTTACCCCAACCAGGATTGCAGTAGTTTGGGGCAGTAACCCGTCTTGAGACAGAAACAAACGCAAAACAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCAGGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGc                                                                                                                                             >  2:1073240/1‑150 (MQ=255)
   aaTTTTACCCCAACCAGGATTGCAGTAGTTTGGGGCAGTAACCCGTCTTGAGACAGAAACAAACGCAAAACAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTGTTTTGTAGGCCGGATGCGGCGTAAACGCCTTATCCGGc                                                                                                                                             <  2:1754973/150‑1 (MQ=255)
       ttACCCCAACCAGGATTGCAGTAGTTTGGGGCAGTAACCCGTCTTGAGACAGAAACAAACGCAAAACAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATGCGGCGTAAACGCCTTATCCGGcctac                                                                                                                                         <  2:406651/150‑5 (MQ=255)
        tACCCCAACCAGGATTGCAGTAGTTTGGGGCAGTAACCCGTCTTGAGACAGAAACAAACGCAAAACAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATgg                                                                                                                                        <  2:1928644/150‑1 (MQ=255)
             cAACCAGGATTGCAGTAGTTTGGGGCAGTAACCCGTCTTGAGACAGAAACAAACGCAAAACAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTc                                                                                                                                   <  1:746760/150‑1 (MQ=255)
              aaCCAGGATTGCAGTAGTTTGGGGCAGTAACCCGTCTTGAGACAGAAACAAACGCAAACCAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCa                                                                                                                                  <  1:866223/150‑1 (MQ=255)
                ccAGGATTGCAGTAGTTTGGGGCAGTAACCCGTCTTGAGACAGAAACAAACGCAAAACAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAAc                                                                                                                                >  2:462689/1‑150 (MQ=255)
                 cAGGATTGCAGTAGTTTGGGGCAGTAACCCGTCTTGAGACAGAAACAAACGCAAAACAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACg                                                                                                                               <  2:949129/150‑1 (MQ=255)
                   ggATTGCAGTAGTTTGGGGCAGTAACCCGTCTTGAGACAGAAACAAACGCAAAACAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGcc                                                                                                                             >  1:1689789/1‑150 (MQ=255)
                   ggATTGCAGTAGTTTGGGGCAGTAACCCGTCTTGAGACAGAAACAAACGCAAAACAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGcc                                                                                                                             >  1:986358/1‑150 (MQ=255)
                      ttGCAGTAGTTTGGGGCAGTAACCCGTCTTGAGACAGAAACAAACGCAAAACAGGCCAGAAGGATATATTTCAACATTTTTAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGa                                                                                                                          <  2:1643116/150‑1 (MQ=255)
                      ttGCAGTAGTTTCGGGCAGTAACCCGTCTTGAGACAGAAACAAACGCAAAACAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGa                                                                                                                          <  2:1547355/150‑1 (MQ=255)
                        gCAGTAGTTTGGGGCAGTAACCCGTCTTGAGACAGAAACAAACGCAAAACAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATg                                                                                                                        <  2:173635/150‑1 (MQ=255)
                          agtagtTTGGGGCAGTAACCCGTCTTGAGACAGAAACAAACGCAAAACAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCg                                                                                                                      <  2:1750524/150‑1 (MQ=255)
                             agtTTGGGGCAGTAACCCGTCTTGAGACAGAAACAAACGCAAAACAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACg                                                                                                                   <  2:1738727/150‑1 (MQ=255)
                                    ggCAGTAACCCGTCTTGAGACAGAAACAAACGCAAAACAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGgcgc                                                                                                            >  2:901254/1‑150 (MQ=255)
                                          aaCCCGTCTTGAGACAGAAACAAACGCAAAACAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTa                                                                                                      <  2:858534/150‑1 (MQ=255)
                                          aaCCCGTCTTGAGACAGAAACAAACGCAAAACAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTa                                                                                                      <  1:1790854/150‑1 (MQ=255)
                                                 cTTGAGACAGAAACAAACGCAAAACAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGcc                                                                                               <  1:462689/150‑1 (MQ=255)
                                                   tGAGACAGAAACAAACGCAAAACAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATg                                                                                                                                         >  2:349707/1‑106 (MQ=255)
                                                   tGAGACAGAAACAAACGCAAAACAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATg                                                                                                                                         <  1:349707/106‑1 (MQ=255)
                                                   tGAGACAGAAACAAACGCAAAACAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTa                                                                                             >  2:327217/1‑150 (MQ=255)
                                                    gagaCAGAAACAAACGCAAAACAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCt                                                                                              <  2:452558/148‑1 (MQ=255)
                                                    gagaCAGAAACAAACGCAAAACAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCt                                                                                              >  1:452558/1‑148 (MQ=255)
                                                     agaCAGAAACAAACGCAAAACAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACa                                                                                           <  2:959280/150‑1 (MQ=255)
                                                     agaCAGAAACAAACGCAAAACAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACa                                                                                           >  2:108645/1‑150 (MQ=255)
                                                        cAGAAACAAACGCAAAACAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAcc                                                                                        >  1:1292086/1‑150 (MQ=255)
                                                        cAGAAACAAACGCAAAACAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAcc                                                                                        >  2:463938/1‑150 (MQ=255)
                                                             acaaacGCAAAACAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCt                                                                                   <  1:111605/150‑1 (MQ=255)
                                                              caaacGCAAAACAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTc                                                                                  >  1:1116986/1‑150 (MQ=255)
                                                              caaacGCAAAACAGGCAAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTc                                                                                  >  2:1072106/1‑150 (MQ=255)
                                                                aacGCAAAACAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCAt                                                                                >  2:700956/1‑150 (MQ=255)
                                                                      aaaCAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTa                                                                          >  1:1517662/1‑150 (MQ=255)
                                                                          aGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGc                                                                      <  2:1233707/150‑1 (MQ=255)
                                                                           ggCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCGACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCg                                                                     >  1:65808/1‑150 (MQ=255)
                                                                            gCCAGACGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCgg                                                                    <  2:1136963/150‑1 (MQ=255)
                                                                             ccAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGa                                                                   <  1:1025284/150‑1 (MQ=255)
                                                                              cAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTTTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGAt                                                                  <  2:1146257/150‑1 (MQ=255)
                                                                                   ggATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGc                                                             >  1:368592/1‑150 (MQ=255)
                                                                                      tatatTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCAtt                                                          >  2:449447/1‑150 (MQ=255)
                                                                                        tatTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCATTTa                                                        >  1:457306/1‑150 (MQ=255)
                                                                                         atTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCATTTAc                                                       >  2:842093/1‑150 (MQ=255)
                                                                                            tCAACATTTTGAATTTGCAAGTTTTTTGTAGGCCGGATAACGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCATTTACttt                                                    >  2:714419/1‑150 (MQ=255)
                                                                                             cAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTa                                                   <  2:1292086/150‑1 (MQ=255)
                                                                                             cAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTa                                                   <  2:88030/150‑1 (MQ=255)
                                                                                                 aTTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCATTTACATTATCAc                                               >  2:1568210/1‑150 (MQ=255)
                                                                                                   tttGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCATTTACtt                                                     <  2:1634518/142‑1 (MQ=255)
                                                                                                   tttGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCATTTACtt                                                     >  1:1634518/1‑142 (MQ=255)
                                                                                                     tGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCggg                                           <  1:1072106/150‑1 (MQ=255)
                                                                                                       aaTTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTc                                         <  2:1689789/150‑1 (MQ=255)
                                                                                                            gCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAAcc                                    >  2:377252/1‑150 (MQ=255)
                                                                                                            gCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAAcc                                    >  2:1271698/1‑150 (MQ=255)
                                                                                                             cACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACca                                   <  1:1884143/150‑1 (MQ=255)
                                                                                                               cGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACcaac                                 >  2:462172/1‑150 (MQ=255)
                                                                                                                   ttttGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAg                             <  1:40718/150‑1 (MQ=255)
                                                                                                                   ttttGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAg                             <  1:1271698/150‑1 (MQ=255)
                                                                                                                       gTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGAc                         >  1:71264/1‑150 (MQ=255)
                                                                                                                       gTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGAc                         >  2:1624607/1‑150 (MQ=255)
                                                                                                                       cattgCCGGATGCGGCGTAAACGCCTTATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGAc                         >  2:1412808/5‑150 (MQ=38)
                                                                                                                          ggCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGACCCg                      >  1:44666/1‑150 (MQ=255)
                                                                                                                          ggCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGACCCg                      >  2:1346920/1‑150 (MQ=255)
                                                                                                                           gCCGGATGCGGCGTAAACGCCTTATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGACCCGc                     >  2:1484329/1‑150 (MQ=38)
                                                                                                                           gCCGGATGCGGCGTAAACGCCTTATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGACCCGc                     >  2:8591/1‑150 (MQ=38)
                                                                                                                             cGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGACCCGCCt                   <  1:361109/150‑1 (MQ=255)
                                                                                                                              ggATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGACCCGCCtt                  <  1:700956/150‑1 (MQ=255)
                                                                                                                                 tAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGACCCGCCTTGat               <  2:986358/150‑1 (MQ=255)
                                                                                                                                  aaGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGACCCGCCTTGata              <  2:1116986/150‑1 (MQ=255)
                                                                                                                                   aGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGACCCGCCTTGatat             >  1:830638/1‑150 (MQ=255)
                                                                                                                                      cGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGACCCGCCTTGATATaaa          >  2:392797/1‑150 (MQ=255)
                                                                                                                                        ccattgccggATGCGGGAATGGTGCTCAACGCCTGATGCGCCGCTGGCGCGGCTTATAATGCCGACAACCCCCCGCATAACGTAAGGCGGATAAGGGATTTACTTTATAACCCGGTGAAACCAACAAAGCGACCCGCCGTGATATaaaaa        >  2:338458/6‑148 (MQ=255)
                                                                                                                                          tACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGACCCGCCTTGATATAAATccc      >  2:1142331/1‑150 (MQ=255)
                                                                                                                                            cgccgcATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGACCCGCCTTGATATAAATCCCAg    <  1:449447/150‑1 (MQ=255)
                                                                                                                                            cgccgcATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGACCCGCCTTGATATAAATCCCAg    >  2:1774804/1‑150 (MQ=255)
                                                                                                                                              ccgcATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGACCCGCCTTGATATAAATCCCAGGc  >  1:823332/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
GGAAATTTTACCCCAACCAGGATTGCAGTAGTTTGGGGCAGTAACCCGTCTTGAGACAGAAACAAACGCAAAACAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGACCCGCCTTGATATAAATCCCAGGC  >  NC_000913/215936‑216227

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: