Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 216,067 A→C 6.9% intergenic (+88/+112) 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,0670AC 6.9% 226.1 / 7.5 72intergenic (+88/+112)nlpE/yaeFlipoprotein NlpE/peptidase C92 family protein
Reads supporting (aligned to +/- strand):  ref base A (31/36);  new base C (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 = 9.92e-01

TGGCACCCGATACGGCAGGGAAATTTTACCCCAACCAGGATTGCAGTAGTTTGGGGCAGTAACCCGTCTTGAGACAGAAACAAACGCAAAACAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGACCCGCCTTGATAT  >  NC_000913/215918‑216216
                                                                                                                                                     |                                                                                                                                                     
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: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                                                                                                                              <  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:108645/1‑150 (MQ=255)
                                                                       agaCAGAAACAAACGCAAAACAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACa                                                                                <  2:959280/150‑1 (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:88030/150‑1 (MQ=255)
                                                                                                               cAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTa                                        <  2:1292086/150‑1 (MQ=255)
                                                                                                                   aTTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCATTTACATTATCAc                                    >  2:1568210/1‑150 (MQ=255)
                                                                                                                     tttGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCATTTACtt                                          >  1:1634518/1‑142 (MQ=255)
                                                                                                                     tttGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCATTTACtt                                          <  2:1634518/142‑1 (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:1271698/150‑1 (MQ=255)
                                                                                                                                     ttttGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAg                  <  1:40718/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)
                                                                                                                                                     aGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGACCCGCCTTGatat  >  1:830638/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGGCACCCGATACGGCAGGGAAATTTTACCCCAACCAGGATTGCAGTAGTTTGGGGCAGTAACCCGTCTTGAGACAGAAACAAACGCAAAACAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGACCCGCCTTGATAT  >  NC_000913/215918‑216216

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: