Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 216,074 T→A 7.0% intergenic (+95/+105) 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,0740TA 7.0% 222.1 / 5.2 71intergenic (+95/+105)nlpE/yaeFlipoprotein NlpE/peptidase C92 family protein
Reads supporting (aligned to +/- strand):  ref base T (33/33);  new base A (3/2);  total (36/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 = 9.81e-01

CGATACGGCAGGGAAATTTTACCCCAACCAGGATTGCAGTAGTTTGGGGCAGTAACCCGTCTTGAGACAGAAACAAACGCAAAACAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGACCCGCCTTGATATAAATCCC  >  NC_000913/215925‑216223
                                                                                                                                                     |                                                                                                                                                     
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                                                                                                  <  1:1790854/150‑1 (MQ=255)
                                                     aaCCCGTCTTGAGACAGAAACAAACGCAAAACAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTa                                                                                                  <  2:858534/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                                                                                          >  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                                                                                    >  2:463938/1‑150 (MQ=255)
                                                                   cAGAAACAAACGCAAAACAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAAcc                                                                                    >  1:1292086/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                                                 <  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                     >  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)
                                                                                                                                                 cGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGACCCGCCTTGATATaaa      >  2:392797/1‑150 (MQ=255)
                                                                                                                                                     tACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGACCCGCCTTGATATAAATccc  >  2:1142331/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGATACGGCAGGGAAATTTTACCCCAACCAGGATTGCAGTAGTTTGGGGCAGTAACCCGTCTTGAGACAGAAACAAACGCAAAACAGGCCAGAAGGATATATTTCAACATTTTGAATTTGCACGTTTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATGGTGCTCAACGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACAACCCCCCTCATACCTTAAGGCGGATAAGGCATTTACTTTATCACCGGGTCAAACCAACAAAGCGACCCGCCTTGATATAAATCCC  >  NC_000913/215925‑216223

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: