Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 216,066 A→G 6.9% intergenic (+87/+113) 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,0660AG 6.9% 218.0 / 5.4 72intergenic (+87/+113)nlpE/yaeFlipoprotein NlpE/peptidase C92 family protein
Reads supporting (aligned to +/- strand):  ref base A (31/36);  new base G (3/2);  total (34/38)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.76e-01

TGGCACCCGATACGGCAGGGAAATTTTACCCCAACCAGGATTGCAGTAGTTTGGGGCAGTAACCCGTCTTGAGACAGAAACAAACGCAAAACAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGACCCGCCTTGATA  >  NC_000913/215918‑216215
                                                                                                                                                    |                                                                                                                                                     
tGGCACCCGATACGGCAGGGAAATTTTACCCCAACCAGGATTGCAGTAGTTTGGGGCAGTAACCCGTCTTGAGACAGAAACAAACGCAAAACAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATaa                                                                                                                                                      >  1:858534/1‑150 (MQ=255)
  gCACCCGATACGGCAGGGAAATTTTACCCCAACCAGGATTGCAGTAGTTTGGGGCAGTAACCCGTCTTGAGACAGAAACAAACGCAAAACAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAgg                                                                                                                                                    <  1:1373189/150‑1 (MQ=255)
     cccGATACGGCAGGGAAATTTTACCCCAACCAGGATTGCAGTAGTTTGGGGCAGTAACCCGTCTTGAGACAGAAACAAACGCAAAACAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGt                                                                                                                                                 <  1:1188942/150‑1 (MQ=255)
      ccGATACGGCAGGGAAATTTTACCCCAACCAGGATTGCAGTAGTTTGGGGCAGTAACCCGTCTTGAGACAGAAACAAACGCAAAACAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGtt                                                                                                                                                <  1:1491096/150‑1 (MQ=255)
       cGATACGGCAGGGAAATTTTACCCCAACCAGGATTGCAGTAGTTTGGGGCAGTAACCCGTCTTGAGACAGAAACAAACGCAAAACAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGttt                                                                                                                                               <  2:323394/150‑1 (MQ=255)
                  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:986358/1‑150 (MQ=255)
                                     ggATTGCAGTAGTTTGGGGCAGTAACCCGTCTTGAGACAGAAACAAACGCAAAACAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGcc                                                                                                                 >  1:1689789/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                                                                                                                             <  1:349707/106‑1 (MQ=255)
                                                                     tGAGACAGAAACAAACGCAAAACAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATg                                                                                                                             >  2:349707/1‑106 (MQ=255)
                                                                     tGAGACAGAAACAAACGCAAAACAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTa                                                                                 >  2:327217/1‑150 (MQ=255)
                                                                      gagaCAGAAACAAACGCAAAACAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCt                                                                                  >  1:452558/1‑148 (MQ=255)
                                                                      gagaCAGAAACAAACGCAAAACAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCt                                                                                  <  2:452558/148‑1 (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:1271698/1‑150 (MQ=255)
                                                                                                                              gCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAAcc                        >  2:377252/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             >  2:1624607/1‑150 (MQ=255)
                                                                                                                                         gTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGAc             >  1:71264/1‑150 (MQ=255)
                                                                                                                                         cattgCCGGATGCGGCGTAAACGCCTTATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGAc             >  2:1412808/5‑150 (MQ=38)
                                                                                                                                            ggCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGACCCg          >  2:1346920/1‑150 (MQ=255)
                                                                                                                                            ggCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGACCCg          >  1:44666/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)
                                                                                                                                                    |                                                                                                                                                     
TGGCACCCGATACGGCAGGGAAATTTTACCCCAACCAGGATTGCAGTAGTTTGGGGCAGTAACCCGTCTTGAGACAGAAACAAACGCAAAACAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGACCCGCCTTGATA  >  NC_000913/215918‑216215

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: